Interferometric observations of magnetic fields in forming stars

The magnetic field is a key ingredient in the recipe of star formation. Over the past two decades, millimeter and submillimeter interferometers have made major strides in unveiling the role of the magnetic field in star formation at progressively smaller spatial scales. From the kiloparsec scale of molecular clouds down to the inner few hundred au immediately surrounding forming stars, the polarization at millimeter and submillimeter wavelengths is dominated by polarized thermal dust emission, where the dust grains are aligned relative to the magnetic field. Interferometric studies have focused on this dust polarization and occasionally on the polarization of spectral-line emission. We review the current state of the field of magnetized star formation in the context of several questions that continue to motivate the studies of high- and low-mass star formation. By aggregating and analyzing the results from individual studies, we come to several conclusions: (1) Magnetic fields and outflows from low-mass protostellar cores are randomly aligned, suggesting that the magnetic field at ~1000 au scales is not the dominant factor in setting the angular momentum of embedded disks and outflows. (2) Recent measurements of the thermal and dynamic properties in high-mass star-forming regions reveal small virial parameters, challenging the assumption of equilibrium star formation. However, we estimate that a magnetic field strength of a fraction of a mG to several mG in these objects could bring the dense gas close to a state of equilibrium. Finally, (3) We find that the small number of sources with hourglass-shaped magnetic field morphologies at 0.01 -- 0.1 pc scales cannot be explained purely by projection effects, suggesting that while it does occur occasionally, magnetically dominated core collapse is not the predominant mode of low- or high-mass star formation. [Abridged]


INTRODUCTION
Magnetic fields are known to play a critical role in many aspects of both low-and high-mass star formation. Even weakly ionized star-forming material is coupled to the ambient magnetic field, and thus the field can regulate (or prevent) the collapse and fragmentation of star-forming clouds (Hennebelle and Inutsuka, 2019;Krumholz and Federrath, 2019;Teyssier and Commerçon, 2019, in this volume), can influence the formation of protoplanetary disks (Wurster and Li, 2018, in this volume), and can launch bipolar outflows and jets from young protostars (Pudritz and Ray, 2019, in this volume). Mapping the morphology of magnetic fields in low-and high-mass star-forming regions is therefore critical to better understand how magnetic fields affect the star-formation process at early times, and how the role of the field changes relative to other dynamical effects (e.g., turbulence, rotation, thermal and radiation pressure, and gravitational collapse) as a function of spatial scale, source environment, and source mass. Over more than fifty years, studies of magnetic fields in low-and high-mass star-forming regions have been performed across a wide range of spatial scales, from the > 100 pc scale of molecular clouds, to the 1 pc scale of clumps, to the 0.1 pc scale of dense cores, and finally to the 1000-100 au scale of protostellar envelopes surrounding forming protostellar systems. 1 In this review we introduce the typical tracers of the magnetic field in star-forming regions, as well as the methods used to analyze the observations. We then discuss the state of the field of magnetized star formation, focusing primarily on the last three decades of high-resolution polarization studies by millimeter and submillimeter (hereafter combined into "(sub)millimeter") interferometers including the Berkeley-Illinois-Maryland Association (BIMA) millimeter array, the Combined Array for Research in Millimeter-wave Astronomy (CARMA), the Submillimeter Array (SMA), and the Atacama Large Millimeter/submillimeter Array (ALMA). 2

Polarized dust emission
Polarized thermal emission from dust grains is the main tracer of the magnetic field in star-forming regions observed at high resolution and at (sub)millimeter wavelengths. Under most circumstances, oblong interstellar medium (ISM) dust grains with sizes 100 µm are aligned with their long axes perpendicular to magnetic field lines (e.g., Hildebrand, 1988). The currently accepted way to achieve this alignment is via the "radiative torque" (RAT) mechanism, where an anisotropic radiation field (e.g., the external UV field in the ISM, or the radiation from a deeply embedded protostar) causes grains to become aligned relative to the magnetic field (Lazarian, 2007;Hoang and Lazarian, 2009;Andersson et al., 2015). 3 Thus, at the physical scales of star-forming clouds, cores, and envelopes (i.e., scales 100 au), magnetically aligned dust grains emit thermal radiation that is polarized perpendicular to the magnetic field. Observations of dust polarization, which comprise the vast majority of the single-dish and interferometric polarization observations to date, are discussed at length throughout this article.

Polarized molecular-line emission
Polarization of molecular-line emission is another tracer of the magnetic field in star-forming regions. Molecular and atomic lines are sensitive to magnetic fields, which cause their spectral levels to split into magnetic sub-levels. For some molecules, linear polarization can arise when an anisotropy in the radiation and/or velocity field yields a population of magnetic sub-levels that are not in local thermodynamic equilibrium (LTE); this is known as the Goldreich-Kylafis (G-K) effect. Polarization from the G-K effect is most easily detected where the spectral line emission has an optical depth τ ≈ 1, when the ratio of the collision rate to the radiative transition rate (i.e., the spontaneous emission rate) is ∼ 1, and where the gradient in the radiation and/or velocity field is large. The effect is strongest in the lowest rotational transitions of simple molecules such as CO, CS, HCN, SiO, or HCO + . Spectral-line polarization from the G-K effect can be parallel or perpendicular to the magnetic field. Due to the different optical depths of the parallel and perpendicularly polarized components in different locations on the sky, polarization from the G-K effect ultimately traces the plane-of-sky magnetic field orientation with an ambiguity of 90 • Kylafis, 1981, 1982;Kylafis, 1983;Deguchi and Watson, 1984;Lis et al., 1988).
The first detections of the G-K effect were by Glenn et al. (1997), who presented National Radio Astronomy Observatory (NRAO) 12 m telescope observations of the evolved star IRC +10216, and by Greaves et al. (1999), who observed the galactic center and the high-mass star-forming clouds S140 and DR 21 using the James Clerk Maxwell Telescope (JCMT). In addition to these first observations, others have detected the G-K effect in molecular outflows from low-mass protostellar sources (e.g., Girart et al., 1999;Lai et al., 2002;Lee et al., 2014;Ching et al., 2016; as well as in high-mass protostellar sources (e.g., Lai et al., 2003;Cortes et al., 2005;Cortes et al., 2008). The G-K effect is a powerful way to trace the magnetic field in regions where the brightness of the thermal dust emission is too low to detect polarization at the few-percent level, either because the source is too distant or because the column density of dust is too low (e.g., in an outflow cavity where the gas and dust have been driven away).
One way to resolve the 90 • ambiguity in G-K observations is via radiative transfer modeling. This was done in Cortes et al. (2005), who expanded on the original G-K models (which assumed that gradients in the CO optical depth were responsible for the necessary anisotropy in the radiation field) by including bright, central sources as additional causes of anisotropy. They used these models to successfully reproduce the 90 • difference in polarization angles seen toward the high-mass star-forming region DR 21(OH) in their 3 mm BIMA observations of CO (J = 1 → 0) versus earlier 1.3 mm observations of CO (J = 2 → 1) by Lai et al. (2003). While the modeling in Cortes et al. (2005) was successful, in general it is difficult to know the structure of the radiation field, especially in high-mass star-forming regions. Other methods exist to resolve the ambiguity, such as comparing polarization from both dust and spectral lines in the same region. This was done in DR 21(OH) by Lai et al. (2003) and Ching et al. (2016); however, it is not clear how strong the correlation should be between the polarization angles from the two types of emission, as the dust emission traces much denser material than the CO.
Looking to the future, in low-mass sources with well behaved outflows, or in protoplanetary disks, the radiation and velocity fields may be simple enough to allow us to more confidently break the ambiguity in polarization orientation. Given the high quality of ALMA data toward these types of sources and the fact that ALMA's sensitivity will allow us to probe polarization in multiple transitions of many different molecules, spectral-line polarization studies have a bright future.

Frontiers
Finally, one must exercise caution when interpreting polarized spectral-line observations, as linearly polarized spectral-line emission can be converted into circularly polarized emission via anisotropic resonant scattering. Several studies have detected non-Zeeman circularly polarized emission, including Houde et al. (2013) in Orion KL and Hezareh et al. (2013) in the supernova remnant SNR IC 443. These observations were from the polarimeters at the Caltech Submillimeter Observatory (CSO) and the Institut de Radioastronomie Millimétrique (IRAM) 30 m telescope, respectively (see Chamma et al. 2018 for more results from the SMA). The model described in those publications shows that resonant scattering can cause a rotation in the polarization position angle in maps of linear spectral-line polarization. Hezareh et al. (2013) used the detected Stokes V signal toward SNR IC 443 to correct the map of linear spectral-line polarization. After doing so, the polarization was well aligned with dust polarization observations using the PolKA polarimeter at the Atacama Pathfinder Experiment (APEX) telescope. The authors note that this effect is proportional to the square of the magnitude of the plane-of-sky magnetic field, which opens up the possibility of detecting the magnetic field via circular polarization observations of Zeeman-insensitive spectral lines.

The Velocity-gradient technique
Another method for probing the plane-of-sky magnetic field in star-forming regions is the "velocitygradient technique" (VGT), which estimates the magnetic field orientation using the velocity gradients present in spectral-line observations. In turbulent regions that are not gravitationally bound, González-Casanova and  and Yuen and Lazarian (2017b) showed that the magnetic field is expected to be oriented perpendicular to velocity gradients. The authors also extended their study to shocked and gravitationally bound regions (Yuen and Lazarian, 2017a), and argued that the VGT provides an even better view of the plane-of-sky magnetic field morphology in the interstellar medium (particularly in diffuse regions with a low column density of dust) than both the Planck polarization data (e.g., Planck Collaboration et al., 2015b) and observations of neutral HI "fibers," which also trace the interstellar magnetic field (Clark et al., 2014(Clark et al., , 2015.

Magnetic field tracers (line-of-sight component): the Zeeman effect
The Zeeman effect is another important tracer of the magnetic field that has been observed primarily with single-dish radio telescopes to infer the line-of-sight magnetic field strength 4 in molecular clouds (Troland and Heiles, 1986;Crutcher et al., 1993;Crutcher, 1999). When threaded by a magnetic field, atomic hydrogen and molecules with a strong magnetic dipole moment will have the degeneracy in magnetic sub-levels lifted for states with non-zero angular momentum. This will split the radio frequency transitions into a number of linearly and elliptically polarized components separated slightly in frequency. Measuring this Zeeman splitting is the only way to directly measure a component of the magnetic field strength. However, we will not focus more on the thermal Zeeman effect in this review, as there are no reported observations with a (sub)millimeter-wavelength interferometer. For reviews of single-dish observations of magnetic fields in molecular clouds via the Zeeman effect, see Crutcher (2012) and Crutcher and Kemball (2019, in this volume). 4 Note that under certain conditions one can derive the total magnetic field strength from the Zeeman effect Heiles and Crutcher, 2005). This has been seen several times toward galactic OH masers (Hutawarakorn et al., 2002;Fish and Reid, 2006;Caswell et al., 2014). However, in typical Zeeman observations of star-forming molecular clouds, the Zeeman signal is only strong enough to yield the strength of the magnetic field along the line of sight.

Indirect estimates of magnetic field strength
The polarization arising from magnetically aligned dust grains and from the G-K effect yields the planeof-sky magnetic field morphology. However, these observations do not contain information about the magnetic field strength. Since knowing the field strength is critical to our understanding of the importance of magnetic fields, a variety of indirect methods have thus been developed to estimate the field strength from these types of observations. The most longstanding of these methods is the Davis-Chandrasekhar-Fermi (DCF) method (Davis, 1951;Chandrasekhar and Fermi, 1953), which uses the fact that turbulent motions should cause an observable scatter in what would otherwise be a well ordered mean magnetic field. The original, simplest form of the DCF formula yields an estimate of the plane-of-sky magnetic field strength where ρ is the gas density, δV is the one-dimensional velocity dispersion, and δφ is the dispersion in polarization position angles. δV and δφ are both assumed to be caused by turbulent motions in the region being studied. The derivation of this expression also assumes that there is a mean magnetic field in the region, that the turbulence is incompressible and isotropic, and that the turbulent components of the kinetic and magnetic energies are in equipartition. Note that the DCF method also assumes that the dispersion is "small," i.e., that the turbulent magnetic energy is small compared with the mean-field magnetic energy in the system.
The DCF method was originally developed to estimate magnetic field strengths in the diffuse ISM, where perturbations in the magnetic field can be assumed to be caused by turbulent motions in the magnetized medium. Comparisons with MHD simulations of giant molecular clouds (GMCs) by Ostriker et al. (2001) found that at these scales, the original DCF method typically overestimates the magnetic field strength by a factor of ∼ 2 because line-of-sight field tangling is not taken into account. Further simulations by Heitsch et al. (2001) also found that the true field strength is overestimated unless finite telescope resolution and self-gravity within the GMC are accounted for.
At the smaller scales of the individual star-forming clumps and cores that are probed by interferometers, gravity is almost always the dominant dynamical factor, and thus the structure of the background field must be removed before calculating the dispersion in polarization position angles. In two ideal cases with very clean examples of hourglass-shaped magnetic fields, Girart et al. (2006) and Qiu et al. (2014) removed the background hourglass field by subtracting a set of parabolic field lines from the data, after which they calculated the magnetic field dispersion relative to the background structure. In complicated cases, however, a more general method is necessary to remove arbitrary background field structure. This has been achieved by employing a second-order structure function of the polarization angle that allows the separation of the turbulent and mean magnetic-field components, with the underlying assumption that the mean-field component has an intrinsically larger spatial scale than the turbulent component (Falceta- Gonçalves et al., 2008;Hildebrand et al., 2009;Houde et al., 2009Houde et al., , 2011Chitsazzadeh et al., 2012). Further refinements of these studies have also taken into account interferometric filtering effects (Houde et al., 2016).
Falceta- Gonçalves et al. (2008) used the structure-function approach to test the validity of the DCF technique using their MHD simulations. In addition to analyzing the effect of different telescope resolutions Frontiers (their conclusion: lower-resolution observations tend to overestimate the field strength), they also derived a generalized form of the DCF equation, which allows for the separation of the turbulent and mean magnetic field components, and yielded magnetic-field estimates that were accurate to within ∼ 20%. More recently, Juárez et al. (2017) performed SMA observations of a magnetized high-mass star-forming region and used the structure-function method to compare the data with synthetic observations of gravity-dominated MHD simulations. They found the magnetic field strength estimates from both the observations and simulations to be in good agreement. Koch et al. (2012) proposed another method of indirectly measuring the magnetic field strength that is also based on an ideal MHD framework. They argue that a position-dependent magnetic field strength in a polarization map can be calculated using the angle between the magnetic field and the gradient of the total-intensity emission, and a second angle between the local gravity direction and the gradient of the total-intensity emission. This is based on the assumption that the intensity gradient is a proxy for the direction of motion in the MHD force equation. For example, in a centrally concentrated, gravitationally bound core, the infalling motion will be along the intensity gradient (across the iso-intensity contours). This method requires observations of the magnetic-field morphology (i.e., via observations of thermal dust polarization or the G-K effect) in order to produce a spatial distribution of the magnetic field strength. However, unlike the DCF method, this method does not require spectral-line observations in the analysis. Estimations of magnetic field strengths using this method have been found to be consistent with previous estimations from the DCF method (see, e.g., Section 3.1).

The mass-to-flux ratio
Merely measuring the magnetic field strength does not allow us to determine immediately the importance of the magnetic field in a given star-forming region. Therefore, in magnetic field studies it is common to compare the magnetic energy density with that of other dynamical quantities such as gravity, turbulence, and outflow feedback (see, e.g., Hull et al. 2017b for a comparison of magnetic energy with gravitational and outflow energy densities in the intermediate-mass Class 0 protostellar core Serpens SMM1).
Historically, there has been a strong focus on the comparison of gravity (which causes inward motion) and the magnetic field (whose tension provides resistance against infall across the field lines). The common quantity quoted as a metric for the importance of the magnetic field with respect to gravity is the "massto-flux ratio" M/Φ, where M is the mass of the object of interest and Φ is the flux of the magnetic field threading the object. As discussed in Crutcher (2004, and references therein), the maximum mass that can be supported by a given magnetic flux is given by M crit = Φ/2π √ G. However, it is more useful to discuss the dynamical status of an object by measuring the ratio λ of the observed mass-to-flux ratio to the critical mass-to-flux ratio: Clouds that are supported by the magnetic field and are not collapsing are deemed "subcritical" (λ < 1), whereas those where gravity has overcome the resistance of the magnetic field are referred to as "supercritical" (λ > 1). Crutcher (2012) analyzed data across a wide range of spatial scales that trace more than five orders of magnitude in densities, and found that when the hydrogen column density N H > 10 21 cm −2 , all star-forming objects are supercritical (i.e., are collapsing). This value is less than the typical column densities of low-mass (N H ∼ 10 22 − 10 23 cm −2 ; e.g., Girart et al. 2006;Hull et al. 2017b) and high-mass (N H ∼ 10 23 − 10 24 cm −2 ; e.g., Girart et al. 2009) protostellar cores, and thus the types of objects we review in this article are all supercritical. This is reasonable, as most of them have already formed stars, as revealed by the presence of bipolar outflows. Furthermore, due to the sensitivity limits of CARMA and the SMA, most of the sources in previous interferometric surveys of polarization Zhang et al., 2014) were chosen based on their strong millimeter flux, which correlates with the presence of embedded star formation.

Core-mass estimates from dust emission
In order to convert the observed millimeter-wave flux density S ν contained within a given spatial area into a corresponding gas mass M gas , we can use the following relation: where d is the distance to the source, κ ν is the opacity of the dust (Ossenkopf and Henning, 1994), and B ν (T d ) is the Planck function at the frequency of the observations. T d is the temperature of the dust, which is usually ∼ 20-50 K in a low-mass protostellar core (Girart et al., 2006), and as high as (or greater than) 100 K in a high-mass core . Once the dust mass is calculated, a gas-to-dust mass ratio of 100 is usually assumed in order to calculate the total (gas + dust) mass of the protostellar core. Note that such an estimate does not include the mass of the central star(s), which must be obtained by other means, e.g., via direct detection of a Keplerian disk around the source (e.g., Tobin et al., 2012;Ohashi et al., 2014) or via determination of the source's bolometric luminosity.

Motivating questions in low-and high-mass star formation
The primary goal of observing the magnetic field at any spatial scale is to determine the importance (or lack thereof) of the magnetic field in the star-formation process. The steady progress toward this goal over the last two decades began with single-dish submillimeter polarization surveys probing 20 scales using the Viper 2 m telescope at the South Pole (SPARO polarimeter, e.g., Dotson et al. 1998;Renbarger et al. 2004;Li et al. 2006), the JCMT (850 µm SCUBA polarimeter, e.g., Matthews et al. 2009), and the CSO (350 µm SHARP  and Hertz [Dotson et al. 2010] polarimeters). A resurgence of single-dish studies has been brought about by results from the PolKa polarimeter at the APEX telescope (Siringo et al., 2004(Siringo et al., , 2012Hezareh et al., 2013;Wiesemeyer et al., 2014;Alves et al., 2014), the BISTRO survey with the upgraded POL-2 polarimeter at the JCMT (e.g., Ward-Thompson et al., 2017;Pattle et al., 2017;Pattle et al., 2018;Soam et al., 2018), results from the polarimeter on the Balloon-borne Large Aperture Submillimeter Telescope (BLAST; Roy et al., 2011;Gandilo et al., 2016;Fissel et al., 2016), observations from the HAWC+ polarimeter (Vaillancourt et al., 2007) on the Stratospheric Observatory for Infrared Astronomy (SOFIA; e.g., Chuss et al., 2018;Gordon et al., 2018;Lopez-Rodriguez et al., 2018), and the galactic polarization maps from the Planck satellite (e.g., Planck Collaboration et al., 2015bCollaboration et al., ,a,c, 2016a. These new studies will pave the way for future work with even more sensitive instruments such at the next-generation BLAST instrument (BLAST-TNG; Galitzki et al. 2014), the NIKA2 polarimeter at the IRAM 30 m telescope (Ritacco et al., 2017), and the TolTEC polarimeter at the Large Millimeter Telescope (LMT). For reviews on multi-scale/multi-wavelength studies and single-dish observations of magnetic fields, respectively, see Li and Law (2019) and Pattle and Fissel (2019), both in this volume.
There are a number of questions applicable to both low-and high-mass star formation that have been investigated using (sub)millimeter polarimetric observations. These include the (direct or indirect) measurements of the absolute magnetic field strength in star-forming material at different spatial scales, Frontiers as well as the estimation of the dynamical importance of the magnetic field with respect to gravity (i.e., the mass-to-flux ratio; see Section 1.3.2). Observations of magnetic fields across multiple spatial scales toward both low-mass (Li et al., 2009;Hull et al., 2014) and high-mass Li et al., 2015) protostellar sources have also been used to constrain the dynamical importance of the magnetic field based on the morphological consistency (or lack thereof) of the field as a function of scale. Generally, however, these observations have only compared two or three of the relevant spatial scales (i.e., 100 pc cloud scales, 1 pc clump scales, 0.1 pc dense-core scales, 1000 au protostellar envelope scales, and 100 au disk scales). A full characterization of the magnetic field from galactic (Planck) scales down to scales approaching the 100 au size of protoplanetary disks has yet to be accomplished, but will be possible in the near future when upcoming polarization surveys of the full populations of protostars in entire molecular clouds are completed.
In the low-mass regime, single-dish observations probe the magnetic field in star-forming clouds at large scales, revealing the magnetic field from the scale of entire molecular clouds (Planck, BLAST) to the canonical, ∼ 0.1 pc dense core, where one or a few protostars will form (JCMT, CSO). One of the main benefits of single-dish studies is their ability to recover a larger range of spatial scales than interferometers, thus enabling an accurate characterization of the magnetic field in ambient cloud material. However, higher resolution is needed in order to probe the environments of individual stars; this is where results from the BIMA, CARMA, SMA, and ALMA interferometers dominate the discussion, allowing us to characterize the magnetic field from scales of several × 1000 au down to the scales of a few × 10 au accessible to ALMA. 5 The main questions that have been tackled over the last two decades using data from these interferometers include: (1) What is the importance of the magnetic field in regulating the collapse of star-forming cores? (Section 2.1); (2) What is the relationship of bipolar outflows with the magnetic field? (Section 2.2); and (3) What is the role of the magnetic field in the launching and collimation of bipolar outflows in low-mass protostars? (Section 2.3) Moving to the high-mass regime: high-mass stars (M * > 8 M ) form predominantly in clustered environments where a population of stars are born with a range of stellar masses (Lada and Lada, 2003). These high-mass stellar populations form in dense cores that are embedded in parsec-scale, massive molecular clumps. This review will focus on several questions raised in recent studies of these high-mass sources, including (4) What is the dynamical role of magnetic fields in dense cores? (Section 3.1); (5) What is the role of the magnetic field in the formation of disks and the launching of protostellar outflows in high-mass protostars? (Section 3.2); (6) Do magnetic fields play a significant role in the fragmentation of molecular clumps and the formation of dense cores? (Section 3.3); and (7) Does high-mass star formation proceed in virial equilibrium? (Section 3.4) The dense clustering nature of high-mass star formation implies considerable fragmentation within these massive molecular clumps, which distinguishes high-mass star formation from the more isolated process of low-mass star formation. High-mass regions also tend to have much more intense radiation environments, hosting HII regions whose radiative feedback and ionization can impact the ambient magnetic field. Another clear difference arises from the fact that the best-studied low-mass stars are forming at distances ∼ 10 × closer than typical high-mass star forming regions, allowing us to study the formation of individual low-mass protostellar systems and their associated outflows, jets, and disks in much greater detail than is possible in high-mass systems. Ultimately, however, studies of both low-and high-mass star formation use the same observing techniques and confront many of the same questions. In this review we focus on those questions that have been of the most interest to both communities in recent years.

MAGNETIC FIELDS IN LOW-MASS STAR FORMATION
The revolution of high-resolution, interferometric observations of polarization began with BIMA and the Owens Valley Radio Observatory (OVRO). These two sets of antennas were later combined into CARMA (Bock et al., 2006). Early observations with BIMA and OVRO 6 covered a wide range of topics, including polarization observations of dust, SiO, CO, and SiO masers toward iconic regions in Orion (Rao et al., 1998;Plambeck et al., 2003;Girart et al., 2004; as well as observations of individual protostars (Girart et al., 1999;Kwon et al., 2006;. These first observations, combined with the extensive follow-up from CARMA, the SMA, ALMA form the body of work motivating this review. Below we put this work in the context of a narrative addressing several of the major open questions in the field of magnetized low-mass 7 star formation.

The role of the magnetic field in protostellar collapse
2002), might occur if, like their low-speed counterparts in Figure 5, they escape preferentially along the low-density axial direction of the ambient medium. In any case, interaction between the fast-moving jet or wind of an embedded protostar appears u this intera The slo lapse proc region (wh toward th The pseud and magn tional pul falls inwa angular m threading more slow twisting o Alfvén wa plane. Th pseudodis lope mate low-speed where the twisted (se plasma almost ins twisted en show belo outflow gr The ang braking is showing an hourglass-shaped magnetic field configuration in a collapsing, magnetized starforming core. Reproduced with permission from the American Astronomical Society (AAS).
In models of magnetically regulated protostellar collapse (e.g., Mouschovias, 1976aMouschovias, ,b, 1991Mouschovias and Ciolek, 1999), a strong, well ordered magnetic field provides outward pressure support of the infalling material. This is because the field is coupled (or "frozen") to the small fraction of charged particles in the weakly ionized gas. However, in non-turbulent models, the non-ideal MHD effect of ambipolar diffusion (Mestel and Spitzer, 1956) enables the neutral material (which comprises the bulk of the star-forming core) to slip slowly past the magnetic field lines, thus removing magnetic flux and eventually allowing collapse to proceed once the mass-to-flux ratio exceeds the critical value.
One of the predicted signposts of highly magnetized star formation is that at high enough densities ( 10 4 cm −3 ), the collapse of strongly magnetized gas should pinch the magnetic field into an "hourglass" shape with a symmetry axis perpendicular to the major axis of a flattened, ∼ 1000 au "pseudodisk" (Galli and Shu, 1993a,b). The hourglass is expected to persist down to scales < 1000 au (Fiedler and Mouschovias, 1993;Galli and Shu, 1993b;Allen et al., 2003;Gonçalves et al., 2008;Frau et al., 2011;Kataoka et al., 2012;Mocz et al., 2017); see Figure 1. And indeed, the predicted hourglass has now been seen in a number of interferometric observations of low-mass protostellar cores (Girart et al., 1999(Girart et al., , 2006Rao et al., 2009;Stephens 6 The two polarization results from OVRO are observations toward NGC 1333-IRAS 4A and IRAS 16293 (Akeson et al., 1996;Akeson and Carlstrom, 1997); however, OVRO was known to have issues with polarization calibration, which is the most likely explanation for the inconsistency of those results with later observations of the same sources (Girart et al., 1999(Girart et al., , 2006Rao et al., 2009). 7 We do not treat the topic of intermediate-mass star formation separately in this review. As many of the characteristics of the early stages of intermediate-and low-mass star formation are thought to be similar (Beltrán, 2015), we include references to several intermediate-mass objects and regions in this section. Many of these intermediate-mass sources are in Orion (e.g., Takahashi et al., 2006Takahashi et al., , 2018Hull et al., 2014), but there are also objects in other regions, such as Serpens A second signpost of strong-field star formation is the consistency of the magnetic field orientation across multiple spatial scales. If the field is strong relative to other dynamical effects, observations at small scales should reveal a magnetic field whose original orientation is preserved from the parsec scale of the cloud in which the source is embedded. To date, multi-scale studies of the magnetic field in low-and high-mass star-forming regions have compared two or three scales: i.e., ∼ kpc galactic scales to ∼ 0.1 pc dense-core scales in Stephens et al. (2011); 100 pc cloud scales to 0.1 pc dense-core scales in Li et al. (2009Li et al. ( , 2015; 1 pc clump scales to 0.1 pc dense-core scales in Zhang et al. (2014), continuing down to ∼ 0.01 pc scales in Girart et al. (2013) and Ching et al. (2017); 0.1 pc dense core scales to 1000 au protostellar envelope scales in Hull et al. (2014) and Davidson et al. (2014); and 0.1 pc to 1000 au to 100 au scales in Hull et al. (2017a,b).
In the low-mass regime, Li et al. (2009) found striking consistency between the magnetic field orientation in the Orion molecular cloud derived from background starlight polarization at ∼ 100 pc scales versus polarized thermal dust emission at ∼ 0.1 pc scales. Hull et al. (2014) took this one step further, finding consistency in the field orientation in just a few of the low-mass protostars in their sample from scales of ∼ 0.1 pc to ∼ 1000 au. The sample of sources from Hull et al. (2014) that maintained consistency in the magnetic field orientation down to scales of ∼ 1000 au tended to be those sources with a higher polarization fraction, which implies that the magnetic fields in those sources are more ordered, and thus may be more dynamically important.
Among those sources with consistent large-to-small scales magnetic fields are several with known hourglass morphologies, including OMC3 MMS6 , NGC 1333-IRAS 4A (Girart et al., 2006), L1157 (Stephens et al., 2013), and L1448 IRS 2 . The magnetic field strengths have been estimated toward the latter three objects, an are all relatively high, on the order of 1 mG, which is similar to values obtained in high-mass regions (see Section 3.1). However, while the values are high (5 mG in IRAS 4A [Girart et al. 2006]; 1.3-3.5 mG in L1157 [Stephens et al. 2013]; and 750 µG in L1448 IRS 2 ), the mass-to-flux ratios calculated for the two most magnetized sources (IRAS 4A and L1157) are both slightly greater than the critical value (1.7 and 1.1 for IRAS 4A and L1157, respectively), which is reasonable, considering that the objects have already collapsed to form protostars.
Follow-up studies with ALMA of individual sources from previous surveys (e.g., Hull et al. 2014) have suggested that consistency in the magnetic field orientation across spatial scales may be the exception rather than the rule at scales smaller than ∼ 0.1 pc. Indeed, detailed multi-scale follow-up studies by Hull et al. (2017a,b) of Ser-emb 8 and Serpens SMM1, two Class 0 protostellar sources in the Serpens Main molecular cloud, found significant deviations in the magnetic field morphology across spatial scales ranging from ∼ 0.1 pc down to ∼ 100 au.
In the case of Serpens SMM1, it appears that the deviations at the ∼ 100 au scales probed by ALMA are due to shaping of the magnetic field by the bipolar outflow. This is in spite of the fact that the magnetic field strength is estimated to be quite high, ∼ 5.7 mG (Hull et al., 2017b). The fact that dynamics are shaping the magnetic field morphology in this source is not surprising, however, as the magnetic is comparable to the kinetic energy in the outflow. This is in contrast to sources like NGC 1333-IRAS 4A, where estimates by Girart et al. (2006) show that the magnetic energy is clearly dominant over turbulent motions in the source.
In the case of Ser-emb 8, the outflow does not have a clear effect on the magnetic field, and yet the field morphology is still not consistent across scales. A comparison of the Ser-emb 8 data with moving-mesh AREPO gravoturbulent MHD simulations (Mocz et al., 2017) suggests that Ser-emb 8 may have formed in an environment where dynamical effects such as turbulence and infall dominate the magnetic field, in contrast to the strong-field examples described above; see Figure 3. These results are qualitatively consistent with other simulations studying star formation in weakly magnetized turbulent environments (e.g., Seifried et al., 2015;Offner and Chaban, 2017), suggesting that the importance of the magnetic field at the smallest scales of star formation lies along a continuum, from the bright, highly polarized examples of strong-field star formation to the more complex examples of star formation in regions dominated by dynamical processes. Large, high-resolution polarization surveys by ALMA will reveal the distribution of low-mass protostars across this continuum of magnetic importance, and will allow us to compare with results from high-mass regions such as those by Zhang et al. (2014), who found that the magnetic field is dynamically important at the larger spatial scales probed by their observations (see Section 3.3).

Misalignment of outflows and magnetic fields
For more than a decade, one of the primary ways that astronomers have tested the importance of the magnetic field in star-forming regions has been via comparisons of the orientations of bipolar outflows/jets and the ambient magnetic field. This is because if a protostellar core is very strongly magnetized, the magnetic field has the ability to align all of the relevant axes: the axis of the (well ordered and poloidal) magnetic field, the symmetry axis of the pseudodisk, the rotational axis of the protostellar disk, and the axis of the outflow and/or jet emanating from the central source. This happens because of the "magnetic  Hull et al. 2014), and ALMA at 870 µm (right). The red and blue arrows indicate the red-and blueshifted lobes of the bipolar outflow from Ser-emb 8 traced in CO (J = 2 → 1). The data show the inconsistency of the magnetic field across spatial scales, and are consistent with moving-mesh AREPO MHD simulations (bottom) whose initial conditions include a magnetic field that is weak compared with turbulence at large (5 pc) scales. Adapted from Hull et al. (2017a). Reproduced with permission from the AAS.
braking" phenomenon (Allen et al., 2003), where a strong magnetic field removes angular momentum from the central source and causes the angular momentum (and thus disk/outflow) and magnetic axes to align (e.g., Machida et al., 2006). This same magnetic braking phenomenon can potentially lead to what has come to be known as the "magnetic braking catastrophe," where a strong magnetic field aligned with the core rotation axis can suppress the formation of a Keplerian disk in MHD simulations (e.g., Galli et al., 2006;Mellon and Li, 2008;Li et al., 2011). This may lead to the formation of sources like L1157 (Stephens et al., 2013) and B335 , which have outflows aligned with the magnetic field, and which have as-of-yet undetectably small disks ( 10 au). However, since it is known that large Keplerian disks form around many protostellar sources, a variety of methods have been proposed to overcome this problem, including the introduction of an initial misalignment between the rotation axis and the magnetic field, which enhances disk formation (e.g., Hennebelle and Ciardi, 2009;Joos et al., 2012;Krumholz et al., 2013;Li et al., 2013)), and the consideration of non-ideal MHD effects such as ambipolar diffusion (Dapp et al., 2012;Masson et al., 2016;Tsukamoto et al., 2018), Ohmic dissipation (Dapp et al., 2012;Tomida et al., 2015;Tsukamoto et al., 2018), the Hall effect (Tsukamoto et al., 2015a(Tsukamoto et al., ,b, 2017, and magnetic reconnection (Santos-Lima et al., 2012;Li et al., 2014).
Observationally, if very strong magnetic fields were the norm, then the rotational axes of protostellar disks, and the jets and outflows that emanate from them, should all be parallel with the ambient magnetic field. A study of seven low-mass protostellar cores by Chapman et al. (2013) found a correlation between outflows and magnetic fields at ∼ 0.1 pc scales. However, the majority of the studies of this type have come to the opposite conclusion. For example, Ménard and Duchêne (2004) found that the optical jets from classical T Tauri stars in the Taurus-Auriga molecular cloud are randomly oriented with respect to the parsec-scale magnetic field observed via background-starlight polarization observations. Targon et al. (2011) obtained a similar result for 28 regions spread over the Galaxy, finding no strong correlations between protostellar jets and the ambient magnetic field. On the ∼ 0.1 pc scales of high-mass star-forming cores, Curran and Chrysostomou (2007) and Zhang et al. (2014) used thermal dust polarization observations by the JCMT and the SMA, respectively, to determine that outflows and inferred magnetic fields are randomly aligned. Finally, Hull et al. (2013Hull et al. ( , 2014 used the 1.3 mm polarization system at CARMA (Hull and Plambeck, 2015) to observe dust polarization toward a sample of low-and high-mass sources, and found that the outflows and ∼ 1000 au-scale magnetic fields in the low-mass sources were randomly aligned. In Figure 4, we compile all of the outflow-versus-magnetic-field angles derivable to date from interferometric observations of low-mass protostellar cores (Girart et al., 1999(Girart et al., , 2006Rao et al., 2009;Hull et al., 2013;Stephens et al., 2013;Hull et al., 2014Hull et al., , 2017aCox et al., 2018;Galametz et al., 2018;Maury et al., 2018;Sadavoy et al., 2018a;Harris et al., 2018), and come to the same conclusion: while a few sources have well aligned outflows and magnetic fields (e.g., those on the very bottom-left of the plot in Figure 4 that are climbing the 0-20 • curve, several of which have hourglass-shaped field morphologies; see Section 4), overall protostellar outflows and magnetic fields measured at 1000 au-scales are randomly aligned.
The finding from several studies that magnetic fields and outflows are randomly oriented suggests that most protostars form out of material with a magnetic field that is too weak to maintain a consistent orientation all the way down to the 0.1-10 au scales where outflows are launched. Rather, it seems likely that dynamical effects such as turbulence and infall are more important than the magnetic field when it comes to dictating the ultimate angular-momentum direction at the small (< 1000 au) scales relevant for the formation of protostellar disks and outflows/jets. This claim is corroborated by the analysis by Lee et al. The plot shows that outflows appear to be randomly aligned with magnetic fields in the sample of low-mass sources whose magnetic fields have been observed with ∼ 500-1000 au resolution.
(2017) of the synthetic observations (produced using the Dustpol radiative transfer software Padovani et al. 2012) of two MHD simulations with different initial mass-to-flux ratios (Offner and Chaban, 2017). They concluded that while the protostar that formed in the strong-field case exhibited a correlation between the outflow and the magnetic field, the weaker-field case showed a random orientation. This is most likely because asymmetric accretion from the turbulent envelope stochastically reoriented the disk/outflow during the earliest (Class 0) formation stage, ultimately decoupling the source from the natal magnetic field (see, e.g., simulations by Chen and Ostriker 2018).
Two non-magnetic studies focused on the alignment of protostellar outflows with respect to the natal filamentary structure in which the sources formed, and on the relative alignment of the outflows from wide (∼ 1000 au) and tight (< 250 au) binary/multiple systems. Both studies used data from the MASSES survey at the SMA (Lee et al., 2015;Stephens et al., 2018). Regarding outflows versus filaments, Stephens et al. (2017a) studied the relative orientation of 57 protostellar outflows in the Perseus molecular cloud (derived from the SMA MASSES data) with the local filamentary structure (derived from Herschel observations), and found that the orientations are consistent with a random distribution. Their results held regardless of protostellar age, multiplicity, or the opacity of the dense core, suggesting that the angular momentum of the protostellar cores and outflow-launching disks are independent of the large-scale structure out of which the protostars are forming. Regarding the orientation of outflows from multiple systems, Lee et al. (2016) used SMA MASSES observations to determine that the outflows from proto-multiple systems in the Perseus molecular cloud are randomly aligned. Offner et al. (2016) followed up with MHD simulations that are consistent with the SMA observations, arguing that multiple systems with randomly oriented angular momenta are likely to have formed from turbulent fragmentation.
The turbulent fragmentation scenario is consistent with recent ALMA observations of multiple systems with misaligned protoplanetary disks (Jensen and Akeson, 2014;. However, observations by Tobin et al. (2016a) have found evidence for an alternative theory of multiple formation: disk fragmentation (Kratter et al., 2010). The VLA Nascent Disk and Multiplicity Survey of Perseus Protostars (VANDAM) by Tobin et al. (2016b) found a bimodality in the multiplicity of systems; they argued that the wide multiples are the result of turbulent fragmentation, whereas the tight multiples are the result of disk fragmentation. This latter conclusion was strengthened in a follow-up study by Tobin et al. (2018), who observed a sample of tight binaries from the VANDAM survey with ALMA in both continuum and spectral lines.
Recent work by Galametz et al. (2018) focused on the role of the magnetic field in the formation of multiple systems using an SMA survey of magnetic fields in protostellar cores. They found tantalizing evidence that a large misalignment between the outflow and the magnetic field is found preferentially in protostars with higher rotational energies. This ∼ 90 • misalignment observed in some objects could thus be due to the winding of the magnetic field lines in the equatorial plane by strong rotation. Furthermore, they found hints that many of those same sources (i.e., those with approximately perpendicular outflows and magnetic fields) are wide multiple sources and/or have large disks, whereas the sources with well aligned magnetic fields and outflows tend to be single objects with small (or undetected) disks at the ∼ 500 -1500 au resolution of their observations. The trend of large disks being associated with sources that have misaligned magnetic fields and outflows was also tentatively seen in an analysis of Karl G. Jansky Very Large Array (VLA) observations of Class 0 and I protostars by Segura- Cox et al. (2018). These results suggest that the morphology and dynamical importance of magnetic fields at the scale of the protostellar envelope may significantly impact the outcome of protostellar collapse as well as the formation of disks and multiple systems. Large ALMA surveys in polarization toward binary/multiple systems will shed light on the impact of the magnetic field on multiple formation via both turbulent and disk fragmentation.

The importance of the magnetic field in launching jets and outflows
Magnetic fields play a critical role in launching and collimating both bipolar outflows and jets from young forming stars (Frank et al., 2014). Several theories exist to explain how outflows and jets are generated, including the "disk wind" theory where an outflow is launched from the magnetized, rotating surface of a disk (Konigl and Pudritz, 2000), and the "X-wind" theory where jets are launched close to the central protostar itself (Shu et al., 2000). Both of these theories require a magnetic field to function, and that magnetic field is expected to have both poloidal (i.e., along the outflow) and toroidal (perpendicular to the outflow) components due to the combination of infall, outflow, and rotational motions present near a forming star. Characterizing the magnetic field in outflows and jets can thus allow us to investigate the origin of outflows in the context of these different theories.
Historically, observations of dust polarization have been used mainly to study the magnetic field morphology in the optically thin dense cores of dust and gas surrounding embedded protostars. With the sensitivity and resolution of ALMA, it is now possible to detect polarized dust emission along the edges of outflow cavities (Hull et al., 2017b;Maury et al., 2018). However, there is not enough dust in the cavity itself (where the outflow has evacuated most of the material) to allow for a detection of the polarized emission. Therefore, in order to probe the magnetic field in the outflowing material, one must turn to observations of spectral-line polarization (see Section 1.1.2).
While several studies have focused on SiO maser polarization (e.g., Plambeck et al., 2003), the majority of spectral-line polarization studies toward low-mass forming stars have targeted thermal CO and SiO emission polarized due to the G-K effect. Pioneering CO (J = 2 → 1) polarization observations were performed with BIMA by Girart et al. (1999Girart et al. ( , 2004; Kwon et al. (2006); and . So far there has been one detection of SiO (J = 8 → 7) polarization toward the low-mass protostar HH 211, tentatively detected by Lee et al. (2014) using the SMA, and recently confirmed by ALMA observations reported in : see the bottom panel of Figure 5.
The BIMA observations of CO (J = 2 → 1) polarization toward the iconic Class 0 protostar NGC 1333-IRAS 4A by Girart et al. (1999) were the first interferometric detection of the G-K effect. Toward the central core, these observations are consistent with the magnetic field inferred from polarized dust emission. First was the fact that the data are consistent with a magnetic field in IRAS 4A that is poloidal at the base of the outflows (there are two outflows, each launched by a member of the embedded binary) and toroidally wrapped up further out in the outflow cavity. This observation, combined with the coexistence of a low-velocity outflow and a high-velocity jet in the source, led the authors to conclude that the outflows in IRAS 4A are most likely driven by MHD winds from the surface of a rotating disk.

Magnetic field measurements at core scales
The first pioneering high-resolution observations of linearly polarized continuum and spectral-line emission toward high-mass star-forming regions were made with the BIMA interferometer. Rao et al. (1998) reported the first interferometric polarization observations of a high-mass star-forming region, toward Orion KL. They detected linear polarization at both 3.3 mm and 1.3 mm in the BIMA data at a resolution of 1000 -3000 au, revealing abrupt changes in the magnetic field orientations among the continuum emission peaks. This chaotic distribution is in contrast to the uniform magnetic field topology in the lower resolution polarization maps revealed by single dish telescopes (e.g., Schleuning, 1998;Houde et al., 2004;Pattle et al., 2017;Ward-Thompson et al., 2017).
Shortly after the Orion KL study, Lai et al. (2001) reported polarization observations of W51 e1 and e2 in the 1.3 mm continuum emission using BIMA. Later, Tang et al. (2009b) and Koch et al. (2018) published results toward the same source using the SMA and ALMA, respectively. W51 is a cloud complex that harbors massive star formation at various evolutionary stages (Ginsburg et al., 2015;Saral et al., 2017). W51 east hosts an active star-forming molecular clump with as many as 10 compact radio continuum sources over the 0.2 pc projected area of sources e1 and e2 (Ginsburg et al., 2016). Zhang and Ho (1997) reported inverse P-Cygni profile in the e2 core in the NH 3 (J,K) = (1,1), (2,2) and (3,3) spectral lines, consistent with infall motions of the dense gas. Lai et al. (2001) detected linear polarization in the continuum emission in the e2 and e8 cores at a resolution of 14,000 au using BIMA and found that the inferred plane-of-sky components of the magnetic fields are mostly uniform, with an average position angle of 113 • in the e2 et al. 2000). In the wind-driven model, molecular outflows are the ambient material sweptup by the wide-angle wind from a young star. The intersection of the wind and the ambient gas determines the shape and velocity of the outflow. For a constant-velocity wind with its mass-loss rate ∝sin q -2 blowing into an ambient medium with density profile ∝sin q 2 r 2 , the outflow has an approximately parabolic shape in which the velocity increases with the distance from the star (Li & Shu 1996). As a result, the shell of the parabolic outflow gives rise to ellipse or parabolic shapes (depending on the inclination angle of outflow) in channel maps and parabolic shapes in PV diagrams.
Because the observed cone-like and parabolic features in Figures 1 and 2 are similar to the structures predicted by the wind-driven model, we adopted the simplified analytical model in Lee et al. (2000) to study the driving mechanism of IRAS 4A outflows. In the cylindrical coordinate system, the structure and velocity of the outflow can be written as follows: where z is the distance along the outflow axis,R is the radial size of the outflow perpendicular to z,andC and v 0 are free parameters that describe the spatial and velocity distributions of the outflow shell, respectively. The observed features of 4A1 and 4A2 outflows are reproduced by the models of C = 3 arcsec −1 , v 0 = 5 km s −1 arcsec −1 , and an inclination of 14°(with respect to the plane of the sky) for the 4A1 outflow, and C = 2 arcsec −1 , v 0 = 3 km s −1 arcsec −1 , and an inclination of 20°for the 4A2 outflow. The predicted PV structures are shown in Figure 2. In the 4A1 and 4A2 PV diagrams, our models show that the blueshifted parabolic structures are associated with the wind-driven outflows. The displacement between the redshifted emission and the models is probably caused by the distortion of merging the 4A1 and 4A2 redshifted outflows. In the 4A1 PV diagram, the fast linear structure cannot be modeled by the wind-driven outflows. The

Frontiers
and 105 • in the e8 core. Using the DCF method, the authors estimated a magnetic field strength of 0.8 mG and 1.3 mG in the e2 and e8 cores, respectively.
More sensitive observations of continuum emission at 870 µm toward W51 using the SMA revealed a non-uniform magnetic field morphology at a higher resolution of 3300 au (Tang et al., 2009b). Tang et al. explored the possible reasons for the different distributions between the BIMA and the SMA images, finding that interferometric spatial filtering is the most likely cause. This spatial filtering by interferometers can be an advantage, as it allows us to probe magnetic fields at different spatial scales, thus revealing the dynamical role of the magnetic field across the many spatial scales relevant to the star-formation process. Figure 6 shows the magnetic field maps obtained by the SMA and later by ALMA at various spatial resolutions toward W51 e2 and e8. The figure shows a pinched, hourglass morphology in the e2 core at a resolution of 3300 au (the top-middle panel), and significantly more substructures at a resolution of 1500 au (top-right panel). Koch et al. (2018) speculate that the additional substructure in the magnetic field is the result of gravitational collapse at high densities that pulls and/or bends the field lines.
One of the key questions in the fields of magnetized low-and high-mass star formation is whether the magnetic field is dynamically important relative to turbulence and gravity. Koch et al. (2010) applied the two-point correlation function technique  and found a ratio of turbulent to magnetic energy ranging from 0.7 to 1.27 at scales ranging from ∼ 15,000 au down to ∼ 500 au. Koch et al. (2012) proposed a polarization-intensity-gradient technique that can be used to derive a distribution of magnetic field strength in a given source (see Section 1.3.1). Applying the technique to the polarization map of W51 e2, Koch et al. (2012) found a field strength of 7.7 mG. This value is in agreement with the estimate by Tang et al. (2009b), who derived the field strength from the SMA polarization observations. However, it is higher than that found by Lai et al. (2001), who used lower-resolution BIMA polarization observations probing larger spatial scales. The different values reported by Lai et al. (2001) and Koch et al. (2012) can be explained by the fact that the magnetic field strength increases toward smaller spatial scales where the gas density is higher.
Besides W51, another region of massive star formation whose magnetic field has been well studied is DR 21(OH) and its neighboring regions along the DR 21 filament (Lai et al., 2003;Cortes et al., 2005;Girart et al., 2013;Hull et al., 2014;Houde et al., 2016;Ching et al., 2017Ching et al., , 2018. DR 21(OH) was first imaged in linear polarization using BIMA in the CO (J = 2 → 1) transition and 1.3 mm continuum emission (Lai et al., 2003). These results were followed up in the CO (J = 1 → 0) transition and 3.4 mm continuum emission using BIMA (Cortes et al., 2005); in the CO (J = 3 → 2) transition and 870 µm continuum emission using the SMA Zhang et al., 2014); and in the CO (J = 2 → 1) transition and 1.3 mm continuum emission with CARMA Houde et al., 2016). Strong polarization was detected in both 1.3 mm and 870 µm continuum emission. The field strength derived using the two-point correlation function is 2.1 mG, yielding a mass-to-flux ratio of 6 × the critical value . The field strength is consistent with the value reported by Lai et al. (2003) using the DCF method. Ching et al. (2017) imaged five additional regions along the DR 21 filament using the SMA. The magnetic fields in all six cores display large dispersions, in contrast to the ordered magnetic field obtained with the JCMT at lower angular resolution (Vallée and Fiege, 2006;Matthews et al., 2009). The field strength derived from the DCF method ranges from 0.4 to 1.7 mG among the five cores, with mass-to-flux ratios from 1 to 4.3 × the critical value.
As was the case in many low-mass star formation studies (see Section 2.1), significant effort has been devoted to the search for pinched, hourglass-shaped magnetic field morphologies in massive dense cores. So far, the best examples of hourglass-shaped magnetic fields in massive dense cores are G31.41+0.31 and G240.31+0.07 Qiu et al., 2014). G31.41+0.31 is a 500 M hot molecular core with a luminosity of 3 × 10 5 L . Observed at resolutions as high as 2400 au, the dense core does not appear to fragment in the dust continuum emission. However, sensitive observations with the VLA at 1.3 and 0.7 cm wavelength revealed two compact continuum objects with a projected separation of 1300 au (Cesaroni et al., 2010). The dense core exhibits infall motions as well as rotation over a scale of 14,000 au. Magnetic fields inferred from the 870 µm continuum emission reveal a distribution that is pinched along the major axis of the flattened core . No molecular outflows have been definitively identified in this region. Observations in the CO (J = 2 → 1) transition found a velocity gradient along the major axis of the flattened core (Cesaroni et al., 2011); however, it is not clear if the gradient represents an outflow, or if it is due to core rotation. The strength of the plane-of-sky component of the magnetic field is 9.7 mG, implying a turbulence-to-magnetic-energy ratio of 0.35. The rotational velocity within the core inferred Frontiers from spectral-line observations of high-density tracers indicates significant magnetic braking. G31.41+0.31 is a case where the magnetic field dominates the turbulence and the dynamics in the system. Similar to G31.41+0.31, G240.31+0.07 is a massive star-forming region with an hourglass magnetic field morphology . As shown in Figure 7, the dust continuum emission reveals a flattened structure extended along the northeast-southwest direction that has fragmented into three cores, each harboring at least one massive young star. A wide-angle, bipolar outflow is seen in the CO emission (Qiu et al., 2009), with the outflow axis parallel to the minor axis of the flattened dense core. Polarization observations at 870 µm reveal a magnetic field topology pinched along the major axis of the core. The magnetic field strength estimated from the DCF method is 1.2 mG, with a mass-to-flux ratio of 1.2 × the critical value, and a turbulent-to-magnetic-energy ratio of 0.4. G240.31+0.07 is another clear example of a massive star-forming core in which the magnetic field dominates the turbulence and the dynamics in the system.  (Rao et al., 1998;Lai et al., 2001Lai et al., , 2003Cortes et al., , 2008Girart et al., 2009;Tang et al., 2009a,b;Beuther et al., 2010;Tang et al., 2010;Chen et al., 2012;Girart et al., 2013;Liu et al., 2013;Qiu et al., 2013;Frau et al., 2014;Hull et al., 2014;Qiu et al., 2014;Sridharan et al., 2014;Wright et al., 2014;Zhang et al., 2014;Li et al., 2015;Cortes et al., 2016;Houde et al., 2016;Ching et al., 2017;Juárez et al., 2017;Koch et al., 2018) Magnetic fields toward these sources display diverse topologies. In Sections 3.2, 3.3, 3.4, and 4 we discuss findings and interesting trends from the statistical analysis of this sample. While the analysis represents a significant advance in the study of the role of the magnetic field in high-mass star formation, the limitations and biases in the sample used in the analysis cannot be overlooked. One of the most significant limitations is that non-detections are not included in the statistics. The advent of ALMA will increase the size of the sample drastically, enabling significantly more robust analysis within the next decade.

Outflow-core magnetic field connection
Linear polarization from CO rotational transitions (i.e., from the G-K effect) probes magnetic fields in molecular gas with densities from 10 2 -10 3 cm −3 . This can be an effective tool for tracing magnetic fields in protostellar outflows. In the context of high-mass star formation, Lai et al. (2003) reported the first detection of linear polarization in the CO (J = 2 → 1) transition in DR 21(OH). The inferred magnetic field orientation is in the east-west direction, aligned with the major axis of the CO outflows. Subsequent polarimetric observations in the CO (J = 1 → 0) transition by Cortes et al. (2005) found linear polarization perpendicular to that of the CO (J = 2 → 1) transition. While this disagreement may be reconciled by anisotropy in the CO optical depth and external radiation field (see Section 1.1.2), it highlights the difficulty of interpreting spectral-line polarization from the G-K effect. Finally, Beuther et al. (2010) reported detections of spectral-line polarization in the CO (J = 3 → 2) transition in IRAS 18089-1732. Using the DCF method, they derived a magnetic field strength of 28 µG. A similar field strength of 10 µG is reported in DR 21(OH) by Cortes et al. (2005).
Despite the early success in detecting spectral-line polarization primarily in DR 21(OH), very few protostellar sources have shown detectable spectral-line linear polarization since those early results. In the survey of 14 high-mass protostellar clumps by Zhang et al. (2014), only DR 21(OH) had detectable polarization in the CO (J = 3 → 2) transition . The limited sensitivity of the SMA may be a contributing factor to these non-detections. Furthermore, the CO emission is typically spatially extended, which presents an additional challenge when imaging linear polarization, since Stokes Q and U can be either positive or negative, and a lack of short-spacing information in sparsely sampled interferometric data can produce negative emission that may be confused with the polarization signal in Stokes Q and U . Both factors are mitigated by ALMA, which provides significant improvements in sensitivity and imaging fidelity over BIMA, CARMA, and the SMA.
Probing magnetic fields in accretion disks around high-mass protostars is challenging in two regards. First, extremely high angular resolution observations are required to achieve the ∼ 100 au linear scales of a disk at a typical source distance of several kpc. Second, at the typical high densities and optical depths in these disks, the polarized emission may be contaminated (and sometimes dominated) by self-scattering of large dust grains (see Section 5). For example, Girart et al. (2018) observed polarized emission from the disk associated with the radio jet HH 80-81. These deep ALMA observations in the 1.14 mm continuum emission spatially resolved a flattened disk with a radius of 291 au. The disk is partially optically thick and its polarized emission appears to be dominated by self-scattering of large dust grains.
In the absence of magnetic field information in protostellar outflows and disks, the comparison of magnetic fields in cores with outflow axes offers an alternative to assess the importance of magnetic fields in the formation of disks and outflows in high-mass systems. As discussed in Section 2.1, when collapsing to form stars, strongly magnetized cores are expected to form a disk and outflow system with the major axis of the outflow parallel to the magnetic field. Zhang et al. (2014) compared outflow axes with the plane-of-the-sky component of magnetic fields in cores from which the outflows originate. They found a slight preference around 0 • and 90 • in angles between the magnetic field and the outflow axis. However, due to the small sample size, the data are consistent with a random orientation of magnetic fields and outflows. This lack of correlation, if confirmed by more statistically significant samples, suggests that angular momentum in accretion disks is not dictated by the orientation of the magnetic field in the cores. The dynamical importance of the magnetic field thus appears to weaken relative to gravity and angular momentum from the core to the disk scales.

Fragmentation and cluster formation
The clustering nature of massive stars implies that there must be considerable fragmentation in parsecscale molecular clumps in order to produce a cluster of stars. Gravity, turbulence, magnetic fields, and stellar feedback all influence fragmentation and the outcome of cluster formation. The dynamical role of the magnetic field during the fragmentation of molecular clumps can be assessed if the magnetic field strength is measured directly. As mentioned in Section 1.2, this can be achieved by observing circular polarization from the Zeeman effect. While observations of the Zeeman effect have been carried in HI and OH line emission using the VLA (Crutcher, 2012;Crutcher and Kemball, 2019), these observations probe a low-density medium that may not be directly involved in gravitational collapse. There are no reported interferometric observations of the Zeeman effect in dense molecular gas, although ALMA is likely to reshape this field once precise measurements of circular polarization become available to the user community. Future interferometric observations of the Zeeman effect (measuring the line-of-sight magnetic field), when combined with linear dust and spectral-line polarization observations (measuring the plane-of-sky field), will allow much more robust estimates of the total magnetic field strengths in star-forming regions.
A powerful indirect method that can be used to assess the dynamical role of magnetic fields in protocluster formation is the analysis of the distribution of magnetic fields within cluster-forming molecular clumps. Numerical simulations of turbulent, magnetized molecular clouds offer clues about magnetic field topologies in strong and weak field regimes. When the magnetic field is strong relative to turbulence, the field is less disturbed and appears to be ordered (e.g. Ostriker et al., 2001).
Observationally, such a study becomes meaningful only when a statistically significant sample is involved. Zhang et al. (2014) compared dust polarization in dense cores probed by the SMA with the polarization in the parental molecular clumps observed by single-dish telescopes in a sample of 14 high-mass star forming regions; the results show a bimodal distribution in polarization angles. As shown in Figure 8, magnetic fields on dense-core scales are not randomly distributed, but are either parallel or perpendicular to the field orientations in their parental clumps. A later study of a larger sample of 50 primarily high-mass sources by Koch et al. (2014) compared the magnetic field orientation with the gradient of the total dust emission and came to the same conclusion. These findings indicate that the magnetic fields are dynamically important in cluster-forming clumps, and that the field is strong enough on the clump scale to channel the material along the field lines into dense cores during the gravitational collapse. While Hull et al. (2014) found hints of consistency in the magnetic field orientation from ∼ 0.1 pc to ∼ 1000 au scales in a few low-mass sources (see Section 2.1), the results at the ∼ 1 -0.1 pc scales from Koch et al. (2014) and Zhang et al. (2014) are much more significant, suggesting that the magnetic field may be more dynamically important at parsec scales.

Virial equilibrium in massive cores and cluster formation
An essential component of the classical view of star formation is that dense cores should be in a state of an approximate virial equilibrium (Larson, 1981;Shu et al., 1987;McKee and Tan, 2003). However, there is increasing evidence from recent observations that cores forming high-mass protostars may not be in virial equilibrium. Pillai et al. (2011) carried out a stability analysis in two infrared dark clouds (IRDCs) and found that the virial parameter α vir , defined as the ratio between the virial mass M Virial and the gas mass M gas , is typically between 0.1 to 0.3. These results are corroborated by Kauffmann et al. (2013), who compiled a large sample of massive clumps from surveys and found that a significant fraction of them have virial parameters 2, where 2 is the value expected from a pressure confined, self-gravitating Bonnor-Ebert sphere. More recently, Traficante et al. (2018) pointed out that the virial mass can be underestimated because the spectral line emission of a tracer molecule preferentially probes sections of a cloud with gas densities above the critical density of the line transition. This effect may result in an underestimate of the observed line width, and hence lead to a small virial parameter. However, since the effective excitation density of a line transition is typically lower than the critical density due to radiative trapping of photons when the optical depth is non-negligible, this effect may not be as significant. For example, many of the virial analyses in the literature use the NH 3 and N 2 H + lines, which often exhibit consistent line widths despite having critical densities that differ by more than a factor of 5.
These small virial parameters appear to challenge the picture of star formation in which gas evolves in a state of equilibrium. However, the virial analyses discussed above do not include contributions from magnetic fields, which may provide significant support in dense cores. Frau et al. (2014) carried out continuum polarization measurements of NGC 7538 with the SMA and found a magnetic field strength of 2.3-2.6 mG in the region. Frau et al. also performed a detailed energetics comparison of the gravitational potential energy, turbulent support, thermal pressure, and magnetic support. Among the 13 dense cores analyzed, the magnetic support amounts to 0.2 to 2.4 × the combined effect from turbulent and thermal pressure. More than half (eight out of 14) of the cores have magnetic support comparable to the turbulent and thermal support. In addition, the virial parameters including the magnetic support vary from 0.1 in the central region where star formation takes place to 8 in cores that are influenced significantly by molecular outflows.
The studies of NGC 7538 by Frau et al. (2014) and Wright et al. (2014) demonstrate that magnetic fields can indeed be a significant contributor to the support of dense cores. Such an analysis should be extended to high-mass star-forming regions at earlier evolutionary stages when stellar feedback has not significantly altered the initial physical conditions. To that end, Pillai et al. (2011);Zhang et al. (2015); Ohashi et al. (2016); Sanhueza et al. (2017);and Lu et al. (2018) performed virial analyses in high-mass star forming regions prior to the development of strong HII regions. The interferometric observations allowed the identification of structures < 0.1 pc in size within massive clumps. Figure 9 presents a comparison of the virial mass and the gas mass for structures identified in 12 high-mass star forming regions (Pillai et al., 2011;Zhang et al., 2015;Ohashi et al., 2016;Lu et al., 2018). The virial mass is computed as where σ v is the line-of-sight velocity dispersion due to both non-thermal and thermal broadening, R is the radius, and G the gravitational constant. k = 5−2a 3−a is a correction factor related to the density profile ρ ∝ r −a . We assume a constant density in the structure, i.e., a = 0, which leads to the maximum value in the correction factor (k = 5/3) and thus the maximum virial mass. The gas mass is derived from the dust continuum emission (see Section 1.4). Figure 9 reveals that the dense gas structures in these high-mass star forming regions have virial parameters < 2. Furthermore, a large fraction of the dense gas structures have virial parameters 2. Since there are no direct measurements of magnetic fields for this sample, we compute the magnetic field required to increase the virial parameter from 0.5 to 1.0. The magnetic virial mass (Hennebelle and Chabrier, 2008) is computed as where σ A = B √ 4πρ is the Alfvénic velocity corresponding to a magnetic field strength B and density ρ. We compute the magnetic field strengths when M Mag = M Virial = 0.5M gas . Assuming a representative gas mass M gas and radius R for the clumps (M gas = 1000 M , R = 0.5 pc), cores (M gas = 50 M , R = 0.05 pc), and condensations (M gas = 5 M , R = 0.005 pc), we find required magnetic field strengths of 0.29 mG, 1.46 mG and 14.4 mG, respectively. Under these conditions, the virial parameters would be 0.5 without the contribution of magnetic fields, and would increase to 1 after the inclusion of magnetic fields. The required field strengths are in broad agreement with typical literature values of magnetic fields that were derived using the DCF method based on polarization observations of massive star-forming regions.
Recent measurements of the thermodynamic properties in high-mass star-forming regions reveal small virial parameters that appear to challenge the assumption of equilibrium star formation. However, a lack of magnetic field measurements in these same regions leaves open the possibility of virialized star formation, since the field strength of a fraction of a mG to several mG could bring the dense gas close to a state of equilibrium. With the advent of ALMA, we expect significant progress to be made on this vital question in high-mass star formation as more observations of both spectral lines and polarization are carried out over the coming years. The data were taken from Pillai et al. (2011);Zhang et al. (2015); Ohashi et al. (2016);and Lu et al. (2018). The virial mass includes contributions from thermal and non-thermal line widths. The gas mass was derived from dust continuum emission. The dashed, solid, and dotted lines denote virial parameters α vir = M Virial /M gas of 2, 1 and 0.1, respectively. The majority of the data points have α vir < 1, indicating that clumps and cores are far from virial equilibrium when contributions from magnetic fields are not included. The red arrows indicate the strength of the magnetic field (in mG) that are required to increase α vir from 0.5 to 1.0 (see Section 3.4).
A pinched/hourglass-shaped magnetic field configuration can yield a range of magnetic-field morphologies when projected in the plane of the sky (Frau et al., 2011). When the main axis of the system is aligned with the line of sight, the plane-of-sky component of the field is oriented radially outward. This scenario was proposed by Girart et al. (2013) to explain the magnetic field topology in DR 21(OH). Among the low-and high-mass objects with robust detections of polarization (and thus inferred magnetic fields), 28% of the low-mass sources and 21% of the high-mass sources exhibit hourglass configurations. Assuming that the sample is randomly oriented in three dimensions, the probability than an hourglass axis will be oriented within 30 • of the line of sight is ∼ 13% (Frau et al. 2011 show that an angle closer to the line of sight than ∼ 30 • yields a mostly radial pattern). Therefore, the majority of the objects (∼ 87%) should display an hourglass shape if the variations in magnetic field morphologies were solely due to projection effects.
The fact that only a small fraction of the objects exhibit hourglass morphologies suggests that the scenario of magnetically dominant core collapse is not the predominant mode of low-or high-mass star formation. However, we should note that when taking into account the projection effect, the fraction of the detected hourglass reaches 28%/87% ≈ 1/3 for low-mass and 21%/87% ≈ 1/4 for high-mass sources. Such fractions are not negligible, and are even more significant considering the various environmental and dynamical differences between sources, as well as the observational limitations that may hinder the detection of an hourglass-in particular the sensitivity and dynamic range limitations of the pre-ALMA interferometers.
Of the sources observed by ALMA that clearly lack hourglass magnetic field morphologies, only one so far (Ser-emb 8 in Hull et al. 2017a; see Figure 3) shows a truly chaotic morphology, presumably having been influenced by turbulence and infall. However, the few ALMA polarization observations published to date have shown more examples of magnetic fields that seem to have been shaped by bipolar outflows. These observations are challenging to interpret, as magnetic fields aligned with the walls of an outflow cavity can look quite similar to an hourglass when observed with low resolution. Thus, future studies will need to carefully consider how to determine which "hourglasses" are preserved from the natal collapse phase, such as in NGC 1333-IRAS 4A, where at ∼ 1000 au scales the hourglass appears to be unperturbed by the two outflows emanating at different angles from the embedded binary (Girart et al., 2006); versus which are caused by outflow/jet feedback, such as in Serpens SMM1 (Hull et al., 2017b) and B335 .
We further caution that the above values are lower limits that are limited by telescope resolution: i.e., in the high-mass case, we consider each source as unique, despite the fact that each clump is likely to have more than one protostellar object. This is even true in the low-mass case: for example, SMA observations were able to distinguish the magnetic fields in SVS 13A and B, whereas CARMA observations were not; and ALMA observations were able to map the magnetic fields in Serpens SMM1-a and b separately, whereas CARMA and SMA observations were not. Finally, and perhaps most important, the number of objects observed with interferometers remains small: observations by ALMA are likely to increase the sample size significantly in the near future, enabling studies that will establish far better statistical constraints.

POLARIZATION FROM DUST SCATTERING
While polarized emission from dust scattering is independent of the magnetic field, we briefly mention it here because this young, quickly growing subfield emerged from studies that were attempting to make resolved maps of the magnetic field in protoplanetary disks, which is one of the longstanding goals of diskand star-formation studies. Rao et al. (2014) reported a tentative detection of a toroidally wrapped magnetic field in the Class 0 protostar IRAS 16293B at ∼ 75 au resolution using the SMA. However, on the whole, pre-ALMA full polarization observations of disks did not have the sensitivity or the resolution necessary to make well resolved maps of low-level ( 1%) polarized dust emission (Hughes et al., 2009(Hughes et al., , 2013. Around the same time that ALMA started producing polarization results, several theoretical studies suggested that polarized (sub)millimeter-wave emission from protoplanetary disks could be caused by the self-scattering of emission from (sub)millimeter-sized dust grains (Kataoka et al., 2015(Kataoka et al., , 2016aPohl et al., 2016;Yang et al., 2016a,b), consistent with previous work by Cho and Lazarian (2007). Kataoka et al. (2016a) and Yang et al. (2016b) explained the 1.3 mm CARMA polarization observations by Stephens et al. (2014) of the Class I/II source HL Tau in the context of this self-scattering mechanism. Self-scattering can also explain the polarization pattern observed in the Class II protoplanetary disk IM Lup in 870 µm ALMA observations by Hull et al. (2018b), and similar observations of HD 163296 by Dent et al. (2019).
The situation became even more complex when Tazaki et al. (2017) proposed yet another mechanism that can cause polarization in disks: dust grains whose minor axes are aligned with the gradient in the dust emission (this work is rooted in the radiative torque model by Lazarian 2007). This mechanism, which is sometimes referred to as "k-RATs" (i.e., where k is the orientation of the gradient in the dust emission; see also Footnote 3), is different from self-scattering by dust grains. However, similar to self-scattering, it is independent of any magnetic field that might be in the disk. Polarization from k-RATs, which has a morphology that is azimuthal, is broadly consistent with 3 mm ALMA observations of polarization in the HL Tau disk (Kataoka et al., 2017), although Yang et al. (2019) found that their k-RAT model did not reproduce the data when examined in detail. The transition from possible k-RAT alignment at long (3 mm) wavelengths to scattering at short (870 µm) wavelengths was shown clearly by Stephens et al. (2017b), who reported high-resolution ALMA polarization results at 870 µm and 1.3 mm, complementing the 3 mm observations reported by Kataoka et al. (2017). The intermediate-wavelength 1.3 mm data exhibit a morphology consistent with roughly equal contributions of self-scattering and k-RATs.
Other recent results have interpretations that are not as clear. These include observations of the edge-on HH 111 and HH 212 disks , the high-mass protostellar disk Cepheus A HW2 (Fernández-López et al., 2016), and the transition disk HD 142527 (Kataoka et al., 2016b;Ohashi et al., 2018); the polarized emission from these objects may be caused by either scattering or magnetically aligned grains (B-RATs). Observations of the low-mass disks CW Tau and DG Tau (Bacciotti et al., 2018) and the massive GGD27 MM1 disk  associated with the iconic HH 80-81 radio jet show polarization that may be caused by dust self-scattering and/or k-RATs. Polarization observations of a number of embedded Class 0 protostars by Segura- Cox et al. (2015); Cox et al. (2018); Harris et al. (2018); Sadavoy et al. (2018a,b) and Takahashi et al. (2018) sometimes show signs of scattering in the inner regions of sources and magnetic alignment in the outer regions. Finally, Alves et al. (2018) reported ∼ 30 au-resolution observations of the Class I binary source BHB07-11, which exhibits consistent polarization patterns at 3 mm, 1.3 mm, and 850 µm. They argue that the polarized emission most likely arises from magnetically aligned dust grains.
We emphasize that studies of (sub)millimeter polarization from disks are in their infancy. More observations and more mature, detailed models are necessary before we will be able to determine whether there is, in fact, any contribution from magnetic fields in the polarized emission from protostellar disks.

SUMMARY
The steady improvement in telescope sensitivity at the (sub)millimeter-wave bands where dust continuum emission is strong has led to a revolution in interferometric observations of linear polarization over the past two decades. In recent years, both CARMA and the SMA conducted surveys of magnetic fields in samples of > 10 star-forming cores/clumps. These studies have proven insightful in revealing the statistical behavior of magnetic fields in star formation. However, while individual case studies continue to be fruitful, large, less biased surveys that target statistically significant samples of sources are the only way to further constrain the dynamical role of magnetic fields in the star formation process.
Studies of low-mass star formation have found that while there are a few cases where the magnetic field seems to regulate the collapse of star-forming material across the scales from clouds/clumps to cores to envelopes, there are also many cases where the multi-scale magnetic field morphology shows no consistency. Furthermore, the resolution and sensitivity of ALMA have yielded new observations that show evidence of shaping of the magnetic field by bipolar outflows at the ∼ 100 au scale in some (but not all) sources. A few recent observations of spectral-line polarization by the SMA and ALMA have shed light on the importance of magnetic fields in the launching of outflows. More observations of this type with ALMA will enhance this subfield significantly in the coming years.
Studies of high-mass star formation found that magnetic fields are not the dominant force within the parsec scale of molecular clumps. However, fields appear to be dynamically important in the fragmentation of clumps into 0.1 pc-scale cores. Furthermore, molecular gas with densities > 10 4 cm −3 in high-mass star formation appears to be far from virial equilibrium if magnetic fields do not contribute significantly to the dynamical process. This significant departure from the state of virialization appears to challenge the basic assumption of equilibrium star formation. Direct measurements of magnetic fields will shed more insight on this important issue.
In both the low-and high-mass regimes, comparisons of outflow versus magnetic field orientations have yielded random alignment. The outflow is a probe of the angular momentum at the smallest scales in the source, and thus the limited data currently available point to a scenario where angular momentum is far more important than the magnetic field at the small (< 1000 au) scales of individual forming protostars.
The role of magnetic fields in the formation of high-and low-mass disks is less clear due to a small number of observations, and the issue of scattering by dust. Near-future studies targeting the Zeeman and G-K effects may finally be able to access information about the magnetic field in disks (e.g., Brauer et al., 2017). Furthermore, high resolution studies at longer wavelengths in regions that are optically thick (and thus dominated by scattering) in the ALMA data will be made possible by future telescopes such as the Next Generation Very Large Array (ngVLA; Carilli et al. 2015). The ngVLA will be a powerful tool for studying magnetized star formation Hull et al., 2018a), as it will offer dramatic improvements in the sensitivity, resolution, and image fidelity over the current VLA, which has enabled just a few studies of polarization in the very brightest low-and high-mass star-forming sources (Carrasco-González et al., 2010;Cox et al., 2015;Liu et al., 2016Liu et al., , 2018. Despite major advances in recent years, the studies presented in this review-in particular the survey studies-are biased toward objects with detectable polarization. Nearly all low-and high-mass star forming objects have embedded protostars (and sometimes HII regions, in the case of high-mass sources). A crucial, under-explored area of star formation involves studying the initial conditions of the magnetic field before feedback (e.g., infall, outflows, and radiation) has altered its morphology. However, this has rarely been achieved due to a lack of sensitivity. In the coming years, large, multi-scale, multi-telescope surveys of magnetic fields in star-forming clouds covering the full range of evolutionary stages will allow us to break new ground in our studies of magnetized star formation.

AUTHOR CONTRIBUTIONS
CH led the writing of the sections about low-mass star formation and scattering. QZ led the section about high-mass star formation. Both authors wrote the abstract, introduction, hourglass-statistics section, and summary.

FUNDING
CH acknowledges the support of both the NAOJ Fellowship as well as JSPS KAKENHI grant 18K13586. QZ acknowledges the support of the the Smithsonian Institute SSA grant, "Are Magnetic Fields Dynamically Important in Massive Star Formation?" ACKNOWLEDGMENTS CH and QZ acknowledge the two anonymous referees for their careful reading of the manuscript and for the insightful comments, which led to a more clear and thorough presentation of this review. CH acknowledges Paulo Cortés, Martin Houde, and Josep Miquel Girart for the useful discussion. The BIMA array was operated by the Berkeley-Illinois-Maryland Association with support from the National Science Foundation. Support for CARMA construction was derived from the states of California, Illinois, and Maryland, the James S. McDonnell Foundation, the Gordon and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the University of Chicago, the Associates of the California Institute of Technology, and the National Science Foundation. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. Software: APLpy, an open-source plotting package for Python hosted at aplpy.github.com (Robitaille and Bressert, 2012). Common Astronomy Software Applications (CASA, McMullin et al. (2007)). Astropy (Astropy Collaboration et al., 2018).