We have developed a numerical MHD model of the propeller candidate star AE Aqr using axisymmetric magneto-hydrodynamic (MHD) simulations. We suggest that AE Aqr is an intermediate polar-type star, where the magnetic field is relatively weak and an accretion disc may form around the white dwarf. The star is in the propeller regime, and many of its observational properties are determined by the disc-magnetosphere interaction. Comparisons of the characteristics of the observed versus modelled AE Aqr star show that the model can explain many observational properties of AE Aqr. In a representative model, the magnetic field of the star is B\approx 3.3x10^5 G and the time averaged accretion rate in the disc is 5.5\times 10^{16} g/s. Most of this matter is ejected into conically-shaped winds. The numerical model explains the rapid spin-down of AE Aqr through the outflow of angular momentum from the surface of the star to the wind, corona and disc. The energy budget in the outflows, 9x10^{33} erg/s, is sufficient for explaining the observed flaring radiation in different wavebands. The time scale of ejections into the wind matches the short time scale variability in the light curves of AE Aqr.

Inward migration of low-mass planets and embryos of giant planets can be stopped at the disc-cavity boundaries due to co-orbital corotation torque. We performed the first global three-dimensional (3D) simulations of planet migration at the disc-cavity boundary, and have shown that the boundary is a robust trap for low-mass planets and embryos. A protoplanetary disc may have several such trapping regions at various distances from the star, such as at the edge of the stellar magnetosphere, the inner edge of the dead zone, the dust-sublimation radius and the snow lines. Corotation traps located at different distances from a star, and moving outward during the disc dispersal phase, may possibly explain the observed homogeneous distribution of low-mass planets with distance from their host stars.

Many pulsars propagate through the interstellar medium (ISM) with supersonic velocities, and their pulsar winds interact with the interstellar medium (ISM), forming bow shocks and magnetotails (PWN). We model the propagation of pulsars through the inhomogeneous ISM using non-relativistic axisymmetric magneto-hydrodynamic (MHD) simulations. We take into account the wind from the star, and the azimuthal and poloidal components of the magnetic field, and investigate the PWN at different levels of magnetization (the ratio of magnetic to matter energy-densities) in the wind. We consider the interaction of PWN with small-scale and large-scale imhomogeneities in the ISM at different values of magnetization. We conclude that the inhomogeneities in the ISM can change the shapes of the bow shocks and magnetotails at different values of the magnetization.We compare the results of our simulations with the images of the Guitar Nebula and other PWN that show irregularities in the shapes of their bow shocks and magnetotails. We conclude that these irregularities may be caused by the interaction of PWN with the inhomogeneities in the ISM.

We investigate the properties of magnetized stars in the propeller regime using axisymmetric numerical simulations. We modelled the propeller regime for stars with realistically large magnetospheres (5-7 stellar radii), so that our results could be applied to different types of magnetized stars, including Classical T Tauri stars, cataclysmic variables, and accreting millisecond pulsars. A wide range of propeller strengths has been studied, from very strong propellers (where the magnetosphere rotates much more rapidly than the inner disk) to very weak propellers (where the magnetosphere rotates only slightly faster than the inner disc. In both propellers, matter is accumulated at the inner disc for the majority of the time, while episodes of accretion onto the star and ejection into the wind are relatively brief. The efficiency of the propeller, which characterizes the part of inner disk matter flowing into the wind, strongly depends on the fastness parameter omega_s: propeller efficiency increases with omega_s. The properties of the winds are different in strong and weak propellers. In the strong propellers, matter is accelerated rapidly above the escape velocity and flows at a relatively small opening angle of 40-45 degrees. In the weak propellers, matter may flow faster or slower than the escape velocity and at a large opening angle of 60-70 degrees. A star-disk system loses energy and angular momentum. A part of the rotational energy of the star is ejected to the magnetically-dominated (Poynting flux) jet, which is only present in the strong propellers. The other part of the energy flows from the inner disk into a propeller-driven wind. A star spins down partly due to the flow of angular momentum from the star to the corona or to the Poynting flux jet along the open field lines, and partly due to the flow of angular momentum to the inner disk along the closed field lines.

We present a study of the Rayleigh-Taylor unstable regime of accretion onto rotating magnetized stars in a set of high grid resolution three-dimensional (3D) magnetohydrodynamic (MHD) simulations performed in low-viscosity discs. We find that the boundary between the stable and unstable regimes is determined almost entirely by the fastness parameter omega_s=Omega_star/Omega_K(r_m), where Omega_star is the angular velocity of the star and Omega_K(r_m) is the angular velocity of the Keplerian disc at the disc-magnetosphere boundary r=r_m. We found that accretion is unstable if omega_s < 0.6. Accretion through instabilities is present in stars with different magnetospheric sizes. However, only in stars with relatively small magnetospheres, r_m/R_star < 7, do the unstable tongues produce chaotic hot spots on the stellar surface and irregular light-curves. At even smaller values of the fastness parameter, omega_s < 0.45, multiple irregular tongues merge, forming one or two ordered unstable tongues that rotate with the angular frequency of the inner disc. This transition occurs in stars with even smaller magnetospheres, r_m/R_star < 4.2. Most of our simulations were performed at a small tilt of the dipole magnetosphere, Theta=5 degrees, and a small viscosity parameter alpha=0.02. Test simulations at higher alpha values show that many more cases become unstable, and the light-curves become even more irregular. Test simulations at larger tilts of the dipole Theta show that instability is present, however, accretion in two funnel streams dominates if Theta > 15 degrees. The results of these simulations can be applied to accreting magnetized stars with relatively small magnetospheres: Classical T Tauri stars, accreting millisecond X-ray pulsars, and cataclysmics variables.

Many types of stars have strong magnetic fields that can dynamically influence the flow of circumstellar matter. In stars with accretion disks, the stellar magnetic field can truncate the inner disk and determine the paths that matter can take to flow onto the star. These paths are different in stars with different magnetospheres and periods of rotation. External field lines of the magnetosphere may inflate and produce favorable conditions for outflows from the disk-magnetosphere boundary. Outflows can be particularly strong in the propeller regime, wherein a star rotates more rapidly than the inner disk. Outflows may also form at the disk-magnetosphere boundary of slowly rotating stars, if the magnetosphere is compressed by the accreting matter. In isolated, strongly magnetized stars, the magnetic field can influence formation and/or propagation of stellar wind outflows. Winds from low-mass, solar-type stars may be either thermally or magnetically driven, while winds from massive, luminous O and B type stars are radiatively driven. In all of these cases, the magnetic field influences matter flow from the stars and determines many observational properties. In this chapter we review recent studies of accretion, outflows, and winds of magnetized stars with a focus on three main topics: (1) accretion onto magnetized stars; (2) outflows from the disk-magnetosphere boundary; and (3) winds from isolated massive magnetized stars. We show results obtained from global magnetohydrodynamic simulations and, in a number of cases compare global simulations with observations.

We investigate the migration of low-mass planets ($5 M_{\oplus}$ and $20 M_{\oplus}$) in accretion discs threaded with a magnetic field using 2D MHD code in polar coordinates. We observed that, in the case of a strong azimuthal magnetic field where the plasma parameter is $\beta\sim 1-2$, density waves at the magnetic resonances exert a positive torque on the planet and may slow down or reverse its migration. However, when the magnetic field is weaker (i.e., the plasma parameter $\beta$ is relatively large), then non-axisymmetric density waves excited by the planet lead to growth of the radial component of the field and, subsequently, to development of the magneto-rotational instability, such that the disc becomes turbulent. Migration in a turbulent disc is stochastic, and the migration direction may change as such. To understand migration in a turbulent disc, both the interaction between a planet and individual turbulent cells, as well as the interaction between a planet and ordered density waves, have been investigated.

We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.

Observations of jets from young stellar objects reveal the asymmetric outflows from some sources. A large set of $2.5$D MHD simulations has been carried out for axisymmetric viscous/diffusive disc accretion to rotating magnetized stars for the purpose of assessing the conditions where the outflows or jets are asymmetric relative to the equatorial plane. We consider initial magnetic fields that are symmetric about the equatorial plane and consist of a radially distributed field threading the disc (disc-field) and a stellar dipole field.({\bf 1}). For pure disc-fields the symmetry or asymmetry of the outflows is affected by the midplane plasma $\beta$ of the disc (where $\beta$ is the ratio of the plasma pressure to the magnetic pressure). For the low density discs with small plasma $\beta$ values, outflows are observed to be symmetric to within $10\%$ over timescales of hundreds of inner disc orbits. For the denser higher $\beta$ discs, the coupling of the upper and lower coronal plasmas is broken, and quasi-periodic field motion in the two hemispheres becomes different. This asymmetry leads to asymmetric episodic outflows. ({\bf 2.}) Accreting stars with a stellar dipole field and no disc-field exhibit episodic, two component outflows - a magnetospheric wind and an inner disc wind from somewhat larger radial distances. Both are characterized by similar velocity profiles but the magnetospheric wind has densities $\gtrsim 10$ times that of the disc wind. ({\bf 3}.)Adding a disc-field parallel to the stellar dipole field acts to enhance the magnetospheric winds but suppress the disc wind. ({\bf 4}.) In contrast, adding a disc-field which is anti-parallel to the stellar dipole field in the disc acts to suppress the magnetospheric and disc winds. Our simulations reproduce some key features of observations of asymmetric outflows of T Tauri stars.

An analytic model is developed for the erosion of protoplanetary gas discs by high velocity magnetized stellar winds. The winds are centrifugally driven from the surface of rapidly rotating, strongly magnetized young stars. The presence of the magnetic field in the wind leads to Reynolds numbers sufficiently large to cause a strongly turbulent wind/disk boundary layer which entrains and carries away the disc gas. The model uses the conservation of mass and momentum in the turbulent boundary layer. The time-scale for significant erosion depends on the disc accretion speed and accretion rate and on the wind mass loss rate and the wind velocity. The time-scale is estimated to be 2E6 yr. A significant contribution to the disc erosion can come from frequent powerful coronal mass ejections (CMEs) where the average mass loss rate in CMEs and velocities have values comparable to those for the steady wind.

Counter-rotating discs can arise from the accretion of a counter-rotating gas cloud onto the surface of an existing co-rotating disc or from the counter-rotating gas moving radially inward to the outer edge of an existing disc. At the interface, the two components mix to produce gas or plasma with zero net angular momentum which tends to free-fall towards the disc center. We discuss high-resolution axisymmetric hydrodynamic simulations of a viscous counter-rotating disc for cases where the two components are vertically separated and radially separated. The viscosity is described by an isotropic $\alpha-$viscosity including all terms in the viscous stress tensor. For the vertically separated components a shear layer forms between them. The middle of this layer free-falls to the disk center. The accretion rates are increased by factors $\sim 10^2-10^4$ over that of a conventional disc rotating in one direction with the same viscosity. The vertical width of the shear layer and the accretion rate are strongly dependent on the viscosity and the mass fraction of the counter-rotating gas. In the case of radially separated components where the inner disc co-rotates and the outer disc rotates in the opposite direction, a gap between the two components opens and closes quasi-periodically. The accretion rates are $\gtrsim 25$ times larger than those for a disc rotating in one direction with the same viscosity.

The mass outflows in T Tauri stars (TTS) are thought to be an effective mechanism to remove angular momentum during the pre-main-sequence contraction of a low-mass star. The most powerful winds are observed at the FUor stage of stellar evolution. V1331 Cyg has been considered as a TTS at the pre-FUor stage. We analyse high-resolution spectra of V1331 Cyg collected in 1998-2007 and 20-d series of spectra taken in 2012. For the first time the photospheric spectrum of the star is detected and stellar parameters are derived: spectral type G7-K0 IV, mass 2.8 Msun, radius 5 Rsun, vsini < 6 km/s. The photospheric spectrum is highly veiled, but the amount of veiling is not the same in different spectral lines, being lower in weak transitions and much higher in strong transitions. The Fe II 5018, Mg I 5183, K I 7699 and some other lines of metals are accompanied by a `shell' absorption at radial velocity of about -240 km/s. We show that these absorptions form in the post-shock gas in the jet, i.e. the star is seen though its jet. The P Cyg profiles of H-alpha and H-beta indicate the terminal wind velocity of about 500 km/s, which vary on time-scales from several days to years. A model of the stellar wind is developed to interpret the observations. The model is based on calculation of hydrogen spectral lines using the radiative transfer code TORUS. The observed H-alpha and H-beta line profiles and their variability can be well reproduced with a stellar wind model, where the mass-loss rate and collimation (opening angle) of the wind are variable. The changes of the opening angle may be induced by small variability in magetization of the inner disc wind. The mass-loss rate is found to vary within (6-11)x10^{-8} Msun/yr, with the accretion rate of 2.0x10^{-6} Msun/yr.

A brief review is given of the Rossby wave instability (RWI) in astrophysical discs. In non-self-gravitating discs, around for example a newly forming stars, the instability can be triggered by an axisymmetric bump at some radius $r_0$ in the disc surface mass-density. It gives rise to exponentially growing non-axisymmetric perturbation [$\propto \exp({\rm i}m\phi)$, $m=1,2..$] in the vicinity of $r_0$ consisting of {\it anticyclonic} vortices. These vortices are regions of high pressure and consequently act to trap dust particles which in turn can facilitate planetesimal growth in proto-planetary discs. The Rossby vortices in the discs around stars and black holes may cause the observed quasi-periodic modulations of the disc's thermal emission.

We discuss the results of modelling of young magnetized stars, where matter flow is calculated using the three-dimensional (3D) magneto-hydrodynamic (MHD) "Cubed Sphere" code, and the spectra are calculated using the 3D Monte Carlo radiative transfer code "TORUS". Two examples of modelling are shown: (1) accretion onto stars in stable and unstable regimes, and (2) accretion to a young star V2129 Oph, modelled with realistic parameters.

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

We review recent axisymmetric and three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetospheric accretion, plasma-field interaction and outflows from the disk-magnetosphere boundary.

In this work, we consider the physics of the interaction between a planet and a magnetized gaseous protoplanetary disc. We investigate the migration of a planet in a disc that is threaded with an azimuthal magnetic field. We find that, for a larger magnetic field amplitude, there is an increasingly large positive torque on the planet from the disc, resulting in slowed and even outward migration. Our results indicate that magnetic resonances due to a purely azimuthal, ordered magnetic field can slow or stop the inward migration of Jupiter-mass, Saturn-mass, and $5 M_{\oplus}$ planets.

Recent spectropolarimetric observations suggest that young low-mass stars such as classical T Tauri stars (CTTSs) possess relatively strong (~kG) magnetic field. This supports a scenario in which the final accretion onto the stellar surface proceeds through a magnetosphere, and the winds are formed in magnetohydrodynamics (MHD) processes. We examine recent numerical simulations of magnetospheric accretions via an inclined dipole and a complex magnetic fields. The difference between a stable accretion regime, in which accretion occurs in ordered funnel streams, and an unstable regime, in which gas penetrates through the magnetosphere in several unstable streams due to the magnetic Rayleigh-Taylor instability, will be discussed. We describe how MHD simulation results can be used in separate radiative transfer (RT) models to predict observable quantiles such as line profiles and light curves. The plausibility of the accretion flows and outflows predicted by MHD simulations (via RT models) can be tested against observations. We also address the issue of outflows/winds that arise from the innermost part of CTTSs. First, we discuss the line formations in a simple disk wind and a stellar wind models. We then discuss the formation of the conically shaped magnetically driven outflow that arises from the disk-magnetosphere boundary when the magnetosphere is compressed into an X-type configuration.

We observed persistent high-frequency oscillations of the boundary layer near an accreting, weakly-magnetized star in global 3D MHD simulations. The tilted dipole magnetic field is not strong enough to open a gap between the star and the disk. Instead, it forms a highly-wrapped azimuthal field near the surface of the star which slows down rotation of the disk matter, while a small tilt of the field excites oscillations of the boundary layer with a frequency below the Keplerian frequency. This mechanism may be responsible for the high-frequency oscillations in accreting neutron stars, white dwarfs and classical T Tauri stars.

We investigated the boundary between stable and unstable regimes of accretion and its dependence on different parameters. Simulations were performed using a "cubed sphere" code with high grid resolution (244 grid points in the azimuthal direction), which is twice as high as that used in our earlier studies. We chose a very low viscosity value, with alpha-parameter alpha=0.02. We observed from the simulations that the boundary strongly depends on the ratio between magnetospheric radius r_m (where the magnetic stress in the magnetosphere matches the matter stress in the disk) and corotation radius r_cor (where the Keplerian velocity in the disk is equal to the angular velocity of the star). For a small misalignment angle of the dipole field, Theta=5 degrees, accretion is unstable if r_cor/r_m>1.35, and is stable otherwise. In cases of a larger misalignment angle of the dipole, Theta=20 degrees, instability occurs at slightly larger values, r_cor/r_m>1.41.

We report here results of spectropolarimetric observations of the classical T Tauri star DN Tau carried out (at 2 epochs) with ESPaDOnS at the Canada-France-Hawaii Telescope within the `Magnetic Protostars and Planets' programme. We infer that DN Tau, with a photospheric temperature of 3,950+-50 K, a luminosity of 0.8+-0.2 Lsun and a rotation period of 6.32 d, is a ~2Myr-old fully-convective 0.65+-0.05 Msun star with a radius of 1.9+-0.2 Dsun, viewed at an inclination of 35+-10degr. Clear circularly-polarized Zeeman signatures are detected in both photospheric and accretion-powered emission lines, probing longitudinal fields of up to 1.8 kG (in the He1 D3 accretion proxy). Rotational modulation of Zeeman signatures, detected both in photospheric and accretion lines, is different between our 2 runs, providing further evidence that fields of cTTSs are generated by non-stationary dynamos. Using tomographic imaging, we reconstruct maps of the large-scale field, of the photospheric brightness and of the accretion-powered emission at the surface of DN Tau at both epochs. We find that the magnetic topology is mostly poloidal, and largely axisymmetric, with an octupolar component (of polar strength 0.6-0.8 kG) 1.5-2.0x larger than the dipolar component (of polar strength 0.3-0.5 kG). DN Tau features dominantly poleward accretion at both epochs. The large-scale dipole component of DN Tau is however too weak to disrupt the surrounding accretion disc further than 65-90% of the corotation radius (at which the disc Keplerian period matches the stellar rotation period), suggesting that DN Tau is already spinning up despite being fully convective.

We discuss key observational signatures of Classical T Tauri stars (CTTSs) accreting through Rayleigh-Taylor instability, which occurs at the interface between an accretion disk and a stellar magnetosphere. In this study, the results of global 3-D MHD simulations of accretion flows, in both stable and unstable regimes, are used to predict the variability of hydrogen emission lines and light curves associated with those two distinctive accretion flow patterns. In the stable regime, a redshifted absorption component (RAC) periodically appears in some hydrogen lines, but only during a fraction of a stellar rotation period. In the unstable regime, the RAC is present rather persistently during a whole stellar rotation period, and its strength varies non-periodically. The latter is caused by multiple accreting streams, formed randomly due to the instability, passing across the line of sight to an observer during one stellar rotation. This results in the quasi-stationarity appearance of the RAC because at least one of the accretion stream is almost always in the line of sight to an observer. In the stable regime, two stable hot spots produce a smooth and periodic light curve that shows only one or two peaks per stellar rotation. In the unstable regime, multiple hot spots formed on the surface of the star, produce the stochastic light curve with several peaks per rotation period.

We discuss recent progress in understanding the launching of outflows/jets from the disc-magnetosphere boundary of slowly and rapidly rotating magnetized stars. In most of the discussed models the interior of the disc is assumed to have a turbulent viscosity and magnetic diffusivity (as described by two ``alpha'' parameters), whereas the coronal region outside of the disc is treated using ideal magnetohydrodynamics (MHD). Extensive MHD simulations have established the occurrence of long-lasting outflows in both the cases of slowly and rapidly rotating stars. In the case of {\it slowly rotating stars}, a new type of outflow, {\it a conical wind}, is found and studied in simulations. The conical winds appear in cases where the magnetic flux of the star is bunched up by the inward motion of the accretion disc. Near their region of origin, the winds have the shape of a thin conical shell with a half-opening angle of $ \sim 30^\circ$. At large distances these outflows become magnetically collimated by their toroidal magnetic field and form matter dominated jets. That is, the jets are current carrying. The predominant driving force for the conical winds is the magnetic force proportional to the negative gradient of the square of the toroidal magnetic field and not the centrifugal force. In the case of {\it rapidly rotating stars} in the ``propeller regime,'' a two-component outflow is observed. The first component is similar to the matter dominated conical winds. A large fraction of the disc matter may be ejected into the winds in this regime. The second component is a high-velocity, low-density magnetically dominated {\it axial jet} where matter flows along the opened polar field lines of the star.

We present the results of axisymmetric simulations of MRI-driven accretion onto a rapidly rotating, magnetized star accreting in the propeller regime. The stellar magnetosphere corotates with the star, forming a centrifugal barrier at the disc-magnetosphere boundary which inhibits matter accretion onto the star. Instead, the disc matter accumulates at the disc-magnetosphere interface and slowly diffuses into the inner magnetosphere where it picks up angular momentum and is quickly ejected from the system as an outflow. Due to the interaction of the matter with the magnetosphere, this wind is discontinuous and is launched as discrete plasmoids. If the ejection rate is lower than the disc accretion rate, the matter accumulates at the disc-magnetosphere boundary faster than it can be ejected. In this case, accretion onto the star proceeds through the episodic accretion instability in which episodes of matter accumulation are followed by simultaneous accretion and ejection. During the accretion phase of this instability in which matter flows onto the star in funnel streams, we observe a corresponding rise in the outflow rate. Both the accretion and ejection processes observed in our simulations are highly non-stationary. The stars undergo strong spin-down due to the coupling of the stellar field with the disc and corona and we measure the spin-down timescales of around 1 Myr for a typical CTTS in the propeller regime.

We have carried out and analyzed a set of axisymmetric MHD simulations of the evolution of a turbulent/diffusive accretion disc around an initially unmagnetized star. The disc is initially threaded by a weak magnetic field where the magnetic pressure is significantly less than the kinetic pressure in the disc. The viscosity and magnetic diffusivity are modelled by two "alpha" parameters, while the coronal region above the disc is treated using ideal MHD. The initial magnetic field is taken to consist of three poloidal field loops threading the disc. The motivation for this study is to understand the advection of disc matter and magnetic field by the turbulent/diffusive disc. At early times the innermost field loop twists and its field lines become open. The twisting of the opened field lines leads to the formation of both an inner collimated, magnetically-dominated jet, and at larger distances from the axis a matter dominated uncollimated wind. For later times, the strength of the magnetic field decreases owing to field reconnection and annihilation in the disc. For the early times, we have derived from the simulations both the matter accretion speed in the disc and the accretion speed of the magnetic field. We show that the derived matter accretion speed agrees approximately with the predictions of a model where the accretion speed is the sum of two terms, one due to the disc's viscosity (which gives a radial outflow of angular momentum in the disc), and a second due to the twisted magnetic field at the disc's surface (which gives a vertical outflow of angular momentum). For early times we find that the magnetic contribution is roughly twice the viscous contribution for the case where the alpha parameters are both equal to 0.1. At later times the magnetic contribution to the matter speed becomes small compared to the viscous contribution.

We report results of the first global three-dimensional magnetohydrodynamic simulations of the waves excited in an accretion disc by a rotating star with a dipole magnetic field misaligned from the star's rotation axis (which is aligned with the disc axis). The main results are the following: (1) If the magnetosphere of the star corotates approximately with the inner disc, then we observe a strong one-armed bending wave (a warp). This warp corotates with the star and has a maximum amplitude between corotation radius and the radius of the vertical resonance. The disc's center of mass can deviate from the equatorial plane up to the distance of z_w\approx 0.1 r. However, the effective height of the warp can be larger, h_w \approx 0.3 r due to the finite thickness of the disc. Stars with a range of misalignment angles excite warps. However, the amplitude of the warps is larger for misalignment angles between 15 and 60 degrees. (2) If the magnetosphere rotates slower, than the inner disc, then a bending wave is excited at the disc-magnetosphere boundary, but does not form a large-scale warp. Instead, high-frequency oscillations become strong at the inner region of the disc. These are (a) trapped density waves which form inside the radius where the disc angular velocity has a maximum, and (b) inner bending waves which appear in the case of accretion through magnetic Raleigh-Taylor instability. These two types of waves are connected with the inner disc and their frequencies will vary with accretion rate. Bending oscillations at lower frequencies are also excited including global oscillations of the disc. In cases where the simulation region is small, slowly-precessing warp forms. Simulations are applicable to young stars, cataclysmic variables, and accreting millisecond pulsars.

We present an analytical formula for the position and shape of the spots on the surface of accreting magnetized stars in cases where a star has a dipole magnetic field tilted at a small misalignment angle Theta < 30 degrees about the rotational axis, and the magnetosphere is 2.5-5 times the radius of the star. We observed that the azimuthal position of the spots varies significantly when the position of the inner disc varies. In contrast, the polar position of the spots varies only slightly because of the compression of the magnetosphere. The azimuthal width of the spots strongly varies with Theta: spots have the shape of an arc at larger misalignment angles, and resemble a ring at very small misalignment angles. The polar width of the spots varies only slightly with changes in parameters. The motion of the spots in the azimuthal direction can provide phase-shifts in accreting millisecond pulsars, and the "drift" of the period in Classical T Tauri stars. The position and shape of the spots are determined by three parameters: misalignment angle Theta; normalized corotation radius, r_c/R_* and normalized magnetospheric radius, r_m/R_*, where R_* is the stellar radius. We also use our data to check the formula for the Alfv\'en radius, where the main dependencies are r_m \sim (\mu^2/\dot M)^{2/7}, where \mu is the magnetic moment of the star, and \dot M is the accretion rate. We found that the dependence is more gradual, r_m \sim (\mu^2/\dot M)^{1/5}, which can be explained by the compression of the magnetosphere by the disc matter and by the non-dipole shape of the magnetic field lines of the external magnetosphere.

Classical T Tauri stars (CTTSs) are variable in different time-scales. One type of variability is possibly connected with the accretion of matter through the Rayleigh-Taylor instability that occurs at the interface between an accretion disc and a stellar magnetosphere. In this regime, matter accretes in several temporarily formed accretion streams or `tongues' which appear in random locations, and produce stochastic photometric and line variability. We use the results of global three-dimensional magnetohydrodynamic simulations of matter flows in both stable and unstable accretion regimes to calculate time-dependent hydrogen line profiles and study their variability behaviours. In the stable regime, some hydrogen lines (e.g. H-beta, H-gamma, H-delta, Pa-beta and Br-gamma) show a redshifted absorption component only during a fraction of a stellar rotation period, and its occurrence is periodic. However, in the unstable regime, the redshifted absorption component is present rather persistently during a whole stellar rotation cycle, and its strength varies non-periodically. In the stable regime, an ordered accretion funnel stream passes across the line of sight to an observer only once per stellar rotation period while in the unstable regime, several accreting streams/tongues, which are formed randomly, pass across the line of sight to an observer. The latter results in the quasi-stationarity appearance of the redshifted absorption despite the strongly unstable nature of the accretion. In the unstable regime, multiple hot spots form on the surface of the star, producing the stochastic light curve with several peaks per rotation period. This study suggests a CTTS that exhibits a stochastic light curve and a stochastic line variability, with a rather persistent redshifted absorption component, may be accreting in the unstable accretion regime.

We present the emission line profile models of hydrogen and helium based on the results from axisymmetric magnetohydrodynamics (MHD) simulations of the wind formed near the disk-magnetosphere boundary of classical T Tauri stars (CTTSs). We extend the previous outflow models of `the conical-shell wind' by Romanova et al. to include a well defined magnetospheric accretion funnel flow which is essential for modelling the optical and near-infrared hydrogen and helium lines of CTTSs. The MHD model with an intermediate mass-accretion rate shows outflows in conical-shell shape with a half opening angle about 35 degrees. The flow properties such as the maximum outflow speed in the conical-shell wind, maximum inflow speed in the accretion funnel, mass-accretion and mass-loss rates are comparable to those found in a typical CTTS. The density, velocity and modified temperature from the MHD simulations are used in a separate radiative transfer model to predict the line profiles and test the consistency of the MHD models with observations. The line profiles are computed with various combinations of X-ray luminosities, temperatures of X-ray emitting plasma, and inclination angles. A rich diversity of line profile morphology is found, and many of the model profiles are very similar to those found in observations. We find that the conical-shell wind may contribute to the emission in some hydrogen lines (e.g. H-alpha, H-beta, Pa-beta and Pa-gamma) significantly when the temperature in the wind is relatively high (e.g. \sim 10^{4} K); however, the wind contribution decreases rapidly when a lower wind temperature is adopted. The model well reproduces a relatively narrow and low-velocity blueshifted absorption component in He I (10830), which are often seen in observations.

We report here results of spectropolarimetric observations of the classical T Tauri star (cTTS) GQ Lup carried out with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) in the framework of the "Magnetic Protostars and Planets" (MaPP) programme, and obtained at 2 different epochs (2009 July & 2011 June). From these observations, we first infer that GQ Lup has a photospheric temperature of 4,300+-50\^A K and a rotation period of 8.4+-0.3 d; it implies that it is a 1.05+-0.07 Msun star viewed at an inclination of ~30deg, with an age of 2-5 Myr, a radius of 1.7+-0.2 Rsun, and has just started to develop a radiative core. Large Zeeman signatures are clearly detected at all times, both in photospheric lines & in accretion-powered emission lines, probing longitudinal fields of up to 6 kG and hence making GQ Lup the cTTS with the strongest large-scale fields known as of today. Rotational modulation of Zeeman signatures is clearly different between our 2 runs, demonstrating that large-scale fields of cTTSs are evolving with time and are likely produced by non-stationary dynamo processes. Using tomographic imaging, we reconstruct maps of the large-scale field, of the photospheric brightness & of the accretion-powered emission of GQ Lup. We find that the magnetic topology is mostly poloidal & axisymmetric; moreover, the octupolar component of the large-scale field (of strength 2.4 & 1.6 kG in 2009 & 2011) dominates the dipolar component (of strength ~1 kG) by a factor of ~2, consistent with the fact that GQ Lup is no longer fully-convective. GQ Lup also features dominantly poleward magnetospheric accretion at both epochs. The large-scale dipole of GQ Lup is however not strong enough to disrupt the surrounding accretion disc further than about half-way to the corotation radius, suggesting that GQ Lup should rapidly spin up like other similar partly-convective cTTSs (abridged).

We analyze the photometric and spectroscopic variability of the classical T Tauri star V2129 Oph over several rotational cycles to test the dynamical predictions of magnetospheric accretion models. The photometric variability and the radial velocity variations in the photospheric lines can be explained by rotational modulation due to cold spots, while the radial velocity variations of the He I (5876 \AA) line and the veiling variability are due to hot spot rotational modulation. The hot and cold spots are located at high latitudes and about the same phase, but the hot spot is expected to sit at the chromospheric level, while the cold spot is at the photospheric level. Using the dipole+octupole magnetic-field configuration previously proposed in the literature for the system, we compute 3D MHD magnetospheric simulations of the star-disk system. We use the simulation's density, velocity and scaled temperature structures as input to a radiative transfer code, from which we calculate theoretical line profiles at all rotational phases. The theoretical profiles tend to be narrower than the observed ones, but the qualitative behavior and the observed rotational modulation of the H\alpha and H\beta emission lines are well reproduced by the theoretical profiles. The spectroscopic and photometric variability observed in V2129 Oph support the general predictions of complex magnetospheric accretion models with non-axisymmetric, multipolar fields.

We show results of global 3D MHD simulations of accretion onto stars with superposition of the dipole and octupole fields, where we vary the ratio between components. Simulations show that if octupolar field strongly dominates at the disc-magnetosphere boundary, then matter flows into the ring-like octupolar poles, forming ring-shape spots at the surface of the star above and below equator. The light-curves are complex and may have two peaks per period. In case where the dipole field dominates, matter accretes in two ordered funnel streams towards poles, however the polar spots are meridionally-elongated due to the action of the octupolar component. In the case when the fields are of similar strengths, both, polar and belt-like spots are present. In many cases the light-curves show the evidence of complex fields, excluding the cases of small inclinations angles, where sinusoidal light-curve %is observed and `hides' the information about the field complexity. We also propose new mechanisms of phase shift in stars with complex magnetic fields. We suggest that the phase shifts can be connected with: (1) temporal variation of the star's intrinsic magnetic field and subsequent redistribution of main magnetic poles; (2) variation of the accretion rate, which causes the disc to interact with the magnetic fields associated with different magnetic moments. We use our model to demonstrate these phase shift mechanisms, and we discuss possible applications of these mechanisms to accreting millisecond pulsars and young stars.

We discuss results of global 3D MHD simulations of accretion on to a rotating magnetized star with a tilted dipole magnetic field, where the accretion is driven by the magneto-rotational instability (MRI). The simulations show that MRI-driven turbulence develops in the disc, and angular momentum is transported outwards due primarily to the magnetic stress. The turbulent flow is strongly inhomogeneous and the densest matter is in azimuthally-stretched turbulent cells. We investigate two regimes of accretion: a magnetospheric regime and a boundary layer (BL) regime. In the magnetospheric regime, the accretion disc is truncated by the star's magnetic field within a few stellar radii from the star, and matter flows to the star in funnel streams. The funnel streams flowing towards the south and north magnetic poles but are not equal due to the inhomogeneity of the flow. In the BL regime, matter accretes to the surface of the star through the boundary layer. The magnetic field in the inner disc is strongly amplified by the shear of the accretion flow, and the matter and magnetic stresses become comparable. Accreting matter forms a belt-shaped region on the surface of the star. The belt has inhomogeneous density distribution which varies in time due to variable accretion rate. Results of simulations can be applied to classical T Tauri stars, accreting brown dwarfs, millisecond pulsars, dwarf novae cataclysmic variables, and other stars with magnetospheres smaller than several stellar radii.

Axisymmetric MHD simulations are used to investigate the Bondi-Hoyle accretion onto an isolated magnetized neutron star moving supersonically (with Mach number of 3) through the interstellar medium. The star is assumed to have a dipole magnetic field aligned with its motion and a magnetospheric radius R_m less then the accretion radius R_BH, so that the gravitational focusing is important. We find that the accretion rate to a magnetized star is smaller than that to a non-magnetized star for the parameters considered. Close to the star the accreting matter falls to the star's surface along the magnetic poles with a larger mass flow to the leeward pole of the star. In the case of a relatively large stellar magnetic field, the star's magnetic field is stretched in the direction of the matter flow outside of R_m (towards the windward side of the star). For weaker magnetic fields we observed oscillations of the closed magnetosphere frontward and backward. These are accompanied by strong oscillations of the mass accretion rate which varies by factors ~ 3. Old slowly rotating neutron stars with no radio emission may be visible in the X-ray band due to accretion of interstellar matter. In general, the star's velocity, magnetic moment, and angular velocity vectors may all be in different directions so that the accretion luminosity will be modulated at the star's rotation rate.

We use axisymmetric magnetohydrodynamic (MHD) simulations to investigate the launching and collimation of jets emerging from the disk-magnetosphere boundary of accreting magnetized stars. Our analysis shows that the matter flows into the jet from the inner edge of the accretion disk. It is magnetically accelerated along field lines extending up from the disk and simultaneously collimated by the magnetic pinch force. In the reference run which we use for analysis, the matter in the jet crosses the Alfv\'en surface a few R_* above the disk and the fast magnetosonic surface ~13 R_* above the disk. At larger distances, the magnetic pressure is a few times smaller than the total matter pressure but the magnetic force continues to accelerate and collimate the jet. In steady state, we observe a matter ejection-to-accretion ratio of ~0.2. Across different simulation runs, we measure a range of half-opening angles between {\Theta} \approx 4{\deg} and 20{\deg} at the top of the simulation region, depending on the degree of magnetization in the outflow. We consider the case of stars undergoing epochs of high accretion (such as EXors, FUORs, and CTTSs) where the stellar magnetosphere is strongly compressed by the incoming accretion disk. For a typical EXor (mass 0.8 M_sun, radius 2 R_sun) accreting at ~10^-5 M_sun/yr, we measure poloidal velocities in the jet ranging from 30 km/s on the outer edge of the jet to more than 260 km/s on the inner edge. In general, the models can be applied to a variety of magnetized stars--white dwarfs, neutron stars, and brown dwarfs--which exhibit periods of high accretion.

We report here results of spectropolarimetric observations of the ~8Myr classical TTauri star (cTTS) TWHya carried out with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) in the framework of the `Magnetic Protostars and Planets' (MaPP) programme, and obtained at 2 different epochs (2008 March and 2010 March). Obvious Zeeman signatures are detected at all times, both in photospheric lines and in accretion-powered emission lines. Significant intrinsic variability and moderate rotational modulation is observed in both photospheric and accretion proxies. Using tomographic imaging, we reconstruct maps of the large-scale field, of the photospheric brightness and of the accretion-powered emission at the surface of TWHya at both epochs. We find that the magnetic topology is mostly poloidal and axisymmetric with respect to the rotation axis of the star, and that the octupolar component of the large-scale field (2.5-2.8kG at the pole) largely dominates the dipolar component. This large-scale field topology is characteristic of partly-convective stars, supporting the conclusion (from evolutionary models) that TWHya already hosts a radiative core. We also show that TWHya features a high-latitude photospheric cool spot overlapping with the main magnetic pole (and producing the observed radial velocity fluctuations); this is also where accretion concentrates most of the time, although accretion at lower latitudes is found to occur episodically. We propose that the relatively rapid rotation of TWHya (with respect to AATau-like cTTSs) directly reflects the weakness of the large-scale dipole, no longer capable of magnetically disrupting the accretion disc up to the corotation radius (at which the Keplerian period equals the stellar rotation period). We therefore conclude that TWHya is in a phase of rapid spin-up as its large-scale dipole field progressively vanishes.

FU Orionis (FUOR) outbursts are major optical brightening episodes in low-mass protostars that correspond to rapid mass-accretion events in the innermost region of a protostellar disc. The outbursts are accompanied by strong outflows, with the inferred mass outflow rates reaching ~10% of the mass inflow rates. Shu et al. proposed that the outflows represent accreted disc material that is driven centrifugally from the spun-up surface layers of the protostar by the stellar magnetic field. This model was critiqued by Calvet et al., who argued that it cannot reproduce the photospheric absorption-line shifts observed in the prototype object FU Ori. Calvet et al. proposed that the wind is launched, instead, from the surface of the disc on scales of a few stellar radii by a non-stellar magnetic field. In this paper we present results from numerical simulations of disc accretion on to a slowly rotating star with an aligned magnetic dipole moment that gives rise to a kilogauss-strength surface field. We demonstrate that, for parameters appropriate to FU Ori, such a system can develop a strong, collimated disc outflow of the type previously identified by Romanova et al. in simulations of protostars with low and moderate accretion rates. At the high accretion rate that characterizes FUOR outbursts, the radius at which the disc is truncated by the stellar magnetic field moves much closer to the stellar surface, but the basic properties of the outflow, which is launched from the vicinity of the truncation radius along opened-up stellar magnetic field lines, remain the same, and are distinct from those of the mechanism proposed by Shu et al. We show that the simulated outflow can in principle account for the main observed characteristics of FUOR winds, including the photospheric line shifts measured in FU Ori.

We present multidimensional non-LTE radiative transfer models of hydrogen and helium line profiles formed in the accretion flows and the outflows near the star-disk interaction regions of classical T Tauri stars (CTTSs). The statistical equilibrium calculations, performed under the assumption of the Sobolev approximation using the radiative transfer code TORUS, has been improved to include He I and He II energy levels. This allows us to probe the physical conditions of the inner wind of CTTSs by simultaneously modelling the robust wind diagnostic line He I (10830) and the accretion diagnostic lines such as Pa-beta, Br-gamma and He I (5876). The code has been tested in 1 and 2-D problems, and we have shown that the results are in agreement with established codes. We apply the model to the complex flow geometries of CTTSs. Example model profiles are computed using the combinations of (1) magnetospheric accretion and disc wind, and (2) magnetospheric accretion and the stellar wind. In both cases, the model produces line profiles which are qualitatively similar to those found in observations. Our models are consistent with the scenario in which the narrow blueshifted absorption component of He I (10830) seen in observations is caused by a disc wind, and the wider blueshifted absorption component (the P-Cygni profile) is caused by a bipolar stellar wind. However, we do not have a strong constraint on the relative importance of the wind and the magnetosphere for the `emission' component. Our preliminary calculations suggest that the temperatures of the disc wind, stellar wind and the magnetosphere cannot be much higher than ~10,000 K, on the basis of the strengths of hydrogen lines. With these low temperatures, we find that the photoionzation by high energy photons (e.g. X-rays) is necessary to produce He I (10830) in emission and to produce the blueshifted absorption components.

We present the first results of a global axisymmetric simulation of accretion onto rotating magnetized stars from a turbulent, MRI-driven disk. The angular momentum is transported outward by the magnetic stress of the turbulent flow with a rate corresponding to a Shakura-Sunyaev viscosity parameter alpha\approx 0.01-0.04. The result of the disk-magnetosphere interaction depends on the orientation of the poloidal field in the disk relative to that of the star at the disk-magnetosphere boundary. If fields have the same polarity, then the magnetic flux is accumulated at the boundary and blocks the accretion which leads to the accumulation of matter at the boundary. Subsequently, this matter accretes to the star in outburst before accumulating again. Hence, the cycling, `bursty' accretion is observed. If the disc and stellar fields have opposite polarity, then the field reconnection enhances the penetration of the disk matter towards the deeper field lines of the magnetosphere. However, the magnetic stress at the boundary is lower due to the field reconnection. This decreases the accretion rate and leads to smoother accretion at a lower rate. Test simulations show that in the case of higher accretion rate corresponding to alpha=0.05-0.1, accretion is bursty in cases of both polarities. On the other hand, at much lower accretion rates corresponding to alpha < 0.01, accretion is not bursty in any of these cases. We conclude that the episodic, bursty accretion is expected during periods of higher accretion rates in the disc, and in some cases it may alternate between bursty and smooth accretion, if the disk brings the poloidal field of alternating polarity. We find that a rotating, magnetically-dominated corona forms above and below the disk, and that it slowly expands outward, driven by the magnetic force.

We report here the first results of a multi-wavelength campaign focussing on magnetospheric accretion processes of the classical TTauri star (cTTS) V2129Oph. In this paper, we present spectropolarimetric observations collected in 2009 July with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT). Circularly polarised Zeeman signatures are clearly detected, both in photospheric absorption and accretion-powered emission lines, from time-series of which we reconstruct new maps of the magnetic field, photospheric brightness and accretion-powered emission at the surface of V2129Oph using our newest tomographic imaging tool - to be compared with those derived from our old 2005 June data set, reanalyzed in the exact same way. We find that in 2009 July, V2129Oph hosts octupolar & dipolar field components of about 2.1 & 0.9kG respectively, both tilted by about 20deg with respect to the rotation axis; we conclude that the large-scale magnetic topology changed significantly since 2005 June (when the octupole and dipole components were about 1.5 and 3 times weaker respectively), demonstrating that the field of V2129Oph is generated by a non-stationary dynamo. We also show that V2129Oph features a dark photospheric spot and a localised area of accretion-powered emission, both close to the main surface magnetic region (hosting fields of up to about 4kG in 2009 July). We finally obtain that the surface shear of V2129Oph is about half as strong as solar. From the fluxes of accretion-powered emission lines, we estimate that the observed average logarithmic accretion rate (in Msun/yr) at the surface of V2129Oph is -9.2+-0.3 at both epochs, peaking at -9.0 at magnetic maximum. It implies in particular that the radius at which the magnetic field of V2129Oph truncates the inner accretion disc is 0.93x and 0.50x the corotation radius in 2009 July and 2005 June respectively.

We present recent results of 3D magnetohydrodynamic simulations of neutron stars with small misalignment angles, as regards the features in lightcurves produced by regular movements of the hot spots during accretion onto the star. In particular, we show that the variation of position of the hot spot created by the infalling matter, as observed in 3D simulations, can produce high frequency Quasi Periodic Oscillations with frequencies associated with the inner zone of the disk. Previously reported simulations showed that the usual assumption of a fixed hot spot near the polar region is valid only for misalignment angles relatively large. Otherwise, two phenomena challenge the assumption: one is the presence of Rayleigh-Taylor instabilities at the disk-magnetospheric boundary, which produce tongues of accreting matter that can reach the star almost anywhere between the equator and the polar region; the other one is the motion of the hot spot around the magnetic pole during stable accretion. In this paper we start by showing that both phenomena are capable of producing short-term oscillations in the lightcurves. We then use Monte Carlo techniques to produce model lightcurves based on the features of the movements observed, and we show that the main features of kHz QPOs can be reproduced. Finally, we show the behavior of the frequencies of the moving spots as the mass accretion rate changes, and propose a mechanism for the production of double QPO peaks.

3D MHD simulation of accretion onto neutron stars have shown in the last few years that the footprint (hotspot) of the accretion flow changes with time. Two different kinds of accretion, namely the funnel flow and the equatorial accretion produced by instabilities at the inner disk, produce different kinds of motion of the hotspot. The funnel flow produces hotspots that move around the magnetic pole, while instabilities produce other hotspots that appear randomly and move along the equator or slightly above. The angular velocities of the two hotspots are different, the equatorial one being higher and both close to the Keplerian velocity in the inner region. Modeling of the lightcurves of these hotspots with Monte Carlo simulations show that the signatures produced in power specra by them, if observed, are QPOs plus low frequency components. Their frequencies, general behavior and features describe correctly most of the properties of kHz QPOs, if we assume the funnel flow hotspots as the origin of the lower kHz QPO and instabilities as the origin of the upper kHz QPO.

The magnetic field of the classical T Tauri star V2129 Oph can be modeled approximately by superposing slightly tilted dipole and octupole moments, with polar magnetic field strengths of 0.35kG and 1.2kG respectively (Donati et al. 2007). Here we construct a numerical model of V2129 Oph incorporating this result and simulate accretion onto the star. Simulations show that the disk is truncated by the dipole component and matter flows towards the star in two funnel streams. Closer to the star, the flow is redirected by the octupolar component, with some of the matter flowing towards the high-latitude poles, and the rest into the octupolar belts. The shape and position of the spots differ from those in a pure dipole case, where crescent-shaped spots are observed at the intermediate latitudes. Simulations show that if the disk is truncated at the distance of 6.2 R_* which is comparable with the co-rotation radius, 6.8 R_*, then the high-latitude polar spots dominate, but the accretion rate obtained from the simulations is about an order of magnitude lower than the observed one. The accretion rate matches the observed one if the disk is disrupted much closer to the star, at 3.4 R_*. However, the octupolar belt spots strongly dominate. Better match has been obtained in experiments with a dipole field twice as strong. The torque on the star from the disk-magnetosphere interaction is small, and the time-scale of spin evolution, 2 x10^7-10^9 years is longer than the 2x10^6 years age of V2129 Oph. The external magnetic flux of the star is strongly influenced by the disk: the field lines connecting the disk and the star inflate and form magnetic towers above and below the disk. The potential (vacuum) approximation is still valid inside the Alfv\'en (magnetospheric) surface where the magnetic stress dominates over the matter stress.

The magnetic field of the classical T Tauri star BP Tau can be approximated as a superposition of dipole and octupole moments with respective strengths of the polar magnetic fields of 1.2 kG and 1.6 kG (Donati et al. 2008). We adopt the measured properties of BP Tau and model the disc accretion onto the star. We observed in simulations that the disc is disrupted by the dipole component and matter flows towards the star in two funnel streams which form two accretion spots below the dipole magnetic poles. The octupolar component becomes dynamically important very close to the star and it redirects the matter flow to higher latitudes. The spots are meridionally elongated and are located at higher latitudes, compared with the pure dipole case, where crescent-shaped, latitudinally elongated spots form at lower latitudes. The position and shape of the spots are in good agreement with observations. The disk-magnetosphere interaction leads to the inflation of the field lines and to the formation of magnetic towers above and below the disk. The magnetic field of BP Tau is close to the potential only near the star, inside the magnetospheric surface, where magnetic stress dominates over the matter stress. A series of simulation runs were performed for different accretion rates. They show that an accretion rate is lower than obtained in many observations, unless the disc is truncated close to the star. The torque acting on the star is about an order of magnitude lower than that which is required for the rotational equilibrium. We suggest that a star could lose most of its angular momentum at earlier stages of its evolution.

We investigate the generation of intrinsically asymmetric or {\it one-sided} outflows or jets from disk accretion onto rotating stars with complex magnetic fields using axisymmetric (2.5D) magnetohydrodynamic simulations. The intrinsic magnetic field of the star is assumed to consist of a superposition of an aligned dipole and an aligned quadrupole in different proportions. The star is assumed to be rapidly rotating in the sense that the star's magnetosphere is in the propeller regime where strong outflows occur. Our simulations show that for conditions where there is a significant quadrupole component in addition to the dipole component, then a dominantly {\it one-sided} conical wind tends to form on the side of the equatorial plane with the larger value of the intrinsic axial magnetic field at a given distance. For cases where the quadrupole component is absent or very small, we find that dominantly one-sided outflows also form, but the direction of the flow "flip-flops" between upward and downward on a time-scale of $\sim 30$ days for a protostar. The average outflow will thus be symmetrical. In the case of a pure quadrupole field we find symmetric outflows in the upward and downward directions.

3D Magnetohydrodynamic simulations show that when matter accretes onto neutron stars, in particular if the misalignment angle is small, it does not constantly fall at a fixed spot. Instead, the location at which matter reaches the star moves. These moving hot spots can be produced both during stable accretion, where matter falls near the magnetic poles of the star, and unstable accretion, characterized by the presence of several tongues of matter which fall on the star near the equator, due to Rayleigh-Taylor instabilities. Precise modeling with Monte Carlo simulations shows that those movements could be observed as high frequency Quasi Periodic Oscillations. We performed a number of new simulation runs with a much wider set of parameters, focusing on neutron stars with a small misalignment angle. In most cases we observe oscillations whose frequency is correlated with the mass accretion rate $\dot{M}$. Moreover, in some cases double QPOs appear, each of them showing the same correlation with $\dot{M}$.

This work investigates the short wavelength stability of the magnetopause between a rapidly-rotating, supersonic, dense accretion disc and a slowly-rotating low-density magnetosphere of a magnetized star. The magnetopause is a strong shear layer with rapid changes in the azimuthal velocity, the density, and the magnetic field over a short radial distance and thus the Kelvin-Helmholtz (KH) instability may be important. The plasma dynamics is treated using non-relativistic, compressible (isentropic) magnetohydrodynamics. It is necessary to include the displacement current in order that plasma wave velocities remain less than the speed of light. We focus mainly on the case of a star with an aligned dipole magnetic field so that the magnetic field is axial in the disc midplane and perpendicular to the disc flow velocity. However, we also give results for cases where the magnetic field is at an arbitrary angle to the flow velocity. For the aligned dipole case the magnetopause is most unstable for KH waves propagating in the azimuthal direction perpendicular to the magnetic field which tends to stabilize waves propagating parallel to it. The wave phase velocity is that of the disc matter. A quasi-linear theory of the saturation of the instability leads to a wavenumber ($k$) power spectrum $\propto k^{-1}$ of the density and temperature fluctuations of the magnetopause, and it gives the mass accretion and angular momentum inflow rates across the magnetopause. For self-consistent conditions this mass accretion rate will be equal to the disc accretion rate at large distances from the magnetopause.

We investigate the launching of outflows from the disk-magnetosphere boundary of slowly and rapidly rotating magnetized stars using axisymmetric and exploratory 3D magnetohydrodynamic (MHD) simulations. We find long-lasting outflows in both cases. (1) In the case of slowly rotating stars, a new type of outflow, a conical wind, is found and studied in simulations. The conical winds appear in cases where the magnetic flux of the star is bunched up by the disk into an X-type configuration. The winds have the shape of a thin conical shell with a half-opening angle 30-40 degrees. The conical winds may be responsible for episodic as well as long-lasting outflows in different types of stars. (2) In the case of rapidly rotating stars (the "propeller regime"), a two-component outflow is observed. One component is similar to the conical winds. A significant fraction of the disk matter may be ejected into the winds. A second component is a high-velocity, low-density magnetically dominated axial jet where matter flows along the opened polar field lines of the star. The jet has a mass flux about 10% that of the conical wind, but its energy flux (dominantly magnetic) can be larger than the energy flux of the conical wind. The jet's angular momentum flux (also dominantly magnetic) causes the star to spin-down rapidly. Propeller-driven outflows may be responsible for the jets in protostars and for their rapid spin-down. The jet is collimated by the magnetic force while the conical winds are only weakly collimated in the simulation region.

We investigate the photometric variability of magnetized stars, particularly neutron stars, accreting through a magnetic Rayleigh-Taylor-type instability at the disk-magnetosphere interface, and compare it with the variability during stable accretion, with the goal of looking for possible quasi-periodic oscillations. The lightcurves during stable accretion show periodicity at the star's frequency and sometimes twice that, due to the presence of two funnel streams that produce antipodal hotspots near the magnetic poles. On the other hand, lightcurves during unstable accretion through tongues penetrating the magnetosphere are more chaotic due to the stochastic behaviour of the tongues, and produce noisier power spectra. However, the power spectra do show some signs of quasi-periodic variability. Most importantly, the rotation frequency of the tongues and the resulting hotspots is close to the inner-disk orbital frequency, except in the most strongly unstable cases. There is therefore a high probability of observing QPOs at that frequency in longer simulations. In addition, the lightcurves in the unstable regime show periodicity at the star's rotation frequency in many of the cases investigated here, again except in the most strongly unstable cases which lack funnel flows and the resulting antipodal hotspots. The noisier power spectra result in the fractional rms amplitudes of the Fourier peaks being smaller. We also study in detail the effect of the misalignment angle between the rotation and magnetic axes of the star on the variability, and find that at misalignment angles $\gtrsim 25^\circ$, the star's period always appears in the lightcurves.

We report on the numerical discovery of quasi-periodic oscillations (QPOs) associated with accretion through a non-axisymmetric magnetic boundary layer in the unstable regime, when two ordered equatorial streams form and rotate synchronously at approximately the angular velocity of the inner disk The streams hit the star's surface producing hot spots. Rotation of the spots leads to high-frequency QPOs. We performed a number of simulation runs for different magnetospheric sizes from small to tiny, and observed a definite correlation between the inner disk radius and the QPO frequency: the frequency is higher when the magnetosphere is smaller. In the stable regime a small magnetosphere forms and accretion through the usual funnel streams is observed, and the frequency of the star is expected to dominate the lightcurve. We performed exploratory investigations of the case in which the magnetosphere becomes negligibly small and the disk interacts with the star through an equatorial belt. We also performed investigation of somewhat larger magnetospheres where one or two ordered tongues may dominate over other chaotic tongues. In application to millisecond pulsars we obtain QPO frequencies in the range of 350 Hz to 990 Hz for one spot. The frequency associated with rotation of one spot may dominate if spots are not identical or antipodal. If the spots are similar and antipodal then the frequencies are twice as high. We show that variation of the accretion rate leads to drift of the QPO peak.

We discuss three modes of oscillation of accretion disks around rotating magnetized neutron stars which may explain the separations of the kilo-Hertz quasi periodic oscillations (QPO) seen in low mass X-ray binaries. The existence of these compressible, non-barotropic magnetohydrodynamic (MHD) modes requires that there be a maximum in the angular velocity $\Omega_\phi(r)$ of the accreting material larger than the angular velocity of the star $\Omega_*$, and that the fluid is in approximately circular motion near this maximum rather than moving rapidly towards the star or out of the disk plane into funnel flows. Our MHD simulations show this type of flow and $\Omega_\phi(r)$ profile. The first mode is a Rossby wave instability (RWI) mode which is radially trapped in the vicinity of the maximum of a key function $g(r){\cal F}(r)$ at $r_{R}$. The real part of the angular frequency of the mode is $\omega_r=m\Omega_\phi(r_{R})$, where $m=1,2...$ is the azimuthal mode number. The second mode, is a mode driven by the rotating, non-axisymmetric component of the star's magnetic field. It has an angular frequency equal to the star's angular rotation rate $\Omega_*$. This mode is strongly excited near the radius of the Lindblad resonance which is slightly outside of $r_R$. The third mode arises naturally from the interaction of flow perturbation with the rotating non-axisymmetric component of the star's magnetic field. It has an angular frequency $\Omega_*/2$. We suggest that the first mode with $m=1$ is associated with the upper QPO frequency, $\nu_u$; that the nonlinear interaction of the first and second modes gives the lower QPO frequency, $\nu_\ell =\nu_u-\nu_*$; and that the nonlinear interaction of the first and third modes gives the lower QPO frequency $\nu_\ell=\nu_u-\nu_*/2$, where $\nu_*=\Omega_*/2\pi$.

The jets observed to emanate from many compact accreting objects may arise from the twisting of the magnetic field threading a differentially rotating accretion disk which acts to magnetically extract angular momentum and energy from the disk. Two main regimes have been discussed, hydromagnetic outflows, which have a significant mass flux and have energy and angular momentum carried by both matter and electromagnetic field and, Poynting outflows, where the mass flux is negligible and energy and angular momentum are carried predominantly by the electromagnetic field. We describe recent theoretical work on the formation of relativistic Poynting jets from magnetized accretion disks and new relativistic, fully-electromagnetic, particle-in-cell simulations of the formation of jets from accretion disks.

We present results of global 3D MHD simulations of disk accretion to a rotating star with dipole and more complex magnetic fields using a Godunov-type code based on the "cubed sphere" grid developed earlier in our group. We describe the code and the grid and show examples of simulation results.

We carried out axisymmetric simulations of disk accretion to a rapidly rotating magnetized star in the "propeller" regime. Simulations show that propellers may be "weak" (with no outflows), and "strong" (with outflows). Investigation of the difference between these two regimes have shown that outflows appear only in the case where the "friction" between the disk and magnetosphere is sufficiently large, and when accreting matter flux is not very small. Matter outflows in a wide cone and is magneto-centrifugally ejected from the inner regions of the disk. Closer to the axis there is a strong, collimated, magnetically dominated outflow of energy and angular momentum carried by the open magnetic field lines from the star. The "efficiency" of the propeller may be very high in the respect that most of the incoming disk matter is expelled from the system in winds. The star spins-down rapidly due to the magnetic interaction with the disk through closed field lines and with corona through open field lines. This mechanism may act in a variety of situations where magnetized star rotates with super-Keplerian velocity at the magnetospheric boundary. We speculate that in general any object rotating with super-Keplerian velocity may drive outflows from accreting disk, if the friction between them is sufficiently large.

A new type of wind - a conical wind - has been discovered in axisymmetric magnetohydrodynamic simulations of the disk-magnetosphere interaction in cases where the magnetic field of the star is bunched into an X-type configuration. Such a configuration arises if the effective viscosity of the disk is larger than the effective diffusivity, or if the accretion rate in the disk is enhanced. Conical outflows flow from the inner edge of the disk into a narrow shell with half-opening angle of 30-45 degrees. The outflow carries about 0.1-0.3 of the disk mass accretion rate and part of the disk's angular momentum. The conical winds are driven by the gradient of the magnetic pressure which exists above the disk due to the winding of the stellar magnetic field. Exploratory 3D simulations show that conical winds are symmetric about rotation axis of the disk even if the magnetic dipole is significantly misaligned with the disk's rotation axis. Conical winds appear around stars of different periods. However, in the case of a star in the "propeller" regime, an additional - much faster component appears: an axial jet, where matter is accelerated up to very high velocities at small distances from the star by magnetic pressure force above the surface of the star. The simulations are done in dimensionless units and are applicable to a variety of the disk-accreting magnetized stars: young stars, white dwarfs, neutron stars, and possibly black holes. For the case of young stars, conical winds and axial jets may appear in different cases, including Class I young stars, classical T Tauri stars, and EXors. In EXors periods of enhanced accretion may lead to the formation of conical winds which correspond to the outflows observed from these stars.

We investigate outflows from the disk-magnetosphere boundary of rotating magnetized stars in cases where the magnetic field of a star is bunched into an X-type configuration using axisymmetric and full 3D MHD simulations. Such configuration appears if viscosity in the disk is larger than diffusivity, or if the accretion rate in the disk is enhanced. Conical outflows flow from the inner edge of the disk to a narrow shell with an opening angle 30-45 degrees. Outflows carry 0.1-0.3 of the disk mass and part of the disk's angular momentum outward. Conical outflows appear around stars of different periods, however in case of stars in the "propeller" regime, an additional - much faster component appears: an axial jet, where matter is accelerated up to very high velocities at small distances from the star by magnetic pressure force above the surface of the star. Exploratory 3D simulations show that conical outflows are symmetric about rotational axis of the disk even if magnetic dipole is significantly misaligned. Conical outflows and axial jets may appear in different types of young stars including Class I young stars, classical T Tauri stars, and EXors.

The appearance and time variability of accreting millisecond X-ray pulsars (hereafter AMXPs, e.g. Wijnands & van der Klis 1998) depends strongly on the accretion rate, the effective viscosity and the effective magnetic diffusivity of the disk-magnetosphere boundary. The accretion rate is the main parameter which determines the location of the magnetospheric radius of the star for a given stellar magnetic field. We introduce a classification of accreting neutron stars as a function of the accretion rate and show the corresponding stages obtained from our global 3D magnetohydrodynamic (MHD) simulations and from our axisymmetric MHD simulations. We discuss the expected variability features in each stage of accretion, both periodic and quasi-periodic (QPOs). We conclude that the periodicity may be suppressed at both very high and very low accretion rates. In addition the periodicity may disappear when ordered funnel flow accretion is replaced by disordered accretion through the interchange instability.

We present results of 3D simulations of MHD instabilities at the accretion disk-magnetosphere boundary. The instability is Rayleigh-Taylor, and develops for a fairly broad range of accretion rates and stellar rotation rates and magnetic fields. It produces tall, thin tongues of plasma that penetrate the magnetosphere in the equatorial plane. The shape and number of the tongues changes with time on the inner-disk dynamical timescale. In contrast with funnel flows, which deposit matter mainly in the polar region, the tongues deposit matter much closer to the stellar equator. The instability appears for relatively small misalignment angles, $\Theta\lesssim30^\circ$, between the star's rotation and magnetic axes, and is associated with higher accretion rates. The hot spots and light curves during accretion through instability are generally much more chaotic than during stable accretion. The unstable state of accretion has possible implications for quasi-periodic oscillations and intermittent pulsations from accreting systems.

This paper investigates the influence of magneto-centrifugally driven or simply magnetic winds of rapidly-rotating, strongly-magnetized T Tauri stars in causing the inward or outward migration of close-in giant planets. The azimuthal ram pressure of the magnetized wind acting on the planet tends to increase the planet's angular momentum and cause outward migration if the star's rotation period $P_*$ is less than the planet's orbital period $P_p$. In the opposite case, $P_* > P_p$, the planet migrates inward. Thus, planets orbiting at distances larger (smaller) than $0.06 {\rm AU}(P_*/5{\rm d})^{2/3}$ tend to be pushed outward (inward), where $P_*$ is the rotation period of the star assumed to have the mass of the sun. The magnetic winds are likely to occur in T Tauri stars where the thermal speed of the gas close to the star is small, where the star's magnetic field is strong, and where the star rotates rapidly. The time-scale for appreciable radial motion of the planet is estimated as $\sim 2 - 20$ Myr. A sufficiently massive close-in planet may cause tidal locking and once this happens the radial migration due to the magnetic wind ceases. The magnetic winds are expected to be important for planet migration for the case of a multipolar magnetic field rather than a dipole field where the wind is directed away from the equatorial plane and where a magnetospheric cavity forms. The influence of the magnetic wind in eroding and eventually destroying the accretion disk is analyzed. A momentum integral is derived for the turbulent wind/disk boundary layer and this is used to estimate the disk erosion time-scale as $\sim 1-10^2$ Myr, with the lower value favored.

We discuss a number of topics relevant to disk-magnetosphere interaction and how numerical simulations illuminate them. The topics include: (1) disk-magnetosphere interaction and the problem of disk-locking; (2) the wind problem; (3) structure of the magnetospheric flow, hot spots at the star's surface, and the inner disk region; (4) modeling of spectra from 3D funnel streams; (5) accretion to a star with a complex magnetic field; (6) accretion through 3D instabilities; (7) magnetospheric gap and survival of protoplanets. Results of both 2D and 3D simulations are discussed.

This work treats the matter deceleration in a magnetohydrodynamics radiative shock wave at the surface of a star. The problem is relevant to classical T Tauri stars where infalling matter is channeled along the star's magnetic field and stopped in the dense layers of photosphere. A significant new aspect of the present work is that the magnetic field has an arbitrary angle with respect to the normal to the star's surface. We consider the limit where the magnetic field at the surface of the star is not very strong in the sense that the inflow is super Alfv\'enic. In this limit the initial deceleration and heating of plasma (at the entrance to the cooling zone) occurs in a fast magnetohydrodynamic shock wave. To calculate the intensity of radiative losses we use "real" and "power-law" radiative functions. We determine the stability/instability of the radiative shock wave as a function of parameters of the incoming flow: velocity, strength of the magnetic field, and its inclination to the surface of the star. In a number of simulation runs with the "real" radiative function, we find a simple criterion for stability of the radiative shock wave. For a wide range of parameters, the periods of oscillation of the shock wave are of the order 0.02-0.2 sec.

We present results of 3D simulations of MHD instabilities at the accretion disk-magnetosphere boundary. The instability is Rayleigh-Taylor, and develops for a fairly broad range of accretion rates and stellar rotation rates and magnetic fields. It manifests itself in the form of tall, thin tongues of plasma that penetrate the magnetosphere in the equatorial plane. The shape and number of the tongues changes with time on the inner-disk dynamical timescale. In contrast with funnel flows, which deposit matter mainly in the polar region, the tongues deposit matter much closer to the stellar equator. The instability appears for relatively small misalignment angles, $\Theta\lesssim30^\circ$, between the star's rotation and magnetic axes, and is associated with higher accretion rates. The hot spots and light curves during accretion through instability are generally much more chaotic than during stable accretion. The unstable state of accretion has possible implications for quasi-periodic oscillations and intermittent pulsations from accreting systems, as well as planet migration.

Disk accretion to rotating stars with complex magnetic fields is investigated using full three-dimensional magnetohydrodynamic (MHD) simulations. The studied magnetic configurations include superpositions of misaligned dipole and quadrupole fields and off-centre dipoles. The simulations show that when the quadrupole component is comparable to the dipole component, the magnetic field has a complex structure with three major magnetic poles on the surface of the star and three sets of loops of field lines connecting them. A significant amount of matter flows to the quadrupole "belt", forming a ring-like hot spot on the star. If the maximum strength of the magnetic field on the star is fixed, then we observe that the mass accretion rate, the torque on the star, and the area covered by hot spots are several times smaller in the quadrupole-dominant cases than in the pure dipole cases. The influence of the quadrupole component on the shape of the hot spots becomes noticeable when the ratio of the quadrupole and dipole field strengths $B_q/B_d\gtrsim0.5$, and becomes dominant when $B_q/B_d\gtrsim1$. In the case of an off-centre dipole field, most of the matter flows through a one-armed accretion stream, forming a large hot spot on the surface, with a second much smaller secondary spot. The light curves may have simple, sinusoidal shapes, thus mimicking stars with pure dipole fields. Or, they may be complex and unusual. In some cases the light curves may be indicators of a complex field, in particular if the inclination angle is known independently. We also note that in the case of complex fields, magnetospheric gaps are often not empty, and this may be important for the survival of close-in exosolar planets.

We present three-dimensional (3-D) simulations of rotationally induced line variability arising from complex circumstellar environment of classical T Tauri stars (CTTS) using the results of the 3-D magnetohydrodynamic (MHD) simulations of Romanova et al., who considered accretion onto a CTTS with a misaligned dipole magnetic axis with respect to the rotational axis. The density, velocity and temperature structures of the MHD simulations are mapped on to the radiative transfer grid, and corresponding line source function and the observed profiles of neutral hydrogen lines (H-beta, Pa-beta and Br-gamma) are computed using the Sobolev escape probability method. We study the dependency of line variability on inclination angles (i) and magnetic axis misalignment angles (Theta). By comparing our models with the Pa-beta profiles of 42 CTTS observed by Folha & Emerson, we find that models with a smaller misaligngment angle (Theta<~15 deg.) are more consistent with the observations which show that majority of Pa-beta are rather symmetric around the line centre. For a high inclination system with a small dipole misalignment angle (Theta ~ 15 deg.), only one accretion funnel (on the upper hemisphere) is visible to an observer at any given rotational phase. This can cause an anti-correlation of the line equivalent width in the blue wing (v<0) and that in the red wing (v>0) over a half of a rotational period, and a positive correlation over other half. We find a good overall agreement of the line variability behaviour predicted by our model and those from observations. (Abridged)

We report on the first global three-dimensional (3D) MHD simulations of disk accretion onto a rotating magnetized star through the Rayleigh-Taylor instability. The star has a dipole field misaligned relative to the rotation axis by a small angle Theta. Simulations show that, depending on the accretion rate, a star may be in the stable or unstable regime of accretion. In the unstable regime, matter penetrates deep into the magnetosphere through several elongated "tongues" which deposit matter at random places on the surface of the star, leading to stochastic light-curves. In the stable regime, matter accretes in ordered funnel streams and the light-curves are almost periodic. A star may switch between these two regimes depending on the accretion rate and may thus show alternate episodes of ordered pulsations and stochastic light-curves. In the intermediate regime, both stochastic and ordered pulsations are observed. For Theta > 30 degrees, the instability is suppressed and stable accretion through funnel streams dominates.

We calculated profiles of CIV 1550, Si IV 1400, NV 1240 and OVI 1035 doublet lines using results of 3D MHD simulations of disc accretion onto young stars with dipole magnetic field. It appeared that our calculations can not reproduce profiles of these lines observed (HST/GHRS-STIS and FUSE) in CTTSs's spectra. We also found that the theory predicts much larger C IV 1550 line flux than observed (up to two orders of magnitude in some cases) and argue that the main portion of accretion energy in CTTSs is liberated outside accretion shock. We conclude that the reason of disagreement between the theory and observation is strongly non-dipole character of CTTS's magnetic field near its surface.

We discuss three specific modes of accretion disks around rotating magnetized neutron stars which may explain the separations of the kilo Hertz quasi periodic oscillations (QPO) seen in low mass X-ray binaries. The existence of these modes requires that there be a maximum in the angular velocity of the accreting material, and that the fluid is in stable, nearly circular motion near this maximum rather than moving rapidly towards the star or out of the disk plane into funnel flows. It is presently not known if these conditions occur, but we are exploring this with 3D magnetohydrodynamic simulations and will report the results elsewhere. The first mode is a corotation mode which is radially trapped in the vicinity of the maximum of the disk rotation rate and is unstable. The second mode, relevant to relatively slowly rotating stars, is a magnetically driven eccentric ($m=1$) oscillation of the disk excited at a Lindblad radius in the vicinity of the maximum of the disk rotation. The third mode, relevant to rapidly rotating stars, is a magnetically coupled eccentric ($m=1$) and an axisymmetric ($m=0$) radial disk perturbation which has an inner Lindblad radius also in the vicinity of the maximum of the disk rotation. We suggest that the first mode is associated with the upper QPO frequency, $\nu_u$, the second with the lower QPO frequency, $\nu_\ell =\nu_u-\nu_*$, and the third with the lower QPO frequency, $\nu_\ell=\nu_u-\nu_*/2$, where $\nu_*$ is the star's rotation rate.

Disc accretion to a rotating star with a non-dipole magnetic field is investigated for the first time in full three-dimensional (3D) magnetohydrodynamic (MHD) simulations. We investigated the cases of (1) pure dipole, (2) pure quadrupole, and (3) dipole plus quadrupole fields. Simulations have shown that in each case the structure of the funnel streams and associated hot spots on the surface of the star have specific features connected with the magnetic field configuration. In the case (1), matter accretes in two funnel streams which form two arch-like spots near the magnetic poles. In the case (3), most of the matter flows through the quadrupole "belt" forming a ring-shaped hot region on the magnetic equator. In the case (2), magnetic flux in the northern magnetic hemisphere is larger than that in the southern, and the quadrupole belt and the ring are displaced to the south. The stronger the quadrupole, the closer the ring is to the magnetic equator. At sufficiently large misalignment angle $\Theta$, matter also flows to the south pole, forming a hot spot near the pole. The light curves have a variety of different features which makes it difficult to derive the magnetic field configuration from the light curves. There are specific features which are different in cases of dipole and quadrupole dominated magnetic field: (1) Angular momentum flow between the star and disc is more efficient in the case of the dipole field; (2) Hot spots are hotter and brighter in case of the dipole field because the matter accelerates over a longer distance compared with the flow in a quadrupole case.

We discuss axisymmetric force-free pulsar magnetospheres with magnetically collimated jets and a disk-wind obtained by numerical solution of the pulsar equation. This solution represents an alternative to the quasi-spherical wind solutions where a major part of the current flow is in a current sheet which is unstable to magnetic field annihilation.

We discuss a self-consistent unified model of the matter outflows from AGNs based on a theoretical approach and involving data on AGN evolution and structure. The model includes a unified geometry, two-phase gas dynamics, radiation transfer, and absorption spectrum calculations in the UV and X-ray bands. We briefly discuss several questions about the mass sources of the flows, the covering factors, and the stability of the narrow absorption details.

Thermonuclear burning on the surface of a neutron star causes the expansion of a thin outer layer of the star, $\Delta R(t)$. The layer rotates slower than the star due to angular momentum conservation. The shear between the star and the layer acts to twist the star's dipole magnetic field giving at first a trailing spiral field. The twist of the field acts in turn to `torque up' the layer increasing its specific angular momentum. As the layer cools and contracts, its excess specific angular momentum causes it to {\it rotate faster} than the star which gives a leading spiral magnetic field. The process repeats, giving rise to torsional oscillations. We derive equations for the angular velocity and magnetic field of the layer taking into account the diffusivity and viscosity which are probably due to turbulence. The magnetic field causes a nonuniformity of the star's photosphere (at the top of the heated layer), and this gives rise to the observed X-ray oscillations. The fact that the layer periodically rotates faster than the star means that the X-ray oscillation frequency may ``overshoot'' the star's rotation frequency. Comparison of the theory is made with observations of Chakrabarty et al. (2003) of an X-ray burst of SAX J1808.4-3658.

We investigate spherical accretion to a rotating magnetized star in the "propeller" regime using axisymmetric resistive magnetohydrodynamic simulations. The regime is predicted to occur if the magnetospheric radius is larger than the corotation radius and smaller than the light cylinder radius. The simulations show that accreting matter is expelled from the equatorial region of the magnetosphere and that it moves away from the star in a supersonic, disk-shaped outflow. At larger radial distances the outflow slows down and becomes subsonic. The equatorial matter outflow is initially driven by the centrifugal force, but at larger distances the pressure gradient force becomes significant. We find the fraction of the Bondi accretion rate which accretes to the surface of the star.

We use axisymmetric magnetohydrodynamic simulations to investigate the spinning-down of magnetars rotating in the propeller regime and moving supersonically through the interstellar medium. The simulations indicate that magnetars spin-down rapidly due to this interaction, faster than for the case of a non-moving star. From many simulation runs we have derived an approximate scaling laws for the angular momentum loss rate, \dot{L} \propto \~\eta_m^{0.3}\mu^{0.6}\rho^{0.8}{\cal M}^{-0.4} \Omega_*^{1.5}, where \rho is the density of the interstellar medium, \cal M is Mach number, \mu is the star's magnetic moment, \Omega_* is its angular velocity, and \eta_m is magnetic diffusivity. A magnetar with a surface magnetic field of 10^{13} - 10^{15} G is found to spin-down to a period P > 10^5-10^6 s in \sim 10^4 - 10^5 years. There is however uncertainty about the value of the magnetic diffusivity so that the time-scale may be longer. We discuss this model in respect of Soft Gamma Repeaters (SGRs) and the isolated neutron star candidate RXJ1856.5-3754.

The bunching of giant planets at a distance of several stellar radii may be explained by the disruption of the inner part of the disk by the magnetosphere of the star during the T Tauri stage of evolution. The rotating magnetic field of the star gives rise to a low density magnetospheric gap where stellar migration is strongly suppressed. We performed full 3D magnetohydrodynamic simulations of the disk-magnetosphere interaction and examined conditions for which the magnetospheric gap is "empty", by changing the misalignment angle between magnetic and rotational axes of the star, Theta, and by lowering the adiabatic index gamma, which mocks up the effect of heat conductivity and cooling. Our simulations show that for a wide range of plausible conditions the gap is essentially empty. However, in the case of large misalignment angles Theta, part of the funnel stream is located in the equatorial plane and the gap is not empty. Furthermore, if the adiabatic index is small (gamma=1.1) and the rotational and magnetic axes are almost aligned, then matter penetrates through the magnetosphere due to 3D instabilities forming high-density equatorial funnels. For these two limits there is appreciable matter density in the equatorial plane of the disk so that a planet may migrate into the star.

The inner 0.1 AU around accreting T Tauri stars hold clues to many physical processes that characterize the early evolution of solar-type stars. The accretion-ejection connection takes place at least in part in this compact magnetized region around the central star, with the inner disk edge interacting with the star's magnetosphere thus leading simultaneously to magnetically channeled accretion flows and to high velocity winds and outflows. The magnetic star-disk interaction is thought to have strong implications for the angular momentum evolution of the central system, the inner structure of the disk, and possibly for halting the migration of young planets close to the stellar surface. We review here the current status of magnetic field measurements in T Tauri stars, the recent modeling efforts of the magnetospheric accretion process, including both radiative transfer and multi-D numerical simulations, and summarize current evidence supporting the concept of magnetically-channeled accretion in young stars. We also discuss the limits of the models and highlight observational results which suggest that the star-disk interaction is a highly dynamical and time variable process in young stars.

We performed 2D and full 3D magnetohydrodynamic simulations of disk accretion to a rotating star with an aligned or misaligned dipole magnetic field. We investigated the rotational equilibrium state and derived from simulations the ratio between two main frequencies: the spin frequency of the star and the orbital frequency at the inner radius of the disk. In 3D simulations we observed different features related to the non-axisymmetry of the magnetospheric flow. These features may be responsible for high-frequency quasi-periodic oscillations (QPOs). Variability at much lower frequencies may be connected with restructuring of the magnetic flux threading the inner regions of the disk. Such variability is specifically strong at the propeller stage of evolution.

We present results of axisymmetic magnetohydrodynamic simulations of the interaction of a rapidly-rotating, magnetized star with an accretion disk. The disk is considered to have a finite viscosity and magnetic diffusivity. The main parameters of the system are the star's angular velocity and magnetic moment, and the disk's viscosity, diffusivity. We focus on the "propeller" regime where the inner radius of the disk is larger than the corotation radius. Two types of magnetohydrodynamic flows have been found as a result of simulations: "weak" and "strong" propellers. The strong propeller is characterized by a powerful disk wind and a collimated magnetically dominated outflow or jet from the star. The weak propeller have only weak outflows. We investigated the time-averaged characteristics of the interaction between the main elements of the system, the star, the disk, the wind from the disk, and the jet. Rates of exchange of mass and angular momentum between the elements of the system are derived as a function of the main parameters. The propeller mechanism may be responsible for the fast spinning-down of the classical T Tauri stars in the initial stages of their evolution, and for the spinning-down of accreting millisecond pulsars.

We report the discovery of propeller-driven outflows in axisymmetric magnetohydrodynamic simulations of disk accretion to rapidly rotating magnetized stars. Matter outflows in a wide cone and is centrifugally ejected from the inner regions of the disk. Closer to the axis there is a strong, collimated, magnetically dominated outflow of energy and angular momentum carried by the open magnetic field lines from the star. The ``efficiency'' of the propeller may be very high in the respect that most of the incoming disk matter is expelled from the system in winds. The star spins-down rapidly due to the magnetic interaction with the disk through closed field lines and with corona through open field lines. Diffusive and viscous interaction between magnetosphere and the disk are important: no outflows were observed for very small values of the diffusivity and viscosity. These simulation results are applicable to the early stages of evolution of classical T Tauri stars and to different stages of evolution of cataclysmic variables and neutron stars in binary systems. As an example, we have shown that young rapidly rotating magnetized CTTSs spin-down to their present slow rotation in less than 10^6 years.

We investigate the rotational equilibrium state of a disk accreting magnetized stars using axisymmetric magnetohydrodynamic (MHD) simulations. In this ``locked'' state, the spin-up torque balances the spin-down torque so that the net average torque on the star is zero. We investigated two types of initial conditions, one with a relatively weak stellar magnetic field and a high coronal density, and the other with a stronger stellar field and a lower coronal density. We observed that for both initial conditions the rotation of the star is locked to the rotation of the disk. In the second case, the radial field lines carry significant angular momentum out of the star. However, this did not appreciably change the condition for locking of the rotation of the star. We find that in the equilibrium state the corotation radius $r_{co}$ is related to the magnetospheric radius $r_A$ as $r_{co}/r_A\approx 1.2-1.3$ for case (1) and $r_{co}/r_A\approx 1.4-1.5$ for case (2). We estimated periods of rotation in the equilibrium state for classical T Tauri stars, dwarf novae and X-ray millisecond pulsars.

An analytical model is developed for the screening of the external magnetic field of a rotating, axisymmetric neutron star due to the accretion of plasma from a disk. The decrease of the field occurs due to the electric current in the infalling plasma. The deposition of this current carrying plasma on the star's surface creates an induced magnetic moment with a sign opposite to that of the original magnetic dipole. The field decreases independent of whether the star spins-up or spins-down. The time-scale for an appreciable decrease (factor of $>100$) of the field is found to be $\sim 1.6 \times 10^7$ yr, for a mass accretion rate $\dot{M}=10^{-9} M_\odot/$yr and an initial magnetic moment $\mu_i = 10^{30}{\rm G cm}^3$ which corresponds to a surface field of $10^{12} $G if the star's radius is $10^6$ cm. The time-scale varies approximately as $\mu_i^{3.8}/\dot{M}^{1.9}$. The decrease of the magnetic field does not have a simple relation to the accreted mass. Once the accretion stops the field leaks out on an Ohmic diffusion time scale which is estimated to be $ > 10^9$ yr.

We model the variability profiles of millisecond period X-ray pulsars. We performed three-dimensional magnetohydrodynamic simulations of disk accretion to millisecond period neutron stars with a misaligned magnetic dipole moment, using the pseudo-Newtonian Paczynski-Wiita potential to model general relativistic effects. We found that the shapes of the resulting funnel streams of accreting matter and the hot spots on the surface of the star are quite similar to those for more slowly rotating stars obtained from earlier simulations using the Newtonian potential. The funnel streams and hot spots rotate approximately with the same angular velocity as the star. The spots are bow-shaped (bar-shaped) for small (large) misalignment angles. We found that the matter falling on the star has a higher Mach number when we use the Paczynski-Wiita potential than in the Newtonian case. Having obtained the surface distribution of the emitted flux, we calculated the variability curves of the star, taking into account general relativistic, Doppler and light-travel-time effects. We found that general relativistic effects decrease the pulse fraction (flatten the light curve), while Doppler and light-travel-time effects increase it and distort the light curve. We also found that the light curves from our hot spots are reproduced reasonably well by spots with a gaussian flux distribution centered at the magnetic poles. We also calculated the observed image of the star in a few cases, and saw that for certain orientations, both the antipodal hot spots are simultaneously visible, as noted by earlier authors.

A model is developed for the time dependent electromagnetic - radio to gamma-ray - emission of active galactic nuclei, specifically, the blazars, based on the acceleration and creation of leptons at a propagating discontinuity or {\it front} of a Poynting flux jet. The front corresponds to a discrete relativistic jet component as observed with very-long-baseline-interferometry (VLBI). Equations are derived for the number, momentum, and energy of particles in the front taking into account synchrotron, synchrotron-self-Compton (SSC), and inverse-Compton processes as well as photon-photon pair production. The apparent synchrotron, SSC, and inverse-Compton luminosities as functions of time are determined. Predictions of the model are compared with observations in the gamma, optical and radio bands. The delay between the high-energy gamma-ray flare and the onset of the radio is explained by self-absorption and/or free-free absorption by external plasma. Two types of gamma-ray flares are predicted depending on pair creation in the front.

We investigate the emission of rotating magnetized neutron stars due to the acceleration and radiation of particles in the relativistic wind and in the magnetotail of the star. We consider that the charged particles are accelerated by driven collisionless reconnection. Outside of the light cylinder, the star's rotation acts to wind up the magnetic field to form a predominantly azimuthal, slowly decreasing with distance, magnetic field of opposite polarity on either side of the equatorial plane normal to the star's rotation axis. The magnetic field annihilates across the equatorial plane with the magnetic energy going to accelerate the charged particles to relativistic energies. For a typical supersonically moving pulsar, the star's wind extends outward to the standoff distance with the interstellar medium. At larger distances, the power output of pulsar's wind $\dot{E}_w$ of electromagnetic field and relativistic particles is {\it redirected and collimated into the magnetotail} of the star. In the magnetotail it is proposed that equipartition is reached between the magnetic energy and the relativistic particle energy. For such conditions, synchrotron radiation from the magnetotails may be a significant fraction of $\dot{E}_w$ for high velocity pulsars. An equation is derived for the radius of the magnetotail $r_m(z^\prime)$ as a function of distance $z^\prime$ from the star. For large distances $z^\prime$, of the order of the distance travelled by the star, we argue that the magnetotail has a `trumpet' shape owing to the slowing down of the magnetotail flow.

It is shown that the existence of hot accretion spots at the surface of T Tau practically has no effect on the accuracy of estimation of its magnetic field strength at photospheric level. We also found that one can interpret results of T Tau's photospheric magnetic field strength measurements carried out via different methods in the frame of the following alternative: 1) if T Tau's inclination angle $i\le 10^\circ,$ then magnetic field of the star may be dipolar with the angle between rotational and magnetic axes is near $85^\circ;$ 2) if it will be found (e.g. from interferometric observations) that $i>10^\circ,$ then magnetic field of T Tau is essentially non-dipolar or/and non-stationary.

The propeller regime of disk accretion to a rapidly rotating magnetized star is investigated here for the first time by axisymmetric 2.5D magnetohydrodynamic simulations. An expanded, closed magnetosphere forms in which the magnetic field is predominantly toroidal. A smaller fraction of the star's poloidal magnetic flux inflates vertically, forming a magnetically dominated tower. Matter accumulates in the equatorial region outside magnetosphere and accretes to the star quasi-periodically through elongated funnel streams which cause the magnetic field to reconnect. The star spins-down owing to the interaction of the closed magnetosphere with the disk. For the considered conditions, the spin-down torque varies with the angular velocity of the star omega* as omega*^1.3 for fixed mass accretion rate. The propeller stage may be important in the evolution of X-ray pulsars, cataclysmic variables and young stars. In particular, it may explain the present slow rotation of the classical T Tauri stars.

The jets observed to emanate from many compact accreting objects may arise from the twisting of a magnetic field threading a differentially rotating accretion disk which acts to magnetically extract angular momentum and energy from the disk. Two main regimes have been discussed, hydromagnetic jets, which have a significant mass flux and have energy and angular momentum carried by both matter and electromagnetic field and, Poynting jets, where the mass flux is small and energy and angular momentum are carried predominantly by the electromagnetic field. Here, we describe recent theoretical work on the formation of relativistic Poynting jets from magnetized accretion disks. Further, we describe new relativistic, fully-electromagnetic, particle-in-cell simulations of the formation of jets from accretion disks. Analog Z-pinch experiments may help to understand the origin of astrophysical jets.

The physics of the "hot spots" on stellar surfaces and the associated variability of accreting magnetized rotating stars is investigated for the first time using fully three-dimensional magnetohydrodynamic simulations. The magnetic moment of the star is inclined relative to its rotation axis by an angle Theta. A sequence of misalignment angles was investigated, between Theta=0 and 90 degrees.Typically at small Theta the spots are observed to have the shape of a bow which is curved around the magnetic axis, while at largest Theta the spots have a shape of a bar, crossing the magnetic pole. The physical parameters (density, temperature, etc.) increase toward the central regions of the spots. At relatively low density and temperature, the spots occupy approximately 10-20 % of the stellar surface, while at the highest values of these parameters this area may be less than 1 % of the area of the star. The light curves were calculated for different Theta and inclination angles of the disk i. They show a range of variability patterns, including one maximum-per-period curves (at most of angles Theta and i), and two maximum-per-period curves (at large Theta and i). At small Theta, the funnel streams may rotate faster/slower than the star, and this may lead to quasi-periodic variability of the star. The results are of interest for understanding the variability and quasi-variability of Classical T Tauri Stars, millisecond pulsars and cataclysmic variables.

A model is developed for relativistic Poynting jets from the inner region of a disk around a rotating black hole. The disk is initially threaded by a dipole-like magnetic field. The model is derived from the special relativistic equation for a force-free electromagnetic field. The ``head'' of the Poynting jet is found to propagate outward with a velocity which may be relativistic. The Lorentz factor of the head (Gamma) is found to be dependent on the magnetic field strength close to the black hole, B_0, the density of the external medium n_ext, and on the ratio R=r_0/r_g >1, where r_g is the gravitational radius of the black hole, and r_0 is the radius of the O-point of the initial dipole field threading the disk. For conditions pertinent to an active galactic nuclei, Gamma is approximately equal to 8 (10/R)^(1/3) (B_0/10^3 Gauss)^(1/3) (1/cm^3/n_ext)^(1/6). This model offers an explanation for the observed Lorentz factors which are of the order of 10 for the parsec-scale radio jets measured with very long baseline interferometry.

We present results of fully three-dimensional MHD simulations of disk accretion to a rotating magnetized star with its dipole moment inclined at an angle Theta to the rotation axis of the disk. We observed that matter accretes from the disk to a star in two or several streams depending on Theta. Streams may precess around the star at small Theta. The inner regions of the disk are warped. The warping is due to the tendency of matter to co-rotate with inclined magnetosphere. The accreting matter brings positive angular momentum to the (slowly rotating) star tending to spin it up. The corresponding torque N_z depends only weakly on Theta. The angular momentum flux to the star is transported predominantly by the magnetic field; the matter component contributes < 1 % of the total flux. Results of simulations are important for understanding the nature of classical T Tauri stars, cataclysmic variables, and X-ray pulsars.

We give further considerations on the problem of the evolution of a coronal, force-free magnetic field which threads a differentially rotating, conducting Keplerian disk, extending the work of Li {\it et al.} (2001). This situation is described by the force-free Grad-Shafranov (GS) equation for the flux function $\Psi(r,z)$ which labels the poloidal field lines (in cylindrical coordinates). The GS equation involves a function $H(\Psi)$ describing the distribution of poloidal current which is determined by the differential rotation or {\it twist} of the disk which increases linearly with time. We numerically solve the GS equation in a sequence of volumes of increasing size corresponding to the expansion of the outer perfectly conducting boundaries at ($R_{m}, Z_{m}$). The outer boundaries model the influence of an external non-magnetized plasma. The sequence of GS solutions provides a model for the dynamical evolution of the magnetic field in response to (1) the increasing twist of the disk and (2) the pressure of external plasma. We find solutions with {\it magnetically collimated} Poynting jets where there is a {\it continuous} outflow of energy, angular momentum, and toroidal magnetic flux from the disk into the external space.

We have calculated self-consistent equilibria of a collisionless relativistic electron-positron gas in the vicinity of a magnetic $X-$point. For the considered conditions, pertinent to extra-galactic jets, we find that leptons are accelerated up to Lorentz factors $\Gamma_0 = \kappa e B_0L {\cal E}^2 /mc^2 \gg 1$, where $B_0$ is the typical magnetic field strength, ${\cal E}\equiv E_0/B_0$, with $E_0$ the reconnection electric field, $L$ is the length scale of the magnetic field, and $\kappa \approx 12$. The acceleration is due to the dominance of the electric field over the magnetic field in a region around the $X-$point. The distribution function of the accelerated leptons is found to be approximately $dn/d\gamma \propto \gamma^{-1}$ for $\gamma \lesssim \Gamma_0$. The apparent distribution function may be steeper than $\gamma^{-1}$ due to the distribution of $\Gamma_0$ values and/or the radiative losses. Self-consistent equilibria are found only for plasma inflow rates to the $X-$point less than a critical value.

We describe a three-dimensional, Godunov-type numerical magnetohydrodynamics (MHD) method designed for studying disk accretion to a rotating magnetized star in the general case where the star's rotation axis, its magnetic moment, and the normal to the disk all have different directions. The numerical method uses a "cubed sphere" coordinate system which has advantages of Cartesian and spherical coordinate systems but does not have the singular axis ofthe spherical system. The grid is formed by a sequence of concentric spheres of radii $R_j \propto q^j$ with $j=1..N_R$ and $q={\rm const}>1$. The grid on the surface of the sphere consists of six sectors with the grid on each sector topologically equivalent to the equidistant grid on the face of a cube. Simulation results are discussed for the funnel flows (FF) to a star with dipole moment $\rvecmu$ at an angle $\Theta=30^\circ$ to the star's rotation axis ${\bf \Omega}$ which is aligned with the normal to the disk. Two important new 3D features are found in these simulations: (1) The funnel flow to the stellar surface is mainly in two streams which approach the star from opposite directions. (2) In the $x-z$ cross section of the flow containing $\rvecmu$ and ${\bf \Omega}$, the funnel flow often takes the longer of the two possible paths along magnetic field lines to the surface of the star. A subsequent paper will give a detailed description of the method and results on 3D funnel flows at different inclination angles $\Theta$.

This work investigates Bondi accretion to a rotating magnetized star in the "propeller" regime using axisymmetric resistive, magnetohydrodynamic simulations. In this regime accreting matter tends to be expelled from the equatorial region of the magnetosphere where the centrifugal force on matter rotating with the star exceeds the gravitational force. The regime is predicted to occur if the magnetospheric radius larger than the corotation radius and less than the light cylinder radius. The simulations show that accreting matters is expelled from the equatorial region of the magnetosphere and that it moves away from the star in a supersonic, disk-shaped outflow. At larger radial distances the outflow slows down and becomes subsonic. The equatorial matter outflow is initially driven by the centrifugal force, but at larger distances the pressure gradient becomes significant. We find that the star is spun-down mainly by the magnetic torques at its surface with the rate of loss of angular momentum $\dot{L}$ proportional to $-\Omega_*^{1.3}\mu^{0.8}$, where $\Omega_*$ is the star's rotation rate and $\mu$ is its magnetic moment. Further, we find that $\dot{L}$ is approximately independent of the magnetic diffusivity of the plasma $\eta_m$. The fraction of the Bondi accretion rate which accretes to the surface of the star is found to be $\propto\Omega_*^{-1.0}\mu^{-1.7}\eta_m^{0.4}$. Predictions of this work are important for the observability of isolated old neutron stars and for wind fed pulsars in X-ray binaries.

Magnetohydrodynamic (MHD) simulations have been used to study disk accretion to a rotating magnetized star with an aligned dipole moment. Quiescent initial conditions were developed in order to avoid the fast initial evolution seen in earlier studies. A set of simulations was performed for different stellar magnetic moments and rotation rates. Simulations have shown that the disk structure is significantly changed insidea radius r_{br} where magnetic braking is significant. In this region the disk is strongly inhomogeneous. Radial accretion of matter slows as it approaches the area of strong magnetic field and a dense ring and funnel flow form at the magnetospheric radius r_m where the magnetic pressure is equal to the total, kinetic plus thermal, pressure of the matter.

This work investigates the propagation of magnetized, isolated old neutron stars through the interstellar medium. We performed axisymmetric non relativistic magnetohydrodynamic (MHD) simulations of the propagation of a non-rotating star with dipole magnetic field aligned with its velocity through the interstellar medium (ISM). Effects of rotation will be discussed in a subsequent work. We consider two cases: (1) where theaccretion radius $R_{acc}$ is comparable to the magnetic standoff distance or Alfv\'en radius $R_A$ and gravitational focusing is important; and (2) where $R_{acc}<< R_A$ and the magnetized star interacts with the ISM as a ``magnetic plow'', without significant gravitational focusing. For the first case simulations were done at a low Mach number ${\cal M}=3$ for a range of values of the magnetic field $B_*$. For the second case, simulations were done for higher Mach numbers, ${\cal M}=10,$ 30, and 50. In both cases, the magnetosphere of the star represents an obstacle for the flow, and a shock wavestands in front of the star. Magnetic field lines are stretched downwind from the star and form a hollow elongated magnetotail. Reconnection of the magnetic field is observed in the tail which may lead to acceleration of particles. Similar powers are estimated to be released in the bow shock wave and in the magnetotail. The estimated powers are, however, below present detection limits. Results of our simulations maybe applied to other strongly magnetized stars, for example, white dwarfs and magnetic Ap stars. Future more sensitive observations may reveal long magnetotails of magnetized stars moving through the ISM.

This work considers flows from an accretion disk corotating with the aligned dipole magnetic field of a rotating star. Ideal magnetohydrodynamics (MHD) is assumed with the pressure and density related as $p \propto \rho^\gamma$ and with $\rho {\bf v}^2 \ll {\bf B}^2/4\pi$, where ${\bf v}$ is the flow velocity. This limit corresponds to the Alfv\'en radius for the disk accretion larger than the corotation radius. Transonic flows, which go from subsonic motion near the disk to supersonic inflow near the star, are shown to be possible only for a narrow range of $R_d \sim r_c \equiv (GM/\Omega^2)^{1/3}$, where $R_d$ is the radius at which the dipole field line intersects the disk, $r_c$ the corotation radius, $M$ the mass of the star, and $\Omega$ its angular rotation rate. The transonic flows have very different behaviors for $\gamma > 7/5$ and $< 7/5$. In both cases, the plasma flow velocity ${\bf v}$ (which is parallel to ${\bf B}$) increases with decreasing distance $R$ from the star. However, for $\gamma >7/5$, the Mach number ${\cal M} \equiv |{\bf v}|/c_s$ (with $c_s$ the sound speed) initially increases with $R$ decreasing from $R_d$, but for $R$ decreasing from $\approx 0.22 R_d$ (for $\gamma=5/3$)the Mach number surprisingly {\it decreases}. In the other limit, $\gamma < 7/5$, ${\cal M}$ increases monotonically with decreasing $R$. Application of these results is made to funnel flows to rotating magnetized neutron stars and young stellar objects.

The powerful narrow jets observed to emanate from many compact accreting objects may arise from the twisting of a magnetic field threading a differentially rotating accretion disk which acts to magnetically extract angular momentum and energy from the disk. Two main regimes have been discussed, {\it hydromagnetic outflows}, which have a significant mass flux and have energy and angular momentum carried by both the matter and the electromagnetic field and, Poynting outflows, where the mass flux is negligible and energy and angular momentum are carried predominantly by the electromagnetic field. Here we consider a Keplerian disk initially threaded by a dipole-like magnetic field and we present solutions of the force-free Grad-Shafranov equation for the coronal plasma. We find solutions with Poynting jets where there is a continuous outflow of energy and toroidal magnetic flux from the disk into the external space. This behavior contradicts the commonly accepted ``theorem'' of Solar plasma physics that the motion of the footpoints of a magnetic loop structure leads to a stationary magnetic field configuration with zero power and flux outflows. In addition we discuss recent magnetohydrodynamic (MHD) simulations which establish that quasi-stationary collimated Poynting jets similar to our Grad-Shafranov solutions arise from the inner part of a disk threaded by a dipole-like magnetic field. At the same time we find that there is a steady uncollimated hydromagnetic outflow from the outer part of the disk. The Poynting jets represent a likely model for the jets from active galactic nuclei, microquasars, and gamma ray burst sources.

Models of Blazars based on the propagation of finite discontinuities or fronts in the Poynting flux jet from the innermost regions of an accretion disk around a black hole are discussed. Such fronts may be responsible for short time-scale (from less than hours to days) flares in different wavebands from high frequency radioband to TeV, with delay in low radio frequencies as a result of synchrotron self-absorption. The cases of magnetic fields of one and opposite polarities across the front are investigated. We find that annihilation of magnetic field in the front leads to higher energy spectrum of leptons and possibility of strong TeV flares. Electron-positron pairs form in most cases as a result of interaction between numerous synchrotron photons and SSC photons, and constitute the majority species, compared with the ions at subparsec scales. Frequent weak outbursts may be responsible for flickering core radiation in all wavebands, while the stronger outbursts may be observed as short time-scale flares.

Stationary magnetohydrodynamic (MHD) outflows from a rotating, conducting Keplerian accretion disk threaded by B-field are investigated numerically by time-dependent, axisymmetric (2.5D) simulations using a Godunov-type code. A large class of stationary magneto-centrifugally driven winds are found where matter is accelerated from a thermal speed at the disk to much larger velocity, greater than the fast magnetosonic speed and larger than the escape speed. The flows are approximately spherical outflows with only small collimation within the simulation region. Numerical results are shown to coincide with the theoretical predictions of ideal, axisymmetric MHD to high accuracy. Investigation of the influence of outer boundary conditions, particularly that on the toroidal component of magnetic field shows that the commonly used ``free'' boundary condition leads to artificial magnetic forces which can act to give spurious collimation. New boundary conditions are proposed which do not generate artificial forces. Artificial results may also arise for cases where the Mach cones on the outer boundaries are partially directed into the simulation region.

A model is developed for magnetic `propeller'-driven outflows which cause a rapidly rotating magnetized star accreting from a disk to spin-down. Energy and angular momentum lost by the star goes into expelling most of the accreting disk matter. The theory gives an expression for the effective Alfven radius $R_A$ (where the inflowing matter is effectively stopped) which depends on the mass accretion rate, the star's mass and magnetic moment, and the star's rotation rate. The model points to a mechanism for `jumps' between spin-down and spin-up evolution and for the reverse transition, which are changes between two possible equilibrium configurations of the system. In for example the transistion from spin-down to spin-up states the Alfven radius $R_A$ decreases from a value larger than the corotation radius to one which is smaller. In this transistion the `propeller' goes from being ``on'' to being ``off.'' The ratio of the spin-down to spin-up torque (or the ratio for the reverse change) in a jump is shown to be of order unity.

Quasi-spherical supersonic (Bondi-type) accretion to a star with a dipole magnetic field is investigated using resistive magnetohydrodynamic simulations. A systematic study is made of accretion to a non-rotating star, while sample results for a rotating star are also presented. A new stationary subsonic accretion flow is found with a steady rate of accretion to the magnetized star smaller than the Bondi accretion rate. Dependences of the accretion rate and the flow pattern on the magnetic momentum of the star and the magnetic diffusivity are presented. For slow star's rotation the accretion flow is similar to that in non-rotating case, but in the case of fast rotation the structure of the subsonic accretion flow is fundamentally different and includes a region of ``propeller'' outflow. The methods and results described here are of general interest and can be applied to systems where matter accretes with low angular momentum.

Hydrodynamic simulations have been used to study accretion disks consisting of counterrotating components with an intervening shear layer(s). Configurations of this type can arise from the accretion of newly supplied counterrotating matter onto an existing corotating disk. The grid-dependent numerical viscosity of our hydro code is used to simulate the influence of a turbulent viscosity of the disk. Firstly, we consider the case where the gas well above the disk midplane rotates with angular rate +\Omega(r) and that well below has the same properties but rotates with rate -\Omega(r). We find that there is angular momentum annihilation in a narrow equatorial boundary layer in which matter accretes supersonically with a velocity which approaches the free-fall velocity and the average accretion speed of the disk can be enormously larger than that for a conventional \alpha-disk rotating in one direction. Secondly, we consider the case of a corotating accretion disk for r<r_t and a counterrotating disk for r>r_t. In this case we observed, that matter from the annihilation layer lost its stability and propagated inward pushing matter of inner regions of the disk to accrete. Thirdly, we investigated the case where counterrotating matter inflowing from large radial distances encounters an existing corotating disk. Friction between the inflowing matter and the existing disk is found to lead to fast boundary layer accretion along the disk surfaces and to enhanced accretion in the main disk. These models are pertinent to the formation of counterrotating disks in galaxies and possibly in Active Galactic Nuclei and in X-ray pulsars in binary systems.

An investigation is made of disk accretion of matter onto a rotating star with an aligned dipole magnetic field. A new aspect of this work is that when the angular velocity of the star and disk differ substantially we argue that the $\bf B$ field linking the star and disk rapidly inflates to give regions of open field lines extending from the polar caps of the star and from the disk. The open field line region of the disk leads to the possibility of magnetically driven outflows. An analysis is made of the outflows and their back affect on the disk structure assuming an ``$\ap$" turbulent viscosity model for the disk and a magnetic diffusivity comparable to this viscosity. The outflows are found to extend over a range of radial distances inward to a distance close to $r_{to}$, which is the distance of the maximum of the angular rotation rate of the disk. We find that $r_{to}$ depends on the star's magnetic moment, the accretion rate, and the disk's magnetic diffusivity. The outflow regime is accompanied in general by a spin-up of the rotation rate of the star. When $r_{to}$ exceeds the star's corotation radius $r_{cr} = (GM/\om_*^2)^{1\ov 3}$, we argue that outflow solutions do not occur, but instead that ``magnetic braking" of the star by the disk due to field-line twisting occurs in the vicinity of $r_{cr}$. The magnetic braking solutions can give spin-up or spin-down (or no spin change) of the star depending mainly on the star's magnetic moment and the mass accretion rate. For a system with $r_{to}$ comparable to $r_{cr}$, bimodal behavior is possible where extraneous perturbations cause the system to flip between spin-up and spin-down.