Skip to main content

On the origin of jets from disc-accreting magnetized stars

Abstract

A brief review of the origin of jets from disc-accreting rotating magnetized stars is given. In most models, the interior of the disc is characterized by a turbulent viscosity and magnetic diffusivity (‘alpha’ discs) whereas the coronal region outside the disc is treated using ideal magnetohydrodynamics (MHD). Extensive MHD simulations have established the occurrence of long-lasting outflows in the case of both slowly and rapidly rotating stars. (1) Slowly rotating stars exhibit a new type of outflow, conical winds. Conical winds are generated when stellar magnetic flux is bunched up by the inward motion of the accretion disc. Near their region of origin, the winds have a thin conical shell shape with half opening angle of 30. At large distances, their toroidal magnetic field collimates the outflow forming current carrying, matter dominated jets. These winds are predominantly magnetically and not centrifugally driven. About 10-30% of the disc matter from the inner disc is launched in the conical wind. Conical winds may be responsible for episodic as well as long lasting outflows in different types of stars. (2) Rapidly rotating stars in the ‘propeller regime’ exhibit twocomponent outflows. One component is similar to the matter dominated conical wind, where a large fraction of the disc matter may be ejected in this regime. The second component is a high-velocity, low-density magnetically dominated axial jet where matter flows along the open polar field lines of the star. The axial jet has a mass flux of about 10% that of the conical wind, but its energy flux, due to the Poynting flux, can be as large as for the conical wind. The jet’s magnetically dominated angular momentum flux causes the star to spin down rapidly. Propeller-driven outflows may be responsible for protostellar jets and their rapid spin-down.

When the artificial requirement of symmetry about the equatorial plane is dropped, the conical winds are found to come alternately from one side of the disc and then the other, even for the case where the stellar magnetic field is a centered axisymmetric dipole.

Recent MHD simulations of disc accretion to rotating stars in the propeller regime have been done with no turbulent viscosity and no diffusivity. The strong turbulence observed is due to the magneto-rotational instability. This turbulence drives accretion in the disc and leads to episodic conical winds and jets.

1 Review

Outflows in the form of jets and winds are observed from many disc accreting objects ranging from young stars to systems with white dwarfs, neutron stars and black holes. A large body of observations exists for outflows from young stars at different stages of their evolution, ranging from protostars, where powerful collimated outflows - jets - are observed, to classical T Tauri stars (CTTSs) where the outflows are weaker and often less collimated (see review by Ray et al. [2007]). Correlation between the disc’s radiated power and the jet power has been found in many CTTSs (Cabrit et al. [1990]; Hartigan et al. [1995]). A significant number of CTTSs show signs of outflows in spectral lines, in particular in He I where two distinct components of outflows had been found (Edwards et al. [2003], [2006], [2009]; Kwan et al. [2007]). Outflows are also observed from accreting compact stars such as accreting white dwarfs in symbiotic binaries (Sokoloski and Kenyon [2003]), or from the vicinity of neutron stars, such as from Circinus X-1 (Heinz et al. [2007]).

Different theoretical models have been proposed to explain the outflows from protostars and CTTSs (see review by Ferreira et al. [2006]). The commonly favored model for the origin of protostellar jets and outflows are the radially distributed magneto-centrifugal disc winds which originate from discs threaded by a poloidal magnetic field (Blandford and Payne [1982]; Königl and Pudritz [2000]). MHD simulations disc winds were pioneered by Shibata and Uchida ([1985]) and Uchida and Shibata ([1985]) who used a Lax-Wendroff method to solve the axisymmetric MHD equations for a sub-Keplerian disc initially threaded by a vertical magnetic field. Subsequently, a large number of MHD simulation studies of the disc winds have been carried out with different codes and different assumptions (e.g., Ustyugova et al. [1995], [1999]; Ouyed and Pudritz [1997]; Romanova et al. [1997]; Krasnopolsky et al. [1999]; Casse and Keppens [2004]; Ferreira et al. [2006]; Matt and Pudritz [2008a], [2008b]; Tzeferacos et al. [2009]).

A less favored model for the origin of protostellar jets discussed in this review is one where the jets originate from the innermost region of the accretion disc (Lovelace et al. [1991]) or the disc/magnetosphere boundary. This model is related to the X-wind model (Shu et al. [1994], [2007]; Najita and Shu [1994]; Cai et al. [2008]) where the outflow originates from the vicinity of the disc-magnetosphere boundary. Progress in understanding the theoretical models has come from MHD simulations of accretion discs around rotating magnetized stars as discussed below. Laboratory experiments are also providing insights into jet formation processes (Hsu and Bellan [2002]; Lebedev et al. [2005]) but these are not discussed here.

Outflows or jets from the disc-magnetosphere boundary were found in early axisymmetric MHD simulations by Hayashi et al. ([1996]) and Miller and Stone ([1997]). A one-time episode of outflows from the inner disc and inflation of the innermost field lines connecting the star and the disc were observed for a few dynamical time-scales. Somewhat longer simulation runs were performed by Goodson et al. ([1997], [1999]), Hirose et al. ([1997]), Matt et al. ([2002]) and Küker et al. ([2003]) where several episodes of field inflation and outflows were observed. These simulations hinted at a possible long-term nature for the outflows. However, the simulations were not sufficiently long to establish the behavior of the outflows. MHD simulations showing long-lasting (thousands of orbits of the inner disc) outflows from the disc-magnetosphere have been obtained by our group (Romanova et al. [2009]; Lii et al. [2012], [2014]) and independently by Fendt ([2009]). We obtained these outflows/jets in two main cases: (1) where the star rotates slowly but the field lines are bunched up into an X-type configuration, and (2) where the star rotates rapidly, in the ‘propeller regime’ (Illarionov and Sunyaev [1975]; Alpar and Shaham [1985]; Lovelace et al. [1999]). Field bunching occurs for conditions where the viscosity is larger than the magnetic diffusivity. Figure 1 shows sketches of the equatorial angular rotation rate Ω(r,z=0) of the plasma in the two cases. Here, r is the radius of the star; r m is the magnetospheric radius where the kinetic energy density of the disc matter is about equal to the energy density of the magnetic field; and r cr = ( G M / Ω 2 ) 1 / 3 is the co-rotation radius where the angular rotation rate of the star Ω equals that of the Keplerian disc Ω K = ( G M / r 3 ) 1 / 2 . For a slowly rotating star r m < r cr whereas for a rapidly rotating star in the propeller regime r m > r cr .

Figure 1
figure 1

Slowly/rapidly rotating stars. Schematic profiles of the midplane angular velocity of the plasma for the case of a slowly rotating star (left-hand panel) and a rapidly rotating star (right-hand panel) which is in the propeller regime. Here, Ω is the angular rotation rate of the star, Ω K is the Keplerian rotation rate of the disc, r is the star’s radius, r m is the radius of the magnetosphere, and r cr is the co-rotation radius.

Figure 2 shows examples of the outflows in the two cases. In both cases, two-component outflows are observed: One component originates at the inner edge of the disc near r m and has a narrow-shell conical shape close to the disc and therefore is termed a ‘conical wind’. It is matter dominated but can become collimated at large distances due to its toroidal magnetic field. The other component is a magnetically dominated high-velocity ‘axial jet’ which flows along the open stellar magnetic field lines. The axial jet may be very strong in the propeller regime. A full discussion of the simulations and analysis can be found in Romanova et al. ([2009]) and Lii et al. ([2012]).

Figure 2
figure 2

Two-component outflows observed in slowly (left) and rapidly (right) rotating magnetized stars. The background shows the poloidal matter flux F m =ρ v p , the arrows are the poloidal velocity vectors, and the lines are sample magnetic field lines. The labels point to the main outflow components. In the left-hand panel α ν =0.3 and α η =0.3 while in the right-hand panel α ν =0.1= α η .

The simulation codes used by our US/Russia group have been extensively tested and refined in many respects over the past fifteen years. The tests include the different well-known shock problems described for example by Mignone et al. ([2007]) in regard to the testing of the PLUTO code as well as the magnetic rotor tests described by Romanova et al. ([2009]). More importantly, our group has pioneered the detailed comparison of MHD simulation results (Ustyugova et al. [1999]) with the analytic theory of stationary axisymmetric MHD flows (Lovelace et al. [1986]). Furthermore, detailed analysis of the simulations have been made to evaluate the different forces acting to drive outflows and jets (e.g., Lii et al. [2012]). A major effort by our group has been to implement physically consistent boundary conditions at the outer boundaries of the simulation regions. We were the first to point out the necessity of having the fast-magnetosonic Mach cone of an outflow pointing outwards from the simulation region (Ustyugova et al. [1999]). A number of published MHD simulations of jets in long axial cylindrical regions violate this requirement and are therefore unphysical.

Section 2 describes the simulations. Section 3 discusses the conical winds and axial jets, the driving and collimation forces, and the variability of the winds and jets. Section 4 discusses lopsided jets. Section 5 gives the conclusions.

2 MHD simulations

We simulate the outflows resulting from disc-magnetosphere interaction by solving the equations of axisymmetric MHD on grids using a Godunov type method. Outside of the disc the flow is described by the equations of ideal MHD. Inside the disc the flow is described by the equations of viscous, resistive MHD. In an inertial reference frame the equations are:

ρ t +(ρv)=0,
(1)
( ρ v ) t +T=ρg,
(2)
B t ×(v×B)+×( η t ×B)=0,
(3)
( ρ S ) t +(ρSv)=Q.
(4)

Here, ρ is the density, S is the specific entropy, v is the flow velocity, B is the magnetic field, η t is the magnetic diffusivity, T is the momentum flux-density tensor, Q is the rate of change of entropy per unit volume, and g=(GM/ r 2 ) r ˆ is the gravitational acceleration due to the star which has mass M. In the simulations reviewed here it is assumed that the viscous plus Ohmic heating is balanced by radiative cooling so that Q=0. Most of the volume of the simulated flows does not have shocks and there is no shock heating; however, at the surface of the star where the funnel flows impact the star’s surface there are strong shocks and the shock heating is included (Koldoba et al. [2008]). The total mass of the disc is assumed to be negligible compared to M. Here, T is the sum of the ideal plasma terms and the α-viscosity terms discussed in the next paragraph. The plasma is considered to be an ideal gas with adiabatic index γ=5/3, and S=ln(p/ ρ γ ). We use spherical coordinates (r,θ,ϕ) with θ measured from the symmetry axis. The equations in spherical coordinates are given in Ustyugova et al. ([2006]).

Both the viscosity and the magnetic diffusivity of the disc plasma are considered to be due to turbulent fluctuations of the velocity and the magnetic field. Both effects are non-zero only inside the disc as determined by a density threshold. The microscopic transport coefficients are replaced by turbulent coefficients. The values of these coefficients are assumed to be given by the α-model of Shakura and Sunyaev ([1973]), where the coefficient of the turbulent kinematic viscosity is ν t = α ν c s 2 / Ω K , where c s is the isothermal sound speed and Ω K (r) is the Keplerian angular velocity. We take into account the viscous stress terms T r ϕ vis and T θ ϕ vis (Lii et al. [2012]). Similarly, the coefficient of the turbulent magnetic diffusivity η t = α η c s 2 / Ω K . Here, α ν and α η are dimensionless coefficients which are treated as parameters of the model. The inward advection of matter and large-scale magnetic field in accretion discs with different ( α ν , α η ) values has been studied by Dyda et al. ([2013]). Note that shearing box simulations by Guan and Gammie ([2009]) suggest that α ν α η . In the simulation studies of our group we have studied cases with ( α ν , α η ) in the ranges 0.03-0.3. For these values the viscosity and diffusivity are much larger than the numerical values due to the finite grids.

The MHD equations are solved in dimensionless form so that the results can be readily applied to different accreting stars (see Table 1). Equations (1)-(4) have been integrated numerically in spherical (r,θ,ϕ) coordinates using a Godunov-type numerical scheme. The flux densities of the different quantities are calculated using an eight-wave Roe-type approximate Riemann solver analogous to one described by Powell et al. ([1999]). The calculations were done in the region R in r R out , 0θπ/2. Matter flowing into the star is absorbed. The grid is uniform in the θ-direction with N θ cells. The N r cells in the radial direction have d r j + 1 =(1+0.0523)d r j (j=1,, N r ) so that the poloidal-plane cells are curvilinear rectangles with approximately equal sides. This choice results in high spatial resolution near the star where the disc-magnetosphere interaction takes place while also permitting a large simulation region. We have used a range of resolutions going from N r × N θ =51×31 to 121×51 in order to establish the numerical convergence of our results.

Table 1 Reference values for different types of stars

3 Conical winds and axial jets

A large number of simulations were done in order to understand the origin and nature of conical winds. All of the key parameters were varied in order to ensure that there is no special dependence on any parameter. We observed that the formation of conical winds is a common phenomenon for a wide range of parameters. They are most persistent and strong in cases where the viscosity and diffusivity coefficients are not very small, α ν 0.03, α η 0.03. Another important condition is that α ν α η ; that is, the magnetic Prandtl number of the turbulence, P m = α ν / α η 1. This condition favors the bunching of the stellar magnetic field by the accretion flow.

The velocities in the conical wind component are similar to those in conical winds around slowly rotating stars. Matter launched from the disc-magnetosphere region initially has an approximately Keplerian azimuthal velocity, v K = G M / r . It is gradually accelerated to poloidal velocities v p (0.3-0.5) v K and the azimuthal velocity decreases. The flow has a high density and carries most of the disc mass into the outflows. The situation is the opposite in the axial jet component where the density is 102-103 times lower, while the poloidal and total velocities are significantly higher. Thus we find a two-component outflow: a matter dominated conical wind and a magnetically dominated axial jet.

We observe conical winds in both slowly and rapidly rotating stars. In both cases, matter in the conical winds passes through the Alfvén surface (and shortly thereafter through the fast magnetosonic point), beyond which the flow is matter-dominated in the sense that the energy flow is carried mainly by the matter. The situation is different for the axial jet component where the flow is sub-Alfvénic within the simulation region. For this component the energy flow is carried by the Poynting flux and the angular moment flow is carried by the magnetic field. X-ray observations of the jet from the protostellar object L1551 IRS 5 suggest a high-velocity, highly collimated inner jet and a lower-velocity, less-collimated outer outflow component (Schneider et al. [2011]).

Collimation and driving of the outflows Figure 3 shows the long-distance development of a conical wind from a slowly rotating star. At large distances the conical wind becomes collimated. To understand the collimation we analyzed total force (per unit mass) perpendicular to a poloidal magnetic field line (Lii et al. [2012]). For distances beyond the Alfvén surface of the flow this force is approximately

f tot , = v p 2 Θ s 1 8 π ρ B p 2 n 1 8 π ρ ( r sin θ ) 2 ( r sin θ B ϕ ) 2 n + v ϕ 2 r cos Θ sin θ
(5)

(Ustyugova et al. [1999]). Here, Θ is the angle between the poloidal magnetic field and the symmetry axis, s is the arc length along the poloidal field line, n is a coordinate normal to the poloidal field, and the p-subscripts indicate the poloidal component of a vector. Once the jet begins to collimate, the curvature term v p 2 Θ/s also becomes negligible. The magnetic force may act to either collimate or decollimate the jet, depending on the relative magnitudes of the toroidal ( r sin θ B ϕ ) 2 gradient (which collimates the outflow) and poloidal B p 2 gradient (which ‘decollimates’). In our simulations, the collimation of the matter implies that the magnetic hoop stress is larger than the poloidal field gradient. Thus the main perpendicular forces acting in the jet are the collimating effect of the toroidal magnetic field and the decollimating effect of the centrifugal force and the gradient of B p 2 . The collimated effect of B ϕ dominates. Note that in MKS units 2πrsinθ B ϕ / μ 0 is the poloidal current flowing through a surface of radius r from colatitude zero to θ. For the jets from young stars this current is of the order of 2× 10 13  A.

Figure 3
figure 3

The conical wind/jet from a slowly rotating star at time t=860 . The background shows the poloidal matter flux density ρ v p and the lines show the poloidal projections of the magnetic field. The red vectors show the poloidal matter velocity v p . Dimensional values can be obtained from Table 1. For example for a CTTS, t 0 =0.366 days, R 0 =2 R , t=860 corresponds to 315 days, and the simulation region is 0.39 AU in radius. The horizontal axis shows the distance from the star in units of the reference radius R 0 . For this case α ν =0.3 and α η =0.1.

The driving force for the outflow is simply the force parallel to the poloidal magnetic field of the flow f tot , . This is obtained by taking the dot product of the Euler equation with the b ˆ unit vector which is parallel to the poloidal magnetic field line B p . The derivation by Ustyugova et al. ([1999]) gives

f tot , = 1 ρ P s Φ s + v ϕ 2 r sin θ sin Θ + 1 4 π ρ b ˆ [ ( × B ) × B ] .
(6)

Here, the terms on the right-hand side correspond to the pressure, gravitational, centrifugal and magnetic forces, respectively denoted f P , G , C , M . The pressure gradient force, f P , dominates within the disk. The matter in the disk is approximately in Keplerian rotation such that the sum of the gravitational and centrifugal forces roughly cancel ( f G + C 0). Near the slowly rotating star, however, the matter is strongly coupled to the stellar magnetic field and the disk orbits at sub-Keplerian speeds, giving f G + C 0. The magnetic driving force (the last term of Eq. (6)) can be expanded as

f M , = 1 8 π ρ ( r sin θ ) 2 s ( r sin θ B ϕ ) 2
(7)

(Lovelace et al. [1991]). Figure 4 shows the variation of the total force f tot , the gravitational plus centrifugal force, and the magnetic force along a representative field line. This analysis establishes that the predominant driving force for the outflow is the magnetic force (Eq. (7)) and not the centrifugal force. This in agreement with the analysis of Lovelace et al. ([1991]).

Figure 4
figure 4

Forces along a field line in the jet. Panel (a) shows the poloidal matter flux density ρ| v p | as a background overplotted with poloidal magnetic field lines. The vectors show the total force f tot along a representative field line originating from the disk at r=1.24. Panel (b) plots the angular velocity Ω as the background. The vectors show the sum of the gravitational + centrifugal forces f G + C along the representative field line. Panel (c) shows the poloidal current I p as the background. The vectors show the total magnetic force f M along the representative field line.

Variability For both rapidly and slowly rotating stars the magnetic field lines connecting the disc and the star have the tendency to inflate and open (Lovelace et al. [1995]). Quasi-periodic reconstruction of the magnetosphere due to inflation and reconnection has been discussed theoretically (Aly and Kuijpers [1990]) and has been observed in a number of axisymmetric simulations (Hirose et al. [1997]; Goodson et al. [1997], [1999]; Matt et al. [2002]; Romanova et al. [2002]). Goodson and Winglee ([1999]) discuss the physics of inflation cycles. They have shown that each cycle of inflation consists of a period of matter accumulation near the magnetosphere, diffusion of this matter through the magnetospheric field, inflation of the corresponding field lines, accretion of some matter onto the star, and outflow of some matter as winds, with subsequent expansion of the magnetosphere. There simulations show 5-6 cycles of inflation and reconnection. Our simulations often show 30-50 cycles of inflation and reconnection.

Kurosawa and Romanova ([2012]) have calculated spectra from modeled conical winds and accretion funnels combining the 3D MHD simulations with 3D radiative transfer code TORUS. They have shown that conical winds may explain different features in the hydrogen spectral lines, in the He I line and also a relatively narrow, low-velocity blue-shifted absorption components in the He I λ 10830 which is often seen in observations (Kurosawa et al. [2011]). Further, the 3D MHD+3D radiative transfer codes have been used to model the young star V2129 Oph, where the parameters of the star including the surface magnetic field distribution are known (Alencar et al. [2012]). The spectrum in several Hydrogen lines was calculated and compared it with observed spectrum. A good match was obtained between the modeled and observed spectra (Alencar et al. [2012]).

4 Lopsided jets and outflows from discs

There is clear evidence, mainly from Hubble Space Telescope (HST) observations, of the asymmetry between the approaching and receding jets from a number of young stars. The objects include the jets in HH 30 (Bacciotti et al. [1999]), RW Aur (Woitas et al. [2002]), TH 28 (Coffey et al. [2004]), and LkHα 233 (Perrin and Graham [2007]). Specifically, the radial speed of the approaching jet may differ by a factor of two from that of the receding jet. For example, for RW Aur the radial redshifted speed is 100 km/s whereas the blueshifted radial speed is 175 km/s. The mass and momentum fluxes are also significantly different for the approaching and receding jets in a number of cases. It is possible that the observed asymmetry of the jets could be due to differences in the gas densities on the two sides of the source. However, it is more likely that the asymmetry of the outflows arises from the asymmetry of the star’s magnetic field. Substantial observational evidence points to the fact that young stars often have complex magnetic fields consisting of dipole, quadrupole, and higher order poles misaligned with respect to each other and the rotation axis (Jardine et al. [2002]; Donati et al. [2008]). Analysis of the plasma flow around stars with realistic fields have shown that a significant fraction of the star’s magnetic field lines are open and may carry outflows (Gregory et al. [2006]).

The complex magnetic field of a star will destroy the commonly assumed symmetry of the magnetic field and the plasma about the equatorial plane. MHD simulations by Lovelace et al. ([2010]) fully support the qualitative picture suggested in the sketch in Figure 1 of Lovelace et al. ([2010]). The idea of mixing of even and odd symmetry magnetic fields about z=0 to get lopsided outflows was proposed earlier by Wang et al. ([1992]). The time-scale during which the jet comes from the upper hemisphere is set by the evolution time-scale for the stellar magnetic field. This is determined by the dynamo processes responsible for the generation of the field. Remarkably, once the assumption of symmetry about the equatorial plane is dropped, the conical winds alternately come from one hemisphere and then the other even when the stellar magnetic field is a centered axisymmetric dipole (Lovelace et al. [2010]). Fendt and Sheikhnezami ([2013]) likewise found that symmetric magnetic field configurations produced asymmetric outflows if there were thermal asymetries in the disc. The time-scale for the ‘flipping’ is the accretion time-scale of the inner part of the disc which is expected to be much less than the evolution time of the star’s magnetic field.

We have revisited the problem of the asymmetry of the jets and outflows using a new axisymmetric code with a high-resolution stretched-grid (Dyda et al.: Bipolar MHD Outflows from T Tauri Stars, in preparation). The star has a radius of 1 in our simulation units, and the first 30 grid cells have lengths δR=0.1. At larger R, the cell lengths are given recursively by δ R i + 1 =1.025δ R i . Similarly, in the Z-direction, the first 30 grid cells above and below the equatorial plane have lengths δZ=0.1. At larger |Z|, the cell lengths are given recursively by δ Z j + 1 =1.025δ Z j . This grid gives high resolution in the region occupied by the disc and by the jet. Here, we present sample results with as simulation region of 135×216 (R×Z) cells. Figure 5 shows a sparse version of this grid.

Figure 5
figure 5

Diluted picture of new high-resolution mesh used to study lopsided jets and outflows. The actual mesh has 5-times as many cells in the R-direction and 6-times as many cells in the Z-direction.

The initial magnetic field is taken to be a superposition of a dipole field centered in the star described by the flux function

Ψ D = μ R 2 ( R 2 + Z 2 ) 3 / 2 ,
(8)

and a Zanni-type distributed field in the disc,

Ψ Z = 4 3 B 0 r 0 2 ( r r 0 ) 3 / 4 m 5 / 4 ( m 2 + Z 2 / R 2 ) 5 / 8
(9)

(Zanni et al. [2007]), where μ is the magnetic moment of the star, B 0 is a reference value for the disc field, and m is a dimensionless parameter which controls the initial disc field geometry.

Figure 6 shows a zoomed-in view of results from the new stretched-grid code for an episode of lopsided jet formation for a case where the dipole field of the rotating star is parallel to the disc field in the disc midplane. The outflow from the top side of the disc is at super escape speed velocities with the result that the mass outflow is predominantly from the top side of the disc. During the duration of the run (203 rotation periods of the disc at the corotation radius), the poloidal flux of the disc field advects inward and accumulates in the low-density axial blue region in the figure. In this region there are magnetically collimated Poynting flux outflows of energy and angular momentum in the ±z directions.

Figure 6
figure 6

Lopsided outflow from the disc-magnetosphere boundary at a time corresponding to 203 rotation periods of the disc at the corotation radius r cr , which is twice the star’s radius  r . The initial magnetic field, which is symmetric about z=0, is the sum of the stellar dipole (Eq. (8)) and the disc field (Eq. (9)), where the dipole and disc fields are parallel and equal in magnitude at r=4.4 r and z=0. The color background shows the density and the white lines the poloidal projection of the magnetic field lines. The arrows inside the disc (the orange region) are the poloidal velocity vectors, while the arrows external to the disc show the momentum density ρ v.

5 Conclusions

Detailed magnetohydrodynamic simulations have established that long-lasting outflows of cold disc matter are ejected into a hot, low-density corona from the disc-magnetosphere boundary in the case of both slowly and rapidly rotating stars. The main results are the following:

For slowly rotating stars a new type of outflow - a conical wind - has been discovered. Matter flows out forming a conical wind which has the shape of a thin conical shell with a half-opening angle θ 30 . The outflows appear in cases where the magnetic flux of the star is bunched up by the inward accretion flow of the disc. We find that this occurs when the turbulent magnetic Prandtl number (the ratio of viscosity to diffusivity) P m >1, and when the viscosity is sufficiently high, α ν 0.03.

Winds from the disc-magnetosphere boundary have been proposed earlier by Shu and collaborators and referred to as X-winds (Shu et al. [1994]). In this model, the wind originates from a small region near the corotation radius r cr , while the disc truncation radius r t (or, the magnetospheric radius r m ) is only slightly smaller than r cr ( r m 0.7 r cr , Shu et al. [1994]). It is suggested that excess angular momentum flows from the star to the disc and from there into the X-winds. The model aims to explain the slow rotation of the star and the formation of jets. In the simulations discussed here we have obtained outflows from both slowly and rapidly rotating stars. Both have conical wind components which are reminiscent of X-winds. In some respects the conical winds are similar to X-winds: They both require bunching of the poloidal field lines and show outflows from the inner disc; and they both have high rotation and show gradual poloidal acceleration (e.g., Najita and Shu [1994]).

The main differences are the following: (1) The conical/propeller outflows have two components: a slow high-density conical wind (which can be considered as an analogue of the X-wind), and a fast low-density jet. No jet component is discussed in the X-wind model. (2) Conical winds form around stars with any rotation rate including very slowly rotating stars. They do not require fine tuning of the corotation and truncation radii. For example, bunching of field lines is often expected during periods of enhanced or unstable accretion when the disc comes closer to the surface of the star and r m r cr . Under this condition conical winds will form. In contrast, X-winds require r m r cr . (3) The base of the conical wind component in both slowly and rapidly rotating stars is associated with the region where the field lines are bunched up, and not with the corotation radius. (4) X-winds are driven by the centrifugal force, and as a result matter flows over a wide range of directions below the ‘dead zone’ (Shu et al. [1994]; Ostriker and Shu [1995]). In conical winds the matter is driven by the magnetic force (Lovelace et al. [1991]) which acts such that the matter flows into a thin shell with a cone half-angle θ 30 . The same force tends to collimate the flow.

For rapidly rotating stars in the propeller regime where r m > r cr and where the condition for bunching, P m >1, is satisfied we find two distinct outflow components (1) a relatively low-velocity conical wind and (2) a high-velocity axial jet. A significant part of the disc matter and angular momentum flows into the conical winds. At the same time a significant part of the rotational energy of the star flows into the magnetically-dominated axial jet. This regime is particularly relevant to protostars, where the star rotates rapidly and has a high accretion rate. The star spins down rapidly due to the angular momentum flow into the axial jet along the field lines connecting the star and the corona. For typical parameters a protostar spins down in 3× 10 5 years. The axial jet is powered by the spin-down of the star rather than by disc accretion. The matter fluxes into both components (wind and jet) strongly oscillate due to events of inflation and reconnection. Most powerful outbursts occur every 1-2 months. The interval between outbursts is expected to be longer for smaller diffusivities in the disc. Outbursts are accompanied by higher outflow velocities and stronger self-collimation of both components. Such outbursts may explain the ejection of knots in CTTSs every few months.

When the artificial requirement of symmetry about the equatorial plane is dropped, MHD simulations reveal that the conical winds may alternately come from one side of the disc and then the other even for the case where the stellar magnetic field is a centered axisymmetric dipole (Lovelace et al. [2010]; Fendt and Sheikhnezami [2013]; Dyda et al.: Bipolar MHD Outflows from T Tauri Stars, in preparation).

In recent work we have studied the disc accretion to rotating magnetized stars in the propeller regime using a new code with very high resolution in the region of the disc. In this code no turbulent viscosity or diffusivity is incorporated, but instead strong turbulence occurs due to the magneto-rotational instability. This turbulence drives the accretion and it leads to episodic outflows. The effective ( α ν , α η ) values due to the turbulence arising from the magneto-rotational instability (MRI) are found to be 0.1. These values are much larger than the numerical viscosity and diffusivity values due to the finite grids used which are 0.01. Note however that the characterization of the turbulence by α-values is a rough approximation.

References

  1. Alencar SHP, Bouvier J, Walter FM, Dougados C, Donati J-F, Kurosawa R, Romanova M, Bonfils X, Lima GHRA, Massaro S, Ibrahimov M, Poretti E: Accretion dynamics in the classical T Tauri star V2129 Ophiuchi. Astron. Astrophys. 2012., 541: Article ID A116 10.1051/0004-6361/201118395

    Google Scholar 

  2. Alpar MA, Shaham J: Is GX5–1 a millisecond pulsar? Nature 1985, 316: 239–241. 10.1038/316239a0

    Article  ADS  Google Scholar 

  3. Aly JJ, Kuijpers J: Flaring interactions between accretion disk and neutron star magnetosphere. Astron. Astrophys. 1990, 227: 473–482.

    ADS  Google Scholar 

  4. Bacciotti F, Eislöffel J, Ray TP: The physical properties of the HH 30 jet from HST and ground-based data. Astron. Astrophys. 1999, 350: 917–927.

    ADS  Google Scholar 

  5. Blandford RD, Payne DG: Hydromagnetic flows from accretion discs and the production of radio jets. Mon. Not. R. Astron. Soc. 1982, 199: 883–903. 10.1093/mnras/199.4.883

    Article  ADS  Google Scholar 

  6. Cabrit S, Edwards S, Strom SE, Strom KM: Forbidden-line emission and infrared excesses in T Tauri stars - evidence for accretion-driven mass loss? Astrophys. J. 1990, 354: 687–700. 10.1086/168725

    Article  ADS  Google Scholar 

  7. Cai MJ, Shang H, Lin H-H, Shu FH: X-winds in action. Astrophys. J. 2008, 672: 489–503. 10.1086/523788

    Article  ADS  Google Scholar 

  8. Casse F, Keppens R: Radiatively inefficient magnetohydrodynamic accretion-ejection structures. Astrophys. J. 2004, 601: 90–103. 10.1086/380441

    Article  ADS  Google Scholar 

  9. Coffey D, Bacciotti F, Woitas J, Ray TP, Eislöffel J: Rotation of jets from young stars: new clues from the Hubble Space Telescope imaging spectrograph. Astrophys. J. 2004, 604: 758–765. 10.1086/382019

    Article  ADS  Google Scholar 

  10. Donati J-F, Jardine MM, Gregory SG, Petit P, Paletou F, Bouvier J, Dougados C, Mènard F, Cameron AC, Harries TJ, Hussain GAJ, Unruh Y, Morin J, Marsden SC, Manset N, Aurière M, Catala C, Alecian E: Magnetospheric accretion on the T Tauri star BP Tauri. Mon. Not. R. Astron. Soc. 2008, 386: 1234–1251. 10.1111/j.1365-2966.2008.13111.x

    Article  ADS  Google Scholar 

  11. Dyda S, Lovelace RVE, Ustyugova GV, Lii PS, Romanova MM, Koldoba AV: Advection of matter and B -fields in alpha-discs. Mon. Not. R. Astron. Soc. 2013, 432: 127–137. 10.1093/mnras/stt429

    Article  ADS  Google Scholar 

  12. Edwards S, Fischer W, Hillenbrand L, Kwan J: Probing T Tauri accretion and outflow with 1 micron spectroscopy. Astrophys. J. 2006, 646: 319–341. 10.1086/504832

    Article  ADS  Google Scholar 

  13. Edwards S, Fischer W, Kwan J, Hillenbrandt L, Durpee AK: He I λ 10830 as a probe of winds in accreting young stars. Astrophys. J. 2003, 599: L41-L44. 10.1086/381077

    Article  ADS  Google Scholar 

  14. Edwards S: Winds and accretion in young stars. In Proceedings of the 15th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun. AIP Publishing, Melville; 2009:29–38.

    Google Scholar 

  15. Fendt C: Formation of protostellar jets as two-component outflows from star-disk magnetospheres. Astrophys. J. 2009, 692: 346–363. 10.1088/0004-637X/692/1/346

    Article  ADS  Google Scholar 

  16. Fendt C, Sheikhnezami S: Bipolar jets launched from accretion disks. II. The formation of asymmetric jets and counter jets. Astrophys. J. 2013., 774: Article ID 12 10.1088/0004-637X/774/1/12

    Google Scholar 

  17. Ferreira J, Dougados C, Cabrit S: Which jet launching mechanism(s) in T Tauri stars? Astron. Astrophys. 2006, 453: 785–796. 10.1051/0004-6361:20054231

    Article  ADS  Google Scholar 

  18. Goodson AP, Winglee RM: Jets from accreting magnetic young stellar objects. II. Mechanism physics. Astrophys. J. 1999, 524: 159–168. 10.1086/307780

    Article  ADS  Google Scholar 

  19. Goodson AP, Winglee RM, Böhm K-H: Time-dependent accretion by magnetic young stellar objects as a launching mechanism for stellar jets. Astrophys. J. 1997, 489: 199–209. 10.1086/304774

    Article  ADS  Google Scholar 

  20. Goodson AP, Böhm K-H, Winglee RM: Jets from accreting magnetic young stellar objects. I. Comparison of observations and high-resolution simulation results. Astrophys. J. 1999, 524: 142–158. 10.1086/307779

    Article  ADS  Google Scholar 

  21. Gregory SG, Jardine M, Simpson I, Donati J-F: Mass accretion on to T Tauri stars. Mon. Not. R. Astron. Soc. 2006, 371: 999–1013. 10.1111/j.1365-2966.2006.10734.x

    Article  ADS  Google Scholar 

  22. Guan X, Gammie CF: The turbulent magnetic Prandtl number of MHD turbulence in disks. Astrophys. J. 2009, 607: 1901–1906. 10.1088/0004-637X/697/2/1901

    Article  ADS  Google Scholar 

  23. Hartigan P, Edwards S, Gandhour L: Disk accretion and mass loss from young stars. Astrophys. J. 1995, 452: 736–768. 10.1086/176344

    Article  ADS  Google Scholar 

  24. Hayashi MR, Shibata K, Matsumoto R: X-ray flares and mass outflows driven by magnetic interaction between a protostar and its surrounding disk. Astrophys. J. 1996, 468: L37-L40. 10.1086/310222

    Article  ADS  Google Scholar 

  25. Heinz S, Schulz NS, Brandt WN, Galloway DK: Evidence of a parsec-scale X-ray jet from the accreting neutron star Circinus X-1. Astrophys. J. 2007, 663: L93-L96. 10.1086/519950

    Article  ADS  Google Scholar 

  26. Hirose S, Uchida Y, Shibata K, Matsumoto R: Disk accretion onto a magnetized young star and associated jet formation. Publ. Astron. Soc. Jpn. 1997, 49: 193–205.

    Article  ADS  Google Scholar 

  27. Hsu SC, Bellan PM: A laboratory plasma experiment for studying magnetic dynamics of accretion discs and jets. Mon. Not. R. Astron. Soc. 2002, 334: 257–261. 10.1046/j.1365-8711.2002.05422.x

    Article  ADS  Google Scholar 

  28. Illarionov AF, Sunyaev RA: Why the number of galactic X-ray stars is so small? Astron. Astrophys. 1975, 39: 185–195.

    ADS  Google Scholar 

  29. Jardine M, Collier Cameron A, Donati J-F: The global magnetic topology of AB Doradus. Mon. Not. R. Astron. Soc. 2002, 333: 339–346. 10.1046/j.1365-8711.2002.05394.x

    Article  ADS  Google Scholar 

  30. Koldoba AV, Ustyugova GV, Romanova MM, Lovelace RVE: Oscillations of magnetohydrodynamic shock waves on the surfaces of T Tauri stars. Mon. Not. R. Astron. Soc. 2008, 388: 357–366. 10.1111/j.1365-2966.2008.13394.x

    Article  ADS  Google Scholar 

  31. Königl A, Pudritz RE: Disk winds and the accretion-outflow connection. In Protostars and Planets IV. Edited by: Mannings V, Boss AP, Russell SS. University of Arizona Press, Tucson; 2000:759–791.

    Google Scholar 

  32. Krasnopolsky R, Li Z-Y, Blandford R: Magnetocentrifugal launchinng of jets from accretion disks: I. Cold axisymmetric flows. Astrophys. J. 1999, 526: 631–641. 10.1086/308023

    Article  ADS  Google Scholar 

  33. Küker M, Henning T, Rüdiger G: Magnetic star-disk coupling in classical T Tauri systems. Astrophys. J. 2003, 589: 397–409. 10.1086/374408

    Article  ADS  Google Scholar 

  34. Kurosawa R, Romanova MM, Harries T: Multidimensional models of hydrogen and helium emission line profiles for classical T-Tauri stars: method, tests and examples. Mon. Not. R. Astron. Soc. 2011, 416: 2623–2639. 10.1111/j.1365-2966.2011.19216.x

    Article  ADS  Google Scholar 

  35. Kurosawa R, Romanova MM: Line formation in the inner winds of classical T Tauri stars: testing the conical-shell wind solution. Mon. Not. R. Astron. Soc. 2012, 426: 2901–2916. 10.1111/j.1365-2966.2012.21853.x

    Article  ADS  Google Scholar 

  36. Kwan J, Edwards S, Fischer W: Modeling T Tauri winds from He I λ 10830 profiles. Astrophys. J. 2007, 657: 897–915. 10.1086/511057

    Article  ADS  Google Scholar 

  37. Lebedev SV, Ciardi A, Ampleford DJ, Bland SN, Bott SC, Chittenden JP, Hall GN, Rapley J, Jennings CA, Frank A, Blackman EG, Lery T: Magnetic tower outflows from a radial wire array Z-pinch. Mon. Not. R. Astron. Soc. 2005, 361: 97–108. 10.1111/j.1365-2966.2005.09132.x

    Article  ADS  Google Scholar 

  38. Lii PS, Romanova M, Lovelace R: Magnetic launching and collimation of jets from the disc-magnetosphere boundary: 2.5D MHD simulations. Mon. Not. R. Astron. Soc. 2012, 420: 2020–2033. 10.1111/j.1365-2966.2011.20133.x

    Article  ADS  Google Scholar 

  39. Lii PS, Romanova MM, Ustyugova GV, Koldoba AV, Lovelace RVE: Propeller outflows from an MRI disc. Mon. Not. R. Astron. Soc. 2014, 441: 86–100. 10.1093/mnras/stu495

    Article  ADS  Google Scholar 

  40. Lovelace RVE, Mehanian C, Mobarry CM, Sulkanen ME: Theory of axisymmetric magnetohydrodynamic flows disks. Astrophys. J. Suppl. Ser. 1986, 62: 1–37. 10.1086/191132

    Article  ADS  Google Scholar 

  41. Lovelace RVE, Berk HL, Contopoulos J: Magnetically driven jets and winds. Astrophys. J. 1991, 379: 696–705. 10.1086/170544

    Article  ADS  Google Scholar 

  42. Lovelace RVE, Romanova MM, Bisnovatyi-Kogan GS: Spin-up/spin-down of magnetized stars with accretion discs and outflows. Mon. Not. R. Astron. Soc. 1995, 275: 244–254.

    ADS  Google Scholar 

  43. Lovelace RVE, Romanova MM, Bisnovatyi-Kogan GS: Magnetic propeller outflows. Astrophys. J. 1999, 514: 368–372. 10.1086/306945

    Article  ADS  Google Scholar 

  44. Lovelace RVE, Romanova MM, Ustyugova GV, Koldoba AV: One-sided outflows/jets from rotating stars with complex magnetic fields. Mon. Not. R. Astron. Soc. 2010, 408: 2083–2091. 10.1111/j.1365-2966.2010.17284.x

    Article  ADS  Google Scholar 

  45. Matt S, Goodson AP, Winglee RM, Böhm K-H: Simulation-based investigation of a model for the interaction between stellar magnetospheres and circumstellar accretion disks. Astrophys. J. 2002, 574: 232–245. 10.1086/340896

    Article  ADS  Google Scholar 

  46. Matt S, Pudritz RE: Accretion-powered stellar winds. II. Numerical solutions for stellar wind torques. Astrophys. J. 2008, 678: 1109–1118. 10.1086/533428

    Article  ADS  Google Scholar 

  47. Matt S, Pudritz RE: Accretion-powered stellar winds. III. Spin-equilibrium solutions. Astrophys. J. 2008, 681: 391–399. 10.1086/587453

    Article  ADS  Google Scholar 

  48. Mignone A, Bodo G, Massaglia S, Matsakos T, Tesileanu O, Zanni C, Ferrari A: PLUTO: a numerical code for computational astrophysics. Astrophys. J. Suppl. Ser. 2007, 170: 228–242. 10.1086/513316

    Article  ADS  Google Scholar 

  49. Miller KA, Stone JM: Magnetohydrodynamic simulations of stellar magnetosphere-accretion disk interaction. Astrophys. J. 1997, 489: 890–902. 10.1086/304825

    Article  ADS  Google Scholar 

  50. Najita JR, Shu FH: Magnetocentrifugally driven flows from young stars and disks. 3: Numerical solution of the sub-Alfvenic region. Astrophys. J. 1994, 429: 808–825. 10.1086/174365

    Article  ADS  Google Scholar 

  51. Ostriker EC, Shu FH: Magnetocentrifugally driven flows from young stars and disks. IV. The accretion funnel and dead zone. Astrophys. J. 1995, 447: 813–828. 10.1086/175920

    Article  ADS  Google Scholar 

  52. Ouyed R, Pudritz RE: Numerical simulations of astrophysical jets from Keplerian disks. I. Stationary models. Astrophys. J. 1997, 482: 712–732. 10.1086/304170

    Article  ADS  Google Scholar 

  53. Perrin MD, Graham JR: Laser guide star adaptive optics integral field spectroscopy of a tightly collimated bipolar jet from the Herbig Ae star LkH α 233. Astrophys. J. 2007, 670: 499–508. 10.1086/521643

    Article  ADS  Google Scholar 

  54. Powell KG, Roe PL, Linde TJ, Gombosi TI, De Zeeuw DL: A solution-adaptive upwind scheme for ideal magnetohydrodynamics. J. Comput. Phys. 1999, 154: 284–309. 10.1006/jcph.1999.6299

    Article  MathSciNet  ADS  Google Scholar 

  55. Ray, TP, Dougados, C, Bacciotti, F, Eisloffel, J, Chrysotomuo, A: Protostars and Planets V, 231 (2007)

  56. Romanova MM, Ustyugova GV, Koldoba AV, Chechetkin VM, Lovelace RVE: Formation of stationary magnetohydrodynamic outflows from a disk by time-dependent simulations. Astrophys. J. 1997, 482: 708–711. 10.1086/304199

    Article  ADS  Google Scholar 

  57. Romanova MM, Ustyugova GV, Koldoba AV, Lovelace RVE: Magnetohydrodynamic simulations of disk-magnetized star interactions in the quiescent regime: funnel flows and angular momentum transport. Astrophys. J. 2002, 578: 420–438. 10.1086/342464

    Article  ADS  Google Scholar 

  58. Romanova MM, Ustyugova GV, Koldoba AV, Lovelace RVE: Launching of conical winds and axial jets from the disc-magnetosphere boundary: axisymmetric and 3D simulations. Mon. Not. R. Astron. Soc. 2009, 399: 1802–1828. 10.1111/j.1365-2966.2009.15413.x

    Article  ADS  Google Scholar 

  59. Schneider PC, Günther HM, Schmitt HMM: The X-ray puzzle of the L1551 IRS 5 jet. Astron. Astrophys. 2011., 530: Article ID A123 10.1051/0004-6361/201016305

    Google Scholar 

  60. Shakura NI, Sunyaev RA: Black holes in binary systems: observational appearance. Astron. Astrophys. 1973, 24: 337–355.

    ADS  Google Scholar 

  61. Shibata K, Uchida Y: A magnetodynamic mechanism for the formation of astrophysical jets. I. Dynamical effects of the relaxation of nonlinear magnetic twists. Publ. Astron. Soc. Jpn. 1985, 37: 31–46.

    ADS  Google Scholar 

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

    Article  ADS  Google Scholar 

  63. Shu FH, Galli D, Lizano S, Glassgold AE, Diamond PH: Mean field magnetohydrodynamics of accretion disks. Astrophys. J. 2007, 665: 535–553. 10.1086/519678

    Article  ADS  Google Scholar 

  64. Sokoloski JL, Kenyon SJ: CH Cygni. I: Observational evidence for a disk-jet connection. Astrophys. J. 2003, 584: 1021–1026. 10.1086/345901

    Article  ADS  Google Scholar 

  65. Tzeferacos P, Ferrari A, Mignone A, Zanni C, Bodo G, Massaglia S: On the magnetization of jet-launching discs. Mon. Not. R. Astron. Soc. 2009, 400: 820–834. 10.1111/j.1365-2966.2009.15502.x

    Article  ADS  Google Scholar 

  66. Uchida Y, Shibata K: Magnetodynamical acceleration of CO and optical bipolar flows from the region of star formation. Publ. Astron. Soc. Jpn. 1985, 37: 515–535.

    ADS  Google Scholar 

  67. Ustyugova GV, Koldoba AV, Romanova MM, Chechetkin VM, Lovelace RVE: Magnetohydrodynamic simulations of outflows from accretion disks. Astrophys. J. 1995, 439: L39-L42. 10.1086/187739

    Article  ADS  Google Scholar 

  68. Ustyugova GV, Koldoba AV, Romanova MM, Chechetkin VM, Lovelace RVE: Magnetocentrifugally driven winds: comparison of MHD simulations with theory. Astrophys. J. 1999, 516: 221–235. 10.1086/307093

    Article  ADS  Google Scholar 

  69. Ustyugova GV, Koldoba AV, Romanova MM, Lovelace RVE: ‘Propeller’ regime of disk accretion to rapidly rotating stars. Astrophys. J. 2006, 646: 304–318. 10.1086/503379

    Article  ADS  Google Scholar 

  70. Wang JCL, Sulkanen ME, Lovelace RVE: Intrinsically asymmetric astrophysical jets. Astrophys. J. 1992, 390: 46–65. 10.1086/171258

    Article  ADS  Google Scholar 

  71. Woitas J, Ray TP, Bacciotti F, Davis CJ, Eislöffel J: Hubble Space Telescope space telescope imaging spectrograph observations of the bipolar jet from RW Aurigae: tracing outflow asymmetries close to the source. Astrophys. J. 2002, 580: 336–342. 10.1086/343124

    Article  ADS  Google Scholar 

  72. Zanni C, Ferrari A, Rosner R, Bodo G, Massaglia S: MHD simulations of jet acceleration from Keplerian accretion disks. The effects of disk resistivity. Astron. Astrophys. 2007, 469: 811–828. 10.1051/0004-6361:20066400

    Article  ADS  Google Scholar 

Download references

Acknowledgements

The authors thank GV Ustyugova and AV Koldoba for the development of the codes used in the reviewed simulations. This research was supported in part by NSF grants AST-1008636 and AST-1211318 and by a NASA ATP grant NNX10AF63G; we thank NASA for use of the NASA High Performance Computing Facilities.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Richard VE Lovelace.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RL and MR were principal investigators of this research and drafted this manuscript. PL performed the simulations and analysis of the conical winds and axial jets. SD performed the simulations and analysis of the lopsided jets and disc outflows. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lovelace, R.V., Romanova, M.M., Lii, P. et al. On the origin of jets from disc-accreting magnetized stars. Comput. Astrophys. 1, 3 (2014). https://doi.org/10.1186/s40668-014-0003-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40668-014-0003-5

Keywords