We present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with the circumbinary environment. We focus on studying the impact of the high wind magnetization on the geometrical structure of the shocked flow. This work is an extension of our earlier studies that focused on a purely hydrodynamic interaction and weak wind magnetization. We consider the collision in the axisymmetric approximation, that is, the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star. The increase of the magnetization results in the expansion of the opening cone in which the shocked pulsar wind propagates. This effect is explained in the frameworks of the conventional theory of collimation of magnetized winds. This finding has a direct implication for scenarios that involve Doppler boosting as the primary mechanism behind the GeV flares detected with the Fermi/LAT from PSR B1259-63/LS2883. The maximum enhancement of the apparent emission is determined by the ratio of $4\pi$ to the solid in which the shocked pulsar wind propagates. Our simulations suggest that this enhancement factor is decreased by the impact of the magnetic field.

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.

The very high energy (VHE) gamma-ray emission reported from a number of pulsar wind nebulae (PWNe) is naturally explained by the inverse Compton scattering of multi-TeV electrons. However, the physical dimensions of some gamma-ray-emitting PWNe significantly exceed the scales anticipated by the standard hydrodynamical paradigm of PWN formation. The most "disturbing" case in this regard is HESS J1825-137, which extends to distances $r\approx70\rm\,pc$ from the central pulsar PSR J1826-1334. If the gamma-ray emission is indeed produced inside the PWN, but not by electrons that escaped the nebula and diffuse in the interstellar medium (ISM), the formation of such an anomalously extended plerion could be realized, in a diluted environment with the hydrogen number density $n_{ism}\le10^{-2}\rm\,cm^{-3}$. In this paper, we explore an alternative scenario assuming that the pulsar responsible for the formation of the nebula initially had a very short rotation period. In this case, the sizes of both the PWN and the surrounding supernova remnant depend on the initial pulsar period, the braking index, and the ISM density. To check the feasibility of this scenario, we study the parameter space that would reproduce the size of HESS J1825-137. We show that this demand can be achieved if the braking index is small, $n\leq2$ and the pulsar birth period is short, $P_{\rm b}\simeq1\rm\,ms$. This scenario can reproduce the wind termination position, which is expected at $R_{ts}\simeq0.03\rm\,pc$, only in a dense environment with $n_{ism}\geq\rm1\,cm^{-3}$. The requirement of the dense surrounding gas is supported by the presence of molecular clouds found in the source vicinity.

We numerically investigate the generation of a magnetic field in a protostellar disc via an $\alpha \Omega$-dynamo and the resulting magnetohydrodynamic (MHD) driven outflows. We find that for small values of the dimensionless dynamo parameter $\alpha_d$ the poloidal field grows exponentially at a rate $\sigma \propto \Omega_K \sqrt{\alpha_d}$, before saturating to a value $\propto \sqrt{\alpha_d}$. The dynamo excites dipole and octupole modes, but quadrupole modes are suppressed, because of the symmetries of the seed field. Initial seed fields too weak to launch MHD outflows are found to grow sufficiently to launch winds with observationally relevant mass fluxes of order $10^{-9} M_{\odot}/\rm{yr}$ for T Tauri stars. This suggests $\alpha \Omega$-dynamos may be responsible for generating magnetic fields strong enough to launch observed outflows.

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 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.

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.

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.

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 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.

In this paper, we present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with its environment in binary system. We compare the impact of both the magnetic field and the wind anisotropy to the benchmark case of a purely hydrodynamical (HD) interaction of isotropic winds, which has been studied in detail by Bogovalov et al. (2008). We consider the interaction in axisymmetric approximation, i.e. the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star and the effects related to the pulsar orbiting are neglected. The impact of the magnetic field is studied for the case of weak magnetization (with magnetization parameter $\sigma<0.1$), which is consistent with conventional models of pulsar winds. The effects related to anisotropy in pulsar winds are modeled assuming that the kinetic energy flux in a non-magnetized pulsar wind is strongly anisotropic, with the minimum at the pulsar rotation axis and the maximum in the perpendicular direction. We show that, although both considered effects change the shape of the region occupied by the terminated pulsar wind, their impact appears to be small. In particular, for the magnetization of the pulsar wind below 0.1, the magnetic field pressure remains well below the plasma pressure in the post-shock region. Thus, in the case of interaction of a pulsar with the stellar wind environment (opposite to the case of plerions, i.e. the pulsar interaction with the interstellar medium, when the magnetic field becomes dynamically important independently on the wind magnetization) the HD approach represents a feasible approximation for numerical modelling.

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 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.

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 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.

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.

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 hydrodynamics of the interaction of pulsar and stellar winds in binary systems harboring a pulsar and its impact on the nonthermal radiation of the binary pulsar PSR B1259-63/SS2883 is discussed. The collision of an ultrarelativistic pulsar wind with a nonrelativistic stellar outflow results in significant bulk acceleration of the shocked material from the pulsar wind. Already at distances comparable to the size of the binary system, the Lorentz factor of the shocked flow can be as large as $\gamma$~4. This results in significant anisotropy of the inverse Compton radiation of accelerated electrons. Because of the Doppler boosting of the produced radiation, one should expect a variable gamma-ray signal from the system. In particular, this effect may naturally explain the reported tendency of a decrease of TeV gamma-ray flux close to the periastron. The modeling of the interaction of pulsar and stellar winds allows self-consistent calculations of adiabatic losses. Our results show that adiabatic losses dominate over the radiative losses. These results have direct impact on the orbital variability of radio, X-ray and gamma-ray signals detected from the binary pulsar PSR 1259-63/SS2883.

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.

In this paper, we present a detailed hydrodynamical study of the properties of the flow produced by the collision of a pulsar wind with the surrounding in a binary system. This work is the first attempt to simulate interaction of the ultrarelativistic flow (pulsar wind) with the nonrelativistic stellar wind. Obtained results show that the wind collision could result in the formation of an "unclosed" (at spatial scales comparable to the binary system size) pulsar wind termination shock even when the stellar wind ram pressure exceeds significantly the pulsar wind kinetical pressure. Moreover, the post-shock flow propagates in a rather narrow region, with very high bulk Lorentz factor ($\gamma\sim100$). This flow acceleration is related to adiabatical losses, which are purely hydrodynamical effects. Interestingly, in this particular case, no magnetic field is required for formation of the ultrarelativistic bulk outflow. The obtained results provide a new interpretation for the orbital variability of radio, X-ray and gamma-ray signals detected from binary pulsar system PSR 1259-63/SS2883.

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 the results of the numerical simulations of the interaction between a magnetized star and an imperfectly conducting accretion disk. To analyze the ''star-disk'' interaction we numerically investigate the MHD equations used Godunov-type high resolution numerical method. It was found that the ''star-disk'' interaction occurs with quasi-periodic reconnection of the magnetic field coronal loops and plasmoid ejections. In the case of the perfect disk conductivity the evolution of the coronal magnetic field leads to the periodic outflow of angular momentum from the disk. In the case of an imperfectly conducting disk the configuration of the magnetic field is formed such that the disk angular momentum carried by magnetic field gets balanced by angular momentum carried by matter.

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.

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 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.

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 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$.

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 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.

A review is made of recent magnetohydrodynamic (MHD) theory and simulations of origin of jets from accretion disks. Many compact astrophysical objects emit powerful, highly-collimated, oppositely directed jets. Included are the extra galactic radio jets of active galaxies and quasars, and old compact stars in binaries, and emission line jets in young stellar objects. It is widely thought that these different jets arise from rotating, conducting accretion disks threaded by an ordered magnetic field. The twisting of the magnetic field by the rotation of the disk drives the jets by magnetically extracting matter, angular momentum, and energy from the accretion disk. Two main regimes have been discussed theoretically, hydromagnetic winds which have a significant mass flux, and Poynting flux jets where the mass flux is negligible. Over the past several years, exciting new developments on models of jets have come from progress in MHD simulations which now allow the study of the origin - the acceleration and collimation - of jets from accretion disks. Simulation studies in the hydromagnetic wind regime indicate that the outflows are accelerated close to their region of origin whereas the collimation occurs at much larger distances.

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.