Search

Of all the factors that influence star formation, magnetic fields are perhaps the least well understood. The goal of this paper is to characterize the 3D magnetic field properties of nearby molecular clouds through various methods of statistically analysing maps of polarized dust emission. Our study focuses on nine clouds, with data taken from the Planck Sky Survey as well as data from the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry observations of Vela C. We compare the distributions of polarization fraction (p), dispersion in polarization angles ($\mathcal {S}$), and hydrogen column density (NH) for each of our targeted clouds. To broaden the scope of our analysis, we compare the distributions of our clouds’ polarization observables with measurements from synthetic polarization maps generated from numerical simulations. We also use the distribution of polarization fraction measurements to estimate the inclination angle of each cloud’s cloud-scale magnetic field. We obtain a range of inclination angles associated with our clouds, varying from 16○ to 69○. We establish inverse correlations between p and both $\mathcal {S}$ and NH in almost every cloud, but we are unable to establish a statistically robust $\mathcal {S}$ versus NH trend. By comparing the results of these different statistical analysis techniques, we are able to propose a more comprehensive view of each cloud’s 3D magnetic field properties. These detailed cloud analyses will be useful in the continued studies of cloud-scale magnetic fields and the ways in which they affect star formation within these molecular clouds.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud., 18 pages, 11 figures, accepted for publication in MNRAS

The complete three-dimensional structure of the magnetic field within molecular clouds has eluded determination despite its high value in determining controlling factors in the star formation process, as it cannot be directly probed observationally. Considering that inclination of the magnetic field relative to the plane of sky is one of the major sources of depolarization of thermal emission from dust in molecular clouds, we propose here a new method to estimate the inclination angle of the cloud-scale magnetic field based on the statistical properties of the observed polarization fraction. We test this method using a series of Monte Carlo experiments, and find that the method works well provided that deviations of magnetic field direction from the averaged values are small. When applied to synthetic observations of numerical simulations of star-forming clouds, our method gives fairly accurate measurements of the mean magnetic field inclination angle (within $10^\circ-25^\circ$), which can further be improved if we restrict our technique to regions of low dispersion in polarization angles ${\cal S}$. We tested our method on the BLASTPol polarimetric observations of the Vela C molecular cloud complex, which suggests that the magnetic field of Vela C has a high inclination angle ($\sim 60^\circ$), consistent with previous analyses., Comment: 15 pages, 15 figures, version accepted for publication in MNRAS

We have obtained sensitive dust continuum polarization observations at 850 μm in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope as part of the B-fields in STar-forming Region Observations (BISTRO) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution (∼2000 au or ∼0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis–Chandrasekhar–Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38 ± 14, 44 ± 16, and 12 ± 5 μG, respectively. These cores show distinct mean B-field orientations. The B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. The B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field and not well correlated with other axes. In contrast, Miz-8b exhibits a disordered B-field that shows no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux.

TolTEC is a three-band imaging polarimeter for the Large Millimeter Telescope. Simultaneously observing with passbands at 1.1mm, 1.4mm and 2.0mm, TolTEC has diffraction-limited beams with FWHM of 5, 7, and 11 arcsec, respectively. Over the coming decade, TolTEC will perform a combination of PI-led and Open-access Legacy Survey projects. Herein we provide an overview of the instrument and give the first quantitative measures of its performance in the lab prior to shipping to the telescope in 2021.

We report the first high spatial resolution measurement of magnetic fields surrounding LkH$\alpha$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $\mu$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrm{\mu}$G. The mass-to-magnetic-flux ratio in critical value units, $\lambda\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$\alpha$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $\delta B/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $\alpha=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$\alpha$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend., Comment: 25 pages, 19 figures, Accepted for publication in The Astrophysical Journal

In this work we describe upgrades to the Spider balloon-borne telescope in preparation for its second flight, currently planned for December 2021. The Spider instrument is optimized to search for a primordial B-mode polarization signature in the cosmic microwave background at degree angular scales. During its first flight in 2015, Spider mapped ~10% of the sky at 95 and 150 GHz. The payload for the second Antarctic flight will incorporate three new 280 GHz receivers alongside three refurbished 95- and 150 GHz receivers from Spider's first flight. In this work we discuss the design and characterization of these new receivers, which employ over 1500 feedhorn-coupled transition-edge sensors. We describe pre-flight laboratory measurements of detector properties, and the optical performance of completed receivers. These receivers will map a wide area of the sky at 280 GHz, providing new information on polarized Galactic dust emission that will help to separate it from the cosmological signal., 13 pages, 8 figures; as published in the conference proceedings for SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X (2020)

[Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$\mu$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($\theta_{B}$). We derive the power spectrum of the intensity and $\theta_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $\theta_{B}$ varies coherently along the filaments. The observed power spectrum of $\theta_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $\theta_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $\theta_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars., Comment: Accepted for publication in Astronomy & Astrophysics

The BLASTPol observations of Vela C have provided the most detailed characterization of the polarization fraction $p$ and dispersion in polarization angles $S$ for a molecular cloud. We compare the observed distributions of $p$ and $S$ with those obtained in synthetic observations of simulations of molecular clouds, assuming homogeneous grain alignment. We find that the orientation of the mean magnetic field relative to the observer has a significant effect on the $p$ and $S$ distributions. These distributions for Vela C are most consistent with synthetic observations where the mean magnetic field is close to the line-of-sight. Our results point to apparent magnetic disorder in the Vela C molecular cloud, although it can be due to either an inclination effect (i.e., observing close to the mean field direction) or significant field tangling from strong turbulence/low magnetization. The joint correlations of $p$ with column density and of $S$ with column density for the synthetic observations generally agree poorly with the Vela C joint correlations, suggesting that understanding these correlations require a more sophisticated treatment of grain alignment physics., 24 pages, 10 figures, slightly-revised version accepted by MNRAS

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud., Comment: 27 pages, 15 figures, accepted for publication in ApJ

The mm-wavelength sky reveals the initial phase of structure formation, at all spatial scales, over the entire observable history of the Universe. Over the past 20 years, advances in mm-wavelength detectors and camera systems have allowed the field to take enormous strides forward – particularly in the study of the Cosmic Microwave Background – but limitations in mapping speeds, sensitivity and resolution have plagued studies of astrophysical phenomena. In fact, limitations due to inherent biases in the ground-based mm-wavelength surveys conducted over the last 2 decades continue to motivate the need for deeper and wider-area maps made with increased angular resolution. TolTEC is a new camera that will fill the focal plane of the 50m diameter Large Millimeter Telescope (LMT) and provide simultaneous, polarization-sensitive imaging at 2.0, 1.4, and 1.1mm wavelengths. The instrument, now under construction, is a cryogenically cooled receiver housing three separate kilo-pixel arrays of Kinetic Inductance Detectors (KIDs) that are coupled to the telescope through a series of silicon lenses and dichroic splitters. TolTEC will be installed and commissioned on the LMT in early 2019 where it will become both a facility instrument and also perform a series of 100 hour “Legacy Surveys” whose data will be publicly available. The initial four surveys in this series: the Clouds to Cores Legacy Survey, the Fields in Filaments Legacy Survey, the Ultra-Deep Legacy Survey and the Large Scale Structure Survey are currently being defined in public working groups of astronomers coordinated by TolTEC Science Team members. Data collection for these surveys will begin in late 2019 with data releases planned for late 2020 and 2021. Herein we describe the instrument concept, provide performance data for key subsystems, and provide an overview of the science, schedule and plans for the initial four Legacy Survey concepts.

BFORE is a high-altitude ultra-long-duration balloon mission to map the cosmic microwave background (CMB). During a 28-day mid-latitude flight launched from Wanaka, New Zealand, the instrument will map half the sky to improve measurements of the optical depth to reionization tau. This will break parameter degeneracies needed to detect neutrino mass. BFORE will also hunt for the gravitational wave B-mode signal, and map Galactic dust foregrounds. The mission will be the first near-space use of TES/mSQUID multichroic detectors (150/217 GHz and 280/353 GHz bands) with low-power readout electronics., 11 pages, 2 figures, Submitted to SPIE Astronomical Telescopes and Instrumentation

We present the projected Rayleigh statistic (PRS), a modification of the classic Rayleigh statistic, as a test for non-uniform relative orientation between two pseudo-vector fields. In the application here this gives an effective way of investigating whether polarization pseudo-vectors (spin-2 quantities) are preferentially parallel or perpendicular to filaments in the interstellar medium. For example, there are other potential applications in astrophysics, e.g., when comparing small-scale orientations with larger-scale shear patterns. We compare the efficiency of the PRS against histogram binning methods that have previously been used for characterising the relative orientations of gas column density structures with the magnetic field projected on the plane of the sky. We examine data for the Vela C molecular cloud, where the column density is inferred from Herschel submillimetre observations, and the magnetic field from observations by the Balloon-borne Large-Aperture Submillimetre Telescope in the 250-, 350-, and 500-{\mu}m wavelength bands. We find that the PRS has greater statistical power than approaches that bin the relative orientation angles, as it makes more efficient use of the information contained in the data. In particular, the use of the PRS to test for preferential alignment results in a higher statistical significance, in each of the four Vela C regions, with the greatest increase being by a factor 1.3 in the South-Nest region in the 250-{\mu}m band., Comment: 11 pages, 7 figures, Submitted to MNRAS

We statistically evaluate the relative orientation between gas column density structures, inferred from Herschel submillimetre observations, and the magnetic field projected on the plane of sky, inferred from polarized thermal emission of Galactic dust observed by BLASTPol at 250, 350, and 500 micron, towards the Vela C molecular complex. First, we find very good agreement between the polarization orientations in the three wavelength-bands, suggesting that, at the considered common angular resolution of 3.0 arcminutes that corresponds to a physical scale of approximately 0.61 pc, the inferred magnetic field orientation is not significantly affected by temperature or dust grain alignment effects. Second, we find that the relative orientation between gas column density structures and the magnetic field changes progressively with increasing gas column density, from mostly parallel or having no preferred orientation at low column densities to mostly perpendicular at the highest column densities. This observation is in agreement with previous studies by the Planck collaboration towards more nearby molecular clouds. Finally, we find a correspondence between the trends in relative orientation and the shape of the column density probability distribution functions. In the sub-regions of Vela C dominated by one clear filamentary structure, or "ridges", we find a sharp transition from preferentially parallel or having no preferred relative orientation at low column densities to preferentially perpendicular at highest column densities. In the sub-regions of Vela C dominated by several filamentary structures with multiple orientations, or "nests", such a transition is also present, but it is clearly less sharp than in the ridge-like sub-regions. Both of these results suggest that the magnetic field is dynamically important for the formation of density structures in this region., Comment: 16 pages, 17 figures. Submitted to A&A

Thin, magnetically aligned striations of relatively moderate contrast with the background are commonly observed in both atomic and molecular clouds. They are also prominent in MHD simulations with turbulent converging shocks. The simulated striations develop within a dense, stagnated sheet in the mid plane of the post-shock region where magnetically-induced converging flows collide. We show analytically that the secondary flows are an inevitable consequence of the jump conditions of oblique MHD shocks. They produce the stagnated, sheet-like sub-layer through a secondary shock when, roughly speaking, the Alfv\'enic speed in the primary converging flows is supersonic, a condition that is relatively easy to satisfy in interstellar clouds. The dense sub-layer is naturally threaded by a strong magnetic field that lies close to the plane of the sub-layer. The substantial magnetic field makes the sheet highly anisotropic, which is the key to the striation formation. Specifically, perturbations of the primary inflow that vary spatially perpendicular to the magnetic field can easily roll up the sheet around the field lines without bending them, creating corrugations that appear as magnetically-aligned striations in column density maps. On the other hand, perturbations that vary spatially along the field lines curve the sub-layer and alter its orientation relative to the magnetic field locally, seeding special locations that become slanted overdense filaments and prestellar cores through enhanced mass accumulation along field lines. In our scenario, the dense sub-layer unique to magnetized oblique shocks is the birthplace for both magnetically-aligned diffuse striations and massive star-forming structures., Comment: 13 pages, 13 figures, submitted to ApJ

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) is a suborbital mapping experiment designed to study the role magnetic fields play in star formation. BLASTPol has had two science flights from McMurdo Station, Antarctica in 2010 and 2012. These flights have produced thousands of polarization vectors at 250, 350 and 500 microns in several molecular cloud targets. We present the design, specifications, and progress towards the next-generation BLASTPol experiment (BLAST-TNG). BLAST-TNG will fly a 40% larger diameter primary mirror, with almost 8 times the number of polarization-sensitive detectors resulting in a factor of 16 increase in mapping speed. With a spatial resolution of 22 arcseconds and four times the field of view of BLASTPol, BLAST-TNG will bridge the angular scales between Planck's low resolution all-sky maps and ALMA's ultra-high resolution narrow fields. The new receiver has a larger cryogenics volume, allowing for a 28 day hold time. BLAST-TNG employs three arrays of Microwave Kinetic Inductance Detectors (MKIDs) with 30% fractional bandwidth at 250, 350 and 500 microns. In this paper, we will present the new BLAST-TNG instrument and science objectives., 12 pages, 6 figures. Presented at SPIE Ground-based and Airborne Telescopes V, June 23, 2014. To be published in Proceedings of SPIE Volume 9145

We present first results of neutral carbon ([CI], 3P1 - 3P0 at 492 GHz) and carbon monoxide (13CO, J = 1 - 0) mapping in the Vela Molecular Ridge cloud C (VMR-C) and G333 giant molecular cloud complexes with the NANTEN2 and Mopra telescopes. For the four regions mapped in this work, we find that [CI] has very similar spectral emission profiles to 13CO, with comparable line widths. We find that [CI] has opacity of 0.1 - 1.3 across the mapped region while the [CI]/13CO peak brightness temperature ratio is between 0.2 to 0.8. The [CI] column density is an order of magnitude lower than that of 13CO. The H2 column density derived from [CI] is comparable to values obtained from 12CO. Our maps show CI is preferentially detected in gas with low temperatures (below 20 K), which possibly explains the comparable H2 column density calculated from both tracers (both CI and 12CO underestimate column density), as a significant amount of the CI in the warmer gas is likely in the higher energy state transition ([CI], 3P2 - 3P1 at 810 GHz), and thus it is likely that observations of both the above [CI] transitions are needed in order to recover the total H2 column density., Comment: accepted for publication in ApJ Letters

We present the second generation BLASTbus electronics. The primary purposes of this system are detector readout, attitude control, and cryogenic housekeeping, for balloon-borne telescopes. Readout of neutron transmutation doped germanium (NTD-Ge) bolometers requires low noise and parallel acquisition of hundreds of analog signals. Controlling a telescope's attitude requires the capability to interface to a wide variety of sensors and motors, and to use them together in a fast, closed loop. To achieve these different goals, the BLASTbus system employs a flexible motherboard-daughterboard architecture. The programmable motherboard features a digital signal processor (DSP) and field-programmable gate array (FPGA), as well as slots for three daughterboards. The daughterboards provide the interface to the outside world, with versions for analog to digital conversion, and optoisolated digital input/output. With the versatility afforded by this design, the BLASTbus also finds uses in cryogenic, thermometry, and power systems. For accurate timing control to tie everything together, the system operates in a fully synchronous manner. BLASTbus electronics have been successfully deployed to the South Pole, and flown on stratospheric balloons., Presented at SPIE Ground-based and Airborne Telescopes V, June 23, 2014. To be published in Proceedings of SPIE Volume 9145

We describe the cryogenic system for SPIDER, a balloon-borne microwave polarimeter that will map 8% of the sky with degree-scale angular resolution. The system consists of a 1284 L liquid helium cryostat and a 16 L capillary-filled superfluid helium tank, which provide base operating temperatures of 4 K and 1.5 K, respectively. Closed-cycle helium-3 adsorption refrigerators supply sub-Kelvin cooling power to multiple focal planes, which are housed in monochromatic telescope inserts. The main helium tank is suspended inside the vacuum vessel with thermally insulating fiberglass flexures, and shielded from thermal radiation by a combination of two vapor cooled shields and multi-layer insulation. This system allows for an extremely low instrumental background and a hold time in excess of 25 days. The total mass of the cryogenic system, including cryogens, is approximately 1000 kg. This enables conventional long duration balloon flights. We will discuss the design, thermal analysis, and qualification of the cryogenic system., 11 pages, 4 figures; as published in the conference proceedings for SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy V (2010)

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST., 60 page manuscript, inc. 18 figures. Accepted by PASP

We present results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 {\mu}m. In this initial paper, we show our 500 {\mu}m data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations., Comment: 24 pages, 20 figures, and 2 tables. Accepted for publication in ApJ. Updated from originally submitted version