Far-infrared (FIR) dust polarimetry enables the study of interstellar magnetic fields via tracing of the polarized emission from dust grains that are partially aligned with the direction of the field. The advent of high quality polarimetric data has permitted the use of statistical methods to extract both the direction and magnitude of the magnetic field. This work presents an effort to produce maps (integrated along the line of sight) of the magnetic field strength in both the line-of-sight (LOS) direction and in the plane-of-sky (POS) in the Orion Molecular Cloud (OMC-1) using all four HAWC+ polarimetry bands (53, 89, 154 and 214 $\micron$). Using the statistics of the polarization data (angles and fractions) along with previously-reported Zeeman measurements, information about the three-dimensional magnetic field configuration was inferred over the OMC-1 region. POS magnetic field strengths of up to 2 mG were determined in OMC-1. The strongest fields were observed near the BN/KL object, while the OMC-1 bar shows strengths of up to few hundreds of $\mu$G. The LOS magnetic field has been determined to have comparable strengths to the POS components but it is differently distributed across the cloud. These estimates of the magnetic field components were used to produce more reliable maps of the mass-to-magnetic flux ratio ($M/\Phi$) -- a proxy for probing the conditions for star formation in molecular clouds -- and determine regions of sub- and super-criticality in OMC-1. Such maps can provide invaluable input and comparison to MHD simulations of star formation processes in filamentary structures of molecular clouds.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.