The role of outflows in the formation of stars and the protostellar disks that generate them is a central question in astrophysics. Outflows are associated with star formation across the entire stellar mass spectrum. In this review, we describe the observational, theoretical, and computational advances on magnetized outflows, and their role in the formation of disks and stars of all masses in turbulent, magnetized clouds. The ability of torques exerted on disks by magnetized winds to efficiently extract and transport disk angular momentum was developed in early theoretical models and confirmed by a variety of numerical simulations. The recent high resolution Atacama Large Millimeter Array (ALMA) observations of disks and outflows now confirm several key aspects of these ideas, e.g., that jets rotate and originate from large regions of their underlying disks. New insights on accretion disk physics show that magneto-rotational instability (MRI) turbulence is strongly damped, leaving magnetized disk winds as the dominant mechanism for transporting disk angular momentum. This has major consequences for star formation, as well as planet formation. Outflows also play an important role in feedback processes particularly in the birth of low mass stars and cluster formation. Despite being almost certainly fundamental to their production and focusing, magnetic fields in outflows in protostellar systems, and even in the disks, are notoriously difficult to measure. Most methods are indirect and lack precision, as for example, when using optical/near-infrared line ratios. Moreover, in those rare cases where direct measurements are possible—where synchrotron radiation is observed, one has to be very careful in interpreting derived values. Here we also explore what is known about magnetic fields from observations, and take a forward look to the time when facilities such as SPIRou and the SKA are in routine operation.

1. Introduction

Of the many important roles that magnetic fields play in the formation of stars, perhaps none is more dramatic nor as full of consequence as is the launch and collimation of powerful outflows. It is now nearly 70 years since the first evidence for energetic outflows in star formation regions was discovered (Herbig, 1951; Haro, 1952) although not realized as such. In fact it would be many years later before these nebulous clouds (now known as Herbig-Haro objects) were correctly identified as radiative shocks driven by an outflow from a young star (Schwartz, 1977). The availability of new technologies and telescopes over the subsequent decades has yielded discoveries of bipolar molecular outflows at millimeter wavelengths, high speed (hundreds of km s−1) optical jets in permitted and forbidden line emission, and more recently even synchrotron emitting jets associated with low mass star formation. This extensive body of work shows that bulk flows ranging from several to hundreds of km s−1 are observed in these different tracers. Although first discovered to be associated with the birth of low mass T-Tauri stars, observations over the last two decades have shown outflows to be linked with the formation of objects across the entire mass spectrum: from brown dwarfs (Whelan et al., 2005) on up to O stars (Caratti o Garatti et al., 2017). These observations have made it clear that outflows are an essential component of star formation.

The first clues to the origin of protostellar jets emerged from studies of outflows from low mass stars, whose radiation pressures fail by orders of magnitude to drive them at the observed thrusts (see for example Wu et al., 2004; Vaidya et al., 2011 on the decollimating effects of radiation in the case of massive young stars). It was therefore natural to consider the possibility that T Tauri jets and outflows could be magnetized winds from rotating bodies. Two kinds of magnetized rotors are possible—rapidly spinning, magnetized protostars or magnetized accretions disks out of which all stars and their planetary systems form. While magnetic field measurements needed to confirm such a picture have been long in coming, direct observational evidence is now available using several approaches: the spectro-polarimetry of jet sources (Donati et al., 2010), synchrotron emission in outflows from some young stars and most recently the detection of polarization of the SiO line in a jet from a low mass star arising from the Goldreich-Kylafis effect (Lee et al., 2018). The observed fields range in strength from kilogauss close to the star to a fraction of a milligauss in distant parts for the outflow. Such values suggest that while magnetic forces dominate in the vicinity of the young stellar object (YSO), they are no longer dynamically important further out. This scenario is in agreement with numerical simulations (Hartigan et al., 2007).

The basic physics of how magnetized rotating stars drive winds and undergo magnetic braking as a consequence was developed by Chandrasekhar (1956) and Mestel (1961) before being applied by Weber and Davis (1967) who showed that the magnetized solar wind would spin down the Sun. Magnetic field lines threading a rotating star can enforce the co-rotation of gas out to some distance along the field line—the Alfvén radius r A where the outflow speed equals the speed of a transverse Alfvén wave ( B / 4 π ρ ). In so doing, the wind is able to extract the rotor's angular momentum which is then carried out by the rotating, accelerated outflow. The lever arm of the wind torque is, in effect, the Alfvén radius. It can extend to significant distances from the body making the magnetized outflow a highly efficient mechanism for extracting angular momentum. It is this point, more than any other, that makes magnetic fields so important in the physics of protostellar outflows and underpins their role in star formation. Moreover, the combination of magnetic fields and outflows, or more precisely winds, continue to be important in subsequent evolutionary phases as a means of removing angular momentum even when a star is on the Main Sequence (Reiners and Mohanty, 2012).

Linking magnetized winds to the physics of accretion disks was first suggested in the context of accreting black holes at the centers of radio galaxies (Blandford and Payne, 1982). It became clear that magnetized disk-wind torques, exerted by fields threading vertically through the disk and bending outwards beyond, extract some, if not most of the angular momentum from each annulus of the disk. This, in turn, could drive an accretion flow onto the central black hole. When applied to circumstellar disks around young stars, magnetized disk winds (D-winds), extending over most of the disk surface, could also provide the major driver of accretion onto the central forming star (Pudritz and Norman, 1983). In any event this theory led rather quickly to the prediction that outflows and accretion disks must be fundamentally linked in a structure that undergoes both accretion and ejection. In particular the ratio of the mass outflow rate to the wind driven accretion rate depends on the lever arm of the flow: the ratio of r A to the footpoint radius r o of field line at each disk radius. Observations of a wide variety of systems typically observe a wind mass flux to accretion flux ratio, Ṁ w /Ṁ a ≃ 0.1 (Watson et al., 2016) confirming theoretical predictions (Pudritz and Norman, 1986). Another key prediction of the theory is that jets rotate—a consequence of the fact that they carry off the disks's angular momentum. These predictions have been confirmed by a large number of 2 and 3D hydromagnetic simulations using a wide variety of initial setups (isolated cores, turbulent clumps, etc.) and numerical codes (ZEUS, FLASH AMR, RAMSES AMR, SPH MHD, etc.). A small sampling of these works starts from Shibata and Uchida (1985), Ouyed et al. (1997), Krasnopolsky et al. (1999), Banerjee and Pudritz (2006), Hennebelle and Ciardi (2009), Price et al. (2012), Stepanovs and Fendt (2014), Staff et al. (2015), and Tomida et al. (2015), to the appearance of high resolution zoom-in techniques in Kuffmeier et al. (2017).

While D-wind models are attractive, an alternative mechanism for generating outflows from disks was proposed in a series of papers starting with Shu et al. (1994) and which are comprehensively reviewed in Shu et al. (2000). According to the latter model, high velocity outflows arise not from the surface of the disk but instead from the narrow annulus where the star's magnetosphere interacts directly with the inner edge of the disk, at the so-called co-rotation radius . Here the magnetic field of the young star is assumed to be strong and to clear out disk matter as far as the co-rotation radius. At this radius the magnetic field switches from its inner closed configuration (associated with the star's magnetosphere) to an open configuration. This constitutes a magnetic X-point or more precisely an X ring. Matter is then launched centrifugally from the disk along the open field lines in a flow known as an X-wind. It should be emphasized that both X-winds and D-winds are driven magneto-centrifugally along open field embedded within rotating disks. One of the key physical differences between the two models is that the X-wind only removes most of the Keplerian angular momentum from near the disk's inner edge defined by the magnetopause (where the magnetic field of the star truncates the disk). In contrast, the D-wind removes angular momentum from every disk radius. The D-wind model also allows the Keplerian angular momentum at the disk inner edge to be transported to the star via gas infall along the stellar field lines. Several mechanisms have been proposed to remove this from the star, including coupling to the disk, or by an accretion powered stellar wind. In terms of magnetic field geometries, the main differences (Shang et al., 2007) are where the field lines are anchored: near the co-rotation radius, the X-radius, for X-winds—or over a wide range of radii for D-winds starting from the inner disk edge.

One of the basic questions regarding the application of such magneto-hydrodynamic (MHD) wind models to protostellar disks is the degree to which magnetic fields are really coupled to the disks, given that the very low ionization of high column density disks implies that non-ideal effects are important (Königl, 2009). A great deal of effort has been expended in developing computer codes capable of simulating these effects over the last decade (Duffin and Pudritz, 2009; Seifried et al., 2012; Bai and Stone, 2014; Gressel et al., 2015; Hennebelle et al., 2016). The results of this work have now challenged the traditional idea of accretion disks as systems that are dominated by turbulent angular momentum transport, as we shall see.

At the time bipolar outflows were discovered, protostellar disks had not yet been directly detected although their presence had been inferred from the modeling of observed spectral energy distributions using the IRAS satellite (Rucinski, 1985). Moreover, very little was known about the spins of young stars, except for the fact that a few rotate very slowly (periods of approximately a week). The first optical images of protostellar disks, seen against the bright emission of the Orion Nebula Cluster by the Hubble Space Telescope (O'Dell and Wen, 1994), showed them to be at most 1″ in size for the nearest star formation regions. The search for disks associated with outflows, having the signatures of Keplerian rotation, required both high spatial and spectral resolution. Such capability became available with mm interferometry and an early-known jet source, HL Tau, was shown to possess a rotating circumstellar disk, perpendicular to the jet, using CO as a tracer (Sargent and Beckwith, 1987). Such discoveries, in fact became one of the major science drivers for ALMA, i.e., to discover and map protostellar disks and undertake studies of both star and planet formation within them. The connection of outflows with disks is now being directly addressed as part of the “ALMA revolution” wherein disks around young Class I stars, and their associated outflows (observed at millimeter wavelengths), have been resolved down to scales of a few au in nearby sources. Recent efforts have also focused on studying the formation and growth of disks in the most deeply embedded, even earlier phases of star formation (Class 0 sources), see for example Tobin et al. (2013).

While it is generally acknowledged that magnetic fields likely play a central role in observed protostellar outflows, until recently the possibility that outflows might also control the physics of accretion disks has been largely ignored. One reason is observational. Magnetic fields in jets and disks are notoriously difficult to measure and the needed instruments have not been available. A second reason is theoretical. The seminal papers on accretion disk theory assumed that turbulence in disks would result in a viscous torque. This transports angular momentum radially outwards through the disk, driving most of the disk material to accrete onto the star as the outer disk spreads out to carry off the angular momentum, see for example Armitage (1998). The questions here are—what is the origin of turbulence? Is it of sufficient amplitude to support the large observed disk accretion rates? It took nearly two decades for an explanation to emerge. The Magneto-Rotational Instability (MRI), arising from weak magnetic fields, provided a sound physical basis for turbulence in perfectly conducting, so-called “ideal,” disks (Balbus and Hawley, 1998). However, recent ALMA observations show that turbulence is not strong enough to carry away angular momentum in observed disks, at least at tens of au from a young star (Flaherty et al., 2015, 2017, 2018). While the observations are consistent with numerical simulations for the mid-plane, where MRI-driven turbulence will be strongly damped, they are inconsistent with expected turbulent velocities at the disk surface, v turb ≃ 0.1 − 1.0c s , required to drive the observed accretion rates (Simon et al., 2018). At the same time, numerical simulations of protostellar disks, that account for important non-ideal MHD effects, have discovered that the MRI is strongly suppressed in the expected dense, poorly ionized gas. It gives way instead to magnetized disk winds (Bai and Stone, 2014; Gressel et al., 2015). While further developments are needed—such as including the effects of grain evolution on non-ideal MHD process—progress in the field has been enormous. Thus, recent advances in observations and breakthroughs in difficult numerical simulations both seem to point toward the primacy of magnetized outflows rather than disk turbulence as the driver of accretion disk physics (Simon et al., 2018). Before exploring disk winds further, it should be pointed out that alternative ways of transporting angular momentum outwards in an accretion disk have been suggested. In particular, transport can be by spiral density waves although such waves only form in disks that are a significant fraction of the mass of the central star (Kratter and Lodato, 2016). Moreover embedded masses, e.g., newly formed massive planets, can drain angular momentum away from a disk and support accretion onto the young star; this however does not seem to be an important mechanism (Armitage, 2011).

There is also the question of what causes the surprisingly low angular momentum of many young stars. Early work suggested that one effect of the star's magnetosphere interacting with the disk is to spin-down the star, through magnetic braking, so that it is effectively locked to the same rotation period as the inner edge of the accretion disk (Koenigl, 1991). This would explain why young stars, which show evidence of disks, rotate slower on average than their disk-less counterparts of comparable age and mass (Cieza and Baliber, 2007). More recent work suggests that stellar spin may be controlled by magnetized stellar outflows that originate from the protostars themselves. Here, the accretion of material from the inner disk flows onto the star by falling down magnetospheric field lines. The gravitational energy released in this infall impacts the foot points of these field lines. This ultimately gives rise to a powerful flux of Alfvén waves that heats the stellar corona giving rise to an accretion-powered wind (see for example Matt et al., 2012). As already mentioned, outflows have also been suggested to originate at the disk-magnetosphere boundary—the so-called X-wind—which are proposed to carry off the angular momentum of the inner disk before matter is accreted onto the star (Shu et al., 2000). In any event, as soon as the star loses its disk, it starts to spin-up as it radius shrinks while descending its Hayashi track. Observational advances in measuring both the magnetic field and spins of young stellar objects have opened up the possibility of formulating the complete spin history of a star, the evolution of which appears to be particularly important in the early phases of star formation (Bouvier et al., 2014).

Magnetically driven outflows have a large number of other important consequences for star formation. If outflows originate from disks then the process of disk formation during gravitational collapse in a magnetized medium must be deeply connected with magnetic braking in the early collapse phase and subsequent launching of MHD outflows as the disk, or first core, starts to materialize. Disks and outflows in this view are inseparable, both being tied to collapse in turbulent, clouds. Indeed, as has long been known, one of the central features of the observations is that outflows are most powerful during the most deeply embedded Class 0 stage when at least half the mass of the system is still in the envelope which is collapsing into the disk (Jørgensen et al., 2009). Recent analytic and numerical calculations indeed show that small disks of about 20 au in size should form in magnetized collapse wherein the magnetic field is not frozen into the gas but undergoes (ambipolar) diffusion (Hennebelle et al., 2016). This is very much in keeping with the sense of recent IRAM-PdBI observations which find that disks are typically smaller than expected from purely hydrodynamic models (Maury et al., 2019). Thus, magnetic braking and the launch of outflows early in the gravitational collapse and disk formation stages of star formation, respectively, control one of the most fundamental physical properties of the initial state of a disk—its radius.

Feedback processes arising from outflows have long been thought to play an important role in defining the mass of a star and the efficiency of star formation. Magnetized outflows carry significant amounts of angular momentum and thrust, and are powered by the gravitational potential energy released during the collapse. In this way, outflows can act as an important, and even dominant form of protostellar feedback during star formation. This feedback can cut off the supply of infalling gas to the disk, and thereby help to limit or determine stellar mass (Matzner and McKee, 2000). On physical scales beyond the molecular core radius (typically ≈ 0.04 pc), protostellar outflows could stir up the surrounding molecular cloud and drive turbulence. This would, to some degree, stave off the formation of too much dense, star forming gas as first suggested in the pioneering paper of Norman and Silk (1980), and is in agreement with many current MHD simulations (Federrath, 2016). Feedback from protostellar outflows would then help to reduce star formation rates and efficiencies in molecular clouds addressing one of the big questions of star formation—why is the process so inefficient?

This review examines these issues and questions from a modern perspective. We focus on recent rapid progress in the ALMA era that has been made in the physics and observations of outflows and their role in star formation. ALMA has provided a radically new capability to study disks and their outflows. Most of our attention is on observational and theoretical results that have arisen since 2014 post the publication of Protostars and Planets VI. We first review the observational advances in measuring magnetic fields in stars, jets, and disks, and the properties of outflows (shocks, rotation, etc.). We then go on to examine the theoretical and computational results on how outflows are launched, and their significance in various aspects of star formation. These latter angular momentum related topics include outflows and their effects on—gravitational collapse, disk formation, disk evolution, and protostellar spin. We then discuss the feedback properties of outflows including their effects on stellar mass and molecular cloud turbulence. We address these issues in the context of both low and high mass star formation. The latter subject brings up the difficult question of how powerful protostellar radiation combines with magnetic outflows to determine conditions in massive star forming regions. Finally, since stars and planets form in disks, both must be affected by outflow physics and so we address how planet formation could be affected by outflows. The reader may consult reviews of earlier material in Ray et al. (2007), Pudritz et al. (2007), Frank et al. (2014), Li et al. (2014b), and Bally (2016).

2. Observational Overview

It is remarkable that although magnetic fields are thought to collimate YSO jets, very few observations of field strength or direction are known. In contrast, through optical and near-infrared emission line imaging and spectroscopic studies, many jet parameters can be derived, e.g., neutral density, ionization fraction, temperature, jet opening angle, radial velocity, etc. With the addition of multi-epoch imaging, quantities such as the tangential velocity of knots, post shock cooling times, etc., can also be found (Frank et al., 2014), allowing a full 3-D kinematic study of outflows. We begin our overview from the largest outflow scales, and gradually focus down to accretion and outflow from the protostellar surface.

2.1. Do the Magnetic Fields of the Parent Cloud or Cores Determine the Outflow Direction?

It is now known that outflows from young stars extend to parsec and even tens of parsec scales, i.e., lengths comparable to the size of their parent molecular cloud (Bally, 2016). This is not surprising when one considers that it takes approximately 1 million years (Evans et al., 2009) for a solar mass young star to go through its outflow phase (i.e., from Class 0 to Class 2) and that typical outflows have velocities of several hundreds of kms−1. Therefore an obvious starting point is to consider whether outflows are somehow aligned, either directly or indirectly, with the magnetic field of their parent cloud. Conceivably this could be through some sort of guiding action by the cloud's magnetic field or alternatively, since accretion disks and outflows are perpendicular, by the field determining the disk's orientation. We can immediately dismiss the first suggestion since outflows are not only highly supersonic with respect to their internal and external sound speeds but also jet velocities are much higher than their internal and external Alfvén speeds. The second idea however is worth exploring. Despite low ionization levels, neutral and ionized matter in a molecular cloud are strongly coupled through collisions, and thus matter is preferentially accreted along magnetic field lines. In turn one might naively expect disks to form in the perpendicular plane.

As described elsewhere in this volume, magnetic field directions in a cloud can be obtained either through polarized absorption of continuum light by dust grains, when the optical depth is low, or through polarized emission at infrared/mm wavelengths from dust grains, when depths are high. Pioneering studies by Strom and Strom (1987) suggested there might be a relationship in the sense that the directions of outflows could follow that of the magnetic field in the surrounding cloud. Number statistics at the time however were very poor and subsequent studies have suggested there is no such correlation (Ménard and Duchêne, 2004; Targon et al., 2011) or at best present only for the youngest protostars. Of course expecting such a correlation assumes disks form only by the collapse of material along the field lines and subsequent distortion of the field into the classical “hourglass” pattern as gravity draws material inwards. Many studies however have shown that magnetic braking in such a case is far too effective and that disks, at least of the large sizes observed, should not form even when non-ideal damping MHD effects are taken into account (Duffin and Pudritz, 2009; Krasnopolsky et al., 2012). Moreover we also have to be careful in looking for a correlation with the large scale (≈ 1 pc) magnetic field direction in the parent cloud and the outflow when perhaps we should be searching for a link with the magnetic field direction in the smaller parent core and envelope. As polarization in the latter can only been measured when the optical depth is sufficiently high, and the scales we wish to investigate are at most a few thousand au, one has to turn to millimeter interferometry of younger embedded protostars. Here initial results (Hull et al., 2014) do not suggest any correlation and, in fact, if anything a tendency for outflows and magnetic fields in the cores to be orthogonal. We note however that a smaller study that used a more carefully selected sample of disks, came to the opposite conclusion, namely that there is a good correlation between outflows and magnetic field geometry on the scale of cores in class 0 sources (Chapman et al., 2013). Overall, a lack of correlation is understandable in the context of more up-to-date MHD simulations which show that disks, similar in diameter to what is observed, are only formed if there is misalignment between the rotation axis of the core and the ambient magnetic field providing the field is weak (Joos et al., 2012; Krumholz et al., 2013). Moreover the effectiveness of magnetic braking is also in line with the recent study of Galametz et al. (2018) who has shown there is increased likelihood of alignment between the outflow direction and the orientation of the ambient magnetic field if there is no small-scale multiplicity and no large disks on > 100 au scales. In summary then we do not expect, or find, any strong correlation between outflow orientation and the direction of the ambient magnetic field on scales of several hundred au to parsecs.

2.2. What Do Observations Tell Us About Magnetic Fields in Jets Far From the Young Star?

It is well known that most of the radiation emitted by jets from young stars comes from the cooling zone behind shocks. The properties of such shocks, in the absence of a magnetic field, have been modeled for a number of years and depending on the shock velocity various emission lines and relationships between the fluxes of such lines, are expected in the post shock zone. At optical wavelengths for example, we expect to see both lines from neutral, e.g., from H, O, and singly ionized species, e.g., S+. As jets from young stars, or more precisely their atomic/ionized component, typically travel at hundreds of kms−1, we would also predict the presence of highly ionized species, for example O++, in the post shock cooling zone if most of the outflow energy is converted into radiation at a single shock. In reality, except for the tip of the largest bow shock shaped features (e.g., in HH 34) such emission is not observed. This is because the shock fronts seen in jets are so-called “working surfaces” in which the shock velocity is determined by the difference in velocity between consecutive flows from the young star (Raga et al., 2002). The line emission in this case is then consistent with a shock velocity of several tens of kms−1 rather than a few hundred kms−1, i.e., in line with the velocity differences. Recent 3-D simulations, with realistic cooling, show the rich structure that such variations can lead to, see for example Hansen et al. (2017). What is important to emphasize is that the results of such simulations are in very good agreement with what is found through observations, for example, with HST.

Of course in a standard, i.e., adiabatic, strong shock the density of the flow increases by a factor of 4 and the velocity correspondingly decreases by the same factor. Here we are referring to velocities in the shock frame and by the term strong shock, we mean one in which the incoming velocity V s ≫ c s where c s is the sound speed. In other words M s ≫ 1 where M s = V s /c s is the Mach number of the shock. In a radiative shock, due to the loss in energy, gas in the post shock zone is compressed over and above that expected for an adiabatic shock. The increased density leads to further collisional excitation and strong line emission is produced from oxygen, hydrogen, sulfur, etc. Gas density is then increased further in order to maintain pressure and effectively these lines act as thermostats keeping the temperature around 104 K.

Unfortunately the presence of a magnetic field in a jet cannot be determined in a straightforward manner through, for example, observable Zeeman splitting of its emission lines. If a magnetic field is present however, in particular one parallel to the shock front, or with a significant component that is parallel, then clearly the field will resist compression in the cooling zone and, in turn, this will have an impact on the shock's emission properties, e.g., observed line ratios. However, as pointed out many years ago by Hartigan et al. (1994), interpreting the results is complicated. This is because observable line ratios are degenerate with respect to the field's strength as spectra from low velocity shocks, without a magnetic field, resemble those from higher density, higher velocity shocks with a field. Breaking this degeneracy requires additional information which may be present. For example if the shock is well resolved spatially and is bow shaped, then the extent of the [O III] emission near its tip, can be used directly to infer the shock velocity at the apex. In addition the Hα flux from the bow, assuming it is clearly separated from the cooling zone, and the [SII]λ 6716/λ 6731 line ratio can be used to determine both the pre-shock density and the post-shock compression, respectively. This in turn reveals the strength of the pre- and post-shock fields in the shock plane. In this way Morse et al. (1992, 1993) found weak fields, of around 30 μG, in the gas ahead of bow shocks in the HH 34 and HH 111 outflows. Compression of the field then resulted in around 1 mG in the cooling zone. It should be said that this equates to an Alfvén speed of around 10 kms−1 in most of the observed post-shock zone, i.e., comparable to the sound speed but orders of magnitude smaller than the velocity of the jet itself. Put another way, in these outer regions the magnetic field is no longer dynamically important although the situation may be very different close to the source.

A variation on the above method to measure the magnetic field using shock physics has been proposed by Hartigan and Wright (2015). Here the degeneracy is broken by measuring the extent of the cooling zone which is obtained from the distance, allowing for projection effects, between the Balmer emission lines and the forbidden lines. Applying the method to existing data for a bright knot in the HH 111 jet, they obtain a relatively low Alfvénic Mach number, M A = V S /V A , indicative of super-magnetosonic velocity perturbations in the jet.

Another important effect that magnetic fields in jets have is on softening the effects of a shock. In a standard non-magnetic shock, the transition from pre- to post-shock gas conditions is sudden and the width of the shock is a few mean free paths. The presence of a magnetic field effectively allows energy to be transmitted upstream ahead of the shock, accelerating it, and giving rise to a gentler profile in which the heating of molecules and atoms is more gradual (Smith and Mac Low, 1997). Such shocks are known as C (Continuous) shocks as opposed to their non-magnetic J (Jump) shock counterparts. One effect of the more gradual changes in properties in a C-shock is that molecules, e.g., H 2 , may not be dissociated even at shock velocities where this might be expected (Eislöffel et al., 2000)

The presence of magnetic fields can also be indirectly inferred from jet rotation. The theory is that the wind leaving the accretion disk is collimated into a jet and focused on scales comparable to the Alfvén radius. Up to such distances the ionized wind/jet material is forced to corotate with an angular velocity equal to that of its foot-point. Since the wind is expanding away from the star, along magnetic field lines, it gathers angular momentum as it does so. Through the magnetic field there is then a back reaction on the disk, forcing material in the disk to lose angular momentum and continue to spiral inwards. It is therefore possible for the wind/jet to carry away say 10% of the mass that flows through the disk but practically all of the angular momentum of the remainder. If this scenario is correct, then jets must rotate. The amount of rotation depends on where the jet is launched from. For example, if the X-wind model is correct and the jet comes from close to the disk's co-rotation radius (typically at a few stellar radii), then the expected spin of the jet should be small. If, on the other hand, most of the material comes from further out, say at (~ 1 au), then we expect any measured jet rotation to be corresponding higher. Here we are assuming that the magnetic level arm is similar in both cases, i.e., the ratio of mass outflow to mass accretion is the same.

The search for jet rotation initially used the Space Telescope Imaging Spectrograph (STIS) on the HST (Bacciotti et al., 2002; Coffey et al., 2007) either in the optical or UV. Here the high spatial resolution of this instrument, of order 0.1 arcs was particularly important. The reason for this is that if jets are launched within a few au of a young star, then one might expect a good place to look for signatures of rotation is close to the source, say within a few hundred au. In this region the jet has not had time to interact with its surroundings, an effect that might mask, or even mimic, rotation. In this region however, the jet width is at most a fraction of an arcsecond for the nearest young stars (Dougados et al., 2000). Thus, high spatial resolution, such as afforded by HST, is needed. Moreover, if a typical jet is launched with a velocity of a few hundred km s−1 then at most we expect differences in radial velocity of say a few tens of kms−1 from one side of the jet with respect to the other in a slice transverse to the outflow axis. As shown for example by Bacciotti et al. (2000), a transverse cut through a jet does not show a “top-hat″ velocity profile, even in the absence of rotation. Instead such position velocity (PV) diagrams show a smoother decrease toward the jet edges, i.e., the poloidal velocity gradually decreases away from the jet axis. The search for a rotation signature is then a hunt for a lopsidedness, i.e., a lack of axisymmetry, in this profile.

These observations are extremely challenging however, even with STIS, given the angular widths of the nearest jets and the spectral resolution required. Nevertheless, observations of jet rotation in theory could help discriminate between different MHD jet wind launching models. One other obvious constraint that we might also reasonably impose is that any evidence for rotation must be consistent with the sense of rotation of the underlying accretion disk and that the opposing (red and blue-shifted) jets in a bipolar outflow must have opposing helicity. With these constraints in mind, results from optical/UV observations using HST have been something of a mixed bag. Moreover, even when rotation has been claimed, alternative explanations such as asymmetrical shocks (De Colle et al., 2016) and precession (Cerqueira et al., 2006) have been proposed to account for the observed effects.

An alternative observational approach in the search for rotation has however emerged in recent years. Instead of looking for a rotation signature in the atomic/partially ionized jet, studies are now being undertaken of the molecular outflow in the millimeter band using high spatial and spectral resolution interferometers such as ALMA. Although the available spatial resolution is not as good as HST (although resolutions of 0.02" have been achieved for ALMA observations of HH212, Lee et al., 2017), molecular jets tend to be wider than their atomic counterparts close to their source and the spectral resolution of an instrument like ALMA is much higher than STIS. Thus, it is possible to search amongst molecular jets for rotation signatures.

Figure 1 shows an outstanding example of a rotating, large scale molecular outflow observed by ALMA (Hirota et al., 2017)—a source whose rotating outflow was first discovered by SiO maser observations (Matthews et al., 2010). Observed at resolutions of about 50 au, observations were made in Si18O (484 GHz) and H 2 O (463 GHz) lines of an outflow from a region of massive star formation in the famous Orion nebula (KL region). The former traces emission from the outflow while the latter traces emission from a dense rotating disk. The outflow speed is ~18 km s−1 on an outflow scale of 1000 au. We are afforded an almost edge-on view of the system. Using observations of the rotation curve of the gas in the disk, Ginsburg et al. (2018) were able to estimate the mass of the central star (SrC 1) to be ~(15)M ⊙ . Both these maps (Figures 1b,c) show velocity gradients along the disk plane.

FIGURE 1

Figure 1. ALMA observations of rotation in the bipolar outflow driven by Orion Source I, a high mass YSO candidate. White contours are 490 GHz continuum maps. (a) moment 0 (color; integrated intensity) map of 484 GHz emission in Si18O; (b) moment 0 (black contours) and moment 1 (peak velocity; color) map of the 484 GHz Si18O line; (c) moment 0 (black) and moment 1 (color) maps of the 463 GHz H2O line. These latter two plots show rotational velocity gradients of the outflow- water line emission being more compact. From Hirota et al. (2017) reproduced with permission © Springer Nature.

The authors note that other possible explanations of this velocity gradient such as disk warping and other effects can be ruled out for this source making it likely to be an example of a rotating molecular outflow. The resolution is not sufficient to distinguish Keplerian vs. constant angular momentum ( v ϕ ∝ r - 1 ) for the disk. However the specific angular momentum of the outflow is approximately constant. With a measured value of this wind angular momentum of 400–600 km s−1 au, this corresponds to a source radius on the disk (where centrifugal and gravitational forces balance) of 21–47 au. The observations suggest that this is almost certainly evidence for a disk wind being accelerated off the disk far from the star or star-disk magnetosphere, and that the wind material is not swept up through entrainment by a narrow jet from the surrounding gas envelope. There are however a number of caveats worth pointing out about Source I and its immediate vicinity. First, like most massive stars, Source I is in a cluster that may not only be dynamically active, but which can cause confusion when trying to untangle the effects of one outflow from another. In particular, while previous work suggests a very young jet-like outflow from Source I in SiO and CO, this interpretation is less clear in recent ALMA data (Bally et al., 2017). Another point that should be stressed is that while our understanding of how high mass stars form seems, in the main, to resemble our picture for low mass star formation other physical processes must come into play. For example, as shown by simulations, increasing uv radiation and its ionizing effects, must be taken into account (Krumholz, 2015).

Figure 2 shows ALMA observations of an SiO jet in the famous class 0 outflow source HH 212, located in the Orion L1630 cloud at a distance of 400 pc. The importance of this source is that it allows us to probe the very earliest stages of the launch of outflows and their connection with disks. ALMA observations (Lee et al., 2017) reveal knots in the highly collimated jet within 120 au of the rotating disk (shown in red), down to 10 au scales. Of central importance is the 7σ detection of velocity gradients across these knots, perpendicular to the outflow. The sense of the gradient is the same for all of the knots, ruling out random fluctuations in the jet. The mean specific angular momentum in the jet is l j ~ 10.2±1.0 au km s−1. By assuming the conservation of angular momentum along streamlines, and projecting this back to a source radius on the disk (Anderson et al., 2003), one finds that the footpoint of this flow is r o ≈ 0.05 au, originating from the innermost regions of the disk. This is consistent with an origin from the innermost zones of the disk, from either a close-in disk wind, an X-wind or even a stellar magnetospheric wind.

FIGURE 2

Figure 2. ALMA observations of the rotating outflows of the class 0 protostar HH 212. Shown is a zoom-in to the innermost part of the jet in SiO within ≈ 120 au of the central source, at a resolution of ≈ 8 au on top of the continuum map of the disk. The maps show the intensity (in unit of K km s−1) integrated over certain velocity range. (A) A chain of new knots trace the primary jet emanating from the disk, observed in SiO. (B) Blueshifted and redshifted SiO emission of the jet plotted with the continuum emission. Direction of rotation is shown with arrows. From Lee et al. (2017) reproduced with permission © Springer Nature.

A number of other rotating outflows have now been found, including sources such as DG Tau (Agra-Amboage et al., 2011), Orion BN/KL Source I (Chen et al., 2016), TMC 1A6 (Bjerkeli et al., 2016) with flows originating from extended disk scales, and the rotating and possibly precessing molecular outflow in HH30 (Louvet et al., 2018). In the case of the DG Tau source observed with the Plateau de Bure interferometer, the most likely explanation for the highly collimated jet seen in Fe lines is a quasi-steady centrifugal MHD disk wind ejected over 0.25–1.5 au and/or episodic magnetic tower cavities launched from the disk. The rotation seen in SMA observations of the outflow in Orion BN/KL suggests that the jet launching footpoint on the disk has a radius of 7.2–7.7 au. ALMA observations of the CO outflow in the TMC1A6 outflow indicate a source radius of 25 au from the disk, ruling out X-wind or stellar wind sources for this molecular outflow. In a similar vein, Zhang et al. (2018b) have observed the NGC 1333 IRAS 4C outflow in the Perseus Molecular Cloud with ALMA, tracing CCH and CS emission in an outflow cavity. Outflow rotation is detected from 120 to 1400 au above the disk, and its velocity is highly symmetric about the outflow axis. A flat distribution of specific angular momentum profile is observed in the outflow, with a mean value of 100 au km s−1. Projecting this back onto the disk gives a launch radius of 5-15 au, indicating a wind originating from large portions of the disk. The wind may reach out to the centrifugal radius of the disk (r c ≃ 30 au), so further high resolution observations are needed.

These various examples suffice to make the point that ALMA resolves rotating molecular outflows that originate from large sections of their disks, and perhaps out to their edges.

2.3. Outflows and Disk Formation

The observational evidence shows quite clearly that stars are born within surrounding disks or disk-like structures that, of course, are the birthplaces of planets. Nevertheless, there is still considerable debate in the theoretical models as to how and when protostellar disks first appear and what their basic physical properties are in these early phases (see the article by Würster and Li, this volume). Certainly multiplicity, across the stellar mass spectrum, is likely to be important in how disks evolve, and to effect outflow production. For example most class-0 low mass stars are in multiple systems and the level of multiplicity appears to decline as such systems evolve (Reipurth and Mikkola, 2012, 2015). This immediately suggests that dynamical, so-called N-body interaction, must be taken into account as it can lead to disk truncation and warping, precession of outflows, etc. Despite such complications, it is important to examine what effects magnetic fields might have on disk and outflow evolution. How influential magnetic fields are, can be gleamed from examining their structure, i.e., geometry, during the early stages of star formation.

The study of magnetic field structure in dense molecular clouds relies on the fact that elongated and spinning dust grains are aligned perpendicular to magnetic field lines. While radiation effects produce the spin needed for alignments by magnetic field, it is important to note that as asymmetric dust grains scatter or absorb photons they also undergo radiative torques which can be stochastic or regular (Hoang and Lazarian, 2009). These radiative torques leading to grain alignment would be minimal in regions of sufficiently high optical depths. Polarization observations of background star light passing through cores in molecular clouds are not possible at optical or infrared wavelengths because of the high column densities of such regions. However, thermal emission from these spinning grains can be detected at millimeter and submillimeter wavelengths and is weakly polarized in a direction perpendicular to the field lines (Matthews et al., 2009). Using these techniques, the relative orientation of bipolar outflows with magnetic field directions toward 16 class 0 and I outflows were measured (Hull et al., 2013). These high resolution submm observations were carried out at the CARMA. The results showed that magnetic fields on 1,000 AU scales are consistent with models in which outflows and magnetic fields are either randomly aligned, or preferentially perpendicular. These results may also be interpreted as the random alignment of disks with magnetic fields—as would arise if turbulence plays an important role in disk formation (see review by Hull and Zhang in this volume).

Probing the actual geometry of outflows and magnetic fields on disk scales is now possible with polarimetric observations at mm as well as mid-infrared wavelengths (Chuss et al., 2019; Jones et al., 2019). Observations of L1448 IRS 2 in ALMA band 6 were carried out by Kwon et al. (2019). Dust polarization can be affected by radiation anisotropy as well as self-scattering, but these it is argued may be minimal in this source. Figure 3 shows that the magnetic field takes on the classic hour glass form expected for contracting magnetized disks threaded with an ordered field. The field is aligned with the outflow in the vertical direction, but switches to perpendicular alignment at the disk. The latter effect indicates the presence of a strong toroidal field at the disk. If this magnetic field interpretation is correct, then the field has been strongly wound up in the disk. The Davis-Chandrasekhar-Fermi (DCF) method was used to measure the field strength in the disk region of 720 μG, although some caution needs to be exercised in applying this to regions in which turbulence may not be strong. This strong field would brake the disk quickly; in 1,400 years. The ordered outflow and presence of a large disk perpendicular to it might suggest that magnetic braking effects have been minimized—perhaps by ambipolar diffusion. Chuss et al. (2019) present the best magnetic field map of the OMC1 core and outflow region to date. Again a classic hourglass pattern is seen.

FIGURE 3

Figure 3. Magnetic field morphology around L1448 IRS 2. (Left) The green vectors indicate the inferred magnetic field direction. The gray scales and white contours present the total intensity distribution with levels of 2 to 129 times 0.1 mJy beam. The blue and red contours are CO 2-1 intensity distributions integrated in velocity ranges of -8:0 to 2.0 and 7.5 to 16.0 km s−1, respectively. (Right) Zoomed-in central region. The synthesized beams of the CO and the continuum data are marked in the bottom right corner in blue and white, respectively. From Kwon et al. (2019) reproduced with permission © AAS.

Perhaps the largest sets of high resolution images that we now have of protostellar disks is a consequence of the Disk Substructures at High Angular Resolution Project (DSHARP) campaign at ALMA. The initial survey of 20 protostellar disks at a resolution of 5au FWHM provides the first look at the small-scale features in disks that are directly relevant to the planet formation process. It must be noted that this is a diverse collection of disks that are typically much larger than the bulk of disks now turning up in other surveys. They are images of Class II systems (ie disks well beyond the formation phase), and can be regarded as protolanetary disks (disks in earlier phases of star formation—the so called Class 0 and I phases, are denoted as protostellar disks). These images show a variety of rings and gaps, amplitudes, and disk sizes for targets with a range of stellar and disk properties (Andrews et al., 2018). The ubiquitous ring and gap structures have no obvious spacings that are connected with properties of the host star. The ring systems may arise through modulations of disk pressure that trap dust grains. Specifically, the observed system rings have been modeled as arising by dust trapping in axisymmetric pressure bumps (Dullemond et al., 2018). If the turbulence in disks is relatively weak, then dust can be strongly confined in pressure bumps, while strong turbulence compromises the trap. Theoretical modeling of the ALMA ring data suggests a lower limit to the turbulent amplitude expressed in the standard turbulent viscosity parameter (α) which takes a lower bound of α ≥ 5 × 10−4. Without this level of turbulence, the rings would be even sharper features than they already are. It is possible that these pressure bumps themselves are produced by planets. In any case, the model constraints on the turbulence level is of interest in comparing the relative role of turbulent vs. disk wind torques in disks.

2.4. What Radio Emission From Outflows Can Tell Us

Observations of YSO jets at radio wavelengths (around a few centimeters) have shown them to primarily emit thermal radiation (i.e., radiation with a positive spectral index α where the flux at frequency ν, S ν ∝ ν α ). Moreover, typical flux densities are very low and usually less than 1 mJy (Anglada et al., 2018). Such emission is un-polarized and is often elongated in the direction of the known atomic/molecular outflow. Until recently the amount of information that could be obtained from radio maps was limited, primarily due to the lack of sensitivity of radio interferometers, nevertheless they provided upper limits to the diameter of jets and have proven that they must be collimated on scales less than 50 au in agreement with high spatial resolution HST data (Reipurth et al., 2004). In the past few years however radio astronomy has undergone a revolution that is having a major impact on our ability to image outflows from young stars. This revolution is akin to the replacement of photographic plates in optical astronomy (with quantum efficiencies (QE) of a few per cent) by the CCD (QE ≈ 70% or more). In the past traditional radio interferometers, such as the Very Large Array (VLA) in New Mexico and the high spatial resolution Multi-Element Radio Linked Interferometer (MERLIN) centered on Jodrell Bank, used microwave links between the various telescopes to transmit wavelength, amplitude and phase information. These links were of limited bandwidth and correspondingly only a narrow radio continuum could be observed at any one time. Such old links however have now been replaced by fiber, allowing considerably larger observing bandwidths that are used in conjunction with much-improved correlators, lower noise amplifiers, and increased computer processing power. The net result is an enormous increase in radio continuum sensitivity (typically by an order of magnitude or more). For example the VLA in standard continuum observing mode had 4 sub-bands of 50 MHz each but the JVLA has 2 sub-bands 2 GHz wide! Of course this has a dramatic impact on radio studies of jets from young stars. Complete radio surveys of YSOs in nearby star formation regions, for example, have become possible for the first time (Dzib et al., 2015; Tychoniec et al., 2018).

A second strand to the radio astronomy “revolution” is the opening up of the low frequency radio spectrum with the advent of instruments such as the Low Frequency Array LOFAR (van Haarlem et al., 2013) radio telescope. LOFAR is allowing us to explore jets from young stars in a new band (at meter wavelengths) and, as described below, could reveal any non-thermal emission. A particular advantage of the radio band is that it can potentially provide very valuable information on magnetic fields in YSO jets. Theoretically we only expect magnetic fields to dominate the dynamics of an outflow close to the source, i.e., within 100 au, and in high velocity post-shock regions where the field is compressed and amplified (Hartigan and Morse, 2007). Thus we have a higher chance of measuring jet fields by studying the zone surrounding the young star at high angular resolution and the brightest knots (compressed zones) on more extended scales. But how might we do this? One powerful method is hinted at by the fact that a small number of jets/outflows have non-thermal radio spectra as demonstrated in the case of more massive young stars such as Herbig-Haro 80/81 (Carrasco-González et al., 2010) but also in their lower mass counterparts (Ainsworth et al., 2014). Such radiation appears to come from high-energy electrons (i.e., it is gyro-synchrotron or synchrotron emission) and, as discussed below, it may be common at very weak flux levels in “standard” low luminosity outflows but nevertheless detectable with the new suite of radio interferometers such as e-MERLIN, the Jansky VLA (JVLA) and at low frequencies (e.g., using LOFAR).

Figure 4 shows the radio continuum map of HH 80-81 (Carrasco-González et al., 2010). The analysis of the synchotron emission yields a field strength of ~0.2 mG in the jet. The field is parallel to the jet axis, and the increase of polarized emission toward the limb of the jet is evidence for a collimating toroidal field. These characteristics are very similar to those of AGN jets, providing good evidence to support the idea that the jet mechanism is universal.

FIGURE 4

Figure 4. High resolution, 6 cm radio continuum VLA images of the jet HH80-81. The knots, 0.5 pc from the source are linearly polarized indicating non-thermal synchotron emission from the jet. (A) Radio continuum image of the whole jet. (B) Linearly polarized continuum intensity image. (C) Apparent magnetic field directions superimposed on continuum map. From Carrasco-González et al. (2010) reproduced with permission © AAAS.

As stated previously, we primarily expect to detect free-free emission in the cm (1-010 GHz) band. However, non-thermal emission is observed on occasions and this opens up the intriguing possibility of directly measuring both the magnetic field's strength (using its intensity) and direction (using polarization). Such measurements, although difficult, could provide the “missing link” needed to properly evaluate the importance of magnetic fields in jets, as essentially all other parameters are known. The improved sensitivity and high resolution of e-MERLIN and the JVLA may allow us, in a systematic way, to study this type of emission as faint non-thermal radiation (arising from electrons being accelerated to relativistic speeds near shocks) mixed in with free-free emission (Rodríguez-Kamenetzky et al., 2016).

2.5. Magnetic Fields Close to Young Stars

Magnetic fields are thought to play a very important role in accretion onto young stars and the evolution of their outflows. In particular circumstellar magnetic fields are considered to be sufficiently strong to clear the central region of a young star's accretion disk and force the star to co-rotate with the disk's inner edge (Ferreira, 2013). Moreover, at the same time as they funnel material onto the star, they may also help remove excess angular momentum thereby allowing accretion to proceed and the star to contract to the main sequence. Magnetic fields may even play a role in planet formation since the inner gap they generate could allow hot Jupiters, which have already been discovered around T Tauri stars (Donati et al., 2016), to survive without them being accreted by their parent YSO (Adams et al., 2009).

It is only in recent years that we have been able to determine the geometry and strength of the fields in classical T Tauri stars (Johns-Krull, 2007; Donati et al., 2012). Such studies have revealed some surprises, including the fact that the weakest fields have the most complex geometries, i.e., they are not simple dipoles, and that these tend to be associated with the development of a substantial radiative core (Figure 5, Left). Moreover these fields are found to vary, which with the addition of their complex geometry, suggests they are generated through dynamo action and are not fossil fields as has been assumed in the past. Determining the magnetic field direction and strength can be done using Zeeman Doppler Imaging or ZDI (Morin et al., 2011; Gregory et al., 2012), i.e., by spectro-polarimetry using instruments like HARPS-Pol at ESO. This involves measuring the circular polarization signal (Stokes V parameter) in both photospheric absorption lines and accretion dominated emission lines over several rotation cycles. The signal in accretion lines is usually so strong that measurements can be made using individual lines (e.g., He I 5876) although it should be stressed that this represents the field in the vicinity of an accretion “hotspot.” In contrast the signal in magnetically sensitive photospheric lines is so weak that cross-correlation techniques over many lines have to be employed to get an average Zeeman signature (Gregory et al., 2012). Combining both sources of data together, and by monitoring how the Stokes V parameter changes as a function of rotational phase, it is then possible to map the magnetic field. While very few classical T Tauri star magnetic fields have been mapped in detail (Hussain and Alecian, 2014), application of the method has proven very successful including the case of AA Tau (Donati et al., 2010) a known jet source. The magnetic field of AA Tau is shown in Figure 5.

FIGURE 5

Figure 5. Magnetic field map of the surface of the known jet source AA Tau. The field strength, shown by the colored scale, is in gauss for the radial, azimuthal and meridional components which are shown as a function of stellar phase. From Donati et al. (2010) reproduced with permission © OUP.

In addition to analyzing the magnetic field of the young stars, it may also be possible with the next generation near-infrared (NIR) spectro-polarimeters, such as SPIRou (Donati et al., 2018) to detect spectral features from the innermost regions of the accretion disks of high mass accretion rate stars (Class I and Class II). This could allow us to map both the “normal” disk poloidal and toroidal magnetic fields for the first time in an analogous way to what has been achieved for a CTTS disk in the more excited FU Ori state (Donati et al., 2005). A particular advantage of spectro-polarimeters like SPIRou is that not only do they operate in the NIR where disk emission can come from but also Zeeman splitting increases with λ2, all else being equal, and this makes it easier to measure. Note also that this region of the disk is not accessible to ALMA. In addition, as SPIRou can measure radial velocity with exquisite precision it should be possible from the data to detect the presence of hot Jupiters.

2.6. Feedback: Outflow Interaction With Molecular Cores and Beyond

There is no doubt that the final mass of a young star can be influenced by a whole host of environmental factors. As pointed out earlier, stars are frequently born in multiple systems that are dynamically active. Interaction between the companions can, for example, truncate a star's growth at an early phase (Bate, 2018). Radiation pressure from the newly formed star itself, particularly in the high mass regime, can limit accretion (Sartorio et al., 2019) as can photo-evaporation from external sources such as supernovae and neighboring OB-stars (Haworth and Clarke, 2019). Nevertheless it is thought that outflows may be very important in limiting accretion and thus the final mass a star can attain. Clearly they drive away material from the core and the envelope once they start to operate. It is not immediately obvious however what fraction of the final mass is removed and deposited back into the parent cloud. A related problem is the effect of outflows on the molecular cloud itself. It has long been known that the lifetime of a molecular cloud should be shorter than observed if the only force acting against gravity is thermal energy. Instead an additional force must be present. In some cases this might come directly from the magnetic fields threading through a cloud. We also know however from observations that supersonic and super-Alfvénic turbulence are present although the latter may only exist in the denser regions (Heyer and Brunt, 2012). The presence of turbulence alone however does not resolve the issue as the timescale for turbulence to decay in a typical cloud is short: there must be some local source of turbulence. Can this be outflows from young stars?

Certainly outflows create cavities in their surroundings while at the same time ejecting momentum and energy into it. The net effect is to disperse the core (Arce and Sargent, 2006) and ultimately to terminate accretion. How efficiently outflows can do this determines the final stellar mass and also the fraction of the core's mass that is converted into a young star, i.e., the ratio between the so-called Core Mass Function (CMF) and stellar IMF (e.g., Offner and Arce, 2014). As a specific example, the core mass function of low mass, star forming cores in the Pipe Nebula, was measured by using infrared extinction of back ground stars (Alves et al., 2007). This showed that the IMF is the direct product of the dense core mass function and a uniform star formation efficiency of 30 ± 10%. As already pointed out, both functions have similar profiles (see for example Cheng et al., 2018) the difference being the peak in the CMF is at roughly 3 times the mass of peak in the stellar IMF. Obviously if outflows typically remove 2/3 of the core mass, this could explain the observed relationship. A very good overview of the processes that can determine the relationship between the CMF and the IMF can be found in Offner et al. (2014).

In additon the expectation that outflows sweep away such a large fraction of the surrounding gas over their lifetime is supported by recent observations of the HH46/47 system using ALMA (Zhang et al., 2016). In carrying out such observations it is important to use optically thin tracers, e.g., 13CO, to ensure the motion of the densest gas is not overlooked (Arce and Sargent, 2006).

Numerical simulations, for example by Offner and Chaban (2017), of the collapse and evolution of isolated dense cores now include the effects of turbulence, radiative transfer, and outflow feedback. These show that outflows can drive and maintain turbulence in the core environment even when magnetic fields are initially strong. Moreover, the final efficiencies are 15–40% in line with the observed values. Of course these are simulations at the individual star level: it should not be neglected that multiplicity can play a role in feedback not only in multiple systems but also on the larger cluster scale. In multiple systems, interaction between the components can lead to precession of the associated jets (Wu et al., 2009) thus broadening their impact on their surroundings. Moreover many stars are born in clusters, for example in regions such as NGC 1333 containing hundreds of stars and large numbers of criss-crossing outflows (Walawender et al., 2008). This suggests we need to consider the ensemble to get a clearer picture of how outflows feedback into their environment. Whether magnetic fields or turbulence driven by stellar feedback ultimately control how molecular cloud material is converted into stars is still an open question (Li and Nakamura, 2006; Federrath, 2016).

3. Theoretical Overview

The last 5 years have also seen remarkable progress and new developments in the theory and simulations of magnetized outflows. The emergence of new computational tools and codes has played a leading role in these advances. These include time dependent radiation field and chemistry within MHD simulations (Gressel et al., 2015). On the horizon are new capabilities including zoom-in simulations in radiation MHD simulations using DISPATCH (Nordlund et al., 2018). Post-processing tools now available include Monte Carlo radiation transfer methods in RADMC3D. The ability to compute expected polarization maps from general MHD systems has recently appeared in the POLARIS code (Reissl et al., 2016) which has been applied to the polarimetry of outflows (Reissl et al., 2017).

3.1. Drivers for Outflows

There are several sources for magnetized outflows that are launched during the course of star formation, involving centrifugal or magnetic pressure drives from diverse sources such as forming disks in collapsing molecular cloud cores, the coronae of well established disks, magnetospheric boundaries, or the protostars themselves. Before getting into the details, we outline this general theory landscape.

(i) Magnetic braking and early outflows: As magnetized, pre-stellar cores begin to emerge from the filaments within turbulent clouds, magnetic braking begins to extract their angular momentum. The models of isolated, rotating molecular cloud cores with subsonic line widths that figure so prominently in many early theoretical and numerical studies are overly simplified pictures. In reality, hydrodynamic turbulence always contains small regions of non-zero angular momentum which, fed with enough mass, can produce disks—there is no need to assume that any coherently rotating object was present. In magnetized clouds in this earliest stage, turbulent MHD processes play a significant role in diluting the net braking torques on such pre-stellar cores, allowing disks to begin to form. The magnetic field geometry in these chaotic initial conditions likely become more ordered close to the disks as dissipative effects such as ambipolar diffusion take effect in higher density regions. Finally, at some point, collapsing magnetized cores reach their centrifugal balance radius. With the formation of this disk, or disk-like system, the pinching of field lines and/ or the accumulation of toroidal magnetic field results in the launch of the outflow. These are likely to global encompassing the entire eary disk. Class 0 ouflows likely have their origin in this event.

(ii) Magnetized disk winds: These will arise naturally during the evolution of magnetized accretion disks in their post formation phases. This has been shown by many authors, both theoretically and by means of numerical simulations using a variety of computer codes. These originate from disks that are both highly conducting, or in modern treatments, even from disks for which non-ideal effects dominate the physics of how fields are coupled to the gas and dust. Details of the launching mechanism likely involve a thermal component to the disk wind as it is launched from the upper layers of the disk or base of a disk corona. The large scale suppression of the MRI instability discovered in non-ideal MHD processes means that disk physics itself is not controlled so much by viscosity as it is by the outward transport of disk angular momentum by powerful winds.

(iii) Stellar spin and accretion powered stellar winds: The solution of the angular momentum problem for the protostar itself is another key aspect of star formation physics. If stars were to accrete the angular momentum from the inner edge of a Keplerian disk, they would spin up to break up speeds within a few hundred thousand years (the spin up torque of gas accreting from the magnetopause boundary r MP onto the star is then J ∙ ~ Ṁ w Ω K e p , M P r M P 2 , Matt and Pudritz, 2005). An early suggested solution for this dilemma proposed that magnetic field lines from the star, connect back into the disk such that flux penetrating the disk beyond the co-rotating orbit would effectively brake the stellar spin by depositing angular momentum back into the disk (Koenigl, 1991). This “disk locking” picture has the problem that magnetic field lines would quickly shear out and reconnect, severing the connection between the star and disk. Centrifugally driven winds from the protostars can carry off angular momentum that is deposited by accreting gas streams flowing along magnetospheric field lines that connect the disk's inner edge, to the star (Matt and Pudritz, 2005). This work has met with considerable success in being a part of the rotational history of young stars (Bouvier et al., 2014). Earlier models of X-winds proposed that angular momentum from the disk inner edge is never accreted onto the star but is carried off in an X wind originating from the interface between the disk and the protostellar magnetosphere (Shu et al., 2000).

3.2. MHD Outflow Basics

To understand why magnetized outflows are so powerful and efficient in angular momentum extraction, we recall basic concepts of MHD flows from rotating bodies laid out in many papers and reviews (Konigl and Pudritz, 2000; Pudritz et al., 2007; Spruit, 2010). In the case of axisymmetric systems, we consider a rotating body threaded by magnetic field lines, with the footpoint of a field line located at a radial distance r o from the rotation axis, along the rotor. The angular velocity of this footpoint at the base of the flow is Ω o . The basic equations of stationary, axisymmetric MHD flows (e.g., Pelletier and Pudritz, 1992; Bai, 2016) tell us that the mass flux is conserved along streamlines of the outflow (from the continuity equation), as is the magnetic flux Φ. That is, the mass loss rate per unit area, along streamlines passing through an annular section flow of area dA, is dṀ w /dA = ρv p (where v p is the poloidal velocity along the field line) while the amount of poloidal magnetic flux through this same annulus is dΦ = B p dA. This means that the ratio of these two fluxes along a field line—the mass loss rate per unit magnetic flux, otherwise known as the mass loading of the flow—is also preserved: k = dṀ w /dΦ = const.

The mass loading is determined by the physical conditions at the foot point of the flow, which is essentially at the slow magnetosonic point on each field line. As an example, if we consider the case of accretion disks as our basic rotors, the early 2D axisymmetric disk wind simulations of Ouyed et al. (1997) and Ouyed and Pudritz (1997) assumed that the flow originates at the base of a heated disk corona, whose heat source is not specified. Recent work takes the additional step of explicitly including FUV flux from the host star that heats the disk surface, feeding a disk corona (Bai, 2016). In any event, one needs to have a physical theory of how magnetic flux is distributed across the disk, and how mass is loaded onto it, in order to determine the mass loading function k = k(r o ) which is a function of the conditions at the footpoint r o (ie the slow magnetosonic point) of each field line on the rotor.

No theory that we are aware of has determined what this function is from basic first principles and so this remains one of the lingering uncertainties in the theory of MHD disk winds. Authors adopt some sort of “reasonable” scaling for the magnetic field distribution across the disk. As a specific example, consider outflows from magnetized Keplerian accretion disks. Suppose that the threading vertical field component is a power law B z ∝ r o μ - 1 , the outflow speed at the footpoint scales with the local Kepler speed, and the base of the disk corona can be considered a polytrope (e.g., γ = 5/3). The mass load function then scales as k ∝ r o - ( μ + 1 ) . As an example, Blandford and Payne (1982) assume a self-similar model for disk winds which imposes a specific scaling for the field B z , B P ∝ r o - 5 / 4 , so that μ BP = −1/4 and hence k B P = r o - 3 / 4 . The stable minimum energy jet structures investigated by Pelletier and Pudritz (1992) have k P P ∝ r o - 1 / 2 (Pudritz et al., 2006). In general the power law index will also determine the degree to which poloidal magnetic field lines open up away from the footpoints. Recent global MHD models of outflows (Wang et al., 2018) assume that heating of the disk surface by stellar EUV photons provides thermal heating and provides thermal pressure that assists in driving an MHD disk wind. Here too, the initial distribution of magnetic field is assumed to follow a powerlaw distribution whose field lines initial structure is controlled by the choice of the power law exponent.

It is worth emphasizing that the concept of centrifugal acceleration of the flow arises naturally in the co-rotating frame of the footpoint of the flow (e.g., Spruit, 2010). There is no magnetic force driving the outflow in this frame that role being played by the centrifugal force. The flow in this region can be visualized as “beads on a wire.” This is a very useful way to analyze the dynamics as long as the flow co-rotates with the rotor. This ends roughly at the Alfvén radius, where magnetic and inertial forces balance. The collimation of jets and outflows begins after the outflow passes through the Alfvén point on each field line. The reason is that the field lines now become increasingly toroidal. The degree of outflow collimation depends quite sensitively on the mass loading of field (Ouyed and Pudritz, 1999; Anderson et al., 2005). This is because it is the inertia of the material along the field line that causes the rotating field line to fall back from the rate of rotation of the footpoint, creating the toroidal field component. This is much like the behavior of a mass loaded whip. The solution of the induction equation specifically shows that the toroidal field B ϕ in such an outflow depends on the mass loading. Since collimation of such rotationally driven outflows depends on the radial hoop stress exerted by the tension of the toroidal field line, mass loading must then control collimation properties of such winds. Specifically, magnetic fields in disks with power laws μ > −1/2 will collimate to cylinders, while lower values will lead to wide-angle outflows that are parabolic at infinity (Pudritz et al., 2006).

As already noted, a central aspect of magnetized outflows is that they efficiently extract angular momentum from the rotor. How does this happen? Again, we first begin with ideal MHD since it allows us to see the basic physics quite simply. The conservation of angular momentum equation for our magnetized rotating system says that the conserved angular momentum per unit mass l of the fluid along a field line consists of two parts, the bulk rotation of the fluid itself rv ϕ , and a component that is carried by the twisted mass -loaded field line, giving l(r o ) = rv ϕ − (rB ϕ /4πk). For any theory of stationary winds, solutions are obtained by requiring that they flow smoothly through a critical point on each field line (and for 2D flows, a critical line). Critical points are positions at which the outflow speed is equal to a signal speed in the flow. For purely hydrodynamic flows, the critical point occurs where the outflow speed equals the sound speed in the gas, v p, hydro = c s . For MHD flows, there are three MHD waves: the slow and fast magnetosonic waves (compressive modes) and the transverse Alfvén wave. The value of l on any field line is set by the condition that v p = V A , the Alfvén speed. If r A is the radial distance of the field line from the outflow axis where this condition is met, then l = r A 2 Ω o = ( r A / r o ) 2 l o where l o = Ω o r o 2 is the specific angular momentum of the material at the footpoint of the flow.

This result reveals why outflows are so efficient. Each wind particle carries of a factor of ( r A / r o ) 2 more angular momentum per unit mass than its fellows rotating in the disk at the field line base. Thus, if the typical lever arm is r A /r o ≃ 3, then each wind particle can carry the angular momentum of ten particles in the underlying rotor. From this argument, it is also easy to see that in the limit that winds carry away all of the angular momentum of a disk that is accreting at the rate Ṁ a , the wind mass loss rate needed to drive this should simply scale as Ṁ a ≃ ( r A / r o ) 2 Ṁ w . Indeed, this result follows from disk wind theory outlined in the following subsection.

The conservation of energy in ideal MHD, centrifugally driven winds is expressed by the Bernoulli equation for flow along the magnetic stream line. One finds that the terminal speed of the outflow is greater than the escape speed from flow at the base of any field line, specifically v ∞ ≃ (r A /r o )v esc, o . It is this energy relation that ensures that hydromagnetic outflows scale with the depth of the gravitational well of the rotor. This is why star formation, from brown dwarfs to massive stars, is accompanied by jets and outflows. It is inevitably the efficient tapping of gravitational potential energy released during accretion by magnetized wind torques that is the main driver. This also implies that because the terminal speeds of outflows on each field line scale with the Kepler speed of the footpoint, then disk winds will have a wide range of terminal speeds—the highest (hundreds of km per sec) in the interior parts of the jet or outflow, the lowest on largest outflow scales that can be supported on the disk (10 km per sec or less). This is distinctly different than the prediction of X-wind models which really have only one velocity component originating from the innermost radius of the disk. Observations generally support the former velocity structure (recall however, the SiO observations of HH212, which are consistent with an inner disk-wind or an X-wind origin; Lee et al., 2017).

Another aspect of winds accelerated from disks is the launch angle at the footpoint of the flow. An intriguing aspect of such outflows is that they can be launched for completely cold gas—one only needs the field line to be bent by an angle of 30o from vertical (Blandford and Payne, 1982). Of course, any additional heating of the outflow region near the surface of the disk will mean that radial thermal pressure gradients will increase this opening angle. As an example of this effect, polytropic models for the disk corona will have opening angles for field lines that increase with radial distance across the disk (e.g., Pudritz et al., 2006). Vertical pressure gradients can also lead to better acceleration along the field lines. Detailed heating models of disks by FUV irradiation from the star (Bai, 2016) will have their own prescriptions of field line opening.

Early numerical simulations of jets focused on initial axisymmetric setups, often focused on the launch region near to the disk. However, fully 3D simulations are needed to understand the complexity, stability, and asymmetries that can arise in jet dynamics in both the launch region and far from the disk.

In Figure 6 we show a global, high resolution 3D simulation, (using ZEUSMP, Norman, 2000) of a collimating disk wind whose source is in the corona of a Keplerian disk threaded with magnetic field lines (Staff et al., 2015). This figure shows the toroidal velocity and magnetic field lines in a jet launched from the corona of a Keplerian disk. The magnetic geometry on the disk (and hence mass loading) follows that of Pelletier and Pudritz (1992). The jet is heated by shocks throughout its volume and the densities and temperature profiles are used to compute forbidden line emission in various lines in [CI IV], [S II], and Mg II transitions, whose critical densities are n ≈ 1, 104, 105 cm−3, respectively. The red and blue colors indicate Doppler shifts of rotating knots in the flow. The jet field has a strong poloidal field along the axis of the jet which acts as a “backbone” that helps preserve the jet against instabilities. Note that the toroidal field wraps around the jet, especially in the outer regions, giving rise to a highly collimated jet. This as noted, agrees with observations of field structures seen in discussed in the observations of H80-81 (see Figure 4). Other simulations with a less steep gradient of disk magnetic field did not collimate so well—as discussed above.

FIGURE 6

Figure 6. 3D MHD simulations of a disk wind showing the toroidal velocity and magnetic field lines for a disk wind produced in the Pelletier and Pudritz (1992) model. The outflow source is the base of the corona of a Keplerian disk, shown on the left of the figure. The simulations are run on a grid extending out to 90 au along the jet axis, and 27 au on either side of the jet axis. Brightness indicates higher radiation intensity from shock-produced, forbidden line emission in the jet. The red and blue colors of the knots indicate red and blue shift due to jet rotation. Note also the recollimation of the jet half way down the flow. From Staff et al. (2015) reproduced with permission © OUP.

3.3. Early Phase: Gravitational Collapse, Outflows, and Disk Formation

A problem of critical importance for stellar birth is the formation and evolution of protostellar disks and the concomitant launch of outflows. These are highly nonlinear phenomena whose complexity arises from the interplay of gravity, turbulence, angular momentum, magnetic fields, radiation, thermodynamics and even chemical evolution. Despite this, the appearance of disks and outflows in the earliest stages of star formation appears to be universal and this poses several important challenges for theory.

3.3.1. Idealized Models

The first theoretical treatments of star formation in magnetized clouds focused on magnetic braking of idealized, uniformly rotating gas spheres. This is enabled by a flux of a special wave: the torsional Alfven wave—first discussed in mathematical papers by Gillis et al. (1974) and Gillis et al. (1979). Subsequently analytic solutions in ideal MHD were found by Mouschovias and Paleologou (1980). The idea is simple to visualize. Consider the formation of a condensation that is threaded with a magnetic field in a medium of non-zero angular momentum. The ongoing condensation will generate a shear that twists the magnetic field, creating a toroidal field component. The resulting magnetic torque on the condensation, extracts its angular momentum which is transported away by a flux of torsional Alfvén waves. The angular momentum is transferred to the surrounding medium which is gradually brought into co-rotation with the slowing rotor. This braking picture suggests modest stirring of the environment around the forming star by a flux of these torsional waves. Given that the early stages of star formation are observed to be accompanied by vigorous outflows in the class 0 phase, this suggests that the braking phase occurs earlier, and sets the conditions for the first stages of disk formation.

This simple model, however, was later realized to create a problem for disk formation. Recent simulations have shown that braking in this picture is so efficient that disks fail to form. This has been dubbed as the magnetic braking catastrophe. How does this arise? First we recall that the exact condition for collapse to occur in the first place is that the ratio of the cloud's gravitational to magnetic energy densities exceed unity. The condition may be rephrased. The ratio of these energies can be written in terms of the mass to flux ratio of the cloud—which is equivalent to the ratio of the cloud's column density N to its field strength: (M/Φ) ∝ (N/B). The critical condition is that the ratio λ of the cloud's mass to flux, to the critical mass to flux for stability, exceeds unit: λ > 1. The critical value (N/B) crit = 2.6 must be exceeded (ie. becoming supercritical) for collapse to occur (Heiles and Troland, 2005). Li et al. (2015) performed MHD simulations of magnetized turbulent clumps and compared these to a variety of Zeeman and other observations and found that the typical molecular clump is supercritical with λ ~ 3. Thus, there is now a substantial body of observations confirmed by simulations, that places good limits on the origin and initial strength of magnetic fields in star forming cores—and therefore on the field that is available to drive outflows.

A substantial body of work on collapse calculations focuses on the consequences of the Mestel picture of idealized isolated, spherical magnetized clouds, wherein the MHD is treated as either ideal or non-ideal (Masson et al., 2016). One of the key processes that sets in on small scales of disk formation at higher gas density is non-ideal MHD. Magnetic fields will decouple from the gas at high densities, and this, it was hoped, would lead to substantial decoupling of the field from the collapsing gas. Non-ideal MHD processes such as ambipolar diffusion and Ohmic dissipation, which depend on how gas is ionized, and on grain sizes and charges at high densities, could play a role for the formation of protostellar disks. This has been investigated quite extensively by a number of groups.

As an example, Mellon and Li (2008) showed that braking can prevent the build-up of disks during the early phase of star formation for relatively week fields; λ ≤ 10. Subsequent work established that even weakening the initial field by the process of ambipolar diffusion had little effect on this result (Mellon and Li, 2009). The basic point is that by the time densities at which AD becomes effective are reached, magnetic braking is still so effective that only small disks may form. Normal ambipolar diffusion and Ohmic dissipation are ineffective (Li et al., 2011). High spatial resolution RMHD simulations carried out by Tomida et al. (2015) on the other hand, find that nonideal MHD effects can resolve the magnetic braking catastrophe in the early phase of star formation—before a protostellar core is formed. These simulations (using different inner boundary conditions, and refraining from the use of a sink particle) find the formation of a small rotationally supported (≈5 au) disk in the early phase. It is expected but not shown that the disk will grow later as larger angular momenta accrete.

One way out of this impasse for the classical model setup was suggested by Hennebelle and Ciardi (2009) who showed that by offsetting the magnetic and rotation axes of such rotors by even small angles of ≃10o centrifugally supported disks could form. They found that centrifugally supported disks cannot form for values of λ ≤ 3 when the magnetic field and the rotation axis are perpendicular, and smaller than about λ ≃ 5 − 10 under perfect alignment. However, based on realistic distributions of magnetic field strengths and misalignment angles, Krumholz et al. (2013) showed that this would lead to a Keplerian disk fraction in the Class 0 stage of only 10 to at most 50 per cent, noticeably below the observed fraction of discs around Class I/II objects.

3.3.2. Realistic Initial Conditions: Turbulence

The magnetic braking catastrophe is really only a problem for the model of highly idealized initial conditions that it presupposes. With the benefit of Herschel and ALMA observations, we now have a far better picture. Molecular clouds are characterized by supersonic velocity dispersion, attributed to supersonic turbulence (Larson, 1981). Herschel observations show that stars typically form in filaments (André et al., 2010, 2014). This likely arises from some combination of gravitational fragmentation of filaments into discrete cores, or as members of small groups or clusters at the intersection of filament systems. However, examples of the theorist's dream model—fairly isolated, nearly spherical molecular cloud cores that are virtually identical to gravitationally unstable Bonner-Ebert spheres—do exist (Alves et al., 2001)! Magnetic field observations show that star formation occurs in cores in which gravitational dominates magnetic energy (Crutcher, 2012).

The strength of magnetic fields, and hence the role of outflows in star formation, is rooted in how star forming gas becomes supercritical in the first place. The intriguing point is that the cold neutral medium (CNM) out of which molecular clouds may ultimately form is highly subcritical (Heiles and Troland, 2005). For the CNM, with measured total fields of the order B CNM = 6μG (similar in value to the local Galactic component) these authors measure λ = 0.42—deep into the subcritical, magnetically dominated regime. Starting from this initial state of the ISM, an older approach to star formation in magnetized media supposed that the local mass-to-flux ratio is reduced by ambipolar diffusion until individual cloud cores become supercritical (see the review by Lizano and Galli, 2015). This process is far too slow to be important in quiescent regions of molecular clouds (e.g., star forming cores). However, some analytical models (Zweibel, 2002) and numerical simulations (Nakamura and Li, 2005) have suggested that this can be sped up in regions with supersonic turbulence.

In supersonic turbulence however, intersections of supersonic shocks create filamentary systems (see the review by Mac Low and Klessen, 2004). When magnetic fields are included, observations show that they tend to preserve their orientation at all scales that have been probed—from 100-pc scale down to sub-pc scale cloud cores. This suggests that both gravitational contraction and turbulent velocities should be anisotropic, due to the influence of dynamically important magnetic fields and this is in fact observed (Li et al., 2014a). Simulations do show that because gas is free to flow along magnetic field lines—supersonic turbulence can compress gas flows along field lines (Padoan and Nordlund, 1999). While it has been claimed that gas motions in the denser cores are super Alfvénic (the density has increased with little change in the field Padoan, 2018), the observations suggest that the field is strong enough that motions are trans-Alfvénic. The best studied low-mass cores of nearby star forming regions are typically subsonic.

Recent 3D global simulations of the initially subcritical ISM in a magnetized spiral galaxy show that long filamentary clouds form by the draining and compression of flows into long magnetic valleys of buckling field lines. This onset of such global Parker instabilities can create a system of kpc magnetized filamentary clouds of CNM. Subsequent filamentary flow along these long HI structures gather sufficient gas to quickly produce supercritical GMCs (Körtgen et al., 2018). The action of these global modes addresses the long standing question of how it is possible to gather enough gas from a collection region of the order of 1 kpc needed just to build GMCs (see review McKee and Ostriker, 2007).

The origin of the initial angular moment of pre-stellar cores is equally important. Given the fact the molecular clouds are turbulent, it is questionable to assume that angular momentum is coherently distributed. It is often assumed in analyzing the observations that cores have a uniform rotation and follow a rigid-body rotation law. Their angular velocity Ω is deduced from measurements of global velocity gradients in velocity maps (Myers and Benson, 1983; Goodman et al., 1993). By performing numerical simulations of turbulent, magnetized clouds in 3D, it has been found that the observations of the specific angular momentum of cores (which involve 2D projections) overestimate this quantity by a factor of approximately 10 (Dib et al., 2010). Hence, it is clear that protostellar disks are likely the result of locally generated angular momentum by turbulence (Jappsen and Klessen, 2004). In such situations, the link between the local angular momentum vector that will define the disk and the direction of the magnetic field at larger distances can be essentially random. This, as it turns out, is the clue to resolving disk formation.

The first simulations of disk formation in strongly magnetized, turbulent collapsing clouds were carried out using different numerical methods. Simulations using 3D MHD Adaptive Mesh Refinement (FLASH) simulations on the scale of individual turbulent clumps were carried out by Seifried et al. (2012) and Seifried et al. (2013). This work emphasized the inefficiency of magnetic torques in turbulent media. Another group focused on the role of turbulent reconnection of magnetic field lines in reducing the efficiency of magnetic braking (Santos-Lima et al., 2012). The former approach showed that Keplerian disks of extent 50–150 AU can be formed in the collapse of turbulent clouds with the typical initial cloud magnetic fields; λ = 2.6. Disks form even in turbulent clumps with initial zero total angular momentum. They arise as a consequences of the local angular momentum contained in the convergence of several filaments in the turbulent flow. Disks formed for a range of initial core masses and turbulent settings, and for typical Mach numbers of the turbulence of M = 5, although disks formed in quite subsonic turbulence M = 0.5 as well (Seifried et al., 2015). In all cases, disks are formed by accretion out of discrete filaments and this is a fundamentally different picture than the collapse of idealized monoliths. The latter approach suggested that reconnection of tangled magnetic field lines plays a dominant role. Although the interpretations of the results differ, it appears that turbulence can aid disk formation in magnetized media.

Similar results have been found and extended by a number of different groups, using different numerical setups and codes. Joos et al. (2013) find disk formation in magnetized turbulence and interpret this as a consequence of the misalignment effect that turbulence imposes on the forming disks. Gray et al. (2018) ran simulations with similar initial conditions using their Orion 2 AMR, RMHD code, to the Seifried et al. (2013) setups—with a 300M ⊙ initial clump, λ = 2, similar Mach number for turbulence, and most importantly, similar numerical resolution (1.2 AU). These simulations show, moreover, that rotationally supported disks don't form when the initial turbulent velocity field is aligned with the magnetic field.

An analysis of the statistical properties of such disks formed in purely hydrodynamic simulations of forming star clusters can be found in Bate (2018). Although, this study does not include magnetic fields, it confirms the general picture that initial velocity or density fluctuations are sufficient to provide locally the necessary angular momentum for the formation of protostellar disks.

Recent advances in computational techniques now allow, for the first time, multi-scale MHD simulations in turbulent Giant Molecular Cloud spanning scales from the cloud (40 pc) down to a few AU scale on the disk scale (Kuffmeier et al., 2017). The transport of angular momentum in the disks can be followed. The net inward mechanical transport is compensated for mainly by an outward-directed magnetic transport, with a contribution from gravitational torques by spiral waves that is less than the magnetic transport.

Figure 7 shows a selection of disks that form. One sees that disks are often associated with filaments. Because the most powerful parts of the jets are launched from the central regions of the disk, close to the (proto-)star, it is a numerical challenge to probe the efficiency of protostellar outflows in star cluster environments where spatial scale separation between the molecular cloud scale (~10) and the inner disc scale (~1au) spans more than six orders of magnitude. Attempts to probe self-consistently driven outflows out of turbulent cloud cores have been made. Those simulations don't yet probe the innermost disk scales where the most powerful parts of the jet are launched, but the results are most encouraging.

FIGURE 7

Figure 7. Slices in the plane vertical to the mean angular momentum vector calculated for a sphere of 100 au at t = 50 kyr. First row: sink 1 (left), sink 4 (right); second row: sink 5 (left), sink 6 (right); third row: sink 7 (left), sink 9 (right). From Kuffmeier et al. (2017) reproduced with permission © AAS.

Figure 8 shows the evolution of the angle between the total angular momentum vector and the total magnetic field direction at 100 au around 8 different sink particles (stars) in the simulation. One sees that this distribution is essentially random, which is in good agreement with the observations.

FIGURE 8

Figure 8. Evolution of the angle between total angular momentum vector and total magnetic field direction within a sphere of 100 au around the eight different sinks. Note the high degree of random behavior. From Kuffmeier et al. (2017) reproduced with permission © AAS.

3.3.3. Non-ideal MHD: Effects of Grain Evolution

As already noted, for disks to form in ordered, laminar collapse, something must greatly accelerate the normal ambipolar diffusion time scale to reduce the field while gas is still at relatively low density. This may be possible when one considers the fate of small dust grains in magnetized contracting cores.

The large population of very small dust grains (VSGs) of order 10 - a few 100Å dominates the coupling of magnetic field to gas at densities n > 1010 cm−3. The removal of this dust can lead to an increase of the ambipolar diffusion rate of the magnetic field out of the collapsing gas by two orders of magnitude. This allows 50-100 AU disks to form in well ordered collapses (Zhao et al., 2018). The depletion of VSGs arises as a natural consequence of their accretion of mantles (akin to molecular freeze-out) or coagulation. Grain coagulation has been shown to remove small grains ≤ 0.1μm) within a few yr (Ossenkopf, 1993; Hirashita, 2012). The advantage of ambipolar diffusion driven by the loss of VSGs is that it occurs in the envelope of the forming disk which allows the formation of a large scale disk by avoiding catastrophic braking. 3D simulations of this picture show the presence of rings on scales of ~ 40 AU.

Since even subsonic turbulence allows for the formation of large disks in strong B fields, it is unclear at what stage this mechanism of grain removal in cores becomes important or competitive as subsonic conditions in cores are established.

3.3.4. The Earliest Outflows

With the collapse of a turbulent overdense region of a cloud, the magnetic field configuration winds up creating an ever stronger 