<figure><img src="image_1.png" /></figure>

HAL
open science

# Jets from young stars

T. P. Ray, J. Ferreira

► To cite this version:

T. P. Ray, J. Ferreira. Jets from young stars. New Astronomy Reviews, 2021, 93, 10.1016/j.newar.2021.101615. insu-03705255

HAL Id: insu-03705255

https://insu.hal.science/insu-03705255

Submitted on 27 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

# Jets from young stars

T.P. Ray<sup>a,b,*</sup>, J. Ferreira<sup>c</sup>

## ABSTRACT

Jets are ubiquitous in the Universe and are seen from a large number of astrophysical objects including active galactic nuclei, gamma ray bursters, micro-quasars, proto-planetary nebulae, young stars and even brown dwarfs. In every case they seem to be accompanied by an accretion disk and, while the detailed physics may change, it has been suggested that the same basic mechanism is responsible for generating the jet. Although we do not understand what that mechanism is, or even if it is universal, it is thought to involve the centrifugal ejection of matter from the disk along magnetic field lines. For a number of reasons, in particular their proximity and the abundant range of diagnostics to determine their characteristics, jets from young stars and their associated outflows may offer us the best opportunity to discover how jets are generated and the nature of the link between outflows and their accretion disks. Recently it has become clear that jets may be fundamental to the star formation process in removing angular momentum from the surrounding protoplanetary disk thereby allowing accretion to proceed. Moreover, with the realization that planetary formation begins much earlier than previously thought, jets may also help forge planets by determining initial environmental characteristics. This seems to be particularly true within the so-called terrestrial planet forming zone. Here we review observations of jets from young stars which have greatly benefitted from new facilities such as ALMA, space observatories like Spitzer, Herschel and HST, and radio facilities like LOFAR and the VLA. Interferometers such as CHARA and GRAVITY are starting to make inroads into resolving how they are launched, and we can look forward to a bright future in our understanding of this phenomenon when JWST and the SKA come on stream. In addition, we examine the various magnetohydrodynamic models for how jets from young stars are thought to be generated and how observations may help us select between these various options.

## 1. Historical introduction

In the 1950s, Herbig (1951) and Haro (1952) discovered the nebulous patches in Orion that now bear their name. At first these were thought to be actual “cradles” of stellar birth and it was not until the 1960s that Schwartz (1977) realized that Herbig-Haro (HH) objects may be part of an outflow from a young star. This he thought likely because their spectra bore certain similarities to the shock ejecta from supernova remnants.

It was not however until the development of the Charge-Couple Device (CCD), with its vastly superior quantum efficiency in comparison to photographic plates, that the “dots were joined”. The nebulous clouds seen by Herbig and Haro were found to be driven by, or parts of, a highly collimated line emitting jet. At the same time the use of CCDs was coming to the fore, millimeter wave astronomy was undergoing a revolution. The first large millimeter grade dishes, in combination with superheterodyne receivers, started to map star formation regions in single molecular transitions such as the fundamental 2.3 mm $J = 1 \to 0$ rotational line of CO. Hoping to see matter collapsing inwards to form a star, to their surprise the early pioneers discovered matter was actually flowing away: outflow rather than inflow was observed. Often diametrically opposite red-shifted and blue-shifted lobes were seen moving at velocities of tens of kms$^{-1}$ and straddling the young star (Snell et al.,

1980). Since this gas had to be at very low temperatures, otherwise the 2.3 mm CO transition would not be excited, such molecular outflows had to be highly supersonic.

In a number of cases clear associations were seen between the optical outflows and their molecular counterparts. The optical line emission was observed to be strewn along the same axis as the molecular emission and the latter to encase the former (e.g., Snell et al., 1985). As the area that could be surveyed by CCDs grew, with their increasingly larger formats, it became clear that outflows from young stars could extend up to several parsecs, stretching, in some cases, beyond the boundaries of the young star’s associated molecular cloud (Bally et al., 2012; Qiu et al., 2019).

As outflows were observed from a whole host of young stars of different masses and evolutionary phases, certain patterns were seen to emerge. First of all whether an outflow appeared or not depended on whether the young star was surrounded by an accretion disk. No accretion disk implied no jet, although the converse was not necessarily true. It was also found that the strength of an outflow depended on the evolutionary status of the parent star: less evolved stars have more powerful outflows. As we shall see, this is in line with the fact that such stars have stronger accretion rates, all else being equal, and the relationship between outflow and accretion rates (see 5).

<footnote><sup>a</sup> School of Cosmic Physics, Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin D02 XF86, Ireland</footnote>

<footnote><sup>b</sup> School of Physics, Trinity College, Dublin, D02 W272, Ireland</footnote>

<footnote><sup>c</sup> Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote>* Corresponding author at: School of Cosmic Physics, Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin D02 XF86, Ireland.
E-mail addresses: tr@cp.dias.ie, tray@tcd.ie (T.P. Ray).</footnote>
At this point it may be useful to point out how we distinguish between young stars at different evolutionary phases. The standard picture, for solar-like stars, is that their formation begins through the radial collapse of a molecular core, perhaps a parsec across, and that this gives rise to the birth of one or more stars (e.g Krumholz, 2015). In the earliest phase, most of the mass is still contained in the surrounding envelope, the energy output is entirely dominated by accretion onto the central object, and there are many tens of magnitude of visual extinction present. The source at this stage is referred to as a Class 0 (Andre et al., 1993) young stellar object (YSO). Observations shows that such objects already possesses a disk as well as an outflow. The spectral energy distribution (SED) of a Class 0 YSO closely resembles a blackbody with temperatures of at most a few tens of K and a peak flux in the millimeter regime. The SED is dominated by the optically thick envelope.

In the next phase, Class I (Lada and Wilking, 1984), the protostar has either accreted or driven away enough of its envelope that the effects of extinction are considerably reduced. The YSO now emits in the near-infrared (NIR) and its SED is dominated by the contribution from the disk along with one from the envelope. Once the YSO becomes optically visible, its photosphere can be observed, and it is now a T Tauri star (TTS) or, if a few times more massive than the Sun, a Herbig Ae/Be star. T Tauri stars are further sub-divided into those with accretion disks, so-called classical T Tauri stars (cTTS), or those without disks, referred to as weak-line T Tauri Stars (wTTS). These are also known as Class II and Class III YSOs respectively (Lada and Wilking, 1984). Although estimates vary somewhat as to how long an individual young star takes to go through each phase, as a rule of thumb both cTTSs and wTTSs are at most a few million years old. At all stages of the process, the luminosity of the YSO is derived either from accretion or, in more evolved sources, from gravitational contraction. At no time has fusion been ignited in the core. As pointed out previously, jets are only present when a YSO is surrounded by an accretion disk, thus while a number of cTTS have been found to possess jets, jets are not found around bona fide wTTS.

## 2. Physical characteristics of YSO jets

Jets from young stars have been observed at a variety of wavelengths from the X-ray regime right down to the radio band. In this review we will concentrate on the optical and NIR bands as historically this is where much of the effort on understanding their nature has focused.

The radiation from YSO jets is for the most part characterized by line emission, at least from the UV to the sub-millimeter bands. At optical and NIR wavelengths, this is a combination of both permitted, e.g. $H\alpha$, and forbidden lines, e.g [OI]$\lambda$630 nm, [SII]$\lambda$671.6 nm, 673.1 nm, [FeII]$\lambda$1.64 µm, etc. (see Fig. 1). The presence of such forbidden lines immediately implies that electron densities must be low (at most $10^{4-5}cm^{-3}$) in some parts of the jet (Hartigan et al., 1995). Moreover, as also recognized by Schwartz (1977) many years ago, the line emission in Herbig-Haro jets comes from the cooling zone of shocks, rather than from photo-ionization.

This presence of line emission can thus be used to infer fundamental parameters such as a jet's radial velocity, and typical values of around a few hundred $kms^{-1}$ are found. We expect of course tangential velocities to be similar and this is consistent with observed proper motions (e.g., Eisloffel and Mundt, 1992). In fact, the combination of radial velocities, proper motions and known distance to the source, gives us a full 3-D kinematic picture of a jet. Moreover the inclination angle of the outflow to the plane of the sky can be determined in this way.

Now line spectra can also be used to derive a jet's temperature and therefore its local sound speed. As the temperatures found are around $10^4$ K, i.e. the sound speed is approximately $10 kms^{-1}$, YSO jets are thus highly supersonic with Mach numbers, $M_{Jet}$, of around 20–30. Since the opening angle (in radians) of an unconfined supersonic flow $\approx \frac{1}{M_{Jet}}$,

i.e. the sound speed divided by the jet's velocity, we expect values of a few degrees at most. Since these are the values typically found not only in atomic but even in some molecular species (Codella et al., 2007), we can assume jets for the most part are advancing ballistically or at least they are not confined laterally. Note that this only refers to distances of several hundreds of au from the source where the magnetic field is expected to be too weak to play a role in the jet's dynamics (Hartigan et al., 2007). Closer to the source, magnetic fields are expected to dominate the launching and focusing of the jet.

Raga and Kofman (1992) noted the morphological similarity of YSO jets, such as the HH 34 system (Buehrke et al., 1988), to those from Active Galactic Nuclei (AGNs) and in particular they proposed a similar model to explain their knots, as suggested by Rees (1978) for the M87 jet, i.e. the observed knots are so-called “working surfaces” or internal shocks. Such shocks occur as a result of velocity variations in a highly supersonic flow as more rapidly advancing jet material catches up with slower material ahead of it. The effective shock velocity is then determined by the difference in velocity of material rather than by its absolute value. The latter would only be relevant for a light jet flowing into a relatively dense stationary medium, e.g. a clump in the molecular cloud. The existence of working surfaces explains a conundrum: why the indicated shock velocities at knots are often so low (typically a few tens of $kms^{-1}$ in comparison to their bulk speeds (Hartigan et al., 1994).

That is not to say jet shock speeds are always low: in fact, the largest scale structures which resemble giant bows, such as those in HH 34 and HH 212 are associated with very high shock velocities. This is implied, not only through their spectroscopic line profiles but also by the high degree of ionization at their apex, giving rise to, for example, [OIII] emission (Buehrke et al., 1988).

Of course shocks in YSO jets and those in AGN jets differ in a fundamental way apart from the obvious dissimilarities in size and velocity. Shocks in AGN jets are essentially adiabatic, in relative terms the energy losses at their working surfaces are small. In contrast, shocks in YSO jets are radiative, and sometimes a significant fraction of the shock energy can be converted into line emission.

From a diagnostic perspective, and as already mentioned, the presence of line emission offers a distinct advantage over the broad, often bland, continuum seen in AGN jets: using line diagnostics one can readily determine a jet's fundamental parameters such as temperature, velocity, electron density, etc. Moreover, indirectly we can obtain others. For example, using the assumption of charge exchange, i.e. species such as sulfur or oxygen become ionized, or end up in an excited neutral state, through charge exchange with hydrogen, the total density of neutral plus ionized species can be obtained. Modeling then shows (e.g., Bacciotti and Eisloffel, 1999) that YSO jets are mostly neutral with ionization levels typically around 1%–10% (Podio et al., 2009).

If we have the total density then clearly using the jet radius, in combination with its velocity, we can determine the mass flux through the jet. Such fluxes are found to vary widely ($10^{-5}-10^{-8} M_{\odot}/yr$) but in general a number of expected trends are seen. For example, mass loss rates are observed to decrease with advancing evolutionary phase for a young star of fixed mass but also increase with the mass of the protostar (Purser et al., 2016). Since outflow rates roughly correlate with accretion rates, this reflects what is seen for mass accretion rates onto young stars (e.g., Antoniucci et al., 2014).

## 3. What an outflow from a Young star looks like: The case of HH 211

HH 211 was the first outflow from a young star that was found, not in the optical but, in the NIR using a narrow band filter centered on the 2.12 µm line of shocked molecular hydrogen (McCaughrean et al., 1994). Two clear large bow shocked-shaped features are seen pointing to the northwest and southeast of its source (marked with a white cross in Fig. 2). It is close to the young stellar cluster IC 348 and is located
<figure><img src="image_5.png" /><figcaption>**Fig. 1.** The HH 901 and HH 902 outflows as imaged by HST's Wide Field Camera 3. The colors in this composite image correspond to [OI] (blue), Hα (green), and [SII] (red). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)</figcaption></figure>

Source: Credit: NASA, ESA, M. Livio and the Hubble 20th Anniversary Team (STScI).

<figure><img src="image_6.png" /><figcaption>**Fig. 2.** The HH 211 outflow near IC 348 in shocked molecular hydrogen emission at 2.12 µm. The central source is a highly embedded Class 0 YSO the position of which is marked by a white cross. Red shows shocked SiO emission. This is in the form of a highly collimated bipolar jet that appears to be driving the molecular outflow. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)</figcaption></figure>

Source: Credit: Lee et al. (2018) and reproduced with permission of the lead author.

approximately 300 pc away. The source itself is a highly embedded Class 0 YSO and the outflow appears to be very young: a dynamical timescale of only 1,000 years is inferred from dividing its projected length by its tangential velocity of approximately 100 kms⁻¹ (Jhan and Lee, 2016). Following on from its discovery in the NIR, molecular observations at longer wavelengths, using the CO J=1-0 and J=2-1 lines as tracers, also revealed a very collimated bipolar structure (Gueth and Guilloteau, 1999). In fact if we assume the bow features observed in molecular hydrogen emission are where the bipolar jet rams into and shocks the surrounding interstellar medium, then the morphology of the CO emission clearly suggests cavities in the wake of these terminal

shocks. This is in line with shock-entrainment as the formation mechanisms for young, embedded molecular outflows (Gueth and Guilloteau, 1999).

In other molecules, such as SiO, the outflow appears much more collimated than in low rotational transitions of CO (compare Figs. 2 and 3). Species such as SiO are enhanced, by several orders of magnitude, through grain-sputtering in shocks (Nisini et al., 2007). Thus in Fig. 2, where we see SiO emission, we are really observing where the shocks are strong enough to significantly impact on dust grains. Although the HH 211 region is clearly highly embedded, and thus difficult to observe optically, nevertheless the presence of an underlying atomic jet can be inferred from far-infrared observations, e.g. [OI] 63µm and 145.5µm data from Herschel/PACS (Dionatos et al., 2018).

## 4. What radio emission from outflows has to tell us

Early observations of jets from young stars with, for example, the Very Large Array (VLA) showed them to be faint radio emitters with the flux at frequency v, $S_v \lesssim 1$ mJy at around 5 GHz. Their spectra were found to be thermal with $S_v \propto v^\alpha$ and $\alpha$ in the range 0.1 to 1.0 (Anglada et al., 2018). Models approximating the radio emission as a wind of fixed opening angle (typically a few tens of degree) appeared to fit the data well (Reynolds, 1986). Maps of the source often showed the emission to be elongated in the direction of the known optical/near-infrared jet and occasionally more distant radio emission was seen, again along the jet axis. In a number of cases, proper motion measurements could be made and these were found to be close to those derived from line emission at optical/near-infrared wavelengths, i.e. a few hundred kms⁻¹ (e.g., Rodríguez et al., 1989). Assuming total (radial plus tangential) jet velocities to be similar to those of their optical counterparts, observed radio fluxes can be used in combination with
<figure><img src="image_7.png" /><figcaption>**Fig. 3.** A comparison of NIR $H_2$2.12 µm (white contours) emission and CO J=2→0 emission (in green) for the HH211 blueshifted SE lobe. The CO emission, which covers the velocity interval of −11.2 to −1.7 kms⁻<fnref n="1" /> was mapped using the SMA. Some HCN J=3→2 knots (red contours) near the apex are also seen. At more blueshifted velocities the CO emission is closer to the axis of the SiO jet observed in Fig. 2. Effective beam sizes are shown in the bottom right corner. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)</figcaption></figure>

Source: Credit: Tappe et al. (2012). Reproduced with permission of the authors and
©AAS.

known distances to derive mass loss rates. While these were typically
found to be considerably lower than those inferred from optical line
emission, the different values can be reconciled by noting that the latter
represents both the atomic (neutral) and ionized component whereas
the radio emission only comes from the ionized component. Thus mass
loss rates derived from radio emission are effectively lower bounds and
typically suggest 1%–10% of a jet is ionized.

Since the early observations of YSO jets described above were
made, the sensitivity of interferometers such as the VLA in the US or
MERLIN in the UK have increased enormously. This is as a result of
the broader bands now observable using fiber optic links.<fnref n="1" /> For exam-
ple observations carried out in the so-called C-band (from 4–8 GHz)
routinely encompass 1 GHz either side of the band’s central frequency.
Such an increase in sensitivity is leading to a resurgence of interest
in mapping YSO outflows and even carrying out large surveys of star
formation regions (Tychoniec et al., 2018). Moreover the availability
of very low frequency arrays such as LOFAR are beginning to make an
impact in this field. For example in the case of thermally emitting jets,
the expected turnover at low radio frequencies has been found thus
allowing us an independent estimate of the electron density and the
total mass of ionized gas (Coughlan et al., 2017).

Perhaps however the most intriguing finding in the past few years
is that a growing number of outflows have non-thermal emission or,
at least, a non-thermal component mixed in with a thermal one. The
classic example is the HH80/81 outflow associated with the massive
young star IRAS 18162-2048. This jet stretches over 5 pc and while
the radio spectral index towards the center of the outflow suggests
thermal emission, a non-thermal index is derived further out. Moreover,
the jet has been observed to be linearly polarized, a clear signature of
synchrotron emission (Carrasco-González et al., 2010). In addition, soft
thermal X-ray emission has been observed ahead of the non-thermal
radio peak (Rodríguez-Kamenetzky et al., 2019) suggesting a forward
shock that heats gas up to 10⁶ K and a reverse shock (Mach disk) where
the relativistic particles are likely to be accelerated. Other cases of
non-thermal emission have been found (e.g., Ainsworth et al., 2014) in
outflows from solar-like young stars where the velocity of the outflow is
considerably lower than the 1,000 kms⁻<fnref n="1" /> reached in HH80/81. Clearly,

and in conditions very unlike those seen in supernova remnants for
example, it is possible to accelerate particles to high energies. How
this is done is still under debate (e.g., Bosch-Ramon et al., 2010)
but an interesting prediction of such models is that YSO jets should
produce γ-ray emission through the Inverse Compton process. Recent
observations, based on a time base of 10 years with the Fermi-LAT
satellite, show the HH80/81 jet may well be a source of γ-ray photons
up to 1 GeV (Yan et al., 2019).

Finally it should be remarked that low frequency observations
of non-thermal emission can offer us another method of measuring
the magnetic field in an outflow. This technique relies on the fact
that an outflow lobe can consist of a mixture of relativistic particles
in combination with thermal electrons. The latter, depending on the
electron density, de-collimates the forward-beamed radiation from the
non-thermal electrons, giving rise to a decrease in flux or effectively a
turnover in the synchrotron spectrum. If the electron density is known
from line diagnostics (e.g. using the ratio of the red [SII] doublet)
then the magnetic field can be derived. This, so-called Razin Effect,
may have been seen in a YSO outflow for the first time (Feeney-
Johansson et al., 2019). The derived magnetic fields are in line with
what is expected, and as found, by totally independent line diagnostic
techniques (Hartigan and Wright, 2015).

5. Linking accretion with ejection

Accretion onto a young star can be measured in a variety of
ways (Antoniucci et al., 2014) and an obvious link to look for is
between the level of accretion onto a star and the strength of its
outflow. This is perhaps an appropriate juncture to explain in more
precise terms what is meant by the term “outflow”. The typical veloc-
ities, i.e. the sum of the radial and tangential velocities, of material
observed in optical/NIR lines for example, is a few hundred kms⁻<fnref n="1" />,
i.e. comparable to the escape velocity from the young star. In other
words this component of the outflow clearly originated from deep in
the gravitational well. In contrast molecular outflows, at least in low
excitation lines such as the fundamental rotational transition of the
CO molecule at 2.3 mm, are seen to be slowly moving, i.e. roughly
at a few tens of kms⁻<fnref n="1" />. Despite their slow motion, the total mass of
molecular gas set in motion is considerable. This was recognized very
early on (e.g., Stojimirović et al., 2006) with the discovery that while
the molecular outflow might consist of several solar masses of material,
the star at the center of the outflow could be less than a solar mass.
The inference from this is either star formation is very inefficient in
the sense that many solar masses of material must be funneled in and
away from a protostar for the latter to acquire even a mass comparable
to the sun or alternatively that the high velocity outflow from the young
star is very efficient in driving molecular material around it. Modeling
shows the latter statement is more likely to be correct (Richer et al.,
2000). In other words, to borrow a phrase from Theology used first
by Saint Thomas Aquinas, the high velocity jet is, at least to a first
approximation, the “prime mover”.

Returning then to the question of the link between accretion and
outflow: we can now define the outflow in at least a more consistent
way by looking at the atomic line emission. A number of authors have
investigated this problem using both values for the outflow determined
from extended, i.e. distant, spatially resolved sections of the jet (Harti-
gan et al., 1994) as well as using line emission close to the source that
is identified, through its kinematic and other physical properties, as the
spatially unresolved jet component (Hartigan et al., 1995). Moreover,
aside from optical and NIR lines, mid- and far-infrared atomic lines
can also be used as a diagnostic for more embedded sources and
outflows (Watson et al., 2016). Estimating accretion onto a star can
be done in a number of ways. All pre-main sequence stars derive most
of their luminosity from gravitational contraction and, for those with
accretion disks, the depositing of matter onto the surface of the star
(see Fig. 4). In the latter case matter is deposited along the star's

<footnote n="1">¹ In the past microwave links or coaxial cable were used with very limited bandwidth.</footnote>
<figure><img src="image_8.png" /><figcaption>Fig. 4. Accretion and winds around YSOs. The standard paradigm is that accretion, and hence radially inwards motion, occurs in the first instance through the disk at least up until the so-called co-rotation radius. To do so material must lose angular momentum perhaps as a result of the wind possessing spin. From the vicinity of the co-rotation radius accretion then proceeds along magnetically guided columns onto the star. Before impacting on the stellar surface, the accreted material has to be slowed down through a shock. Emission from the columns is responsible for the time-varying, complex emission line profiles observed at the source.</figcaption></figure>

Source: Credit: Hartmann et al. (2016) and reproduced with the permission of the lead author and ©ARA&A.

magnetic field lines which can essentially act as funnels (Edwards, 1997). The gravitational contraction timescale, also known as the Kelvin-Helmholtz timescale is relatively long, i.e. tens of millions of years for a solar mass star, whereas the accretion timescale is much shorter. In other words the accretion luminosity dominates over any contribution from contraction if the protostar is young enough. Thus the stellar luminosity:

$$L_* \approx GM_* \dot{M}_{acc}/R_* \quad (1)$$

is a good measure of accretion, $\dot{M}_{acc}$, assuming the stellar mass, $M_*$, and the radius, $R_*$, is known. In the case of the optically visible classical T Tauri stars, or their more massive counterparts the Herbig Ae/Be stars, the line profiles from their photospheres in combination with a measure of their photospheric luminosities, allowing for extinction, can also be used as proxies to deduce an accurate accretion luminosity (Antoniucci et al., 2014). An alternative approach, which produces similar results in the case of classical T Tauri stars, is to measure their UV excess (Venuti et al., 2015). This excess, not seen in the case of main sequence stars of similar spectral types, comes from radiation produced by material accreted onto the star as it goes through a strong shock close to the stellar surface (Königl, 1991) and see Fig. 4. In addition to the UV excess, there is a continuum from the accretion shock that extends into the optical. While the UV emission appears as a “bump” on the Wien-like part of the spectral energy distribution, the continuum extends into the photospheric spectrum of the star itself. This means that normal photospheric absorption lines are “filled-in”, a phenomenon known as line veiling (Calvet and Gullbring, 1998). Veiling then can be used as an accretion measure. Finally it should be mentioned that a number of emission lines can be used effectively as proxies for accretion once the correlation between line luminosity, allowing of course for any extinction effects, and accretion luminosity is determined (e.g. Alcalá et al., 2017).

Ultimately of course accretion (Fig. 4) is thought to occur through the disk surrounding the star. In order for disk material to move radially inwards, angular momentum must be removed from it first. Disk winds, which can subsequently be collimated into jets, can do this if they are ejected with spin. Alternatively viscosity in a disk can redistribute angular momentum by moving some material inwards, and allowing it to accrete, while at the same time transporting the excess angular momentum outwards. In this way disks actually grow with time, i.e. viscous spreading occurs. There is some evidence that this might occur, e.g. Najita and Bergin (2018) find that most Class II (T Tauri) disks are typically larger than those from an earlier evolutionary phase. There is,

however, a large spread in sizes and many Class II disks are relatively small. Ultimately disk sizes may be controlled through a combination of disks winds, photo-evaporation and the effects of newly formed planets. Models of disk accretion, i.e. the possible means to transport the disk angular momentum away, will be reviewed below, in the section on disk winds.

## 6. The different scenarios for YSO jets

Discussions about the physical origin of YSO jets followed the same pattern as in AGN jets. Indeed, jets could either be powered by the central object (stellar winds) or linked to the circumstellar accretion disk (disk winds). Moreover, besides finding a model for jet production, one should also seek a mechanism regulating the spin evolution of a protostar. Indeed, it was soon realized that Class II T Tauri stars were slow rotators with typical periods of 8 to 10 days, i.e. only about 10% of their breakup speed (Bouvier et al., 1986; Bertout, 1989). This is one of the yet unsolved mysteries of star formation: where does the stellar angular momentum go during the embedded phase (Classes 0 and I)? In addition, it appears that as long as the circumstellar disk is present (Class II), the protostar maintains a constant period, a situation referred to as “disk locking” (see Gallet and Bouvier 2013 and references therein). This is another mystery since the protostar is actually contracting and accreting material with large angular momentum so that it should actually spin up. Somehow, jets might also be an attractive way to get rid off the stellar angular momentum. Below, we review separately each class of models, as they co-existed in the literature until recently.

### 6.1. Stellar winds

The first scenario invoked to explain the collimation of molecular outflows is based on a de Laval nozzle effect, due to a density gradient in the circumstellar medium, acting on a stellar wind (Canto, 1980). This is the YSO equivalent of the Blandford and Rees (1974) twin exhaust model for AGN radio jets. However, it was also soon realized that YSO jets carry away far too much thrust and that they could not be radiatively driven (DeCampli, 1981). The other way to produce a stellar wind is to rely on magnetic acceleration. It was also known, since Chan and Henriksen (1980) and Blandford and Payne (1982), that large scale magnetic fields would not only accelerate loaded material, but also provide a hoop stress leading to self-confinement, conveniently taking care of the jet collimation issue.
<figure><img src="image_9.png" /><figcaption>**Fig. 5.** Sketch of the magnetized star-disk interaction within the Accretion Powered Stellar Wind model. The stellar magnetic field truncates the disk at $R_*$, assumed to be located inside the corotation radius, $R_{cc}$. The disk may have its own magnetic field, leading to the launching of an outer disk wind. Inside $R_*$, disk material is lifted towards the star but a fraction of it is ejected away.</figcaption></figure>

Source: Credit: Matt and Pudritz (2005).

The idea of magnetic braking by a stellar wind is now very familiar (Mestel and Spruit, 1987). Schatzman (1962) already pointed out that if gas emitted from a star is kept corotating with the star by the magnetic torques out to large distances, it will carry off far more angular momentum per unit mass than gas retaining the angular momentum of the stellar surface. This is the reason why using a magnetized stellar wind to spin down a protostar was one of the first ideas to come to mind. Unfortunately, the low rotation rate of T Tauri stars forbids any efficient centrifugal acceleration and some ad-hoc Alfvénic turbulent pressure must be invoked to lift material away from the protostar (Hartmann and MacGregor, 1982). However, such a model had severe difficulties in explaining YSO jet mass losses and this idea has been discarded for quite some time.

The revival of stellar winds occurred with the Accretion Powered Stellar Wind (APSW) model of Matt and Pudritz (2005). In this scenario, the stellar magnetic field (assumed dipolar) threads the accretion disk and leads to the formation of accretion funnel flows (see Fig. 5). Because of magnetospheric stresses, the disk material is deflected away from the equatorial plane and falls towards the protostar along the dipolar field, leading to an accretion shock. The main assumption of the model is that a fraction of the mass and its associated accretion energy is actually deflected away along open polar stellar field lines, giving rise to an enhanced stellar wind. In principle, such a model is consistent with the observed accretion-ejection correlation and provides efficient magnetic braking. Allowing for protostellar contraction, Matt et al. (2012) computed the evolution of stellar spin subject to an exponentially decreasing rate of accretion with time and showed that, for some parameters, it is possible to achieve low rotation rates over typical TTS lifetimes (1 to 3 Myrs).

This would require however mass loss rates in APSWs comparable to those inferred in YSO jets, implying that YSO jets are mostly made of material ejected from the protostar, in contradiction with jet kinematics (Ferreira et al., 2006; Cabrit, 2007). Moreover, the model has a theoretical drawback summarized in the expression “you can’t have your cake and eat it”. APSWs tap accretion energy, so the more massive

<figure><img src="image_10.png" /><figcaption>**Fig. 6.** Sketch of the X Wind in a typical Class II system. In this star-disk interaction model, the stellar magnetic field is the only relevant one, truncating the disk at a radius $R_X$ assumed to be the co-rotation radius. Stellar field lines inside $R_X$ enforce co-rotation for some disk material, leading to infalling funnel flows. Field lines that would have thread the disk beyond $R_X$ are assumed to have either opened up, and become disconnected from the star leading to some highly localized wind from the disk (X-wind), or to have inflated and still remain connected to the star (giving rise to a mass-less magnetized corona).</figcaption></figure>

Source: Credit: Shu et al. (1994).

they are the more energy they need. Each watt an APSW consumes, is a watt that is not being radiated away in the UV in the accretion shock (assuming a 100% conversion efficiency). Since UV luminosity provides a proxy for the disk accretion rate, tapping this energy while maintaining the same UV emission implies that even more mass is actually being accreted onto the protostar. This enhances the spin up torque which, in turn, would require a larger spin down from the APSW. It can be proven that such a situation does not converge for some objects (Zanni and Ferreira, 2011).

APSWs are thus unlikely to be the main mechanism providing the required spin down torque for TTS. Although they cannot represent the main component of YSO jets, they probably fill in the jet axis, leading to some interaction with the surrounding jet medium (see, for example, Meliani et al., 2006).

## 6.2. Outflows from the star-disk interface

In 1988, F. Shu and collaborators proposed the X-celerator mechanism for YSO jets, where the protostar is assumed to rotate at break up speeds (Shu et al., 1988). The double goal sought by this scenario is to explain YSO jets while simultaneously braking the spin of the protostar. The disk material would reach the star at the equatorial plane and enforce the dipolar stellar magnetic field to open up, leading to centrifugal ejection from the disk mid-plane. But this early model was put aside as it became clear that TTS only rotate at typically 10% of their break up speed.

The X-wind scenario, as put forward by Shu et al. (1994) with its various extensions (for example Cai et al. 2008 and references therein), is an enhanced version of the same idea (see Fig. 6). More precisely, X-winds should be the dominant component in YSO jets (in terms of appearance, mass flux and power) while carrying away enough angular momentum so that the star is not being spun up by the accreting material (zero net torque condition, protostellar contraction neglected). In the X-wind paradigm, the stellar magnetic field is the only relevant large scale magnetic field present in the system. It truncates the accretion disk at the corotation radius<fnref n="2" /> leading to the formation of

<footnote n="2">² The corotation radius is the cylindrical distance where the local Keplerian rotation rate is equal to the stellar rotation rate.</footnote>
accretion funnel flows. However, a fraction of the stellar magnetic field is assumed to open-up and diffuse beyond the co-rotation radius. Then, from what is assumed to remain a point-like source around the co-rotation radius, a magnetized wind is expected to form and carry away all the angular momentum of the accreting material and a significant fraction of the incoming mass.

Because of its elegance, this steady-state jet model had a huge impact on the YSO community. However, it turned out to suffer several severe drawbacks. First, YSO jet kinematics were shown to be inconsistent with published X-wind calculations. For instance, the observed range and steep radial decline of poloidal speeds could not be explained with the chosen Alfvén surface shape (Ferreira et al., 2006; Cabrit, 2007). Second, the conditions required on the underlying disk (used as a boundary condition) turned out to be quite extreme. According to this scenario, material accreted onto the star via funnel flows leaves its angular momentum at the base of the funnels, i.e. inside the co-rotation radius. In order for the X-wind to carry away all the excess angular momentum, some unknown outward (viscous) transport process must be invoked. This is rather unrealistic and probably explains why X-winds have never been obtained so far in MHD star–disk numerical simulations (Ferreira and Casse, 2013).

2-D and 3-D numerical simulations of the star–disk interaction show that the disk truncation radius is always found to be smaller than the co-rotation radius when accretion takes place (e.g. Romanova et al. 2002, Bessolaz et al. 2008, Zanni and Ferreira 2009, Romanova et al. 2012 to cite only a few). However, some of these simulations do show the presence of ejecta from the disk truncation radius, although these events cannot be described within a steady-state MHD jet framework (Goodson et al., 1997; Romanova et al., 2009; Zanni and Ferreira, 2013). The latter authors showed that these magnetospheric ejections (referred to as MEs and see Fig. 7) have two major effects. First, they provide a strong spin down torque on both the star and the underlying disk, so that the combination of MEs and the inner stellar wind could indeed lead to the long awaited spin down torque on the star. Second, since these MEs are plasmoids resulting from reconnection episodes, they propagate almost ballistically and could not represent, alone, the YSO jet phenomenology. At best, they would be an additional, intermittent, plasma component flowing inside a narrow channel delimited between two quasi-steady MHD flows, namely the inner (on axis) stellar wind and the outer disk wind.

## 6.3. Disk winds and their link to accretion

For quite a long time after the seminal paper of Shakura and Sunyaev (1973), the theory of accretion disk was solely based on hydrodynamics and magnetic fields played no role. The plasma in a circumstellar accretion disk was assumed to behave like a neutral fluid, subjected to some turbulence that would (hopefully) provide the required anomalous viscosity. On the other hand, people working on astrophysical jets soon realized that self collimation can only be obtained if large scale magnetic fields are present, requiring the use of MHD (Chan and Henriksen, 1980; Blandford and Payne, 1982). This led Pudritz and Norman (1986) to suggest bipolar outflows are magnetically driven winds emitted from the surface of the molecular disk, providing a torque and leading to accretion. Thus, instead of being evacuated radially by some yet unknown turbulence, the disk angular momentum is transported vertically in the bipolar jets.

The first tentative attempt to couple the Blandford and Payne (1982) MHD jet model to the outer regions of a circumstellar accretion disk was made by Königl (1989). Indeed, since jets exert a torque on the underlying disk, both its radial and vertical structures must be affected. Moreover, bipolar jets require the existence of a large scale vertical magnetic field threading the disk everywhere or, more likely, only in the innermost region. This early work was then followed by many more efforts, of increasing complexity as regards disk-jet interrelations, (Ferreira and Pelletier, 1993, 1995; Li, 1995; Ferreira, 1997; Casse and

Ferreira, 2000; Ogilvie and Livio, 2001; Jacquemin-Ide et al., 2019) or in the inclusion of non-ideal MHD effects within the disk (Wardle and Königl, 1993; Li, 1996; Salmeron et al., 2011).

In parallel to these analytical approaches, huge efforts have been made to provide numerical MHD simulations of jets emitted from accretion disks. The first attempts did not take the disk into account and computed the jet acceleration and propagation with some fixed boundary conditions at the base (Ouyed and Pudritz, 1997a,b; Krasnopolsky et al., 1999, 2003; Ouyed et al., 2003). Such a “platform” approach is very promising for studying the jet collimation issue, allowing for instance to relate jet asymptotic properties to the radial distributions of jet mass loss and disk magnetic field that are imposed at the base (Anderson et al., 2005; Pudritz et al., 2006; Fendt, 2006), up to observable scales (Staff et al., 2010; Ramsey and Clarke, 2011; Staff et al., 2015). However, since the disk is not computed there is no feedback from the jet on the underlying disk structure, leaving thereby too much freedom on the imposed radial distributions (and no constraint on any possible inconsistency). The first MHD simulations that computed simultaneously the accretion disk and their jets are due to Casse and Keppens (2002, 2004). Building upon the analytical results on the disk-jet interrelations found by Ferreira (1997), they obtained super fast-magnetosonic jets using a so-called $\alpha$ prescription for the turbulent disk. Their results have then been extended by other groups using different MHD codes (Zanni et al., 2007; Tzeferacos et al., 2009; Murphy et al., 2010; Sheikhnezami et al., 2012; Tzeferacos et al., 2013; Stepanovs and Fendt, 2016). Fig. 8 provides an illustration of such findings. More recently 3-D MHD simulations of accretion and jet launching have been performed not only from single YSOs but also from disks in binary systems (Sheikhnezami and Fendt, 2015, 2018). Here instabilities in the disk and jet are seen such as disk spiral arms and warps or jet precession. In turn this may lead to enhanced accretion rates.

The main caveat of these $\alpha$-disk simulations is their possible inconsistency with MHD turbulence, since all anomalous transport coefficients (viscosity, magnetic diffusivity) are parametrized. Indeed, since the discovery of the Magneto-Rotational Instability or MRI by Balbus and Hawley (1991), it is known that magnetized accretion disks can sustain turbulence that acts effectively as a viscosity (see Hawley et al. 1995, Stone et al. 1996, Balbus 2003 and references therein). It thus became clear that magnetic fields were also required in accretion disks, regardless of the jet phenomenon. MHD turbulence in a vertically stratified disk needs to be properly followed to make sure that simulations converge. Moreover, when a vertical magnetic field is included (instead of a toroidal one), mass loss is systematically observed from the disk surface and an outflow is obtained. As a consequence, the size of the computational domain needs to be large enough so that boundary conditions do not affect (or at least not too much) the outcome of the simulation. These are the main reasons why early shearing box simulations could hardly address MRI with a non zero net magnetic flux (see for instance Lesur et al. 2013, Fromang et al. 2013, Bai and Stone 2013a and references therein). Therefore, for quite a long time, the main focus of MRI studies was only the measurement of the Shakura–Sunayev $\alpha$ parameter describing the turbulent viscosity and the influence of non-ideal MHD effects.

The presence of MRI relies on a small degree of ionization being present in the disk, so that in the outer cold regions the magnetic field is essentially coupled to the neutrals through ion-neutral collisions. It turns out that the MRI becomes quenched in the presence of non-ideal MHD effects, such as ambipolar diffusion (ion-neutral drift) and the Hall effect (ion–electron drift), both terms that are present whenever the plasma is not ionized enough (Gammie, 1996; Bai and Stone, 2013b; Lesur et al., 2014). In addition expected sources of ionization, for example low and moderate energy galactic cosmic rays ($E_{CR} \le 1$ GeV), cannot penetrate very far into a disk and may even be excluded by the strong magnetosphere of a young star (Cleeves et al., 2013). That said high energy particles generated by the young star itself may at least
<figure><img src="image_11.png" /><figcaption>**Fig. 7.** Temporal evolution of the periodic dynamical processes leading to magnetospheric ejections. Time is given in units of the rotation period of the star, background colors are logarithmic density maps, white solid lines are sampled field lines, blue vectors show the poloidal speed. The yellow solid lines allow us to follow how the system evolves with time. As long as these field lines remain connected to both the star and the disk, angular momentum is being transferred from these two sources to plasma that is magneto-centrifugally ejected (as in a slingshot). This leads unavoidably to magnetic reconnection and plasmoid ejection, channeled between the inner stellar wind and the outer disk wind. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

*Source:* Credit: Zanni and Ferreira (2013).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 8.** Close-up views on the disk of magnetized disk winds. **Left:** Self-similar solution of an isothermal accretion-ejection structure (Ferreira, 1997). Colors show the density, black solid lines show the streamlines and dashed lines isocontours of the total velocity. In the resistive MHD disk (violet-blue region), matter is mostly accreting and rotating (rotation slows down vertically). Above the disk, plasma flows along the field lines (ideal MHD regime) and is progressively accelerated. **Right:** Final output of a numerical simulation made by Stepanovs and Fendt (2016). Colors show the density, thin black lines the poloidal magnetic field. The lines mark the disk surface (black), the sonic surface (red) the Alfvén surface (white), and the fast surface (dashed white). It is striking to recover (away from the axis) conical critical surfaces in the final state, as in analytical self-similar studies, from a time dependent simulation. With kind permission of the authors. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)</figcaption></figure>

assist ionization in the inner regions of the disk (Rodgers-Lee et al., 2017). Nevertheless, in regions that are both too cold for collisional ionization and too dense for cosmic ray or radiative UV-X ionization, such as the disk mid-plane, no MRI can be triggered. As a consequence no “viscously” driven accretion can take place, leading to the concept of “dead zones” (Gammie, 1996). Early simulations did not include non-ideal MHD effects, but when included these confirm the difficulty of providing viscosity through MRI generated turbulence (e.g. Bai, 2015), or even directly using multi-fluid simulations (O’Keeffe and Downes, 2014).

However, as already pointed out in Gammie (1996), MRI may remain active within thin layers close to the disk surface, thanks to cosmic rays ionization for instance. The existence of equatorial “dead zones” with upper and lower layers of MRI-driven accretion were thus thought to be characteristic of YSO disks, leading possibly to planet formation. This paradigm however proved unfruitful when it was realized that layered MRI-driven accretion could not properly account for the observed accretion rates in TTS (Bai, 2011). Moreover, high spectral resolution observations, at least of the outer zones (tens of au) of disks with ALMA also show no evidence for the expected turbulence (Flaherty et al., 2018). Thus both simulations and observations suggest one should look for an alternative to MRI generated turbulence as the source of viscosity, and hence accretion, in the outer parts of accretion disks.

One way out of this dilemma is to rely on an additional mechanism, such as magnetized winds, to redistribute angular momentum. Moreover, magnetized winds are naturally obtained when a large scale vertical field is present in the disk (Suzuki and Inutsuka, 2009; Fromang et al., 2013; Bai and Stone, 2013a). Such a situation has been carefully avoided in early MRI numerical studies, because of the unfortunate tendency of these codes to crash. Global 3-D MHD simulations are thus better suited to address MRI evolution with vertical $B_z$ fields, but achieving reliable simulations of this kind is a herculean task, that has only been tackled quite recently (Suzuki and Inutsuka, 2014; Gressel et al., 2015; Béthune et al., 2017; Zhu and Stone, 2018; Wang et al., 2019).

A fair statement on the current consensus regarding YSO disks, is that a large scale vertical magnetic field is now invoked everywhere, not only to explain jet formation but also to explain accretion onto the central star (Bai, 2011; Bai et al., 2016; Scepi et al., 2018). In some sense, magnetized winds and jets are the main agent allowing accretion in YSO circumstellar disks and not, as thought for a long time, a mere epiphenomenon (e.g. Pudritz and Ray, 2019). The presence of such a large scale field could also be responsible for features such as gaps and rings seen with ALMA (Suriano et al., 2018).

The lack of ionization beyond 1 au in YSO accretion disks has triggered a wealth of studies on the effect of non-ideal MHD terms
on the development of the MRI and this is still a very active area of research. For instance, recent stratified models, that include Ohmic and ambipolar diffusion only, have shown that the disk could become essentially laminar (i.e. not turbulent) in the mid-plane, with a strong magnetically driven outflow launched at the disk surface (Bai and Stone, 2013a; Gressel et al., 2015 and references therein). The effects of including the Hall effect (which may become dominant between 1 and 10 au) remain uncertain, with some discrepancies between single fluid (Lesur et al., 2014; Béthune et al., 2016) and multi-fluid MHD simulations (O’Keeffe and Downes, 2014). However, while relevant in the outer regions of protostellar accretion disks, innermost disk regions as well as disks around active galactic nuclei and X-ray binaries are ionized enough to ensure non-ideal MHD effects vanish. Under these circumstances, one needs to rely on another source for diffusion, which can only be MHD turbulence itself.

Although a large scale vertical magnetic field remains the unavoidable, albeit invisible, main ingredient for jet launching, thermal effects could still play a major role. As first proposed by Casse and Ferreira (2000), UV-X illumination from the central regions could heat up the surface layers of the outer parts of the accretion disks, leading to strong magneto-thermal disk winds. Such a situation could be generic to accretion disks threaded by a large scale vertical field, in AGN (Fukumura et al., 2010), YSOs (Bai et al., 2016; Wang et al., 2019) and X-ray binaries (Chakravorty et al., 2016). It is slightly ironic to realize that such a scenario has already been envisioned to explain the ubiquitous low velocity emission of forbidden lines in YSOs (Kwan and Tademaru, 1995; Kwan, 1997). The novelty is however the growing realization of the existence of a magnetic flux threading the whole accretion disk. Such a field could be a relic of the original infall, modified by the complex interplay between advection and diffusion (see Takeuchi and Okuzumi, 2014 and references therein).

### 6.4. Are jets the prominent aspect of a multi-component structure?

In essence then, the current consensus seems to favor magnetized disk winds as the source of YSO jets. Such jets would be generated from the innermost disk regions down to the disk truncation radius, where the disk magnetization is expected to be highest. This launching zone is expected to exhibit distinct features that could be observationally tested (Combet and Ferreira, 2008). There are however several aspects that cannot be easily explained by these winds: jet asymmetries (Podio et al., 2011; Melnikov et al., 2009), the range in jet time scales from months to hundred of years Hartigan and Morse, 2007; Agra-Amboage et al., 2011), the presence of a hot inner flow (Edwards et al., 2003, 2006), a possible steady component seen in X-rays (Güdel et al., 2005; Skinner et al., 2011; Bonito et al., 2011). Moreover, disk winds are unable by nature to spin down the central star, which is also an absolute requirement. The global picture is therefore more complex than just magnetized disk winds.

A simple explanation could be that jets are only one component of a much broader, multi-component, MHD structure harboring:

* i a warm stellar wind on the axis, emitted from the star and partly powered by an accretion shock;

* ii surrounded by a radially stratified disk wind (going from fast inner jets to outer massive and slow magneto-thermal winds),

* iii with time dependent magnetospheric ejections at the star-disk interface, braking down the star and producing internal jet shocks.

## 7. Distinguishing between models for jet launching: The hunt for jet rotation

Of particular interest as a means of testing rival models for jet launching is the observation, or lack thereof, of jet or outflow rotation. As emphasized already in the previous section, if jets are launched from

a disk initially as an MHD wind, then they must carry away angular momentum. Of course this is true even in the case of stellar winds, the crucial difference is that the degree of angular momentum transported is a function of the launch radius: winds launched at larger radii, normally carry away more angular momentum than those launched from smaller ones. This gives rise to an important difference between outflows/winds launched say at the co-rotation radius, e.g. in the X-wind model, in comparison to those generated at distances of several au in a disk. The latter will contain more angular momentum flux per unit mass. In fact for any reasonable density profile, a protoplanetary disk will contain more angular momentum per unit mass in its periphery than towards its inner radius. There is a clear analogy here with the Solar System where the Sun contains the bulk of the mass but the planets contain most of the angular momentum.

Observations of jet rotation can therefore give us clues as to where a jet was launched from. Early searches for rotation were made with the Space Telescope Imaging Spectrograph (STIS) on HST (Bacciotti et al., 2002; Coffey et al., 2007). These were done as close as possible to the young star itself, in order to minimize any chances of a rotation signature being confused with, for example, the “buffeting” effects of its immediate environment. High spatial resolution is considered important for another reason: typically one is trying to distinguish differences in velocity laterally across a jet, the width of which is at most a fraction of an arcsecond (Dougados et al., 2000).

The need to have an unobscured look close to the source implied classical T Tauri stars with jets, such as DG Tau, are optimum targets. In a MHD disk wind model, for example, where the bulk of the jet is launched within approximately 1 au, the jet can achieve a poloidal velocity of a few hundred kms⁻¹. Correspondingly we would expect velocities differences, from one side of the jet to the other, of at most a few tens of kms⁻¹ as a result of rotation. Although a number of claims of rotation have been made using HST, alternative explanations for the findings have been proposed including precession (Cerreira et al., 2006) and asymmetrical jet shocks (De Colle et al., 2016). Moreover, the observations were in some cases confusing with results indicating, for example, that the jet was rotating in the opposite sense to the disk (e.g., Louvet et al., 2016). Such findings suggest that either optical/near UV measurements with current instrumentation are not an appropriate means of detecting rotation in jets or alternatively the assumption of a steady-state wind is not correct.

Another means of achieving high spatial resolution, and even probing the outflow region close to the source where extinction might be large, is to use millimeter interferometry. Using this approach, observations of rotation have been claimed for example in the large-scale molecular outflow from Orion Source I with ALMA (Hirota et al., 2017). Perhaps the best example of rotation observed to date in a molecular outflow is in HH 212 shown in Fig. 9. Here ALMA detected velocity gradients across the outflow within approximately 100 au from the embedded YSO. While the sense of rotation of all of the knots is the same, and mirrors that of the disk itself, the amount of rotation observed is small. In fact by assuming conservation of angular momentum along streamlines (Anderson et al., 2003) one can deduce that the outflow originates from within ≈ 0.05 au. Such values are consistent with either a close-in disk wind or possibly even an X-wind.

Evidence for jet rotation have been found in a number of other outflows. These include DG Tau (Agra-Amboage et al., 2011), TMC 1A6 (Bjerkeli et al., 2016) and HH 30 (Louvet et al., 2018). In a number of cases such observations suggest that the outflow originates on scales of at most a few au.

## 8. Exploring the launch zone using optical/NIR interferometers

As already stated, we still do not know, with a high degree of certainty, if jets from young stars are launched in the form of a very strong stellar wind, from the zone where the magnetosphere of the star interacts with the circumstellar disk, from the disk itself or from
<figure><img src="image_13.png" /><figcaption>**Fig. 9.** ALMA observations of the launching region of the HH 212 outflow. On the left we can see knots traced in SiO line emission while in red the continuum emission due to dust is observed. The latter traces the disk associated with the HH 212 source. The sense of rotation of the disk is shown on the right and the same sense of rotation is seen in the knots to the north (N) and south (S) with blueshifted SiO emission to the right of each knot and redshifted SiO emission to the left. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)</figcaption></figure>

Source: From Lee et al. (2017) reproduced with permission © Springer Nature.

a combination of all these zones. While we can get some clues from observing the region where the jet reaches its terminal velocity (i.e. at tens of au) and, in particular, from measuring its angular momentum content as explained in Section 7, nevertheless imaging a jet on sub-au scales is clearly the best approach if one is to decisively determine how, and from where, they are launched. Moreover, it is not just spatial but velocity information that is required since different models predict very different kinematic behavior in the vicinity of the launch zone. This type of data can only be obtained using optical/near-infrared spectro-interferometry, where not only the continuum from the disk but also line emission is spectrally and spatially resolved.

Until recently a combination of technical challenges ensured that optical/NIR interferometers could only be used if the target is very bright (e.g., Berger et al., 2012). Major improvements in their design however have pushed the limiting magnitudes to the point where it is now possible not only to observe the immediate circumstellar environment of moderately bright Herbig Ae/Be stars (intermediate mass young stars) in the optical/near-infrared (Caratti o Garatti et al., 2015) but also cTTS (Gravity Collaboration et al., 2017). Clearly to probe the jet on high resolution scales an appropriate diagnostic line, or better still a combination of lines, is required. Forbidden transitions, including those in the NIR, e.g., the [FeII] 1.64 µm line, are unsuitable as their critical electron density is exceeded in an outflow when exploring regions very close to the source. In the past permitted HI lines (e.g. $H\alpha$, $Pa\beta$, and $Br\gamma$) were regarded solely as a measure of magnetospheric accretion (e.g., Calvet et al., 2004; Natta et al., 2004, and see Section 5). Nevertheless, we now know that such lines contain a mixture of accretion and outflow signatures (e.g., Lima et al., 2010), i.e. they are composite as demonstrated many years ago using spectro-astronomy (Whelan et al., 2004). Thus permitted lines can be used to produce rudimentary “images” of the launch zone with an interferometer and may even “image” the accretion zone if the interferometer has sufficiently high spatial resolution.

In this way the new generation of spectro-interferometric instruments, such as GRAVITY and CHARA (Eisenhauer et al., 2011; Labdon et al., 2019), given their high spatial (≈ a few milliarcseconds) and moderate spectral resolution ($R = \lambda/\delta\lambda$ up to a few thousand), as well as their higher sensitivity, are beginning to allow us to measure

the extent, position and velocity of the line emitting regions close to the source. In turn this permits us to start to discriminate between the various theories for their origin as the models (stellar, X- and disk-wind) predict the outflow emission to have different initial spatial widths and kinematic properties.

For example, if disk wind theory is correct, then one expects to see line emission, e.g. due to $Br\gamma$, from an extended region (up to a few au) straddling the source in the disk plane, whereas if the X-wind model is closer to reality this emission will be much more confined (several tens of stellar radii). In both case very compact emission is also expected to coincide with the star's magnetospheric accretion zone. Moreover, for lines that are spectrally resolved, it should be possible to determine spatial offsets from the source as a function of velocity on sub-milliarcsecond scales. Potentially this could reveal the changing kinematics of the gas as it is launched from the star or disk. For example disk wind models predict a gradual change in velocity (both in magnitude and direction) as the wind leaves the disk surface at quasi-Keplerian speeds and then accelerates towards the Alfvén radius (Ferreira, 2009).

As mentioned above, optical/NIR interferometry is most easily attempted when the target is bright and so it is no surprise that Herbig Ae/Be stars where the first YSOs to be studied in this way (e.g., Weigelt et al., 2011). Interestingly extended, i.e. well beyond the expected co-rotation radius, $Br\gamma$ emission was found suggesting at least some of the flux comes from a disk wind (Garcia Lopez et al., 2016) and similar results were found in the optical regime (Perraut et al., 2016). While much more challenging in comparison to the bright Herbig Ae/Be stars, studies with interferometers have been made of the much fainter T Tauri stars. Early attempts, using a single baseline, revealed only very compact line emission within the dust sublimation zone (Eisner et al., 2014), and likely showed only the accretion component. More recent studies however show that extended disk wind structures are present along with a magnetospheric accretion component (Gravity Collaboration et al., 2017) and that the disk wind can be dusty (Labdon et al., 2019). The wind component however is only seen at the earliest phases and once the T Tauri star becomes more evolved the accretion component is only observed (Bouvier et al., 2020).

## 9. Future observational prospects

Over the next decade there are a number of new instruments coming online that will no doubt make an impact on our understanding of the launching and propagation of outflows from young stars. It is not possible in this review to cover all of them in any depth. Instead we point out a few that we think might be important.

We begin by mentioning the High Resolution Mid-InfrarEd Spectrometer (HIRMES) which has been developed by NASA Goddard for the Stratospheric Observatory for Infrared Astronomy (SOFIA). HIRMES is optimized to detect neutral atomic oxygen, water, as well as normal and deuterated hydrogen between 28 and 112 µm (Richards et al., 2018). A primarily scientific goal is to investigate the molecular and atomic gas content of protoplanetary disks as well as their ice features, e.g. at 43, 47 and 62 µm. HIRMES, however, also offers the possibility of probing atomic jets launched from embedded systems using the [OI] 63 µm line previously observed, for example, using Herschel/PACS (see 3). Of particular interest is the fact that HIRMES has a medium ($R \approx 12,000$) spectral resolution mode making it ideal to study the kinematics of atomic jets from embedded sources in contrast to Herschel/PACS. It is hoped to commission the instrument in early 2021.

At the time of writing, the launch of the James Webb Space Telescope (JWST), scheduled for October 2021, is eagerly awaited. A number of guaranteed time programs, and no doubt some Cycle 1 time, will be devoted to observing outflows from young stars as well as probing their circumstellar disks. Both JWST's high spatial resolution near-infrared camera and spectrograph, NIRCAM and NIRSPEC respectively, will be invaluable tools for such studies. Of particular
importance will be the work undertaken using the Mid-Infrared In-
strument (MIRI). MIRI has an integral field unit (IFU) spectrograph
covering the range from approximately 5 µm to 28 µm with a spectral
resolution (λ/Δλ) of roughly 2,500. The combination of this and its high
spatial resolution in the mid-infrared, should make it a very valuable
instrument to probe not only conditions in YSO disks (within a few au
from the central object) but also of outflows from the most embedded
sources (e.g. jets from Class 0/I YSOs). Of note is the intention to map
most of the blue lobe of the HH 211 outflow (see 3) with the MIRI
IFU, the source itself, and some of its red lobe. This will require an
extensive mosaic and almost 20 hours of JWST time. Other embedded
outflow sources will also be mapped but with much smaller mosaics
and in addition a large amount of guaranteed time will be devoted to
understanding disk evolution using MIRI.

In order to advance models for wind/jet launching from disks, and even their parent stars, we need a handle on the magnetic field direction and strength in the launch zone. In the case of the least embedded sources, and normally this means weak-line T Tauri stars, we do have an idea of what their surface magnetic field values are like from optical spectro-polarimetry of both absorption and occasionally emission features in their spectrum (Hussain and Alecian, 2014). Here average surface fields in the several kG range are often seen and, somewhat surprisingly, the overall field structure is not necessarily bipolar as sometimes more complex geometries are seen. Unfortunately many of the sources studied to date are relatively evolved and lack an outflow with the exception of AA Tau (Donati et al., 2010; Cox et al., 2013) which possesses a faint jet discovered with HST. This situation however is about to change. A new spectro-polarimeter however has become available on the Canada France Hawaii Telescope (CFHT): known as SPIRou (Spectro-Polarimeter Infra-ROuge) it will utilize the near-infrared part of the spectrum to investigate not only more embedded, and hence less evolved, sources but also the inner zones of their associated disks. Since spectro-polarimetry exploits the Zeeman effect to measure magnetic fields, and the line separation grows with λ² all else being equal, it is easier to measure magnetic field value in the NIR in comparison with the optical. Certainly a number of young stars with outflows are being targeted as part of the SPIRou Legacy Survey (Donati et al., 2018).

In the past few years it is becoming increasingly obvious that radio
observations can give us valuable information on jets from young stars
and, in particular, their magnetic field geometry and even strength in
those cases where the emission is non-thermal. As stated earlier on, in
part this is because of the increased sensitivity of radio interferometers.
Of course the next major advance at these wavelengths will come when
the Square Kilometer Array (SKA) and next generation VLA (ngVLA)
enter service. With these new facilities, we will be able to explore the
jet launching zone with much better spatial resolution than before and
this will be important from the point of view of distinguishing between
various models for the collimation of jets within tens of au from the
source.

The higher sensitivity of the SKA and ngVLA may also allow us to hopefully employ a new diagnostic technique for imaging and kinematic purposes, i.e. using radio recombination lines (RRLs) instead of the broadband continuum that we now observe. RRL allow us not only to peer into the heart of very embedded, and hence otherwise inaccessible YSOs, but also to understand their dynamics (see, for example, Zhang et al., 2017). At present only a few observations of RRL jet emission from massive young stars have been reported and these are almost certainly contaminated with radio emission from an underlying HII region (Jiménez-Serra et al., 2013; Galván-Madrid et al., 2018). The future however is bright!

Funding

TPR would like to acknowledge funding from the European Re-
search Council under Advanced Grant No. 743029, Ejection, Accretion
Structures in YSOs (EASY).

Acknowledgments

We thank the Editors for their invitation to write this review, and
in particular Rob Fender for his patience and ongoing support of this
work. We also wish to acknowledge useful comments from Christian
Fendt, on an early version of this manuscript, and from an anonymous
referee.

References

Agra-Amboage, V., Dougados, C., Cabrit, S., Reunanen, J., 2011. Sub-arcsecond [Fe II] spectro-imaging of the DG Tauri jet. Periodic bubbles and a dusty disk wind? Astron. Astrophys. 532, A59.

Ainsworth, R.E., Scaife, A.M.M., Ray, T.P., Taylor, A.M., Green, D.A., Buckle, J.V., 2014. Tentative evidence for relativistic electrons generated by the jet of the young sun-like star DG Tau. Astrophys. J. 792 (1), L18.

Alcalá, J.M., Manara, C.F., Natta, A., Frasca, A., Testi, L., Nisini, B., Stelzer, B., Williams, J.P., Antoniucci, S., Biazzo, K., Covino, E., Esposito, M., Getman, F., Rigliaco, E., 2017. X-shooter spectroscopy of young stellar objects in Lupus. Accretion properties of class II and transitional objects. Astron. Astrophys. 600, A20.

Anderson, J.M., Li, Z.-Y., Krasnopolsky, R., Blandford, R.D., 2003. Locating the launching region of T Tauri winds: The case of DG tauri. Astrophys. J. 590 (2), L107-L110.

Anderson, J.M., Li, Z.-Y., Krasnopolsky, R., Blandford, R.D., 2005. The structure of magnetocentrifugal winds. I. Steady mass loading. Astrophys. J. 630, 945-957.

Andre, P., Ward-Thompson, D., Barsony, M., 1993. Submillimeter continuum observations of rho Ophiuchi A: The candidate protostar VLA 1623 and prestellar clumps. Astrophys. J. 406, 122.

Anglada, G., Rodríguez, L.F., Carrasco-González, C., 2018. Radio jets from young stellar objects. Astron. Astrophys. Rev. 26, 3.

Antoniucci, S., García López, R., Nisini, B., Caratti o Garatti, A., Giannini, T., Lorenzetti, D., 2014. POISSON project. III. Investigating the evolution of the mass accretion rate. Astron. Astrophys. 572, A62.

Bacciotti, F., Eislöffel, J., 1999. Ionization and density along the beams of Herbig-Haro jets. Astron. Astrophys. 342, 717-735.

Bacciotti, F., Ray, T.P., Mundt, R., Eislöffel, J., Solf, J., 2002. Hubble space telescope/STIS spectroscopy of the optical outflow from DG Tauri: Indications for rotation in the Initial Jet Channel. Astrophys. J. 576 (1), 222-231.

Bai, X.-N., 2011. Magnetorotational-instability-driven accretion in protoplanetary disks. Astrophys. J. 739, 50.

Bai, X.-N., 2015. Hall effect controlled gas dynamics in protoplanetary disks. II. Full 3D simulations toward the outer disk. Astrophys. J. 798 (2), 84.

Bai, X.-N., Stone, J.M., 2013a. Local study of accretion disks with a strong vertical magnetic field: Magnetorotational instability and disk outflow. Astrophys. J. 767, 30.

Bai, X.-N., Stone, J.M., 2013b. Wind-driven accretion in protoplanetary disks. I. Suppression of the magnetorotational instability and launching of the magnetocentrifugal wind. Astrophys. J. 769 (1), 76.

Bai, X.-N., Ye, J., Goodman, J., Yuan, F., 2016. Magneto-thermal disk winds from protoplanetary disks. Astrophys. J. 818, 152.

Balbus, S.A., 2003. Enhanced angular momentum transport in accretion disks. Annu. Rev. Astron. Astrophys. 41, 555-597.

Balbus, S.A., Hawley, J.F., 1991. A powerful local shear instability in weakly magnetized disks. I - Linear analysis. II - Nonlinear evolution. Astrophys. J. 376, 214-233.

Bally, J., Walawender, J., Reipurth, B., 2012. Deep imaging surveys of star-forming clouds. V. New Herbig-Haro shocks and giant outflows in Taurus. Astron. J. 144 (5), 143.

Berger, J.P., Malbet, F., Baron, F., Chiavassa, A., Duvert, G., Elitzur, M., Freytag, B., Gueth, F., Hönig, S., Hron, J., Jang-Condell, H., Le Bouquin, J.B., Monin, J.L., Monnier, J.D., Perrin, G., Plez, B., Ratzka, T., Renard, S., Stefl, S., Thiébaut, E., Tristram, K.R.W., Verhoelst, T., Wolf, S., Young, J., 2012. Imaging the heart of astrophysical objects with optical long-baseline interferometry. Astron. Astrophys. Rev. 20, 53.

Bertout, C., 1989. T Tauri stars: wild as dust. Annu. Rev. Astron. Astrophys. 27, 351-395.

Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., Bouvier, J., 2008. Accretion funnels onto weakly magnetized young stars. Astron. Astrophys. 478, 155-162.

Béthune, W., Lesur, G., Ferreira, J., 2016. Self-organisation in protoplanetary discs. Global, non-stratified Hall-MHD simulations. Astron. Astrophys. 589, A87.

Béthune, W., Lesur, G., Ferreira, J., 2017. Global simulations of protoplanetary disks with net magnetic flux. I. Non-ideal MHD case. Astron. Astrophys. 600, A75.

Bjerkeli, P., van der Wiel, M.H.D., Harsono, D., Ramsey, J.P., Jørgensen, J.K., 2016. Resolved images of a protostellar outflow driven by an extended disk wind. Nature 540, 406-409.

Blandford, R.D., Payne, D.G., 1982. Hydromagnetic flows from accretion discs and the production of radio jets. Mon. Not. R. Astron. Soc. 199, 883-903.
Blandford, R.D., Rees, M.J., 1974. A "twin-exhaust" model for double radio sources. Mon. Not. R. Astron. Soc. 169, 395-415.

Bonito, R., Orlando, S., Miceli, M., Peres, G., Micela, G., Favata, F., 2011. X-ray emission from protostellar jet HH 154: The first evidence of a diamond shock? Astrophys. J. 737, 54.

Bosch-Ramon, V., Romero, G.E., Araudo, A.T., Paredes, J.M., 2010. Massive protostars as gamma-ray sources. Astron. Astrophys. 511, A8.

Bouvier, J., Bertout, C., Benz, W., Mayor, M., 1986. Rotation in T Tauri stars. I - Observations and immediate analysis. Astron. Astrophys. 165, 110-119.

Bouvier, J., Perraut, K., Le Bouquin, J.B., Duvert, G., Dougados, C., Brandner, W., Benisty, M., Berger, J.P., Alécian, E., 2020. Probing the magnetospheric accretion region of the young pre-transitional disk system DoAr 44 using VLTI/GRAVITY. Astron. Astrophys. 636, A108.

Buehrke, T., Mundt, R., Ray, T.P., 1988. A detailed study of HH 34 and its associated jet. Astron. Astrophys. 200, 99-119.

Cabrit, S., 2007. Jets from young stars: The need for MHD collimation and acceleration processes. In: Ferreira, J., Dougados, C., Whelan, E. (Eds.), Jets from Young Stars: Models and Constraints. In: Lecture Notes in Physics, vol. 723, Berlin Springer Verlag, p. 21.

Cai, M.J., Shang, H., Lin, H.-H., Shu, F.H., 2008. X-Winds in action. Astrophys. J. 672, 489-503.

Calvet, N., Gullbring, E., 1998. The structure and emission of the accretion shock in T Tauri stars. Astrophys. J. 509 (2), 802-818.

Calvet, N., Muzerolle, J., Briceno, C., Hernández, J., Hartmann, L., Saucedo, J.L., Gordon, K.D., 2004. The mass accretion rates of intermediate-mass T Tauri stars. Astron. J. 128 (3), 1294-1318.

Canto, J., 1980. A stellar wind model for Herbig-Haro objects. Astron. Astrophys. 86 (3), 327-338.

Caratti o Garatti, A., Tambovtseva, L.V., Garcia Lopez, R., Kraus, S., Schertl, D., Grinin, V.P., Weigelt, G., Hofmann, K.H., Massi, F., Lagarde, S., Vannier, M., Malbet, F., 2015. AMBER/VLTI high spectral resolution observations of the Bry emitting region in HD 98922. A compact disc wind launched from the inner disc region. Astron. Astrophys. 582, A44.

Carrasco-González, C., Rodríguez, L.F., Anglada, G., Martí, J., Torrelles, J.M., Osorio, M., 2010. A magnetized jet from a massive protostar. Science 311, 1209.

Casse, F., Ferreira, J., 2000. Magnetized accretion-ejection structures. V. Effects of entropy generation inside the disc. Astron. Astrophys. 361, 1178-1190.

Casse, F., Keppens, R., 2002. Magnetized accretion-ejection structures: 2.5-dimensional magnetohydrodynamic simulations of continuous ideal jet launching from resistive accretion disks. Astrophys. J. 581, 988-1001.

Casse, F., Keppens, R., 2004. Radiatively inefficient magnetohydrodynamic accretion-ejection structures. Astrophys. J. 601, 90-103.

Cerqueira, A.H., Velázquez, P.F., Raga, A.C., Vasconcelos, M.J., de Colle, F., 2006. Radial velocity asymmetries from jets with variable velocity profiles. In: Herrera-Velázquez, J.J.E. (Ed.), Plasma and Fusion Science: 16th IAEA Technical Meeting on Research using Small Fusion Devices. In: American Institute of Physics Conference Series, vol. 875, pp. 285-288.

Chakravorty, S., Petrucci, P.-O., Ferreira, J., Henri, G., Belmont, R., Clavel, M., Corbel, S., Rodriguez, J., Coriat, M., Drappeau, S., Malzac, J., 2016. Absorption lines from magnetically driven winds in X-ray binaries. Astron. Astrophys. 589, A119.

Chan, K.L., Henriksen, R.N., 1980. On the supersonic dynamics of magnetized jets of thermal gas in radio galaxies. Astrophys. J. 241, 534-551.

Cleeves, L.I., Adams, F.C., Bergin, E.A., 2013. Exclusion of cosmic rays in protoplanetary disks: Stellar and magnetic effects. Astrophys. J. 772 (1), 5.

Codella, C., Cabrit, S., Gueth, F., Cesaroni, R., Bacciotti, F., Lefloch, B., McCaughrean, M.J., 2007. A highly-collimated sio jet in the HH212 protostellar outflow. Astron. Astrophys. 462 (3), L53-L56.

Coffey, D., Bacciotti, F., Ray, T.P., Eisloffel, J., Woitas, J., 2007. Further indications of jet rotation in new ultraviolet and optical hubble space telescope STIS spectra. Astrophys. J. 663 (1), 350-364.

Combet, C., Ferreira, J., 2008. The radial structure of protostellar accretion disks: influence of jets. Astron. Astrophys. 479, 481-491.

Coughlan, C.P., Ainsworth, R.E., Eisloffel, J., Hoeft, M., Drabent, A., Scaife, A.M.M., Ray, T.P., Bell, M.E., Broderick, J.W., Corbel, S., Grießmeier, J.-M., van der Horst, A.J., van Leeuwen, J., Markoff, S., Pietka, M., Stewart, A.J., Wijers, R.A.M.J., Zarka, P., 2017. A LOFAR detection of the low-mass young star T Tau at 149 MHz. Astrophys. J. 834 (2), 206.

Cox, A.W., Grady, C.A., Hammel, H.B., Hornbeck, J., Russell, R.W., Sitko, M.L., Woodgate, B.E., 2013. Imaging the disk and jet of the classical T Tauri star AA Tau. Astrophys. J. 762, 40.

De Colle, F., Cerqueira, A.H., Riera, A., 2016. Transverse velocity shifts in protostellar jets: Rotation or velocity asymmetries? Astrophys. J. 832 (2), 152.

DeCampli, W.M., 1981. T Tauri winds. Astrophys. J. 244, 124-146.

Dionatos, O., Ray, T., Güdel, M., 2018. Herschel spectral-line mapping of the HH211 protostellar system. Astron. Astrophys. 616, A84.

Donati, J.-F., Kouach, D., Lacombe, M., Baratchart, S., Doyon, R., Delfosse, X., Artigau, É., Moutou, C., Hébrard, G., Bouchy, F., Bouvier, J., Alencar, S., Saddlemyer, L., Parès, L., Rabou, P., Micheau, Y., Dolon, F., Barrick, G., Hernández, O., Wang, S.Y., Reshetov, V., Striebig, N., Challita, Z., Carmona, A., Tibault, S., Martioli, E., Figueira, P., Boisse, I., Pepe, F., 2018. Handbook of Exoplanets. Springer International Publishing, p. 107, SPIRou: A NIR Spectropolarimeter/High-Precision Velocimeter for the CFHT.

Donati, J.-F., Skelly, M.B., Bouvier, J., Gregory, S.G., Grankin, K.N., Jardine, M.M., Hussain, G.A.J., Ménard, F., Dougados, C., Unruh, Y., Mohanty, S., Aurière, M., Morin, J., Farès, R., MAPP Collaboration, 2010. Magnetospheric accretion and spin-down of the prototypical classical T Tauri star AA Tau. Mon. Not. R. Astron. Soc. 409, 1347-1361.

Dougados, C., Cabrit, S., Lavalley, C., Ménard, F., 2000. T Tauri stars microjets resolved by adaptive optics. Astron. Astrophys. 357, L61-L64.

Edwards, S., 1997. Magnetospherically mediated accretion in classical T Tauri stars. Paradigm for low mass stars undergoing disk accretion? In: Reipurth, B., Bertout, C. (Eds.), Herbig-Haro Flows and the Birth of Stars. In: IAU Symposium, vol. 182, pp. 433-442.

Edwards, S., Fischer, W., Hillenbrand, L., Kwan, J., 2006. Probing T Tauri accretion and outflow with 1 micron spectroscopy. Astrophys. J. 646, 319-341.

Edwards, S., Fischer, W., Kwan, J., Hillenbrand, L., Dupree, A.K., 2003. He I λ10830 as a probe of winds in accreting young stars. Astrophys. J. 599, L41-L44.

Eisenhauer, F., Perrin, G., Brandner, W., Straubmeier, C., Perraut, K., Amorim, A., Schöller, M., Gillessen, S., Kervella, P., Benisty, M., Araujo-Hauck, C., Jocou, L., Lima, J., Jakob, G., Haug, M., Clénet, Y., Henning, T., Eckart, A., Berger, J.P., Garcia, P., Abuter, R., Kellner, S., Paumard, T., Hippler, S., Fischer, S., Moulin, T., Villate, J., Avila, G., Gräter, A., Lacour, S., Huber, A., Wiest, M., Nolot, A., Carvas, P., Dorn, R., Pfuhl, O., Gendron, E., Kendrew, S., Yazici, S., Anton, S., Jung, Y., Thiel, M., Choquet, É., Klein, R., Teixeira, P., Gitton, P., Moch, D., Vincent, F., Kudryavtseva, N., Ströbele, S., Sturm, S., Fédou, P., Lenzen, R., Jolley, P., Kister, C., Lapeyrère, V., Naranjo, V., Lucuix, C., Hofmann, R., Chapron, F., Neumann, U., Mehrgan, L., Hans, O., Rousset, G., Ramos, J., Suarez, M., Lederer, R., Reess, J.M., Rohloff, R.R., Haguenauer, P., Bartko, H., Sevin, A., Wagner, K., Lizon, J.L., Rabien, S., Collin, C., Finger, G., Davies, R., Rouan, D., Wittkowski, M., Dodds-Eden, K., Ziegler, D., Cassaing, F., Bonnet, H., Casali, M., Genzel, R., Lena, P., 2011. GRAVITY: Observing the universe in motion. The Messenger 143, 16-24.

Eisloffel, J., Mundt, R., 1992. Proper motion measurements in the HH 34 jet and its associated bow shocks. Astron. Astrophys. 263, 292-300.

Eisner, J.A., Hillenbrand, L.A., Stone, J.M., 2014. Constraining the sub-au-scale distribution of hydrogen and carbon monoxide gas around young stars with the Keck Interferometer. Mon. Not. R. Astron. Soc. 443 (3), 1916-1945.

Feeney-Johansson, A., Purser, S.J.D., Ray, T.P., Eisloffel, J., Hoeft, M., Drabent, A., Ainsworth, R.E., 2019. The first detection of a low-frequency turnover in nonthermal emission from the jet of a young star. Astrophys. J. 885 (1), L7.

Fendt, C., 2006. Collimation of astrophysical jets: The role of the accretion disk magnetic field distribution. Astrophys. J. 651, 272-287.

Ferreira, J., 1997. Magnetically-driven jets from Keplerian accretion discs. Astron. Astrophys. 319, 340-359.

Ferreira, J., 2009. Self-collimated jets from accretion discs and star-disc interaction zones. Astrophys. Space Sci. Proc. 13, 99-110.

Ferreira, J., Casse, F., 2013. On fan-shaped cold MHD winds from Keplerian accretion discs. Mon. Not. R. Astron. Soc. 428, 307-320.

Ferreira, J., Dougados, C., Cabrit, S., 2006. Which jet launching mechanism(s) for T Tauri stars? Astron. Astrophys. 453, 785-796.

Ferreira, J., Pelletier, G., 1993. Magnetized accretion-ejection structures. 1. General statements. Astron. Astrophys. 276, 625.

Ferreira, J., Pelletier, G., 1995. Magnetized accretion-ejection structures. III. Stellar and extragalactic jets as weakly dissipative disk outflows. Astron. Astrophys. 295, 807.

Flaherty, K.M., Hughes, A.M., Teague, R., Simon, J.B., Andrews, S.M., Wilner, D.J., 2018. Turbulence in the TW Hya Disk. Astrophys. J. 856 (2), 117.

Fromang, S., Latter, H., Lesur, G., Ogilvie, G.I., 2013. Local outflows from turbulent accretion disks. Astron. Astrophys. 552, A71.

Fukumura, K., Kazanas, D., Contopoulos, I., Behar, E., 2010. Magnetohydrodynamic accretion disk winds as X-ray absorbers in active galactic nuclei. Astrophys. J. 715, 636-650.

Gallet, F., Bouvier, J., 2013. Improved angular momentum evolution model for solar-like stars. Astron. Astrophys. 556, A36.

Galván-Madrid, R., Beltrán, M., Ginsburg, A., Carrasco-González, C., Liu, H.B., Rodríguez, L.F., Kurtz, S., 2018. Resolving the structure and kinematics of the youngest HII regions and radio jets from young stellar objects. In: Murphy, E. (Ed.), Science with a Next Generation Very Large Array. In: Astronomical Society of the Pacific Conference Series, vol. 517, p. 309.

Gammie, C.F., 1996. Layered accretion in T Tauri disks. Astrophys. J. 457, 355.

García Lopez, R., Kurosawa, R., Caratti o Garatti, A., Kreplin, A., Weigelt, G., Tambovtseva, L.V., Grinin, V.P., Ray, T.P., 2016. Investigating the origin and spectroscopic variability of the near-infrared H I lines in the Herbig star VV Ser. Mon. Not. R. Astron. Soc. 456 (1), 156-170.

Goodson, A.P., Winglee, R.M., Boöhm, K.-H., 1997. Time-dependent accretion by magnetic young stellar objects as a launching mechanism for stellar jets. Astrophys. J. 489, 199.
Gravity Collaboration, Garcia Lopez, R., Perraut, K., Caratti O, Garatti, A., Lazareff, B., Sanchez-Bermudez, J., Benisty, M., Dougados, C., Labadie, L., Brandner, W., Garcia, P.J.V., Henning, T., Ray, T.P., Abuter, R., Amorim, A., Anugu, N., Berger, J.P., Bonnet, H., Buron, A., Caselli, P., Clénet, Y., Coudé Du Foresto, V., de Wit, W., Deen, C., Delplancke-Ströbele, F., Dexter, J., Eckart, A., Eisenhauer, F., Garcia Dabo, C.E., Gendron, E., Genzel, R., Gillessen, S., Haubois, X., Haug, M., Haussmann, F., Hippler, S., Hubert, Z., Hummel, C.A., Horrobin, M., Jocou, L., Kellner, S., Kervella, P., Kulas, M., Kolb, J., Lacour, S., Le Bouquin, J.B., Léna, P., Lippa, M., Mérand, A., Müller, E., Ott, T., Panduro, J., Paumard, T., Perrin, G., Pufhl, O., Ramirez, A., Rau, C., Rohloff, R.R., Rousset, G., Scheithauer, S., Schöller, M., Straubmeier, C., Sturm, E., Thi, W.F., van Dishoeck, E., Vincent, F., Waisberg, I., Wank, I., Wieprecht, E., Wiest, M., Wiezorrek, E., Woillez, J., Yazici, S., Zins, G., 2017. The wind and the magnetospheric accretion onto the T Tauri star S Coronae Australis at sub-au resolution. Astron. Astrophys. 608, A78.

Gressel, O., Turner, N.J., Nelson, R.P., McNally, C.P., 2015. Global simulations of protoplanetary disks with Ohmic resistivity and ambipolar diffusion. Astrophys. J. 801, 84.

Güdel, M., Skinner, S.L., Briggs, K.R., Audard, M., Arzner, K., Telleschi, A., 2005. Evidence for an X-ray jet in DG Tauri A? Astrophys. J. 626, L53-L56.

Gueth, F., Guilloteau, S., 1999. The jet-driven molecular outflow of HH 211. Astron. Astrophys. 343, 571-584.

Haro, G., 1952. Herbig's nebulous objects near NGC 1999. Astrophys. J. 115, 572.

Hartigan, P., Edwards, S., Ghandour, L., 1995. Disk accretion and mass loss from young stars. Astrophys. J. 452, 736.

Hartigan, P., Frank, A., Varniére, P., Blackman, E.G., 2007. Magnetic fields in stellar jets. Astrophys. J. 661, 910-918.

Hartigan, P., Morse, J., 2007. Collimation, proper motions, and physical conditions in the HH 30 jet from Hubble Space Telescope Slitless Spectroscopy. Astrophys. J. 660, 426-440.

Hartigan, P., Morse, J.A., Raymond, J., 1994. Mass-loss rates, ionization fractions, shock velocities, and magnetic fields of stellar jets. Astrophys. J. 436, 125.

Hartigan, P., Wright, A., 2015. A new diagnostic of magnetic field strengths in radiatively-cooled shocks. Astrophys. J. 811, 12.

Hartmann, L., Herczeg, G., Calvet, N., 2016. Accretion onto pre-main-sequence stars. Annu. Rev. Astron. Astrophys. 54, 135-180.

Hartmann, L., MacGregor, K.B., 1982. Wave-driven winds from cool stars. I - Some effects of magnetic field geometry. Astrophys. J. 257, 264-268.

Hawley, J.F., Gammie, C.F., Balbus, S.A., 1995. Local three-dimensional magnetohydrodynamical simulations of accretion disks. Astrophys. J. 440, 742.

Herbig, G.H., 1951. The spectra of two nebulous objects near NGC 1999. Astrophys. J. 113, 697-699.

Hirota, T., Machida, M.N., Matsushita, Y., Motogi, K., Matsumoto, N., Kim, M.K., Burns, R.A., Honma, M., 2017. Disk-driven rotating bipolar outflow in Orion Source I. Nat. Astron. 1, 0146.

Hussain, G.A.J., Alecian, E., 2014. The role of magnetic fields in pre-main sequence stars. In: Petit, P., Jardine, M., Spruit, H.C. (Eds.), Magnetic Fields Throughout Stellar Evolution. In: IAU Symposium, vol. 302, pp. 25-37.

Jacquemin-Ide, J., Ferreira, J., Lesur, G., 2019. Magnetically driven jets and winds from weakly magnetized accretion discs. Mon. Not. R. Astron. Soc. 490 (3), 3112-3133.

Jhan, K.-S., Lee, C.-F., 2016. A multi-epoch SMA study of the HH 211 protostellar jet: Jet motion and knot formation. Astrophys. J. 816 (1), 32.

Jiménez-Serra, I., Báez-Rubio, A., Rivilla, V.M., Martín-Pintado, J., Zhang, Q., Dierickx, M., Patel, N., 2013. A new radio recombination line maser object toward the MonR2 H II region. Astrophys. J. 764 (1), L4.

Königl, A., 1989. Self-similar models of magnetized accretion disks. Astrophys. J. 342, 208-223.

Königl, A., 1991. Disk accretion onto magnetic T Tauri stars. Astrophys. J. 370, L39.

Krasnopolsky, R., Li, Z.-Y., Blandford, R., 1999. Magnetocentrifugal launching of jets from accretion disks. I. Cold axisymmetric flows. Astrophys. J. 526, 631-642.

Krasnopolsky, R., Li, Z.-Y., Blandford, R.D., 2003. Magnetocentrifugal launching of jets from accretion disks. II. Inner disk-driven winds. Astrophys. J. 595, 631-642.

Krumholz, M.R., 2015. Notes on star formation. arXiv e-prints, arXiv:1511.03457.

Kwan, J., 1997. Warm disk coronae in classical T Tauri stars. Astrophys. J. 489, 284.

Kwan, J., Tademaru, E., 1995. Disk winds from T Tauri stars. Astrophys. J. 454, 382.

Labdon, A., Kraus, S., Davies, C.L., Kreplin, A., Kluska, J., Harries, T.J., Monnier, J.D., ten Brummelaar, T., Baron, F., Millan-Gabet, R., Kloppenborg, B., Eisner, J., Sturmann, J., Sturmann, L., 2019. Dusty disk winds at the sublimation rim of the highly inclined, low mass young stellar object SU Aurigae. Astron. Astrophys. 627, A36.

Lada, C.J., Wilking, B.A., 1984. The nature of the embedded population in the rho Ophiuchi dark cloud : mid-infrared observations. Astrophys. J. 287, 610-621.

Lee, C.-F., Ho, P.T.P., Li, Z.-Y., Hirano, N., Zhang, Q., Shang, H., 2017. A rotating protostellar jet launched from the innermost disk of HH 212. Nat. Astron. 1, 0152.

Lee, C.-F., Hwang, H.-C., Ching, T.-C., Hirano, N., Lai, S.-P., Rao, R., Ho, P.T.P., 2018. Unveiling a magnetized jet from a low-mass protostar. Nature Commun. 9, 4636.

Lesur, G., Ferreira, J., Ogilvie, G.I., 2013. The magnetorotational instability as a jet launching mechanism. Astron. Astrophys. 550, A61.

Lesur, G., Kunz, M.W., Fromang, S., 2014. Thanatology in protoplanetary discs. The combined influence of Ohmic, Hall, and ambipolar diffusion on dead zones. Astron. Astrophys. 566, A56.

Li, Z.-Y., 1995. Magnetohydrodynamic disk-wind connection: Self-similar solutions. Astrophys. J. 444, 848-860.

Li, Z.-Y., 1996. Magnetohydrodynamic disk-wind connection: Magnetocentrifugal winds from ambipolar diffusion-dominated accretion disks. Astrophys. J. 465, 855.

Lima, G.H.R.A., Alencar, S.H.P., Calvet, N., Hartmann, L., Muzerolle, J., 2010. Modeling the Hα line emission around classical T Tauri stars using magnetospheric accretion and disk wind models. Astron. Astrophys. 522, A104.

Louvet, F., Dougados, C., Cabrit, S., Hales, A., Pinte, C., Ménard, F., Bacciotti, F., Coffey, D., Mardones, D., Bronfman, L., Gueth, F., 2016. ALMA observations of the <ASTROBJ>Th 28<>/ASTROBJ> protostellar disk. A new example of counter-rotation between disk and optical jet. Astron. Astrophys. 596, A88.

Louvet, F., Dougados, C., Cabrit, S., Mardones, D., Ménard, F., Tabone, B., Pinte, C., Dent, W.R.F., 2018. The HH30 edge-on T Tauri star. A rotating and precessing monopolar outflow scrutinized by ALMA. Astron. Astrophys. 618, A120.

Matt, S.P., Pinzón, G., Greene, T.P., Pudritz, R.E., 2012. Spin evolution of accreting young stars. II. Effect of accretion-powered stellar winds. Astrophys. J. 745, 101.

Matt, S., Pudritz, R.E., 2005. Accretion-powered stellar winds as a solution to the stellar angular momentum problem. Astrophys. J. 632, L135-L138.

McCaughrean, M.J., Rayner, J.T., Zinnecker, H., 1994. Discovery of a molecular hydrogen jet near IC 348. Astrophys. J. 436, L189.

Meliani, Z., Casse, F., Sauty, C., 2006. Two-component magnetohydrodynamical outflows around young stellar objects. Interplay between stellar magnetospheric winds and disc-driven jets. Astron. Astrophys. 460, 1-14.

Melnikov, S.Y., Eislöffel, J., Bacciotti, F., Woitas, J., Ray, T.P., 2009. HST/STIS observations of the RW Aurigae bipolar jet: mapping the physical parameters close to the source. Astron. Astrophys. 506, 763-777.

Mestel, L., Spruit, H.C., 1987. On magnetic braking of late-type stars. Mon. Not. R. Astron. Soc. 226, 57-66.

Murphy, G.C., Ferreira, J., Zanni, C., 2010. Large scale magnetic fields in viscous resistive accretion disks. I. Ejection from weakly magnetized disks. Astron. Astrophys. 512, 82.

Najita, J.R., Bergin, E.A., 2018. Protoplanetary disk sizes and angular momentum transport. Astrophys. J. 864 (2), 168.

Natta, A., Testi, L., Muzerolle, J., Randich, S., Comerón, F., Persi, P., 2004. Accretion in brown dwarfs: An infrared view. Astron. Astrophys. 424, 603-612.

Nisini, B., Codella, C., Giannini, T., Santiago Garcia, J., Richer, J.S., Bachiller, R., Tafalla, M., 2007. Warm SiO gas in molecular bullets associated with protostellar outflows. Astron. Astrophys. 462 (1), 163-172.

Ogilvie, G.I., Livio, M., 2001. Launching of jets and the vertical structure of accretion disks. Astrophys. J. 553, 158-173.

O'Keeffe, W., Downes, T.P., 2014. Multifield simulations of the magnetorotational instability in protostellar discs. Mon. Not. R. Astron. Soc. 441, 571-581.

Ouyed, R., Clarke, D.A., Pudritz, R.E., 2003. Three-dimensional simulations of jets from Keplerian disks: Self-regulatory stability. Astrophys. J. 582, 292-319.

Ouyed, R., Pudritz, R.E., 1997a. Numerical simulations of astrophysical jets from Keplerian disks. I. Stationary models. Astrophys. J. 482, 712.

Ouyed, R., Pudritz, R.E., 1997b. Numerical simulations of astrophysical jets from Keplerian disks. II. Episodic outflows. Astrophys. J. 484, 794.

Perraut, K., Dougados, C., Lima, G.H.R.A., Benisty, M., Mourard, D., Ligi, R., Nardetto, N., Tallon-Bosc, I., ten Brummelaar, T., Farrington, C., 2016. A disk wind in AB Aurigae traced with Hα interferometry. Astron. Astrophys. 596, A17.

Podio, L., Eislöffel, J., Melnikov, S., Hodapp, K.W., Bacciotti, F., 2011. Tracing kinematical and physical asymmetries in the jet from DG Tauri B. Astron. Astrophys. 527, A13.

Podio, L., Medves, S., Bacciotti, F., Eislöffel, J., Ray, T., 2009. Physical structure and dust reprocessing in a sample of HH jets. Astron. Astrophys. 506 (2), 779-788.

Pudritz, R.E., Norman, C.A., 1986. Bipolar hydromagnetic winds from disks around protostellar objects. Astrophys. J. 301, 571-586.

Pudritz, R.E., Ray, T.P., 2019. The role of magnetic fields in protostellar outflows and star formation. Front. Astron. Space Sci. 6, 54.

Pudritz, R.E., Rogers, C.S., Ouyed, R., 2006. Controlling the collimation and rotation of hydromagnetic disc winds. Mon. Not. R. Astron. Soc. 365, 1131-1148.

Purser, S.J.D., Lumsden, S.L., Hoare, M.G., Urquhart, J.S., Cunningham, N., Purcell, C.R., Brooks, K.J., Garay, G., Gúzman, A.E., Voronkov, M.A., 2016. A search for ionized jets towards massive young stellar objects. Mon. Not. R. Astron. Soc. 460 (1), 1039-1053.

Qiu, K., Wyrowski, F., Menten, K., Zhang, Q., Güsten, R., 2019. CO multi-line observations of HH 80-81: A two-component molecular outflow associated with the largest protostellar jet in our galaxy. Astrophys. J. 871 (2), 141.

Raga, A.C., Kofman, L., 1992. Knots in stellar jets from time-dependent sources. Astrophys. J. 386, 222.

Ramsey, J.P., Clarke, D.A., 2011. Simulating protostellar jets simultaneously at launching and observational scales. Astrophys. J. 728, L11.

Rees, M.J., 1978. The M87 jet: internal shocks in a plasma beam? Mon. Not. R. Astron. Soc. 184, 61P-65P.

Reynolds, S.P., 1986. Continuum spectra of collimated, ionized stellar winds. Astrophys. J. 304, 713.
Richards, S.N., Moseley, S.H., Stacey, G., Greenhouse, M., Kutyrev, A., Arendt, R., Atanasoff, H., Banks, S., Brekosky, R.P., Brown, A.-D., Bulcha, B., Cazeau, T., Choi, M., Colazo, F., Engler, C., Hadjimichael, T., Hays-Wehle, J., Henderson, C., Hsieh, W.-T., Huang, J., Jenstrom, I., Kellogg, J., Kimball, M., Kovács, A., Leiter, S., Maher, S., McMurray, R., Melnick, G.J., Mentzell, E., Mikula, V., Miller, T.M., Nagler, P., Nikola, T., Oxborrow, J., Pontoppidan, K.M., Rangwala, N., Rhodes, A., Roberge, A., Rosner, S., Rostem, K., Rustemeyer, N., Sharp, E., Sparr, L., Stevanovic, D., Taraschi, P., Temi, P., Vacca, W.D., de Lorenzo, J.V.H., Wohler, B., Wollack, E.J., Wilks, S., 2018. SOFIA-HIRMES: Looking forward to the High-Resolution Mid-infrarEd spectrometer. J. Astron. Instrum. 7 (4), 1840015.

Richer, J.S., Shepherd, D.S., Cabrit, S., Bachiller, R., Churchwell, E., 2000. Molecular outflows from young stellar objects. In: Mannings, V., Boss, A.P., Russell, S.S. (Eds.), Protostars and Planets IV. p. 867.

Rodgers-Lee, D., Taylor, A.M., Ray, T.P., Downes, T.P., 2017. The ionizing effect of low-energy cosmic rays from a class II object on its protoplanetary disc. Mon. Not. R. Astron. Soc. 472 (1), 26–38.

Rodríguez, L.F., Curiel, S., Moran, J.M., Mirabel, I.F., Roth, M., Garay, G., 1989. Large proper motions in the remarkable triple radio source in serpens. Astrophys. J. 346, L85.

Rodríguez-Kamenetzky, A., Carrasco-González, C., González-Martín, O., Araudo, A.T., Rodríguez, L.F., Vig, S., Hofner, P., 2019. Particle acceleration in the herbig-haro objects HH 80 and HH 81. Mon. Not. R. Astron. Soc. 482 (4), 4687–4696.

Romanova, M.M., Ustyugova, G.V., Koldoba, A.V., Lovelace, R.V.E., 2002. Magnetohydrodynamic simulations of disk-magnetized star interactions in the quiescent regime: Funnel flows and angular momentum transport. Astrophys. J. 578, 420–438.

Romanova, M.M., Ustyugova, G.V., Koldoba, A.V., Lovelace, R.V.E., 2009. Launching of conical winds and axial jets from the disc-magnetosphere boundary: axisymmetric and 3D simulations. Mon. Not. R. Astron. Soc. 399, 1802–1828.

Romanova, M.M., Ustyugova, G.V., Koldoba, A.V., Lovelace, R.V.E., 2012. MRI-Driven accretion on to magnetized stars: global 3D MHD simulations of magnetospheric and boundary layer regimes. Mon. Not. R. Astron. Soc. 421, 63–77.

Salmeron, R., Königl, A., Wardle, M., 2011. Wind-driving protostellar accretion discs - II. Numerical method and illustrative solutions. Mon. Not. R. Astron. Soc. 412, 1162–1180.

Scepi, N., Lesur, G., Dubus, G., Flock, M., 2018. Turbulent and wind-driven accretion in dwarf novae threaded by a large-scale magnetic field. Astron. Astrophys. 620, A49.

Schatzman, E., 1962. A theory of the role of magnetic activity during star formation. Ann. Astrophys. 25, 18.

Schwartz, R.D., 1977. Evidence of star formation triggered by expansion of the Gum Nebula. Astrophys. J. 212, L25–L26.

Shakura, N.I., Sunyaev, R.A., 1973. Black holes in binary systems. Observational appearance. Astron. Astrophys. 24, 337–355.

Sheikhnezami, S., Fendt, C., 2015. Wobbling and precessing jets from warped disks in binary systems. Astrophys. J. 814 (2), 113.

Sheikhnezami, S., Fendt, C., 2018. Long-term simulation of MHD jet launching in an orbiting star-disk system. Astrophys. J. 861 (1), 11.

Sheikhnezami, S., Fendt, C., Porth, O., Vaidya, B., Ghanbari, J., 2012. Bipolar jets launched from magnetically diffusive accretion disks. I. Ejection efficiency versus field strength and diffusivity. Astrophys. J. 757, 65.

Shu, F.H., Lizano, S., Ruden, S.P., Najita, J., 1988. Mass loss from rapidly rotating magnetic protostars. Astrophys. J. 328, L19–L23.

Shu, F., Najita, J., Ostriker, E., Wilkin, F., Ruden, S., Lizano, S., 1994. Magnetocentrifugally driven flows from young stars and disks. 1: A generalized model. Astrophys. J. 429, 781–796.

Skinner, S.L., Audard, M., Güdel, M., 2011. Chandra evidence for extended X-ray structure in RY Tau. Astrophys. J. 737, 19.

Snell, R.L., Bally, J., Strom, S.E., Strom, K.M., 1985. Radio and optical observations of the jets L 1551 IRS 5. Astrophys. J. 290, 587–595.

Snell, R.L., Loren, R.B., Plambeck, R.L., 1980. Observations of CO in L 1551: evidence for stellar wind driven shocks. Astrophys. J. 239, L17–L22.

Staff, J.E., Koning, N., Ouyed, R., Thompson, A., Pudritz, R.E., 2015. Hubble space telescope scale 3D simulations of MHD disc winds: a rotating two-component jet structure. Mon. Not. R. Astron. Soc. 446 (4), 3975–3991.

Staff, J.E., Niebergal, B.P., Ouyed, R., Pudritz, R.E., Cai, K., 2010. Confronting three-dimensional time-dependent jet simulations with hubble space telescope observations. Astrophys. J. 722, 1325–1332.

Stepanovs, D., Fendt, C., 2016. An extensive numerical survey of the correlation between outflow dynamics and accretion disk magnetization. Astrophys. J. 825, 14.

Stojimirović, I., Narayanan, G., Snell, R.L., Bally, J., 2006. Entrainment mechanisms for outflows in the L1551 star-forming region. Astrophys. J. 649 (1), 280–298.

Stone, J.M., Hawley, J.F., Gammie, C.F., Balbus, S.A., 1996. Three-dimensional magnetohydrodynamical simulations of vertically stratified accretion disks. Astrophys. J. 463, 656.

Suriano, S.S., Li, Z.-Y., Krasnopolsky, R., Shang, H., 2018. The formation of rings and gaps in magnetically coupled disc-wind systems: ambipolar diffusion and reconnection. Mon. Not. R. Astron. Soc. 477, 1239–1257.

Suzuki, T.K., Inutsuka, S., 2009. Disk winds driven by magnetorotational instability and dispersal of protoplanetary disks. Astrophys. J. 691, L49–L54.

Suzuki, T.K., Inutsuka, S., 2014. Magnetohydrodynamic simulations of global accretion disks with vertical magnetic fields. Astrophys. J. 784, 121.

Takeuchi, T., Okuzumi, S., 2014. Radial transport of large-scale magnetic fields in accretion disks. II. Relaxation to steady states. Astrophys. J. 797, 132.

Tappe, A., Forbrich, J., Martín, S., Yuan, Y., Lada, C.J., 2012. The anatomy of the young protostellar outflow HH 211. Astrophys. J. 751 (1), 9.

Tychoniec, Ł., Tobin, J.J., Karska, A., Chandler, C., Dunham, M.M., Harris, R.J., Kratter, K.M., Li, Z.-Y., Looney, L.W., Melis, C., Pérez, L.M., Sadavoy, S.I., Segura-Cox, D., van Dishoeck, E.F., 2018. The VLA nascent disk and multiplicity survey of perseus protostars (VANDAM). IV. Free-free emission from protostars: Links to infrared properties, outflow tracers, and protostellar disk masses. Astrophys. J. Suppl. 238 (2), 19.

Tzeferacos, P., Ferrari, A., Mignone, A., Zanni, C., Bodo, G., Massaglia, S., 2009. On the magnetization of jet-launching discs. Mon. Not. R. Astron. Soc. 400, 820–834.

Tzeferacos, P., Ferrari, A., Mignone, A., Zanni, C., Bodo, G., Massaglia, S., 2013. Effects of entropy generation in jet-launching discs. Mon. Not. R. Astron. Soc. 428, 3151–3163.

Venuti, L., Bouvier, J., Irwin, J., Stauffer, J.R., Hillenbrand, L.A., Rebull, L.M., Cody, A.M., Alencar, S.H.P., Micela, G., Flaccomio, E., Peres, G., 2015. UV variability and accretion dynamics in the young open cluster NGC 2264. Astron. Astrophys. 581, A66.

Wang, L., Bai, X.-N., Goodman, J., 2019. Global simulations of protoplanetary disk outflows with coupled non-ideal magnetohydrodynamics and consistent thermochemistry. Astrophys. J. 874 (1), 90.

Wardle, M., Königl, A., 1993. The structure of protostellar accretion disks and the origin of bipolar flows. Astrophys. J. 410, 218–238.

Watson, D.M., Calvet, N.P., Fischer, W.J., Forrest, W.J., Manoj, P., Megeath, S.T., Melnick, G.J., Najita, J., Neufeld, D.A., Sheehan, P.D., Stutz, A.M., Tobin, J.J., 2016. Evolution of mass outflow in protostars. Astrophys. J. 828, 52.

Weigelt, G., Grinin, V.P., Groh, J.H., Hofmann, K.H., Kraus, S., Miroshnichenko, A.S., Schertl, D., Tambovtseva, L.V., Benisty, M., Driebe, T., Lagarde, S., Malbet, F., Meilland, A., Petrov, R., Tatulli, E., 2011. VLTI/AMBER spectro-interferometry of the Herbig Be star MWC 297 with spectral resolution 12 000. Astron. Astrophys. 527, A103.

Whelan, E.T., Ray, T.P., Davis, C.J., 2004. Paschen beta emission as a tracer of outflow activity from T-Tauri stars, as compared to optical forbidden emission. Astron. Astrophys. 417, 247–261.

Yan, D., Zhou, J., Zhang, P., 2019. Detection of γ-rays from the protostellar jet in the HH 80-81 system. arXiv e-prints, arXiv:1908.10994.

Zanni, C., Ferrari, A., Rosner, R., Bodo, G., Massaglia, S., 2007. MHD Simulations of jet acceleration from Keplerian accretion disks: the effects of disk resistivity. Astron. Astrophys. 469, 811.

Zanni, C., Ferreira, J., 2009. MHD simulations of accretion onto a dipolar magnetosphere. I. Accretion curtains and the disk-locking paradigm. Astron. Astrophys. 508, 1117–1133.

Zanni, C., Ferreira, J., 2011. Observational limits on the spin-down torque of accretion powered stellar winds. Astrophys. J. 727, L22.

Zanni, C., Ferreira, J., 2013. MHD simulations of accretion onto a dipolar magnetosphere. II. Magnetospheric ejections and stellar spin-down. Astron. Astrophys. 550, A99.

Zhang, Q., Claus, B., Watson, L., Moran, J., 2017. Angular momentum in disk wind revealed in the young star MWC 349A. Astrophys. J. 837 (1), 53.

Zhu, Z., Stone, J.M., 2018. Global evolution of an accretion disk with a net vertical field: Coronal accretion, flux transport, and disk winds. Astrophys. J. 857, 34.

<figure><img src="image_14.png" /></figure>

Tom Ray is Professor of Astrophysics at the Dublin Institute for Advanced Studies. He began his career in Radio Astronomy at Jodrell Bank before working at the University of Sussex and the Max Planck Institute for Astronomy in Heidelberg. His main area of interest is star formation and, in particular, the supersonic jets that stars like our sun produce when less than a million years old. Tom is also Co-Principal Investigator of the Mid-Infrared Instrument on the James Webb Space Telescope (JWST), a Co-Principal Investigator on the ARIEL mission to explore the atmospheres of exoplanets, and is building a new type of astronomical detector known as a Microwave Kinetic Inductance Detector (MKID). His other interests include ancient astronomical sites (such as Newgrange) and historical astronomy.