Polarization observations of the Milky Way and many other spiral galaxies have found a close correspondence between the orientation of spiral arms and magnetic field lines on scales of hundreds of parsecs. This paper presents polarization measurements at 214 $\mu$m toward ten filamentary candidate ``bones" in the Milky Way using the High-resolution Airborne Wide-band Camera (HAWC+) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). These data were taken as part of the Filaments Extremely Long and Dark: A Magnetic Polarization Survey (FIELDMAPS) and represent the first study to resolve the magnetic field in spiral arms at parsec scales. We describe the complex yet well-defined polarization structure of all ten candidate bones, and we find a mean difference and standard deviation of $-74^{\circ} \pm 32^{\circ}$ between their filament axis and the plane-of-sky magnetic field, closer to a field perpendicular to their length rather than parallel. By contrast, the 850 $\mu$m polarization data from \textit{Planck} on scales greater than 10 pc show a nearly parallel mean difference of $3^{\circ} \pm 21^{\circ}$. These findings provide further evidence that magnetic fields can change orientation at the scale of dense molecular clouds, even along spiral arms. Finally, we use a power law to fit the dust polarization fraction as a function of total intensity on a cloud-by-cloud basis and find indices between $-0.6$ and $-0.9$, with a mean and standard deviation of $-0.7 \pm 0.1$. The polarization, dust temperature, and column density data presented in this work are publicly available online.

Stars primarily form in galactic spiral arms within dense, filamentary molecular clouds. The largest and most elongated of these molecular clouds are referred to as ``bones," which are massive, velocity-coherent filaments (lengths ~20 to >100 pc, widths ~1-2 pc) that run approximately parallel and in close proximity to the Galactic plane. While these bones have been generally well characterized, the importance and structure of their magnetic fields (B-fields) remain largely unconstrained. Through the SOFIA Legacy program FIELDMAPS, we mapped the B-fields of 10 bones in the Milky Way. We found that their B-fields are varied, with no single preferred alignment along the entire spine of the bones. At higher column densities, the spines of the bones are more likely to align perpendicularly to the B-fields, although this is not ubiquitous, and the alignment shows no strong correlation with the locations of identified young stellar objects. We estimated the B-field strengths across the bones and found them to be ~30-150 $\mu$G at pc scales. Despite the generally low virial parameters, the B-fields are strong compared to the local gravity, suggesting that B-fields play a significant role in resisting global collapse. Moreover, the B-fields may slow and guide gas flow during dissipation. Recent star formation within the bones may be due to high-density pockets at smaller scales, which could have formed before or simultaneously with the bones.

The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim$ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mass protobinary systems, and no definitive higher-order multiples, have been detected. Here we report the discovery of one quintuple, one quadruple, one triple and four binary protostellar systems simultaneously forming in a single high-mass protocluster, G333.23--0.06, using Atacama Large Millimeter/submillimeter Array high-resolution observations. We present a new example of a group of gravitationally bound binary and higher-order multiples during their early formation phases in a protocluster. This provides the clearest direct measurement of the initial configuration of primordial high-order multiple systems, with implications for the in situ multiplicity and its origin. We find that the binary and higher-order multiple systems, and their parent cores, show no obvious sign of disk-like kinematic structure. We conclude that the observed fragmentation into binary and higher-order multiple systems can be explained by core fragmentation, indicating its crucial role in establishing the multiplicity during high-mass star cluster formation.

Magnetic fields influence the structure and evolution of protostellar systems, thus understanding their role is essential for probing the earliest stages of star formation. We present ALMA Band 3 and 6 polarized continuum observations at $\sim$0.5$^{\prime \prime}$ resolution toward the Class 0 protostellar system HH 211. Three dust filaments ($\sim$4000 au in length) are found in the HH 211 protostellar envelope, two of which are aligned with core-scale ($\sim$10,000 au) magnetic fields detected by previous JCMT observations. This result suggests that the formation of the dust filaments may be influenced by magnetic fields. In the inner envelope ($\sim$1000 au), we detect a clear hourglass-shaped magnetic field morphology near the protostar and toroidal fields along the outflow directions. We also estimate the line-of-sight-averaged temperature and column density distributions in the inner envelope and find that the temperature is higher in the east, while the column density is enhanced in the southern and western regions. The southern dense regions of the inner envelope may trace either outflow cavity walls, due to their alignment with the outflow, or possible infalling channels in the midplane, given the close correspondence between the observed magnetic fields and the predicted infall trajectories.

Recent (sub)millimeter polarization observations of protoplanetary disks reveal toroidally aligned, effectively prolate dust grains large enough (at least ~100 $\mu$m) to efficiently scatter millimeter light. The alignment mechanism for these grains remains unclear. We explore the possibility that gas drag aligns grains through gas-dust relative motion when the grain's center of mass is offset from its geometric center, analogous to a badminton birdie's alignment in flight. A simple grain model of two non-identical spheres illustrates how a grain undergoes damped oscillations from flow-induced restoring torques which align its geometric center in the flow direction relative to its center of mass. Assuming specular reflection and subsonic flow, we derive an analytical equation of motion for spheroids where the center of mass can be shifted away from the spheroid's geometric center. We show that a prolate or an oblate grain can be aligned with the long axis parallel to the gas flow when the center of mass is shifted along that axis. Both scenarios can explain the required effectively prolate grains inferred from observations. Application to a simple disk model shows that the alignment timescales are shorter than or comparable to the orbital time. The grain alignment direction in a disk depends on the disk (sub-)structure and grain Stokes number (St) with azimuthal alignment for large St grains in sub-Keplerian smooth gas disks and for small St grains near the gas pressure extrema, such as rings and gaps.

We aim to develop a new method to infer the sub-beam probability density function (PDF) of H2 column densities and the dense gas mass within molecular clouds using spatially unresolved observations of molecular emission lines in the 3 mm band. We model spatially unresolved line integrated intensity measurements as the average of an emission function weighted by the sub-beam column density PDF. The emission function, which expresses the line integrated intensity as a function of the gas column density, is an empirical fit to high resolution (< 0.05 pc) multi-line observations of the Orion B molecular cloud. The column density PDF is assumed to be parametric, composed of a lognormal distribution at moderate column densities and a power law distribution at higher column densities. To estimate the sub-beam column density PDF, the emission model is combined with a Bayesian inversion algorithm (the Beetroots code), which takes account of thermal noise and calibration errors. We validate our method by demonstrating that it recovers the true column density PDF of the Orion B cloud, reproducing the observed emission line integrated intensities. We apply the method to 12CO(J=1-0), 13CO(J=1-0), C18O(J=1-0), HCN(J=1-0), HCO+(J=1-0) and N2H+(J=1-0) observations of a 700 x 700 pc2 field of view (FoV) in the nearby galaxy M51. On average, the model reproduces the observed intensities within 30%. The column density PDFs obtained for the spiral arm region within our test FoV are dominated by a power-law tail at high column densities, with slopes that are consistent with gravitational collapse. Outside the spiral arm, the column density PDFs are predominantly lognormal, consistent with supersonic isothermal turbulence. We calculate the mass associated with the powerlaw tail of the column density PDFs and observe a strong, linear correlation between this mass and the 24$\mu$m surface brightness.

Polycyclic aromatic hydrocarbons (PAHs) are among the most ubiquitous compounds in the universe, accounting for up to ~25% of all interstellar carbon. Since most unsubstituted PAHs do not possess permanent dipole moments, they are invisible to radio astronomy. Constraining their abundances relies on the detection of polar chemical proxies, such as aromatic nitriles. We report the detection of 2- and 4-cyanopyrene, isomers of the recently detected 1-cyanopyrene. We find that these isomers are present in an abundance ratio of ~2:1:2, which mirrors the number of equivalent sites available for CN addition. We conclude that there is evidence that the cyanopyrene isomers formed by direct CN addition to pyrene under kinetic control in hydrogen-rich gas at 10 K and discuss constraints on the H/CN ratio for PAHs in TMC-1.

We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from North-West to South-East with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS$(J=5 \rightarrow 4)$ emission yields an outflow energy of approximately $3.5\times10^{45}$ ergs, aligning with previous interferometric studies. Utilizing the Davis-Chandrasekhar-Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 $\pm 1$ mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary {\HII} region is also larger than the magnetic energy, suggesting that protostellar feedback maybe the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions.

The vertical settling of dust grains in a circumstellar disk, characterized by their scale height, is a pivotal process in the formation of planets. This study offers in-depth analysis and modeling of the radial scale height profile of dust grains in the HL Tau system, leveraging high-resolution polarization observations. We resolve the inner disk's polarization, revealing a significant near-far side asymmetry, with the near side being markedly brighter than the far side in polarized intensity. This asymmetry is attributed to a geometrically thick inner dust disk, suggesting a large aspect ratio of $H/R \ge 0.15$. The first ring at 20 au exhibits an azimuthal contrast, with polarization enhanced along the minor axis, indicating a moderately thick dust ring with $H/R \approx 0.1$. The absence of the near-far side asymmetry at larger scales implies a thin dust layer, with H/R &lt; 0.05. Taken together, these findings depict a disk with a turbulent inner region and a settled outer disk, requiring a variable turbulence model with $\alpha$ increasing from $10^{-5}$ at 100 au to $10^{-2.5}$ at 20 au. This research sheds light on dust settling and turbulence levels within protoplanetary disks, providing valuable insights into the mechanisms of planet formation.

We present SCUBA-2/POL-2 850 $\mu$m polarimetric observations of the circumstellar envelope (CSE) of the carbon-rich asymptotic giant branch (AGB) star IRC+10216. Both FIR and optical polarization data indicate grains aligned with their long axis in the radial direction relative to the central star. The 850 $\mu$m polarization does not show this simple structure. The 850 $\mu$m data are indicative, albeit not conclusive, of a magnetic dipole geometry. Assuming such a simple dipole geometry, the resulting 850 $\mu$m polarization geometry is consistent with both Zeeman observations and small-scale structure in the CSE. While there is significant spectral line polarization contained within the SCUBA-2 850 $\mu$m pass-band for the source, it is unlikely that our broadband polarization results are dominated by line polarization. To explain the required grain alignment, grain mineralogy effects, due to either fossil silicate grains from the earlier oxygen-rich AGB phase of the star, or due to the incorporation of ferromagnetic inclusions in the largest grains, may play a role. We argue that the most likely explanation is due to a new alignment mechanism \citep{arXiv:2009.11304} wherein a charged grain, moving relative to the magnetic field, precesses around the induced electric field and therefore aligns with the magnetic field. This mechanism is particularly attractive as the optical, FIR, and sub-mm wave polarization of the carbon dust can then be explained in a consistent way, differing simply due to the charge state of the grains.

Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $\mu$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can exhibit different correlations between observing wavelength and polarization morphology and fraction. We explore possible origins for these differences in polarization, including differences in dust populations and protostar properties. For RY Tau and MWC 480, which are consistent with scattering at both wavelengths, we present models of the scattering polarization from several dust grain size distributions. These models aim to reproduce two features of the observational results for these disks: (1) both disks have an observable degree of polarization at both wavelengths and (2) the polarization fraction is higher at 3 mm than at 870 $\mu$m in the centers of the disks. For both disks, these features can be reproduced by a power-law distribution of spherical dust grains with a maximum radius of 200 $\mu$m and high optical depth. In MWC 480, we can also reproduce features (1) and (2) with a model containing large grains ($a_{max}$ = 490 $\mu$m ) near the disk midplane and small grains ($a_{max}$ = 140 $\mu$m) above and below the midplane.

We investigate the crescent-shaped dust trap in the transition disk, Oph IRS 48, using well-resolved (sub)millimeter polarimetric observations at ALMA Band 7 (870 $\mu$m). The dust polarization map reveals patterns consistent with dust scattering-induced polarization. There is a relative displacement between the polarized flux and the total flux, which holds the key to understanding the dust scale heights in this system. We model the polarization observations, focusing on the effects of dust scale heights. We find that the interplay between the inclination-induced polarization and the polarization arising from radiation anisotropy in the crescent determines the observed polarization; the anisotropy is controlled by the dust optical depth along the midplane, which is, in turn, determined by the dust scale height in the vertical direction. We find that the dust grains can neither be completely settled nor well mixed with the gas. The completely settled case produces little radial displacement between the total and polarized flux, while the well-mixed case produces an azimuthal pattern in the outer (radial) edge of the crescent that is not observed. Our best model has a gas-to-dust scale height ratio of 2, and can reproduce both the radial displacement and the azimuthal displacement between the total and polarized flux. We infer an effective turbulence $\alpha$ parameter of approximately $0.0001-0.005$. The scattering-induced polarization provides insight into a turbulent vortex with a moderate level of dust settling in the IRS 48 system, which is hard to achieve otherwise.

The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as part of the Magnetic Fields in Massive Star-forming Regions (MagMaR) survey. The B$_\mathrm{POS}$ morphology found in this region is consistent with a canonical ``hourglass'' which suggest a dynamically important field. This region is fragmented into two protostars separated by $\sim1740$ au. Interestingly, by analysing H$^{13}$CO$^{+}$ ($J=3-2$) line emission, we find no velocity gradient over the extend of the continuum which is consistent with a strong field. We model the B$_\mathrm{POS}$, obtaining a marginally supercritical mass-to-flux ratio of 1.43, suggesting an initially strongly magnetized environment. Based on the Davis-Chandrasekhar-Fermi method, the magnetic field strength towards G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows towards the central region of G333 suggests strong magnetic braking, consistent with a highly magnetized environment. Our study shows that despite being a strong regulator, the magnetic energy fails to prevent the process of fragmentation, as revealed by the formation of the two protostars in the central region.

Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $\mu$m) and by grains aligned with something other than magnetic fields at longer wavelengths (e.g., 3 mm). Nevertheless, the polarization signal is expected to depend on the underlying substructure, and observations hitherto have been unable to resolve polarization in multiple rings and gaps. HL Tau, a protoplanetary disk located 147.3 $\pm$ 0.5 pc away, is the brightest Class I or Class II disk at millimeter/submillimeter wavelengths. Here we show deep, high-resolution 870 $\mu$m polarization observations of HL Tau, resolving polarization in both the rings and gaps. We find that the gaps have polarization angles with a significant azimuthal component and a higher polarization fraction than the rings. Our models show that the disk polarization is due to both scattering and emission from aligned effectively prolate grains. The intrinsic polarization of aligned dust grains is likely over 10%, which is much higher than what was expected in low resolution observations (~1%). Asymmetries and dust features are revealed in the polarization observations that are not seen in non-polarimetric observations.

We present a study of the relative orientation between the magnetic field and elongated cloud structures for the $\rho$ Oph A and $\rho$ Oph E regions in L1688 in the Ophiuchus molecular cloud. Combining inferred magnetic field orientation from HAWC+ 154 $\mu$m observations of polarized thermal emission with column density maps created using Herschel submillimeter observations, we find consistent perpendicular relative alignment at scales of $0.02$ pc ($33.6"$ at $d \approx 137$ pc) using the histogram of relative orientations (HRO) technique. This supports the conclusions of previous work using Planck polarimetry and extends the results to higher column densities. Combining this HAWC+ HRO analysis with a new Planck HRO analysis of L1688, the transition from parallel to perpendicular alignment in L1688 is observed to occur at a molecular hydrogen column density of approximately $10^{21.7}$ cm$^{-2}$. This value for the alignment transition column density agrees well with values found for nearby clouds via previous studies using only Planck observations. Using existing turbulent, magnetohydrodynamic simulations of molecular clouds formed by colliding flows as a model for L1688, we conclude that the molecular hydrogen volume density associated with this transition is approximately $\sim10^{4}$ cm$^{-3}$. We discuss the limitations of our analysis, including incomplete sampling of the dense regions in L1688 by HAWC+.

Using the upGREAT instrument on SOFIA, we have imaged [C II] 157.74 and [O I] 63.18 micron line emission from a bright photodissociation region (PDR) associated with an ionized ``bubble'' located in the Nessie Nebula, a filamentary infrared dark cloud. A comparison with ATCA data reveals a classic PDR structure, with a uniform progression from ionized gas, to photodissociated gas, and on to molecular gas from the bubble's interior to its exterior. [O I] line emission from the bubble's PDR reveals self-absorption features. Toward a FIR-bright protostar, both [O I] and [C II] show an absorption feature at a velocity of $-18$ km/s, the same velocity as an unrelated foreground molecular cloud. Since the gas density in typical molecular clouds is well below the [O I] and [C II] critical densities, the excitation temperatures for both lines are low (~20 K). The Meudon models demonstrate that the surface of a molecular cloud, externally illuminated by a standard G_0 = 1 interstellar radiation field, can produce absorption features in both transitions. Thus, the commonly observed [O I] and [C II] self-absorption and absorption features plausibly arise from the subthermally excited, externally illuminated, photodissociated envelopes of molecular clouds. The luminous young stellar object AGAL337.916-00.477, located precisely where the expanding bubble strikes the Nessie filament, is associated with two shock tracers: NH3 (3,3) maser emission and SiO 2-1 emission, indicating interaction between the bubble and the filament. The interaction of the expanding bubble with its parental dense filament has triggered star formation.

The relationship between B-field orientation and density structure in molecular clouds is often assessed using the Histogram of Relative Orientations (HRO). We perform a plane-of-the-sky geometrical analysis of projected B-fields, by interpreting HROs in dense, spheroidal, prestellar and protostellar cores. We use James Clerk Maxwell Telescope (JCMT) POL-2 850 $\mu$m polarisation maps and Herschel column density maps to study dense cores in the Ophiuchus molecular cloud complex. We construct two-dimensional core models, assuming Plummer column density profiles and modelling both linear and hourglass B-fields. We find high-aspect-ratio ellipsoidal cores produce strong HRO signals, as measured using the shape parameter $\xi$. Cores with linear fields oriented &lt; 45^{\circ} from their minor axis produce constant HROs with -1 &lt; \xi &lt; 0, indicating fields are preferentially parallel to column density gradients. Fields parallel to the core minor axis produce the most negative value of $\xi$. For low-aspect-ratio cores, $\xi \approx 0$ for linear fields. Hourglass fields produce a minimum in $\xi$ at intermediate densities in all cases, converging to the minor-axis-parallel linear field value at high and low column densities. We create HROs for six dense cores in Ophiuchus. $\rho$ Oph A and IRAS 16293 have high aspect ratios and preferentially negative HROs, consistent with moderately strong-field behaviour. $\rho$ Oph C, L1689A and L1689B have low aspect ratios, and $\xi \approx 0$. $\rho$ Oph B is too complex to be modelled using a simple spheroidal field geometry. We see no signature of hourglass fields, agreeing with previous findings that dense cores generally exhibit linear fields on these size scales.

We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, suggesting that grain growth appears to be significant in disks at earlier protostellar phases. For the rest of the protostars, the dust polarization traces the magnetic field, whose morphology can be approximately classified into three categories: standard-hourglass, rotated-hourglass (with its axis perpendicular to outflow), and spiral-like morphology. 40.0% (+-3.0%) of the protostars exhibit a mean magnetic field direction approximately perpendicular to the outflow on several 100--1000 au scales. However, in the remaining sample, this relative orientation appears to be random, probably due to the complex set of morphologies observed. Furthermore, we classify the protostars into three types based on the C17O (3--2) velocity envelope's gradient: perpendicular to outflow, non-perpendicular to outflow, and unresolved gradient (<1.0~km/s/arcsec). In protostars with a velocity gradient perpendicular to outflow, the magnetic field lines are preferentially perpendicular to outflow, most of them exhibit a rotated hourglass morphology, suggesting that the magnetic field has been overwhelmed by gravity and angular momentum. Spiral-like magnetic fields are associated with envelopes having large velocity gradients, indicating that the rotation motions are strong enough to twist the field lines. All of the protostars with a standard-hourglass field morphology show no significant velocity gradient due to the strong magnetic braking.

We present Atacama Large Millimeter/submillimeter Array observations of the $\sim$10 kAU environment surrounding 21 protostars in the Orion A molecular cloud tracing outflows. Our sample is composed of Class 0 to flat-spectrum protostars, spanning the full $\sim$1 Myr lifetime. We derive the angular distribution of outflow momentum and energy profiles and obtain the first two-dimensional instantaneous mass, momentum, and energy ejection rate maps using our new approach: the Pixel Flux-tracing Technique (PFT). Our results indicate that by the end of the protostellar phase, outflows will remove $\sim$2 to 4 M$_\odot$ from the surrounding $\sim$1 M$_\odot$ low-mass core. These high values indicate that outflows remove a significant amount of gas from their parent cores and continuous core accretion from larger scales is needed to replenish core material for star formation. This poses serious challenges to the concept of ``cores as well-defined mass reservoirs", and hence to the simplified core-to-star conversion prescriptions. Furthermore, we show that cavity opening angles, and momentum and energy distributions all increase with the protostar evolutionary stage. This is clear evidence that even garden-variety protostellar outflows: (a) effectively inject energy and momentum into their environments on $10$ kAU scales, and (b) significantly disrupt their natal cores, ejecting a large fraction of the mass that would have otherwise fed the nascent star. Our results support the conclusion that protostellar outflows have a direct impact on how stars get their mass, and that the natal sites of individual low-mass star formation are far more dynamic than commonly accepted theoretical paradigms.