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.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

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}$.

With the aim of characterizing the dynamical processes involved in the formation of young protostars, we present high angular resolution ALMA dust polarization observations of the Class 0 protostellar cores Serpens SMM1, Emb 8(N), and Emb 8. With spatial resolutions ranging from 150 to 40 au at 870 {\mu}m, we find unexpectedly high values of the polarization fraction along the outflow cavity walls in Serpens Emb8(N). We use 3 mm and 1 mm molecular tracers to investigate outflow and dense gas properties and their correlation with the polarization. These observations allow us to investigate the physical processes involved in the Radiative Alignment Torques (RATs) acting on dust grains along the outflow cavity walls, which experience irradiation from accretion processes and outflow shocks. The inner core of SMM1-a presents a polarization pattern with a poloidal magnetic field at the bases of the two lobes of the bipolar outflow. To the south of SMM1-a we see two polarized filaments, one of which seems to trace the redshifted outflow cavity wall. The other may be an accretion streamer of material infalling onto the central protostar. We propose that the polarized emission we see at millimeter wavelengths along the irradiated cavity walls can be reconciled with the expectations of RAT theory if the aligned grains present at '<' 500 au scales in Class 0 envelopes have grown larger than the 0.1 {\mu}m size of ISM dust grains. Our observations allow us to constrain the star-forming sources magnetic field morphologies within the central cores, along the outflow cavity walls, and in possible accretion streamers.

Context. Submillimeter Array (SMA) 870 micron polarization observations of the hot molecular core G31.41+0.31 revealed one of the clearest examples up to date of an hourglass-shaped magnetic field morphology in a high-mass star-forming region. Aims. To better establish the role that the magnetic field plays in the collapse of G31.41+0.31, we carried out Atacama Large Millimeter/submillimeter Array (ALMA) observations of the polarized dust continuum emission at 1.3 mm with an angular resolution four times higher than that of the previous (sub)millimeter observations to achieve an unprecedented image of the magnetic field morphology. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6). The resulting synthesized beam is 0.28"x0"20 which, at the distance of the source, corresponds to a spatial resolution of ~875 au. Results. The observations resolve the structure of the magnetic field in G31.41+0.31 and allow us to study the field in detail. The polarized emission in the Main core of G31.41+0.41is successfully fit with a semi-analytical magnetostatic model of a toroid supported by magnetic fields. The best fit model suggests that the magnetic field is well represented by a poloidal field with a possible contribution of a toroidal component of ~10% of the poloidal component, oriented southeast to northwest at ~ -44 deg and with an inclination of ~-45 degr. The magnetic field is oriented perpendicular to the northeast to southwest velocity gradient detected in this core on scales from 1E3-1E4 au. This supports the hypothesis that the velocity gradient is due to rotation and suggests that such a rotation has little effect on the magnetic field. The strength of the magnetic field estimated in the central region of the core with the Davis-Chandrasekhar-Fermi method is ~8-13 mG and implies that the mass-to-flux ratio in this region is slightly supercritical ...

L1521F is found to be forming multiple cores and it is cited as an example of the densest core with an embedded VeLLO in a highly dynamical environment. We present the core-scale magnetic fields (B-fields) in the near vicinity of the VeLLO L1521F-IRS using submm polarization measurements at 850$~\mu$m using JCMT POL-2. This is the first attempt to use high-sensitivity observations to map the sub-parsec scale B-fields in a core with a VeLLO. The B-fields are ordered and very well connected to the parsec-scale field geometry seen in our earlier optical polarization observations and the large-scale structure seen in Planck dust polarization. The core scale B-field strength estimated using Davis-Chandrasekhar-Fermi relation is $\rm 330\pm100~\mu$G which is more than ten times of the value we obtained in the envelope (envelope in this paper is "core envelope"). This indicates that B-fields are getting stronger on smaller scales. The magnetic energies are found to be 1 to 2 orders of magnitude higher than non-thermal kinetic energies in the envelope and core. This suggests that magnetic fields are more important than turbulence in the energy budget of L1521F. The mass-to-flux ratio of 2.3$\pm$0.7 suggests that the core is magnetically-supercritical. The degree of polarization is steadily decreasing towards the denser part of the core with a power law slope of -0.86.

We present ALMA observations at 233, 678, and 870 GHz of the Galactic Center black hole, Sagittarius A*. These observations reveal a flat spectrum over this frequency range with spectral index $\alpha \approx -0.3$, where the flux density $S \propto \nu^\alpha$. We model the submm and far infrared spectrum with a one zone synchrotron model of thermal electrons. We infer electron densities $n = (2-5) \times 10^6$ cm$^{-3}$, electron temperatures $T_e = (1-3) \times 10^{11}$ K, and magnetic field strength $B = 10-50$ G. The parameter range can be further constrained using the observed quiescent X-ray luminosity. The flat submm spectrum results in a high electron temperature and implies that the emitting electrons are efficiently heated. We also find that the emission is most likely optically thin at 233 GHz. These results indicate that millimeter and submillimeter wavelength very long baseline interferometry of Sgr A* including those of the Event Horizon Telescope should see a transparent emission region down to event horizon scales.

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.

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process.

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.

Millimeter polarimetry of Sgr A* probes the linearly polarized emission region on a scale of $\sim 10$ Schwarzschild radii ($R_S$) as well as the dense, magnetized accretion flow on scales out to the Bondi radius ($\sim 10^5 R_S$) through Faraday rotation. We present here multi-epoch ALMA Band 6 (230 GHz) polarimetry of Sgr A*. The results confirm a mean rotation measure, ${\rm RM} \approx -5 \times 10^5 {\rm\ rad\ m^{-2}}$, consistent with measurements over the past 20 years and support an interpretation of the RM as originating from a radiatively inefficient accretion flow (RIAF) with $\dot{M} \approx 10^{-8} { \rm\ M_{\odot}\ y^{-1} }$. Variability is observed for the first time in the RM on time scales that range from hours to months. The long-term variations may be the result of changes in the line of sight properties in a turbulent accretion flow. Short-term variations in the apparent RM are not necessarily the result of Faraday rotation and may be the result of complex emission and propagatation effects close to the black hole, some of which have been predicted in numerical modeling. We also confirm the detection of circular polarization at a mean value of $-1.1 \pm 0.2 \%$. It is variable in amplitude on time scales from hours to months but the handedness remains unchanged from that observed in past centimeter- and millimeter-wavelength detections. These results provide critical constraints for the analysis and interpretation of Event Horizon Telescope data of Sgr A*, M87, and similar sources.

Accurately measuring the magnetic field in the ISM is essential for understanding star-formation processes. We searched archival data of the Submillimeter Array (SMA) for evidence of circular polarization in common molecular tracers, most notably CO. This circular polarization possibly arises from anisotropic resonant scattering, which would imply that some background linearly polarized radiation is being converted to circular polarization. We find circular polarization in the star-forming regions NGC7538 (in CO) and Orion KL (in CO and SiO), as well as in the carbon star IRC+10216 (in CS, SiS, H$^{13}$CN and CO) at high enough levels to suggest that the presence of circular polarization in these spectral lines is common for such objects. This implies that measuring circular polarization is important when studying magnetic fields through the linear polarization of molecular spectral lines in the interstellar medium. We also provide a simple Python wrapper for the Miriad data reduction package.

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.

Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87mm, in order to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes, including both single protostars and multiple systems, in nearby clouds. Our SMA observations probe the envelope emission at scales of 600-5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. We report the detection of linearly polarized dust continuum emission in all of our targets, with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources where strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at about 100 au scales.

We present SMA and NOEMA observations of the host galaxy of FRB 121102 in the CO 3-2 and 1-0 transitions, respectively. We do not detect emission from either transition. We set $3\sigma$ upper limits to the CO luminosity $L_{CO} < 2.5 \times 10^7\,{\rm K\,km\,s}^{-1} {\, \rm pc^{-2}}$ for CO 3-2 and $L_{CO} < 2.3 \times 10^9\, {\rm K\,km\,s}^{-1} {\, \rm pc^{-2}}$ for CO 1-0. For Milky-Way-like star formation properties, we set a $3\sigma$ upper limit on the $H_2$ mass of $2.5 \times 10^8 \rm\ M_{\odot}$, slightly less than the predictions for the $H_2$ mass based on the star formation rate. The true constraint on the $H_2$ mass may be significantly higher, however, because of the reduction in CO luminosity that is common forlow-metallicity dwarf galaxies like the FRB host galaxy. These results demonstrate the challenge of identifying the nature of FRB progenitors through study of the host galaxy molecular gas. We also place a limit of 42 $\mu$Jy ($3\sigma$) on the continuum flux density of the persistent radio source at 97 GHz, consistent with a power-law extrapolation of the low frequency spectrum, which may arise from an AGN or other nonthermal source.

Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperature, the intensity profile, and the polarization properties clearly indicate that the disk is optically thick for a radius of R<170 au. The linear polarization of the dust emission is detected almost all along the disk and its properties suggest that dust polarization is produced mainly by self-scattering. However, the polarization pattern presents a clear differentiation between the inner (optically thick) part of the disk and the outer (optically thin) region of the disk, with a sharp transition that occurs at a radius of 0.1 arcsec (~170 au). The polarization characteristics of the inner disk suggest that dust settling has not occurred yet with a maximum dust grain size between 50 and 500 mum. The outer part of the disk has a clear azimuthal pattern but with a significantly higher polarization fraction compared to the inner disk. This pattern is broadly consistent with self-scattering of a radiation field that is beamed radially outward, as expected in the optically thin outer region, although contribution from non-spherical grains aligned with respect to the radiative flux cannot be excluded.

The role of the magnetic field during protostellar collapse is poorly constrained from an observational point of view, although it could be significant if we believe state-of-the-art models of protostellar formation. We present polarimetric observations of the 233 GHz thermal dust continuum emission obtained with ALMA in the B335 Class 0 protostar. Linearly polarized dust emission arising from the circumstellar material in the envelope of B335 is detected at all scales probed by our observations, from radii of 50 to 1000 au. The magnetic field structure producing the dust polarization has a very ordered topology in the inner envelope, with a transition from a large-scale poloidal magnetic field, in the outflow direction, to strongly pinched in the equatorial direction. This is probably due to magnetic field lines being dragged along the dominating infall direction since B335 does not exhibit prominent rotation. Our data and their qualitative comparison to a family of magnetized protostellar collapse models show that, during the magnetized collapse in B335, the magnetic field is maintaining a high level of organization from scales 1000 au to 50 au: this suggests the field is dynamically relevant and capable of influencing the typical outcome of protostellar collapse, such as regulating the disk size in B335.

Context: The importance of magnetic fields at the onset of star formation related to the early fragmentation and collapse processes is largely unexplored today. Aims: We want to understand the magnetic field properties at the earliest evolutionary stages of high-mass star formation. Methods: The Atacama Large Millimeter Array is used at 1.3mm wavelength in full polarization mode to study the polarized emission and by that the magnetic field morphologies and strengths of the high-mass starless region IRDC18310-4. Results: The polarized emission is clearly detected in four sub-cores of the region. In general it shows a smooth distribution, also along elongated cores. Estimating the magnetic field strength via the Davis-Chandrasekhar-Fermi method and following a structure function analysis, we find comparably large magnetic field strengths between ~0.6 and 3.7mG. Comparing the data to spectral line observations, the turbulent-to-magnetic energy ratio is low, indicating that turbulence does not significantly contribute to the stability of the gas clump. A mass-to-flux ratio around the critical value 1.0 - depending on column density - indicates that the region starts to collapse which is consistent with the previous spectral line analysis of the region. Conclusions: While this high-mass region is collapsing and thus at the verge of star formation, the high magnetic field values and the smooth spatial structure indicate that the magnetic field is important for the fragmentation and collapse process. This single case study can only be the starting point for larger sample studies of magnetic fields at the onset of star formation.

Millimeter wavelength polarimetry of accreting black hole systems can provide a tomographic probe of the accretion flow on a wide range of linear scales. We searched for linear polarization in two low luminosity active galactic nuclei (LLAGN), M81 and M84, using the Combined Array for Millimeter Astronomy (CARMA) and the Submillimeter Array (SMA). We find upper limits of $\sim 1 - 2\%$ averaging over the full bandwidth and with a rotation measure (RM) synthesis technique. These low polarization fractions, along with similar low values for LLAGN M87 and 3C84, suggest that LLAGN have qualitatively different polarization properties than radio-loud sources and Sgr A*. If the sources are intrinsically polarized and then depolarized by Faraday rotation then we place lower limits on the RM of a few times $10^7\, {\rm rad\, m^{-2}}$ for the full bandwidth case and $\sim 10^9\, {\rm rad\, m^{-2}}$ for the RM synthesis analysis. These limits are inconsistent with or marginally consistent with expected accretion flow properties. Alternatively, the sources may be depolarized by cold electrons within a few Schwarzschild radii from the black hole, as suggested by numerical models.

The mechanism for producing polarized emission from protostellar disks at (sub)millimeter wavelengths is currently uncertain. Classically, polarization is expected from non-spherical grains aligned with the magnetic field. Recently, two alternatives have been suggested. One polarization mechanism is caused by self-scattering from dust grains of sizes comparable to the wavelength while the other mechanism is due to grains aligned with their short axes along the direction of radiation anisotropy. The latter has recently been shown as a likely mechanism for causing the dust polarization detected in HL Tau at 3.1 mm. In this paper, we present ALMA polarization observations of HL Tau for two more wavelengths: 870 $\mu$m and 1.3 mm. The morphology at 870 $\mu$m matches the expectation for self-scattering, while that at 1.3 mm shows a mix between self-scattering and grains aligned with the radiation anisotropy. The observations cast doubt on the ability of (sub)millimeter continuum polarization to probe disk magnetic fields for at least HL Tau. By showing two distinct polarization morphologies at 870 $\mu$m and 3.1 mm and a transition between the two at 1.3 mm, this paper provides definitive evidence that the dominant (sub)millimeter polarization mechanism transitions with wavelength. In addition, if the polarization at 870 $\mu$m is due to scattering, the lack of polarization asymmetry along the minor axis of the inclined disk implies that the large grains responsible for the scattering have already settled into a geometrically thin layer, and the presence of asymmetry along the major axis indicates that the HL Tau disk is not completely axisymmetric.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

We present subarcsecond angular resolution observations carried out with the Submillimeter Array (SMA) at 880 $\mu$m centered at the B0-type protostar GGD27~MM1, the driving source of the parsec scale HH 80-81 jet. We constrain its polarized continuum emission to be $\lesssim0.8\%$ at this wavelength. Its submm spectrum is dominated by sulfur-bearing species tracing a rotating disk--like structure (SO and SO$_2$ isotopologues mainly), but also shows HCN-bearing and CH$_3$OH lines, which trace the disk and the outflow cavity walls excavated by the HH 80-81 jet. The presence of many sulfurated lines could indicate the presence of shocked gas at the disk's centrifugal barrier or that MM1 is a hot core at an evolved stage. The resolved SO$_2$ emission traces very well the disk kinematics and we fit the SMA observations using a thin-disk Keplerian model, which gives the inclination (47$^{\circ}$), the inner ($\lesssim170$ AU) and outer ($\sim950-1300$~AU) radii and the disk's rotation velocity (3.4 km s$^{-1}$ at a putative radius of 1700 AU). We roughly estimate a protostellar dynamical mass of 4-18\msun. MM2 and WMC cores show, comparatively, an almost empty spectra suggesting that they are associated with extended emission detected in previous low-angular resolution observations, and therefore indicating youth (MM2) or the presence of a less massive object (WMC).

We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the James Clerk Maxwell Telescopes (JCMT), we can compare the magnetic field orientations at different spatial scales. We find major changes in the magnetic field orientation between large (~0.1 pc) scales -- where the magnetic field is oriented E-W, perpendicular to the major axis of the dusty filament where SMM1 is embedded -- and the intermediate and small scales probed by CARMA (~1000 AU resolution), the SMA (~350 AU resolution), and ALMA (~140 AU resolution). The ALMA maps reveal that the redshifted lobe of the bipolar outflow is shaping the magnetic field in SMM1 on the southeast side of the source; however, on the northwestern side and elsewhere in the source, low velocity shocks may be causing the observed chaotic magnetic field pattern. High-spatial-resolution continuum and spectral-line observations also reveal a tight (~130 AU) protobinary system in SMM1-b, the eastern component of which is launching an extremely high-velocity, one-sided jet visible in both CO(2-1) and SiO(5-4); however, that jet does not appear to be shaping the magnetic field. These observations show that with the sensitivity and resolution of ALMA, we can now begin to understand the role that feedback (e.g., from protostellar outflows) plays in shaping the magnetic field in very young, star-forming sources like SMM1.

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 report new Karl G. Jansky Very Large Array (JVLA), 0$"$.5 angular resolution observations of linearly polarized continuum emission at 6.9 mm, towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A. This target source is a collapsing dense molecular core, which was resolved at short wavelengths to have hourglass shaped B-field configuration. We compare these 6.9 mm observations with previous polarization Submillimeter Array (SMA) observations at 0.88 mm, which have comparable angular resolution ($\sim$0$"$7). We found that at the same resolution, the observed polarization position angles at 6.9 mm are slightly deviated from those observed at 0.88 mm. Due to the lower optical depth of the emission at 6.9 mm, and the potential effect of dust grain growth, the new JVLA observations are likely probing B-field alignments in regions interior to those sampled by the previous polarization observations at higher frequencies. Our understanding can be improved by more sensitive observations, and observations for the more extended spatial scales.

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.

We analyze the light curves of 413 radio sources at submillimeter wavelengths using data from the Submillimeter Array calibrator database. The database includes more than 20,000 observations at 1.3 and 0.8 mm that span 13 years. We model the light curves as a damped random walk and determine a characteristic time scale $\tau$ at which the variability amplitude saturates. For the vast majority of sources, primarily blazars and BL Lac objects, we find only lower limits on $\tau$. For two nearby low luminosity active galactic nuclei, M81 and M87, however, we measure $\tau=1.6^{+3.0}_{-0.9}$ days and $\tau=45^{+61}_{-24}$ days, respectively ($2\sigma$ errors). Including the previously measured $\tau=0.33\pm 0.16$ days for Sgr A*, we show an approximately linear correlation between $\tau$ and black hole mass for these nearby LLAGN. Other LLAGN with spectra that peak in the submm are expected to follow this correlation. These characteristic time scales are comparable to the minimum time scale for emission processes close to an event horizon, and suggest that the underlying physics may be independent of black hole mass, accretion rate, and jet luminosity.

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.

We report the discovery of Faraday rotation toward radio source 3C 84, the active galactic nucleus in NGC1275 at the core of the Perseus Cluster. The rotation measure (RM), determined from polarization observations at wavelengths of 1.3 and 0.9 mm, is (8.7 +/- 2.3) x 10^5 radians/m^2, among the largest ever measured. The RM remained relatively constant over a 2 year period even as the intrinsic polarization position angle wrapped through a span of 300 degrees. The Faraday rotation is likely to originate either in the boundary layer of the radio jet from the nucleus, or in the accretion flow onto the central black hole. The accretion flow probably is disk-like rather than spherical on scales of less than a parsec, otherwise the RM would be even larger.

Magnetic fields in accretion disks play a dominant role during the star formation process but have hitherto been observationally poorly constrained. Field strengths have been inferred on T Tauri stars themselves and possibly in the innermost part of the accretion disk, but the strength and morphology of the field in the bulk of the disk have not been observed. Unresolved measurements of polarized emission (arising from elongated dust grains aligned perpendicular to the field) imply average fields aligned with the disks. Theoretically, the fields are expected to be largely toroidal, poloidal, or a mixture of the two, which imply different mechanisms for transporting angular momentum in the disks of actively accreting young stars such as HL Tau. Here we report resolved measurements of the polarized 1.25 mm continuum emission from HL Tau's disk. The magnetic field on a scale of 80 AU is coincident with the major axis (~210 AU diameter) of the disk. From this we conclude that the magnetic field inside the disk at this scale cannot be dominated by a vertical component, though a purely toroidal field does not fit the data well either. The unexpected morphology suggests that the magnetic field's role for the accretion of a T Tauri star is more complex than the current theoretical understanding.

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.

Context: NGC 7538 IRS 1-3 is a high-mass star-forming cluster with several detected dust cores, infrared sources, (ultra)compact H$_{\rm II}$ regions, molecular outflows, and masers. In such a complex environment, important interactions and feedback among the embedded objects are expected to play a major role in the evolution of the region. Aims: We study the dust, kinematic, and polarimetric properties of the NGC 7538 IRS 1-3 region to investigate the role of the different forces interplaying in the formation and evolution of high-mass star-forming clusters. Methods: We perform SMA high angular resolution observations at 880 $\mu$m with the compact configuration. We develop the RATPACKS code to generate synthetic velocity cubes from models of choice to be compared to the observational data. We develop the "mass balance" analysis to quantify the stability against gravitational collapse accounting for all the energetics at core scales. Results: We detect 14 dust cores from 3.5 $M_{\odot}$ to 37 $M_{\odot}$ arranged in two larger scale structures: a central bar and a filamentary spiral arm. The spiral arm presents large scale velocity gradients in H$^{13}$CO$^+$ 4-3 and C$^{17}$O 3-2, and magnetic field segments well aligned to the dust main axis. The velocity gradient is well reproduced by a spiral arm expanding at 9 km s$^{-1}$ with respect to the central core MM1, which is known to power a large precessing outflow. The energy of the outflow is comparable with the spiral arm kinetic energy, which is dominant over gravitational and magnetic energies. In addition, the dynamical ages of the outflow and spiral arm are comparable. ... (Full abstract in the pdf version)

We present the first constraint on Faraday rotation measure (RM) at submillimeter wavelengths for the nucleus of M 87. By fitting the polarization position angles ($\chi$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $\chi$ as a result of \emph{external} Faraday rotation associated with accretion flow, we determine the rotation measure of the M 87 core to be between $-$7.5$\times$10$^{5}$ and 3.4$\times$10$^{5}$ rad/m$^{2}$. Assuming a density profile of the accretion flow that follows a power-law distribution and a magnetic field that is ordered, radial, and has equipartition strength, the limit on the rotation measure constrains the mass accretion rate $\dot{M}$ to be below 9.2$\times$10$^{-4}$ M$_{\odot}$~yr$^{-1}$ at a distance of 21 Schwarzchild radii from the central black hole. This value is at least two orders of magnitude smaller than the Bondi accretion rate, suggesting significant suppression of the accretion rate in the inner region of the accretion flow. Consequently, our result disfavors the classical \emph{advection dominated accretion flow} (ADAF) and prefers the \emph{adiabatic inflow-outflow solution} (ADIOS) or \emph{convection-dominated accretion flow} (CDAF) for the hot accretion flow in M 87.

We present submillimeter spectral line and dust continuum polarization observations of a remarkable hot core and multiple outflows in the high-mass star-forming region W43-MM1 (G30.79 FIR 10), obtained using the Submillimeter Array (SMA). A temperature of $\sim$ 400 K is estimated for the hot-core using CH$_3$CN (J=19-18) lines, with detections of 11 K-ladder components. The high temperature and the mass estimates for the outflows indicate high-mass star-formation. The continuum polarization pattern shows an ordered distribution, and its orientation over the main outflow appears aligned to the outflow. The derived magnetic field indicates slightly super-critical conditions. While the magnetic and outflow energies are comparable, the B-field orientation appears to have changed from parsec scales to $\sim$ 0.1 pc scales during the core/star-formation process.

We present CO 3-2, SiO 8-7, C34S 7-6, and 878 mum dust continuum subarcsecond angular resolution observations with the SMA toward IRAS 16293-2422 (I16293). The C34S emission traces the 878 mum dust continuum well, and clearly shows a smooth velocity gradient along the major axis of component I16293A. The CO shows emission at moderate high velocities arising from two bipolar outflows, which appear to be perpendicular with respect to each other. The high sensitivity and higher angular resolution of these observations allows us to pinpoint well the origin of these two outflows at the center of component I16293A. Interestingly, the most compact outflow appears to point toward I16293B. Our data show that the previously reported monopolar blueshifted CO outflow associated with component I16293B seems to be part of the compact outflow arising from component I16293A. In addition, the SiO emission is also tracing this compact outflow: on one hand, the SiO emission appears to have a jet-like morphology along the southern redshifted lobe; on the other hand, the SiO emission associated with the blueshifted northern lobe traces a well defined arc on the border of component I16293B facing I16293A. The blueshifted CO lobe of the compact outflow splits into two lobes around the position of this SiO arc. All these results lead us to propose that the compact outflow from component I16293A is impacting on the circumstellar gas around component I16293B, possibly being diverged as a consequence of the interaction.

We present subarcsecond resolution polarimetric observations of the 878 mum thermal dust continuum emission obtained with the Submillimeter Array (SMA) towards the IRAS 16293-2422 protostellar binary system. We report the detection of linearly polarized dust emission arising from the circumstellar disk associated with the IRAS 16293-2422 B protostar. The fractional polarization of ~1.4$% is only slightly lower than that expected from theoretical calculations in such disks. The magnetic field structure on the plane of the sky derived from the dust polarization suggests a complex magnetic field geometry in the disk, possibly associated with a rotating disk that is wrapping the field lines as expected from the simulations. The polarization around IRAS 16293-2422 A at subarcsecond angular resolution is only marginally detected.

We present spectro-polarimetric observations of several molecular lines obtained with the Submillimeter Array (SMA) toward the carbon rich AGB star IRC+10216. We have detected and mapped the linear polarization of the CO 3-2, SiS 19-18 and CS 7-6 lines. The polarization arises at a distance of ~450 AU from the star and is blueshifted with respect the Stokes I. The SiS 19-18 polarization pattern appears to be consistent with a locally radial magnetic field configuration. However, the CO 3-2 and CS 7-6 line polarization suggests an overall complex magnetic field morphology within the envelope. This work demonstrates the feasibility of using spectro-polarimetric observations to carry out tomographic imaging of the magnetic field in circumstellar envelopes.

We present the first interferometric polarization map of the W3(OH) massive star-forming region observed with the Submillimeter Array (SMA) at 878 mum with an angular resolution of 1.5 (about 3 \times 10 AU). Polarization is detected in the W3(H2O) hot core, an extended emission structure in the north-west of W3(H2O), and part of the W3(OH) ultracompact HII region. The W3(H2O) hot core is known to be associated with a synchrotron jet along the east-west direction. In this core, the inferred magnetic field orientation is well aligned with the synchrotron jet and close to the plane of sky. Using the Chandrasekhar-Fermi method with the observed dispersion in polarization angle, we estimate a plane-of-sky magnetic field strength of 17.0 mG. Combined with water maser Zeeman measurements, the total magnetic field strength is estimated to be 17.1 mG, comparable to the field strength estimated from the synchrotron model. The magnetic field energy dominates over turbulence in this core. In addition, the depolarization effect is discerned in both SMA and JCMT measurements. Despite the great difference in angular resolutions and map extents, the polarization percentage shows a similar power-law dependence with the beam averaged column density. We suggest that the column density may be an important factor to consider when interpreting the depolarization effect.

Molecular line polarization is a unique source of information about the magnetic fields and anisotropies in the circumstellar envelopes of evolved stars. Here we present the first detection of thermal CO(J=2-1) and SiO(J=5-4, v=0) polarization, in the envelope of the asymptotic giant branch star IK Tau. The observed polarization direction does not match predictions for circumstellar envelope polarization induced only by an anisotropic radiation field. Assuming that the polarization is purely due to the Goldreich-Kylafis effect, the linear polarization direction is defined by the magnetic field as even the small Zeeman splitting of CO and SiO dominates the molecular collisional and spontaneous emission rates. The polarization was mapped using the Submillimeter Array (SMA) and is predominantly north-south. There is close agreement between the CO and SiO observations, even though the CO polarization arises in the circumstellar envelope at ~800 AU and the SiO polarization at <250 AU. If the polarization indeed traces the magnetic field, we can thus conclude that it maintains a large-scale structure throughout the circumstellar envelope. We propose that the magnetic field, oriented either east-west or north-south is responsible for the east-west elongation of the CO distribution and asymmetries in the dust envelope. In the future, the Atacama Large Millimeter/submillimeter Array will be able to map the magnetic field using CO polarization for a large number of evolved stars.

We report the first detections of circularly polarized emission at submillimeter wavelengths from the compact radio source and supermassive black hole candidate Sgr A* at a level of 1.2\pm0.3% at 1.3 mm wavelength (230 GHz) and 1.6\pm0.3% at 860 microns (345 GHz) with the same handedness as observed at lower frequencies (1.4-15 GHz). The observations, taken with the Submillimeter Array in multiple epochs, also show simultaneous linear polarization (LP) at both wavelengths of about 6%. These properties differ sharply from those at wavelengths longer than 1 cm (frequencies below 30 GHz), where weak circular polarization (CP) (~ 0.5%) dominates over LP, which is not detected at similar fractional limits. We describe an extensive set of tests to ensure the accuracy of our measurements. We find no CP in any other source, including the bright quasar 1924-292, which traces the same path on the sky as Sgr A* and therefore should be subject to identical systematic errors originating in the instrument frame. Since a relativistic synchrotron plasma is expected to produce little CP, the observed CP is probably generated close to the event horizon by the Faraday conversion process. We use a simple model to show that the phase shift associated with Faraday conversion can be nearly independent of frequency, a sufficient condition to make the handedness of CP independent of frequency. Because the size of the tau=1-surface changes by more than an order of magnitude between 1.4 and 345 GHz, the magnetic field must be coherent over such scales to consistently produce left CP.

We apply our technique on the dispersion of magnetic fields in molecular clouds to high spatial resolution Submillimeter Array polarization data obtained for Orion KL in OMC-1, IRAS 16293, and NGC 1333 IRAS 4A. We show how one can take advantage of such high resolution data to characterize the magnetized turbulence power spectrum in the inertial and dissipation ranges. For Orion KL we determine that in the inertial range the spectrum can be approximately fitted with a power law k^-(2.9\pm0.9) and we report a value of 9.9 mpc for {\lambda}_AD, the high spatial frequency cutoff presumably due to turbulent ambipolar diffusion. For the same parameters we have \sim k^-(1.4\pm0.4) and a tentative value of {\lambda}_AD \simeq 2.2 mpc for NGC 1333 IRAS 4A, and \sim k^-(1.8\pm0.3) with an upper limit of {\lambda}_AD < 1.8 mpc for IRAS 16293. We also discuss the application of the technique to interferometry measurements and the effects of the inherent spatial filtering process on the interpretation of the results.

We used the Submillimeter Array (SMA) to observe the thermal polarized dust emission from the protostellar source NGC 2024 FIR 5. The polarized emission outlines a partial hourglass morphology for the plane-of-sky component of the core magnetic field. Our data are consistent with previous BIMA maps, and the overall magnetic field geometries obtained with both instruments are similar. We resolve the main core into two components, FIR 5A and FIR 5B. A possible explanation for the asymmetrical field lies in depolarization effects due to the lack of internal heating from FIR 5B source, which may be in a prestellar evolutionary state. The field strength was estimated to be 2.2 mG, in agreement with previous BIMA data. We discuss the influence of a nearby H{\sc ii} region over the field lines at scales of $\sim 0.01$ pc. Although the hot component is probably compressing the molecular gas where the dust core is embedded, it is unlikely that the radiation pressure exceeds the magnetic tension. Finally, a complex outflow morphology is observed in CO (3 $\rightarrow$ 2) maps. Unlike previous maps, several features associated with dust condensations other than FIR 5 are detected.

To characterize the magnetic field structure of the outflow and core region within a prototypical high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the Submillimeter Array -- as well as 880mum submm continuum emission from the high-mass outflow/disk system IRAS18089-1732. Both emission features with polarization degrees at a few percent level indicate that the magnetic field structure is largely aligned with the outflow/jet orientation from the small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic field strength increase from the lower-density envelope to the higher-density core.

We present observational results of the thermal dust continuum emission and its linear polarization in one of the nearest massive star-forming sites Orion BN/KL in Orion Molecular Cloud-1. The observations were carried out with the Submillimeter Array. With an angular resolution of 1" (~2 mpc; 480 AU), we have detected and resolved the densest cores near the BN/KL region. At a wavelength of ~870 micron, the polarized dust emission can be used to trace the structure of the magnetic field in this star-forming core. The dust continuum appears to arise from a V-shaped region, with a cavity nearly coincident with the center of the explosive outflows observed on larger scales. The position angles (P.A.s) of the observed polarization vary significantly by a total of about 90 degree but smoothly, i.e., curl-like, across the dust ridges. Such a polarization pattern can be explained with dust grains being magnetically aligned instead of mechanically with outflows, since the latter mechanism would cause the P.A.s to be parallel to the direction of the outflow, i.e., radial-like. The magnetic field projected in the plane of sky is therefore derived by rotating the P.A.s of the polarization by 90 degree. We find an azimuthally symmetric structure in the overall magnetic field morphology, with the field directions pointing toward 2.5" west to the center of the explosive outflows. We also find a preferred symmetry plane at a P.A. of 36 degree, which is perpendicular to the mean magnetic field direction (120 degree) of the 0.5 pc dust ridge. Two possible interpretations of the origin of the observed magnetic field structure are discussed.

We present polarization observations of the dust continuum emission from IRAS 16293 which is a Class 0 protostar and is known to have at least two components, source A and B. These measurements were conducted by the Submillimeter Array (SMA) at a frequency of ~341 GHz and with high angular resolution (~2-3 arcseconds). We find that the large scale global direction of the field, which is perpendicular to the observed polarization, appears to be along the dust ridge where the emission peaks. On smaller scales we find that the field structure is significantly different for the two components of the binary. The first component, source A, shows a magnetic field structure which is "hourglass" shaped as predicted from theoretical models of low mass star formation in the presence of strong magnetic fields. However, the other component, source B, shows a relatively ordered magnetic field with no evidence of any deformation. There is an observed decrease in polarization towards the center. Our calculations show that in IRAS 16293 the magnetic energy is stronger than the turbulent energy but is approximately similar to the centrifugal energy. Our results provide additional evidence to show that the two protostars appear to be in different stages during their evolution.

We present arcsecond-resolution Submillimeter Array (SMA) polarimetric observations of the 880 um continuum emission from the protoplanetary disks around two nearby stars, HD 163296 and TW Hydrae. Although previous observations and theoretical work have suggested that a 2-3% polarization fraction should be common for the millimeter continuum emission from such disks, we detect no polarized continuum emission above a 3-sigma upper limit of 7 mJy in each arcsecond-scale beam, or <1% in integrated continuum emission. We compare the SMA upper limits with the predictions from the exploratory Cho & Lazarian (2007) model of polarized emission from T Tauri disks threaded by toroidal magnetic fields, and rule out their fiducial model at the ~10-sigma level. We explore some potential causes for this discrepancy, focusing on model parameters that describe the shape, magnetic field alignment, and size distribution of grains in the disk. We also investigate related effects like the magnetic field strength and geometry, scattering off of large grains, and the efficiency of grain alignment, including recent advances in grain alignment theory, which are not considered in the fiducial model. We discuss the impact each parameter would have on the data and determine that the suppression of polarized emission plausibly arises from rounding of large grains, reduced efficiency of grain alignment with the magnetic field, and/or some degree of magnetic field tangling (perhaps due to turbulence). A poloidal magnetic field geometry could also reduce the polarization signal, particularly for a face-on viewing geometry like the TW Hya disk. The data provided here offer the most stringent limits to date on the polarized millimeter-wavelength emission from disks around young stars.

We report our observational results of 870 $\mu$m continuum emission and its linear polarization in the massive star formation site W51 e2/e8. Inferred from the linear polarization maps, the magnetic field in the plane of sky (B$_{\bot}$) is traced with an angular resolution of 0$\farcs$7 with the Submillimeter Array (SMA). Whereas previous BIMA observations with an angular resolution of 3$\arcsec$ (0.1 pc) showed a uniform B field, our revealed B$_{\bot}$ morphology is hourglass-like in the collapsing core near the Ultracompact H II region e2 and also possibly in e8. The decrease in polarization near the continuum peak seen at lower angular resolution is apparently due to the more complex structures at smaller scales. In e2, the pinched direction of the hourglass-like B field morphology is parallel to the plane of the ionized accretion flow traced by H53$\alpha$, suggesting that the massive stars are formed via processes similar to the low mass stars, i.e. accretion through a disk, except that the mass involved is much larger. Furthermore, our finding that the resolved collapsing cores in e2 and e8 lie within one subcritical 0.5 pc envelope supports the scenario of \textit{magnetic fragmentation} via ambipolar diffusion. We therefore suggest that magnetic fields control the dynamical evolution of the envelope and cores in W51 e2 and e8.

We report high angular resolution (3") Submillimeter Array (SMA) observations of the molecular cloud associated with the Ultra-Compact HII region G5.89-0.39. Imaged dust continuum emission at 870 micron reveals significant linear polarization. The position angles (PAs) of the polarization vary enormously but smoothly in a region of 2x10^4 AU. Based on the distribution of the PAs and the associated structures, the polarized emission can be separated roughly into two components. The component "x" is associated with a well defined dust ridge at 870 micron, and is likely tracing a compressed B field. The component "o" is located at the periphery of the dust ridge and is probably from the original B field associated with a pre-existing extended structure. The global B field morphology in G5.89, as inferred from the PAs, is clearly disturbed by the expansion of the HII region and the molecular outflows. Using the Chandrasekhar-Fermi method, we estimate from the smoothness of the field structures that the B field strength in the plane of sky can be no more than 2-3 mG. We then compare the energy densities in the radiation, the B field, and the mechanical motions as deduced from the C^17O 3-2 line emission. We conclude that the B field structures are already overwhelmed and dominated by the radiation, outflows, and turbulence from the newly formed massive stars.

We present 1.35 mm SMA observations around the low-mass Class 0 source IRAS 19156+1906, at the the center of the L723 dark cloud. We detected emission from dust as well as emission from H2CO, DCN and CN, which arise from two cores, SMA 1 and SMA 2, separated by 2.9" (880 AU). SMA 2 is associated with VLA 2. SiO 5-4 emission is detected, possibly tracing a region of interaction between the dense envelope and the outflow. We modeled the dust and the H2CO emission from the two cores: they have similar physical properties but SMA 2 has a larger p-H2CO abundance than SMA 1. The p-H2CO abundances found are compatible with the value of the outer part of the circumstellar envelopes associated with Class 0 sources. SMA 2 is likely more evolved than SMA 1. The kinematics of the two sources show marginal evidence of infall and rotation motions. The mass detected by the SMA observation, which trace scales of ~1000 AU, is only a small fraction of the mass contained in the large scale molecular envelope, which suggests that L723 is still in a very early phase of star formation. Despite the apparent quiescent nature of the L723, fragmentation is occurring at the center of the cloud at different scales. Thus, at 1000 AU the cloud has fragmented in two cores, SMA 1 and SMA 2. At the same time, at least one of these cores, SMA 2, has undergone additional fragmentation at scales of 150 AU, forming a multiple stellar system.

We describe the Submillimeter Array (SMA) Polarimeter, a polarization converter and feed multiplexer installed on the SMA. The polarimeter uses narrow-band quarter-wave plates to generate circular polarization sensitivity from the linearly-polarized SMA feeds. The wave plates are mounted in rotation stages under computer control so that the polarization handedness of each antenna is rapidly selectable. Positioning of the wave plates is found to be highly repeatable, better than 0.2 degrees. Although only a single polarization is detected at any time, all four cross correlations of left- and right-circular polarization are efficiently sampled on each baseline through coordinated switching of the antenna polarizations in Walsh function patterns. The initial set of anti-reflection coated quartz and sapphire wave plates allows polarimetry near 345 GHz; these plates have been have been used in observations between 325 and 350 GHz. The frequency-dependent cross-polarization of each antenna, largely due to the variation with frequency of the retardation phase of the single-element wave plates, can be measured precisely through observations of bright point sources. Such measurements indicate that the cross-polarization of each antenna is a few percent or smaller and stable, consistent with the expected frequency dependence and very small alignment errors. The polarimeter is now available for general use as a facility instrument of the SMA.

Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray observations, but clear counterparts at longer wavelengths have not been seen. We present a flare observed over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, the Caltech Submillimeter Observatory, and the Submillimeter Array. All telescopes observed strong flare events, but the submillimeter peak is found to occur nearly 100 minutes after the X-ray peak. Submillimeter polarization data show linear polarization in the excess flare emission, increasing from 9% to 17% as the flare passes through its peak, consistent with a transition from optically thick to thin synchrotron emission. The temporal and spectral behavior of the flare require that the energetic electrons responsible for the emission cool faster than expected from their radiative output. This is consistent with adiabatic cooling in an expanding emission region, with X-rays produced through self-Compton scattering, although not consistent with the simplest model of such expansion. We also present a submillimeter flare that followed a bright IR flare on 2005 July 31. Compared to 2006, this event had a larger peak IR flux and similar submillimeter flux, but it lacked measurable X-ray emission. It also showed a shorter delay between the IR and submillimeter peaks. Based on these events we propose a synchrotron and self-Compton model to relate the submillimeter lag and the variable IR/X-ray luminosity ratio.

We used the Submillimeter Array to map the angular distribution of the H30$\alpha$ recombination line (231.9 GHz) in the circumstellar region of the peculiar star MWC349A. The resolution was $1\farcs2$, but because of high signal-to-noise ratio we measured the positions of all maser components to accuracies better than $0\farcs01$, at a velocity resolution of $1 kms$. The two strongest maser components (called high velocity components) at velocities near -14 and $32 kms$ are separated by $0\farcs048 \pm 0\farcs001$ (60 AU) along a position angle of $102 \pm 1\arcdeg$. The distribution of maser emission at velocities between and beyond these two strongest components were also provided. The continuum emission lies at the center of the maser distribution to within 10 mas. The masers appear to trace a nearly edge-on rotating disk structure, reminiscent of the water masers in Keplerian rotation in the nuclear accretion disk of the galaxy NGC4258. However, the maser components in MWC349A do not follow a simple Keplerian kinematic prescription with $v \sim r^{-1/2}$, but have a larger power law index. We explore the possibility that the high velocity masers trace spiral density or shock waves. We also emphasize caution in the interpretation of relative centroid maser positions where the maser is not clearly resolved in position or velocity, and we present simulations that illustrate the range of applicability of the centroiding method.

The millimeter/submillimeter wavelength polarization of Sgr A* is known to be variable in both magnitude and position angle on time scales down to a few hours. The unstable polarization has prevented measurements made at different frequencies and different epochs from yielding convincing measurements of Faraday rotation in this source. Here we present observations made with the Submillimeter Array polarimeter at 227 and 343 GHz with sufficient sensitivity to determine the rotation measure at each band without comparing position angles measured at separate epochs. We find the 10-epoch mean rotation measure to be (-5.6+/-0.7)x10^5 rad/m^2; the measurements are consistent with a constant value. We conservatively assign a 3sigma upper limit of 2x10^5 rad/m^2 to rotation measure changes, which limits accretion rate fluctuations to 25%. This rotation measure detection limits the accretion rate to less than 2x10^-7 M_sun/yr if the magnetic field is near equipartition, ordered, and largely radial, while a lower limit of 2x10^-9 M_sun/yr holds even for a sub-equipartition, disordered, or toroidal field. The mean intrinsic position angle is 167+/-7 degrees and we detect variations of 31(+18/-9) degrees. These variations must originate in the submillimeter photosphere, rather than arising from rotation measure changes.

We report high-angular-resolution measurements of polarized dust emission toward the low-mass protostellar system NGC 1333 IRAS 4A. We show that in this system the observed magnetic field morphology is in agreement with the standard theoretical models of the formation of Sun-like stars in magnetized molecular clouds at scales of a few hundred astronomical units; gravity has overcome magnetic support, and the magnetic field traces a clear hourglass shape. The magnetic field is substantially more important than turbulence in the evolution of the system, and the initial misalignment of the magnetic and spin axes may have been important in the formation of the binary system.

We report on the submillimeter properties of Sagittarius A* derived from observations with the Submillimeter Array and its polarimeter. We find that the spectrum of Sgr A* between 230 and 690 GHz is slightly decreasing when measured simultaneously, indicating a transition to optically thin emission around 300-400 GHz. We also present very sensitive and well calibrated measurements of the polarization of Sgr A* at 230 and 345 GHz. With these data we are able to show for the first time that the polarization of Sgr A* varies on hour timescales, as has been observed for the total intensity. On one night we find variability that may arise from a polarized "blob" orbiting the black hole. Finally, we use the ensemble of observations to determine the rotation measure. This represents the first statistically significant rotation measure determination and the only one made without resorting to comparing position angles measured at separate epochs. We find a rotation measure of (-5.6+/-0.7)x10^5 rad/m^2, with no evidence for variability on inter-day timescales at the level of the measurement error. The stability constrains interday fluctuations in the accretion rate to 8%. The mean intrinsic polarization position angle is 167+/-7 degrees and we detect variations of 31+18/-9 degrees. This separation of intrinsic polarization changes and possible rotation measure fluctuations is now possible because of the frequency coverage and sensitivity of our data. The observable rotation measure restricts the accretion rate to the range 2x10^{-7} Msun/yr to 2x10^{-9} Msun/yr, if the magnetic field is near equipartition and ordered.

We report new simultaneous near-infrared/sub-millimeter/X-ray observations of the SgrA* counterpart associated with the massive 3-4x10**6 solar mass black hole at the Galactic Center. The main aim is to investigate the physical processes responsible for the variable emission from SgrA*. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope and the ACIS-I instrument aboard the Chandra X-ray Observatory as well as the Submillimeter Array SMA on Mauna Kea, Hawaii, and the Very Large Array in New Mexico. We detected one moderately bright flare event in the X-ray domain and 5 events at infrared wavelengths.

Using the Submillimeter Array, we have made the first high angular resolution measurements of the linear polarization of Sagittarius A* at submillimeter wavelengths, and the first detection of intra-day variability in its linear polarization. We detected linear polarization at 340 GHz (880um) at several epochs. At the typical resolution of 1.4"x2.2", the expected contamination from the surrounding (partially polarized) dust emission is negligible. We found that both the polarization fraction and position angle are variable, with the polarization fraction dropping from 8.5% to 2.3% over three days. This is the first significant measurement of variability in the linear polarization fraction in this source. We also found variability in the polarization and total intensity within single nights, although the relationship between the two is not clear from these data. The simultaneous 332 and 342 GHz position angles are the same, setting a one-sigma rotation measure (RM) upper limit of 7x10^5 rad/m^2. From position angle variations and comparison of "quiescent" position angles observed here and at 230 GHz we infer that the RM is a few times 10^5 rad/m^2, a factor of a few below our direct detection limit. A generalized model of the RM produced in the accretion flow suggests that the accretion rate at small radii must be low, below 10^{-6} to 10^{-7} M_{sun}/year depending on the radial density and temperature profiles, but in all cases below the gas capture rate inferred from X-ray observations.

We present the first 865 mu continuum image with sub-arcsecond resolution obtained with the Submillimeter Array. These data resolve the Orion-KL region into the hot core, the nearby radio source I, the sub-mm counterpart to the infrared source n (radio source L), and new sub-mm continuum sources. The radio to submillimeter emission from source I may be modeled as either the result of proton-electron free-free emission that is optically thick to ~100 GHz plus dust emission that accounts for the majority of the submillimeter flux, or H- free-free emission that gives rise to a power-law spectrum with power-law index of ~1.6. The latter model would indicate similar physical conditions as found in the inner circumstellar environment of Mira variable stars. Future sub-arcsecond observations at shorter sub-mm wavelengths should easily discriminate between these two possibilities. The sub-mm continuum emission toward source n can be interpreted in the framework of emission from an accretion disk.

In an attempt to probe the magnetic field morphology near the massive young star Orion-IRc2, we mapped the linear polarization of its J=2-1 SiO masers, in both the v=0 and v=1 vibrational levels, with 0.5 arcsec resolution. The intense v=1 masers are confined to a narrow zone 40 AU from the star. Their polarization position angles vary significantly on time scales of years. For the v=1 masers the stimulated emission rate R is likely to exceed the Zeeman splitting g\Omega due to any plausible magnetic field; in this case the maser polarization need not correlate with the field direction. The much weaker v=0 masers in the ground vibrational level lie 100-700 AU from IRc2, in what appears to be a flared disk. Their fractional polarizations are as high as 50%. The polarization position angles vary little across the line profile or the emission region, and appear to be stable in time. The position angle, P.A. = 80 degrees, we measure for the J=2-1 masers differs by 70 degrees from that measured for the J=1-0 SiO transition, possibly because of Faraday rotation in the foreground, Orion A, HII region. A rotation measure RM = 3.3 \times 10^4 rad m$^{-2}$ is required to bring the J=2-1 and J=1-0 position angles into concordance. The intrinsic polarization position angle for both transitions is then 57 degrees, parallel to the plane of the putative disk. Probably the magnetic field threads the disk poloidally. There is little evidence for a pinched or twisted field near the star.

We present the first interferometric polarization maps of the NGC2024 FIR5 molecular core obtained with the BIMA array at approximately 2 arcsec resolution. We measure an average position angle of -60+-6 degrees in the main core of FIR5 and 54+-9 degrees in the eastern wing of FIR5. The morphology of the polarization angles in the main core of FIR5 suggests that the field lines are parabolic with a symmetry axis approximately parallel to the major axis of the putative disk in FIR5, which is consistent with the theoretical scenario that the gravitational collapse pulled the field lines into an hour-glass shape. The polarization percentage decreases toward regions with high intensity and close to the center of the core, suggesting that the dust alignment efficiency may decrease at high density. The plane-of-sky field strength can be estimated with the modified Chandrasekhar-Fermi formula, and the small dispersion of the polarization angles in FIR5 suggests that the magnetic field is strong ($\gtrsim$ 2mG) and perhaps dominates the turbulent motions in the core.

We present the first interferometric polarization map of the W51 e1/e2 molecular cores obtained with the BIMA array at 1.3 mm wavelength with approximately 3 arcsecond resolution. The polarization angle varies smoothly across the double cores with an average position angle of 23+-5 degrees for W51 e1 and 15+-7 degrees for W51 e2. The inferred magnetic field direction is parallel to the minor axis of the double cores, which is consistent with the theoretical picture that clouds collapse along the field lines. However, the magnetic field may not determine the axis of angular momentum of these two cores as the field directions of the two cores significantly differ with the previously measured directions of rotational axes. The polarization percentage decreases toward regions with high intensity, suggesting that the dust alignment efficiency decreases toward high density regions. The field directions are highly ordered, and the small dispersion of the polarization angles implies that magnetic fields are strong ($\gtrsim$ 1 mG) and perhaps dominate turbulence in W51 e1/e2.

We present 1.3 and 3.3 mm polarization maps of Orion-KL obtained with the BIMA array at approximately 4 arcsec resolution. Thermal emission from magnetically aligned dust grains produces the polarization. Along the Orion ``ridge'' the polarization position angle varies smoothly from about 10 degrees to 40 degrees, in agreement with previous lower resolution maps. In a small region south of the Orion ``hot core,'' however, the position angle changes by 90 degrees. This abrupt change in polarization direction is not necessarily the signpost of a twisted magnetic field. Rather, in this localized region processes other than the usual Davis-Greenstein mechanism might align the dust grains with their long axes parallel with the field, orthogonal to their normal orientation.