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.

We study the Rossby wave instability (RWI) in a thin accretion disc threaded by an initially toroidal magnetic field using the magnetohydrodynamics (MHD) code PLUTO. We find that for plasma $10^1 < \beta < 10^3$ , the growth rate and late-time density are suppressed, whereas for plasma $\beta > 10^3$, the magnetic field has negligible effect. The initially toroidal field is twisted inside the vortex and expelled at late times. This creates a radially directed field outside the vortex region, which can be understood via a simple kinematic model of a magnetic field in a rotating fluid and may be observable via polarised dust emission.

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

This paper explores the effects of large-scale magnetic fields in hot accretion flows for asymmetric configurations with respect to the equatorial plane. The solutions that we have found show that the large-scale asymmetric magnetic field can significantly affect the dynamics of the flow and also cause notable outflows in the outer parts. Previously, we treated a viscous resistive accreting disc in the presence of an odd symmetric $\textbf{B-}$ field about the equatorial plane. Now we extend our earlier work by taking into account another configuration of large-scale magnetic field which is no longer symmetric. We provide asymmetric field structures with small deviations from even and odd symmetric $\textbf{B-}$field. Our results show that the disc's dynamics and appearance become different above and below the equatorial plane. The set of solutions also predicts that even a small deviation in a symmetric field causes the disc to compress on one side and expand on the other. In some cases, our solution represents a very strong outflow from just one side of the disc. Therefore, the solution may potentially explain the origin of one-sided jets in radio galaxies.

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.

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.

We use an N-body simulation, constructed using GADGET-2, to investigate an accretion flow onto an astrophysical disk that is in the opposite sense to the disk's rotation. In order to separate dynamics intrinsic to the counter-rotating flow from the impact of the flow onto the disk, we consider an initial condition in which the counter-rotating flow is in an annular region immediately exterior the main portion of the astrophysical disk. Such counter-rotating flows are seen in systems such as NGC 4826 (known as the "Evil Eye Galaxy"). Interaction between the rotating and counter-rotating components is due to two-stream instability in the boundary region. A multi-armed spiral density wave is excited in the astrophysical disk and a density distribution with high azimuthal mode number is excited in the counter-rotating flow. Density fluctuations in the counter-rotating flow aggregate into larger clumps and some of the material in the counter-rotating flow is scattered to large radii. Accretion flow processes such as this are increasingly seen to be of importance in the evolution of multi-component galactic disks.

A non-negligible fraction of a Supermassive Black Hole's (SMBH) rest mass energy gets transported into extragalactic space by a remarkable process in jets which are incompletely understood. What are the physical processes which transport this energy? It is likely that the energy flows electromagnetically, rather than via a particle beam flux. The deduced electromagnetic fields may produce particles of energy as high as $\sim 10^{20}$ eV. The energetics of SMBH accretion disk models and the electromagnetic energy transfer imply that a SMBH should generate a $10^{18} - 10^{19}$ Amp\`eres current close to the black hole and its accretion disk. We describe the so far best observation-based estimate of the magnitude of the current flow along the axis of the jet extending from the nucleus of the active galaxy in 3C303. The current is measured to be $I \sim 10^{18}$ Amp\`eres at $\sim 40$ kpc away from the AGN. This indicates that organized current flow remains intact over multi-kpc distances. The electric current $I$ transports electromagnetic power into free space, $P = I^{2}Z$, where $Z \sim 30$ Ohms is related to the impedance of free space, and this points to the existence of cosmic electric circuit. The associated electric potential drop, $V=IZ$, is of the order of that required to generate Ultra High Energy Cosmic Rays (UHECR). We describe the analogy of electromagnetically dominated jets with transmission lines. High powered jets {\it in vacuo} can be understood by approximate analogy with a waveguide. The importance of inductance, impedance, and other laboratory electrical concepts are discussed in this context. To appear in Proc. 18th International Symposium on Very High Energy Cosmic Ray Interactions (ISVHECR2014), CERN, Switzerland

Observations of galaxies and models of accreting systems point to the occurrence of counter-rotating discs where the inner part of the disc ($r<r_0$) is co-rotating and the outer part is counter-rotating. This work analyzes the linear stability of radially separated co- and counter-rotating thin discs. The strong instability found is the supersonic Kelvin-Helmholtz instability. The growth rates are of the order of or larger than the angular rotation rate at the interface. The instability is absent if there is no vertical dependence of the perturbation. That is, the instability is essentially three-dimensional. The nonlinear evolution of the instability is predicted to lead to a mixing of the two components, strong heating of the mixed gas, and vertical expansion of the gas, and annihilation of the angular momenta of the two components. As a result the heated gas will free-fall towards the disc's center over the surface of the inner disc.

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.

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

A method is proposed for deriving the dispersion measure (related to the line-of-sight integrated free electron density) to high-redshift, powerful synchrotron radio sources by analysis of the temporal/spectral properties of the recorded signals. This would allow direct measurements of the distributed density of ionized baryons along the line-of-sights to the sources.

Radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C 303 have shown that one knot of this jet carries a {\it galactic}-scale electric current and that it is magnetically dominated. We develop the theory of magnetically dominated or Poynting-flux jets by making an analogy of a Poynting jet with a transmission line or waveguide carrying a net current and having a potential drop across it (from the jet's axis to its radius) and a definite impedance which we derive. Time-dependent but not necessarily small perturbations of a Poynting-flux jet are described by the "telegrapher's equations." These predict the propagation speed of disturbances and the effective wave impedance for forward and backward propagating wave components. A localized disturbance of a Poynting jet gives rise to localized dissipation in the jet which may explain the enhanced synchrotron radiation in the knots of the 3C 303 jet, and also in the apparently stationary knot HST-1 in the jet near the nucleus of the nearby galaxy M87. For a relativistic Poynting jet on parsec scales, the reflected voltage wave from an inductive termination or load can lead to a backward propagating wave which breaks down the magnetic insulation of the jet giving $|{\bf E}| /|{\bf B}|\geq 1$. At the threshold for breakdown, $|{\bf E}|/|{\bf B}|=1$, positive and negative particles are directly accelerated in the ${\bf E \times B}$ direction which is approximately along the jet axis. Acceleration can occur up to Lorentz factors $\sim 10^7$. This particle acceleration mechanism is distinct from that in shock waves and that in magnetic field reconnection.

The Rossby wave instability, associated with density bumps in differentially rotating discs, may arise in several different astrophysical contexts, such as galactic or protoplanetary discs. While the linear phase of the instability has been well studied, the nonlinear evolution and especially the saturation phase remain poorly understood. In this paper, we test the non-linear saturation mechanism analogous to that derived for wave-particle interaction in plasma physics. To this end we perform global numerical simulations of the evolution of the instability in a two-dimensional disc. We confirm the physical mechanism for the instability saturation and show that the maximum amplitude of vorticity can be estimated as twice the linear growth rate of the instability. We provide an empirical fitting formula for this growth rate for various parameters of the density bump. We also investigate the effects of the azimuthal mode number of the instability and the energy leakage in the spiral density waves. Finally, we show that our results can be extrapolated to 3D discs.

The Rossby wave instability (RWI) in non-self-gravitating discs can be triggered by a bump at a radius $r_0$ in the disc surface mass-density (which is proportional to the inverse potential vorticity). It gives rise to a growing non-axisymmetric perturbation [$\propto \exp(im\phi)$, $m=1,2..$] in the vicinity of $r_0$ consisting of anticyclonic vortices which may facilitate planetesimal growth in protoplanetary discs. Here, we analyze a continuum of thin disc models ranging from self-gravitating to non-selfgravitating. The key quantities determining the stability/instability are: (1) the parameters of the bump (or depression) in the disc surface density, (2) the Toomre $Q$ parameter of the disc (a non-self-gravitating disc has $Q\gg1$), and (3) the dimensionless azimuthal wavenumber of the perturbation $\bar{k}_\phi =mQh/r_0$, where $h$ is the half-thickness of the disc. For discs stable to axisymmetric perturbations ($Q>1$), the self-gravity has a significant role for $\bar{k}_\phi < \pi/2$ or $m<(\pi/2) (r_0/h)Q^{-1}$; instability may occur for a depression or groove in the surface density if $Q\lesssim 2$. For $\bar{k}_\phi > \pi/2$ the self-gravity is not important, and instability may occur at a bump in the surface density. Thus, for all mode numbers $m \ge 1$, the self-gravity is unimportant for $Q > (\pi/2)(r_0/h)$. We suggest that the self-gravity be included in simulations for cases where $Q< (r_0/h)$.

We discuss the problem of the formation of a large-scale magnetic field in the accretion disks around black holes, taking into account the non-uniform vertical structure of the disk. The high electrical conductivity of the outer layers of the disk prevents the outward diffusion of the magnetic field. This implies a stationary state with a strong magnetic field in the inner parts of the accretion disk close to the black hole, and zero radial velocity at the surface of the disk. Structure of advective accretion disks is investigated, and conditions for formation of optically thin regions in central parts of the accretion disk are found. The problem of jet collimation by magneto-torsion oscillations is considered.

Recent radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C 303 have shown that one knot of this jet has a {\it galactic}-scale electric current of $\sim 3\times 10^{18}$ Amp\`ere flowing along the jet axis (Kronberg et al. 2011). We develop the theory of relativistic Poynting-flux jets which are modeled as a transmission line carrying a DC current $I_0$, having a potential drop $V_0$, and a definite impedance ${\cal Z}_0 =90(u_z/c)\Omega$, where $u_z$ is the bulk velocity of the jet plasma. The electromagnetic energy flow in the jet is ${\cal Z}_0 I_0^2$. The observed current in 3C 303 can be used to calculate the electromagnetic energy flow in this magnetically dominated jet. Time-dependent but not necessarily small perturbations of a Poynting-flux jet - possibly triggered by a gas cloud penetrating the jet - are described by "telegrapher's equations," which predict the propagation speed of disturbances and the effective wave impedance ${\cal Z}$. The disturbance of a Poynting jet by the cloud gives rise to localized dissipation in the jet which may explain the enhanced synchrotron radiation in the knots of the 3C 303 jet, and in the apparently stationary knot HST-1 in the jet from the nucleus of the galaxy M87 (Biretta et al. 1999).

We present radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C303. From this data we derive the magnetoplasma and electrodynamic parameters of this 50 kpc long jet. For a {$\sim 2$} kpc segment of this jet we obtain for the first time a direct determination of a {\it galactic}-scale electric current ($\sim 3 \times 10^{18}$ A), and its direction $-$ {\it positive} away from the AGN. Our analysis strongly supports a model where the jet energy flow is mainly electromagnetic.

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.

In earlier works we pointed out that the disk's surface layers are non-turbulent and thus highly conducting (or non-diffusive) because the hydrodynamic and/or magnetorotational (MRI) instabilities are suppressed high in the disk where the magnetic and radiation pressures are larger than the plasma thermal pressure. Here, we calculate the vertical profiles of the {\it stationary} accretion flows (with radial and azimuthal components), and the profiles of the large-scale, magnetic field taking into account the turbulent viscosity and diffusivity and the fact that the turbulence vanishes at the surface of the disk. Also, here we require that the radial accretion speed be zero at the disk's surface and we assume that the ratio of the turbulent viscosity to the turbulent magnetic diffusivity is of order unity. Thus at the disk's surface there are three boundary conditions. As a result, for a fixed dimensionless viscosity $\alpha$-value, we find that there is a definite relation between the ratio ${\cal R}$ of the accretion power going into magnetic disk winds to the viscous power dissipation and the midplane plasma-$\beta$, which is the ratio of the plasma to magnetic pressure in the disk. For a specific disk model with ${\cal R}$ of order unity we find that the critical value required for a stationary solution is $\beta_c \approx 2.4r/(\alpha h)$, where $h$ the disk's half thickness. For weaker magnetic fields, $\beta > \beta_c$, we argue that the poloidal field will advect outward while for $\beta< \beta_c$ it will advect inward. Alternatively, if the disk wind is negligible (${\cal R} \ll 1$), there are stationary solutions with $\beta \gg \beta_c$.

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 discuss the problem of the formation of a large-scale magnetic field in the accretion disks around black holes, taking into account the nonuniform vertical structure of the disk. The high electrical conductivity of the outer layers of the disk prevents the outward diffusion of the magnetic field. This implies a stationary state with a strong magnetic field in the inner parts of the accretion disk close to the black hole, and zero radial velocity at the surface of the disk. Magnetic jet collimation is considered, when the jet radius is hold due to magneto-torsional oscillations. The range of parameters is found where jet radius is oscillating, regularly or chaotically, within restricted values.

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.

Magnetospheric accretion is an important process for a wide range of astrophysical systems, and may play a role in the formation of gas giant planets. Extending the formalism describing stellar magnetospheric accretion into the planetary regime, we demonstrate that magnetospheric processes may govern accretion onto young gas giants in the isolation phase of their development. Planets in the isolation phase have cleared out large gaps in their surrounding circumstellar disks, and settled into a quasi-static equilibrium with radii only modestly larger than their final sizes (i.e., $ r \sim 1.4 r_{\rm final}$). Magnetospheric accretion is less likely to play a role in a young gas giant's main accretion phase, when the planet's envelope is predicted to be much larger than the planet's Alfv\'en radius. For a fiducial 1 M$_J$ gas giant planet with a remnant isolation phase accretion rate of $\dot{M}_{\odot} =$ 10$^{-10} M_{\odot}{\rm yr}^{-1}=10^{-7}M_{J}{\rm yr}^{-1}$, the disk accretion will be truncated at $\sim 2.7r_J$ (with $r_J$ is Jupiter's radius) and drive the planet to rotate with a period of $\sim$7 hours. Thermal emission from planetary magnetospheric accretion will be difficult to observe; the most promising observational signatures may be non-thermal, such as gyrosynchrotron radiation that is clearly modulated at a period much shorter than the rotation period of the host star.

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.

Massive black hole (BH) mergers can result in the merger remnant receiving a "kick", of order 200 km s$^{-1}$ or more, which will cause the remnant to oscillate about the galaxy centre. Here we analyze the case where the BH oscillates through the galaxy centre perpendicular or parallel to the plane of the galaxy for a model galaxy consisting of an exponential disk, a Plummer model bulge, and an isothermal dark matter halo. For the perpendicular motion we find that there is a strong resonant forcing of the disk radial motion near but somewhat less than the "resonant radii" $r_R$ where the BH oscillation frequency is equal one-half, one-fourth, (1/6, etc.) of the radial epicyclic frequency in the plane of the disk. Near the resonant radii there can be a strong enhancement of the radial flow and disk density which can lead to shock formation. In turn the shock may trigger the formation of a ring of stars near $r_R$. As an example, for a BH mass of $10^8 M_\odot$ and a kick velocity of 150 km s$^{-1}$, we find that the resonant radii lie between 0.2 and 1 kpc. For BH motion parallel to the plane of the galaxy we find that the BH leaves behind it a supersonic wake where star formation may be triggered. The shape of the wake is calculated as well as the slow-down time of the BH. The differential rotation of the disk stretches the wake into ring-like segments.

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.

Activity of the nuclei of galaxies and stellar mass systems involving disk accretion to black holes is thought to be due to (1) a small-scale turbulent magnetic field in the disk (due to the magneto-rotational instability or MRI) which gives a large viscosity enhancing accretion, and (2) a large-scale magnetic field which gives rise to matter outflows and/or electromagnetic jets from the disk which also enhances accretion. An important problem with this picture is that the enhanced viscosity is accompanied by an enhanced magnetic diffusivity which acts to prevent the build up of a significant large-scale field. Recent work has pointed out that the disk's surface layers are non-turbulent and thus highly conducting (or non-diffusive) because the MRI is suppressed high in the disk where the magnetic and radiation pressures are larger than the thermal pressure. Here, we calculate the vertical ($z$) profiles of the stationary accretion flows (with radial and azimuthal components), and the profiles of the large-scale, magnetic field taking into account the turbulent viscosity and diffusivity due to the MRI and the fact that the turbulence vanishes at the surface of the disk. The stationary solutions we find indicate that a weak ($\beta > 1$) large-scale field does not diffuse away as suggested by earlier work.

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

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.

Jets and outflows are observed around a wide variety of accreting objects and seem to be a near-ubiquitous feature of accretion disks. Large-scale magnetic fields are thought to be necessary for jet formation in many systems, but a longstanding puzzle is that the turbulence which is responsible for inward disk accretion should be even more efficient at causing a large-scale magnetic field to diffuse outward; it just doesn't seem possible to build up a strong field in the inner disk through advection of a weak one from outside. Here, we report theoretical work which challenges this conventional wisdom and shows that the surface layer of the disk (which is nonturbulent) can easily advect magnetic fields inward. This opens up the possibility for more diverse and robust jet formation than is seen in numerical simulations which assume that the large-scale field must be generated entirely via a dynamo within the inner disk. In addition, strong disk turbulence (as inferred from observations) can naturally be explained by the magnetorotational instability feeding off a large-scale magnetic field anchored in the surface layer.

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.

A study is made of the behavior of massive black holes in disk galaxies that have received an impulsive kick from a merger or a sustained acceleration from an asymmetric jet. The motion of the gas, stars, dark matter, and massive black hole are calculated using the GADGET-2 simulation code. The massive black hole escapes the galaxy for kick velocities above about 600 km/s or accelerations above about 4*10^{-8} cm/s^2 over time-scales of the order of 10^8 yr. For smaller velocity kicks or smaller accelerations, the black hole oscillates about the center of mass with a frequency which decreases as the kick velocity or acceleration increases. The black hole displacements may give rise to observable nonaxisymmetries in the morphology and dynamics of the stellar and gaseous disk of the galaxy. In some cases the dynamical center of the galaxy is seen to be displaced towards the direction of the BH acceleration with a characteristic ``tongue--'' shaped extension of the velocity contours on the side of the galaxy opposite the acceleration.

We show that a large-scale, weak magnetic field threading a turbulent accretion disk tends to be advected inward, contrary to previous suggestions that it will be stopped by outward diffusion. The efficient inward transport is a consequence of the diffuse, magnetically-dominated surface layers of the disk, where the turbulence is suppressed and the conductivity is very high. This structure arises naturally in three-dimensional simulations of magnetorotationally unstable disks, and we demonstrate here that it can easily support inward advection and compression of a weak field. The advected field is anchored in the surface layer but penetrates the main body of the disk, where it can generate strong turbulence and produce values of alpha (i.e., the turbulent stress) large enough to match observational constraints; typical values of the vertical magnetic field merely need to reach a few percent of equipartition for this to occur. Overall, these results have important implications for models of jet formation which require strong, large-scale magnetic fields to exist over a region of the inner accretion disk.

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

We reconsider the problem of the formation of a large-scale magnetic field in the accretion disks around black holes. In contrast with previous work we take into account the nonuniform vertical structure of the disk. The high electrical conductivity of the outer layers of the disk prevents the outward diffusion of the magnetic field. This implies a stationary state with a strong magnetic field in the inner parts of the accretion disk close to the black hole.

We investigate the effect of shear in the flow of charged particle equilibria that are unstable to the Coherent Synchrotron Radiation (CSR) instability. Shear may act to quench this instability because it acts to limit the size of the region with a fixed phase relation between emitters. The results are important for the understanding of astrophysical sources of coherent radiation where shear in the flow is likely.

This work derives the linearized equations of motion, the Lagrangian density, the Hamiltonian density, and the canonical angular momentum density for general perturbations [$\propto \exp(im\phi)$ with $m=0,\pm 1,..$] of a geometrically thin self-gravitating, homentropic fluid disk including the pressure. The theory is applied to ``eccentric,'' $m=\pm 1$ perturbations of a geometrically thin Keplerian disk. We find $m=1$ modes at low frequencies relative to the Keplerian frequency. Further, it shown that these modes can have negative energy and negative angular momentum. The radial propagation of these low frequency $m=1$ modes can transport angular momentum away from the inner region of a disk and thus increase the rate of mass accretion. Depending on the radial boundary conditions there can be discrete low-frequency, negative-energy, $m=1$ modes.

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.

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.

In order to understand star formation it is important to understand the dynamics of atomic and molecular clouds in the interstellar medium (ISM). Nonlinear hydrodynamic flows are a key component to the ISM. One route by which nonlinear flows arise is the onset and evolution of interfacial instabilities. Interfacial instabilities act to modify the interface between gas components at different densities and temperatures. Such an interface may be subject to a host of instabilities, including the Rayleigh-Taylor, Kelvin-Helmholtz, and Richtmyer-Meshkov instabilities. Recently, a new density interface instability was identified. This self-gravity interfacial instability (SGI) causes any displacement of the interface to gr ow on roughly a free-fall time scale, even when the perturbation wavelength is much less than the Jeans length. In previous work, we used numerical simulations to confirm the expectations of linear theory and examine the nonlinear evolution of the SGI. We now continue our study by generalizing our initial conditions to allow the acceleration due to self-gravity to be non-zero across the interface. We also consider the behaviour of the SGI for perturbation wavelengths near the Jeans wavelength. We conclude that the action of self-gravity across a density interface may play a significant role in the ISM either by fueling the growth of new instabilities or modifying the evolution of existing instabilities.

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.

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 analyze the confining effect of cloud collapse on an expanding supernova shockfront. We solve the differential equation for the forces on the shockfront due to ram pressure, supernova energy, and gravity. We find that the expansion of the shockfront is slowed and in fact reversed by the collapsing cloud. Including radiative losses and a potential time lag between supernova explosion and cloud collapse shows that the expansion is reversed at smaller distances as compared to the non-radiative case. We also consider the case of multiple supernova explosions at the center of a collapsing cloud. For instance, if we scale our self-similar solution to a single supernova of energy 10^51 ergs occurring when a cloud of initial density 10^2 H/cm^3 has collapsed by 50%, we find that the shockfront is confined to ~15 pc in ~1 Myrs. Our calculations are pertinent to the observed unusually compact non-thermal radio emission in blue compact dwarf galaxies (BCDs). More generally, we demonstrate the potential of a collapsing cloud to confine supernovae, thereby explaining how dwarf galaxies would exist beyond their first generation of star formation.

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.

We study the stability of a collisionless, relativistic, finite-strength, cylindrical layer of charged particles in free space by solving the linearized Vlasov-Maxwell equations and compute the power of the emitted electromagnetic waves. The layer is rotating in an external magnetic field parallel to the layer. This system is of interest to understanding the high brightness temperature of pulsars which cannot be explained by an incoherent radiation mechanism. Coherent synchrotron radiation has also been observed recently in bunch compressors used in particle accelerators. We consider equilibrium layers with a `thermal' energy spread and therefore a non-zero radial thickness. The particles interact with their retarded electromagnetic self-fields. The effect of the betatron oscillations is retained. A short azimuthal wavelength instability is found which causes a modulation of the charge and current densities. The growth rate is found to be an increasing function of the azimuthal wavenumber, a decreasing function of the Lorentz factor, and proportional to the square root of the total number of electrons. We argue that the growth of the unstable perturbation saturates when the trapping frequency of electrons in the wave becomes comparable to the growth rate. Owing to this saturation we can predict the radiation spectrum for a given set of parameters. Our predicted brightness temperatures are proportional to the square of the number of particles and scale by the inverse five-third power of the azimuthal wavenumber which is in rough accord with the observed spectra of radio pulsars.

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.

Accretion disk theory was first developed as a theory with the local heat balance, where the whole energy produced by a viscous heating was emitted to the sides of the disk. One of the most important new invention of this theory was the phenomenological treatment of the turbulent viscosity, known the `` alpha'' prescription, where the (r phi) component of the stress tensor was approximated by (alpha p) with a unknown constant alpha. This prescription played the role in the accretion disk theory as well important as the mixing-length theory of convection for stellar evolution. Sources of turbulence in the accretion disk are discussed, including hydrodynamical turbulence, convection and magnetic field role. In parallel to the optically thick geometrically thin accretion disk models, a new branch of the optically thin accretion disk models was discovered, with a larger thickness for the same total luminosity. The choice between these solutions should be done of the base of a stability analysis. The ideas underlying the necessity to include advection into the accretion disk theory are presented and first models with advection are reviewed. The present status of the solution for a low-luminous optically thin accretion disk model with advection is discussed and the limits for an advection dominated accretion flows (ADAF) imposed by the presence of magnetic field are analysed. Related problems of mass ejection from accretion disks and jet formation are discussed.

We revisit the problem of magnetic field generation in accretion flows onto black holes owing to the excess radiation force on electrons. This excess force may arise from the Poynting-Robertson effect. Instead of a recent claim of the generation of dynamically important magnetic fields, we establish the validity of earlier results from 1977 which show only small magnetic fields are generated. The radiative force causes the magnetic field to initially grow linearly with time. However, this linear growth holds for only a {\it restricted} time interval which is of the order of the accretion time of the matter. The large magnetic fields recently found result from the fact that the linear growth is unrestricted. A model of the Poynting-Robertson magnetic field generation close to the horizon of a Schwarzschild black hole is solved exactly using General Relativity, and the field is also found to be dynamically insignificant. These weak magnetic fields may however be important as seed fields for dynamos.

We study the evolution of a magnetic arcade that is anchored to an accretion disk and is sheared by the differential rotation of a Keplerian disk. By including an extremely low external plasma pressure at large distances, we obtain a sequence of axisymmetric magnetostatic equilibria and show that there is a fundamental difference between field lines that are affected by the plasma pressure and those are not (i.e., force-free). Force-free fields, while being twisted by the differential rotation of the disk, expand outward at an angle of $\sim 60^\circ$ away from the rotation axis, consistent with the previous studies. These force-free field lines, however, are enclosed by the outer field lines which originate from small disk radii and come back to the disk at large radii. These outer fields experience most of the twist, and they are also affected most by the external plasma pressure. At large cylindrical radial distances, magnetic pressure and plasma pressure are comparable so that any further radial expansion of magnetic fields is prevented or slowed down greatly by this pressure. This hindrance to cylindrical radial expansion causes most of the added twist to be distributed on the ascending portion of the field lines, close to the rotation axis. Since these field lines are twisted most, the increasing ratio of the toroidal $B_{\phi}$ component to the poloidal component $B_{R,z}$ eventually results in the collimation of magnetic energy and flux around the rotation axis. We discuss the difficulty with adding a large number of twists within the limitations of the magnetostatic approximation.

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.

We have examined the morphological and dynamical HI symmetry properties of a sample of moderately inclined Sa galaxies which are classified as morphologically normal. The sample galaxies were known a prior to exhibit kinematic peculiarities ranging from warps to independent, wholly decoupled disks, and are possibly the remnants of minor mergers. We compare the asymmetry of the rotation curves to global kinematic asymmetry, and find a relationship between rotation curve asymmetry and the kinematic n=2 mode. We have also examined the kinematics of these galaxies following the discussion of Briggs (1990) and find that the warps observed in the HI disks of these galaxies deviate significantly from the simple rules for warps that commonly apply.

Observations of lopsided spiral galaxies motivated us to explore whether the rapid passage of a companion galaxy could cause them. We examine whether the center of mass of the visible matter becomes displaced from the center of mass of the dark halo during the intruder's passage, thereby causing an asymmetric response and asymmetric structure. Two dimensional $N$-body simulations indicate that this can happen. We also explore some consequences of this offset. These include the center of mass of the visible disk following a decaying orbit around the halo center of mass and the development of transient one-armed spirals that persist for up to six rotation periods. We then study the results of a variety of initial conditions based on such offsets. We report on the results of several runs in which we initially offset a disk from its halo's center of mass by an amount typical of the above interaction. In some runs the halo is free to move, while in others it is held fixed. We used three different mass distributions for the halo in these runs. We find that the disk's center of mass spiraled inward creating a variety of observed or observable phenomena including one-armed spirals, massive clumps of particles, and counter-rotating waves. The systems settle into relatively axisymmetric configurations. Whether or not the end states included a bar depended on a variety of initial conditions.

New HI synthesis data have been obtained for six face-on galaxies with the Very Large Array. These data and reanalyses of three additional data sets make up a sample of nine face-on galaxies analyzed for deviations from axisymmetry in morphology and dynamics. This sample represents a subsample of galaxies already analyzed for morphological symmetry properties in the R-band. Four quantitative measures of dynamical nonaxisymmetry are compared to one another and to the quantitative measures of morphological asymmetry in HI and R-band to investigate the relationships between nonaxisymmetric morphology and dynamics. We find no significant relationship between asymmetric morphology and most of the dynamical measures in our sample. A possible relationship is found, however, between morphology and dynamical position angle differences between approaching and receding sides of the galaxy.

In earlier work we identified a global, non-axisymmetric instability associated with the presence of an extreme in the radial profile of the key function ${\cal L}(r) \equiv (\Sigma \Omega/\kappa^2) S^{2/\Gamma}$ in a thin, inviscid, nonmagnetized accretion disk. Here, $\Sigma(r)$ is the surface mass density of the disk, $\Omega(r)$ the angular rotation rate, $S(r)$ the specific entropy, $\Gamma$ the adiabatic index, and $\kappa(r)$ the radial epicyclic frequency. The dispersion relation of the instability was shown to be similar to that of Rossby waves in planetary atmospheres. In this paper, we present the detailed linear theory of this Rossby wave instability and show that it exists for a wider range of conditions, specifically, for the case where there is a ``jump'' over some range of $r$ in $\Sigma(r)$ or in the pressure $P(r)$. We elucidate the physical mechanism of this instability and its dependence on various parameters, including the magnitude of the ``bump'' or ``jump,'' the azimuthal mode number, and the sound speed in the disk. We find large parameter range where the disk is stable to axisymmetric perturbations, but unstable to the non-axisymmetric Rossby waves. We find that growth rates of the Rossby wave instability can be high, $\sim 0.2 \Omega_{\rm K}$ for relative small ``jumps'' or ``bumps''. We discuss possible conditions which can lead to this instability and the consequences of the instability.

A theory is developed for the dynamics of eccentric perturbations $[\propto \exp(\pm i\phi)]$ of a disk galaxy residing in a spherical dark matter halo and including a spherical bulge component. The disk is represented as a large number $N$ of rings with shifted centers and with perturbed azimuthal matter distributions. Account is taken of the dynamics of the shift of the matter at the galaxy's center which may include a massive black hole. The gravitational interactions between the rings and between the rings and the center is fully accounted for, but the halo and bulge components are treated as passive gravitational field sources. Equations of motion and a Lagrangian are derived for the ring + center system, and these lead to total energy and total angular momentum constants of the motion. We study the eccentric dynamics of a disk with an exponential surface density distribution represented by a large number of rings. The inner part of the disk is found to be strongly unstable. Angular momentum of the rings is transferred outward and to the central mass if present, and a trailing one-armed spiral wave is formed in the disk. We also analyze a disk with a modified exponential density distribution where the density of the inner part of the disk is reduced. In this case we find much slower, linear growth of the eccentric motion. A trailing one-armed spiral wave forms in the disk and becomes more tightly wrapped as time increases. The motion of the central mass if present is small compared with that of the disk.

It has recently been suggested that the maximum observed quasi-periodic oscillation (QPO) frequencies, $\nu_{max}$, for several low-mass X-ray binaries, particularly 4U 1820-30, correspond to the orbital frequency at the inner-most stable orbit of the accretion disk. This would imply that the neutron stars in these systems have masses $\go 2~M_\odot$, considerably larger than any well-measured neutron star mass. We suggest that the levelling off of $\nu_{QPO}$ may be also understood in terms of a steepening magnetic field which, although possibly dipolar at the stellar surface, is altered substantially by disk accretion, and presents a ``wall'' to the accretion flow that may be outside the innermost stable orbit. General relativistic effects add to the flattening of the $\nu_{QPO}-\mdot$ relation at frequencies below the Kepler frequency at the innermost stable orbit. We offer two other possible ways to reconcile the low value of $\nu_{max}$ ($\approx 1060$ Hz for 4U 1820-30) with a moderate neutron star mass, $\approx 1.4\msun$: at sufficiently large $\mdot$, either (i) the disk terminates in a very thin boundary layer near the neutron star surface, or (ii) $\nu_{QPO}$ is not the orbital frequency right at the inner edge of the disk, but rather at a somewhat larger radius, where the emissivity of the disk peaks.

Recent papers discussing advection dominated accretion flows (ADAF) as a solution for astrophysical accretion problems should be treated with some caution because of their uncertain physical basis. The suggestions underlying ADAF involve ignoring the magnetic field reconnection in heating of the plasma flow, assuming electron heating due only to binary Coulomb collisions with ions. Here, we analyze the physical processes in optically thin accretion flows at low accretion rates including the influence of an equipartition turbulent magnetic field. For these conditions there is continuous destruction of magnetic flux by reconnection. The reconnection is expected to significantly heat the electrons which can efficiently emit magnetobremstrahlung radiation. Because of this electron emission, the radiative efficiency of the ADAF is not small. We suggest that the small luminosities of nearby galactic black holes is due to outflows rather than ADAF accretion.

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.

The lack of our knowledge on how angular momentum is transported in accretion disks around black holes has prevented us from fully understanding their high energy emissions. We briefly highlight some theoretical models, emphasizing the energy flow and electron energization processes. More questions and uncertainties are raised from a plasma physics point of view.

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.

We find a linear instability of non-axisymmetric Rossby waves in a thin non-magnetized Keplerian disk when there is a local maximum in the radial profile of a key function ${\cal L}(r) \equiv {\cal F}(r) S^{2/\Gamma}(r)$, where ${\cal F}^{-1} = \hat {\bf z}\cdot ({\bf \nabla}\times {\bf v}) /\Sigma$ is the potential vorticity, $S = P/\Sigma^\Gamma$ is the entropy, $\Sigma$ is the surface mass density, $P$ is the vertically integrated pressure, and $\Gamma$ is the adiabatic index. We consider in detail the special case where there is a local maximum in the disk entropy profile $S(r)$. This maximum acts to trap the waves in its vicinity if its height to width ratio ${\rm max}(S)/\Delta r$ is larger than a threshold value. The pressure gradient derived from this entropy variation provides the restoring force for the wave growth. We show that the trapped waves act to transport angular momentum outward. A plausible way to produce an entropy variation is when an accretion disk is starting from negligible mass and temperature, therefore negligible entropy. As mass accumulates by either tidal torquing, magnetic torquing, or Roche-lobe overflow, confinement of heat will lead to an entropy maximum at the outer boundary of the disk. Possible nonlinear developments from this instability include the formation of Rossby vortices and the formation of spiral shocks. What remains to be determined from hydrodynamic simulations is whether or not Rossby wave packets (or vortices) ``hold together'' as they propagate radially inward.

A photometric method for quantifying deviations from axisymmetry in optical images of disk galaxies is applied to a sample of 32 face-on and nearly face-on spirals. The method involves comparing the relative fluxes contained within trapezoidal sectors arranged symmetrically about the galaxy center of light, excluding the bulge and/or barred regions. Such a method has several advantages over others, especially when quantifying asymmetry in flocculent galaxies. Specifically, the averaging of large regions improves the signal-to-noise in the measurements; the method is not strongly affected by the presence of spiral arms; and it identifies the kinds of asymmetry that are likely to be dynamically important. Application of this "method of sectors" to R-band images of 32 disk galaxies indicates that about 30% of spirals show deviations from axisymmetry at the 5-sigma level.

The present work represents an attempt to understand the `rules of behavior' of observed warps in the HI disks of spiral galaxies found by Briggs (1990). In contrast with most earlier theoretical work, the present study investigates different initial value problems of a warped disk in an oblate (or prolate) halo potential, and it represents the disk warp in terms of $N$ independently tilted, self-gravitating, concentric rings. This representation gives new insight into the disk warping. A new constant of the motion of $N$ tilted rings is identified (in addition to the energy). The phenomenon of phase-locking of the lines-of-nodes of nearby rings due to self-gravity is demonstrated. We consider the influence of dynamical friction due to ring motion through the halo matter as well as friction between gaseous rings with different vertical motions due to turbulent viscosity. We first consider the dynamics of one, two, and three tilted rings of different radii in a halo potential. We go on to develop dynamical equations for $N-$rings which are most simply expressed in terms of the complex tilt angles $\Theta_j = \theta_j exp(-i\phi_j)$, where $\theta_j$ is the actual tilt angle and $\phi_j$ the line-of-nodes angle for the $j^{th}$ ring ($j=1..N$). We numerically solve the equations for $\Theta_j$ for four different types of initial conditions: (1) warp excitation by a passing satellite, (2) excitation by a sinking compact minor satellite, (3) warp evolution due toa tilted halo potential, and (4) warp evolution resulting from an initially tilted disk plane.

Advection-dominated, high-temperature, quasi-spherical accretion flow onto a compact object, recently considered by a number of authors, assume that the dissipation of turbulent energy of the flow heats the ions and that the dissipated energy is advected inward. It is suggested that the efficiency of conversion of accretion energy to radiation can be very much smaller than unity. However, it is likely that the flows have an equipartition magnetic field with the result that dissipation of magnetic energy at a rate comparable to that for the turbulence must occur by Ohmic heating. We argue that this heating occurs as a result of plasma instabilities and that the relevant instabilities are current driven in response to the strong electric fields parallel to the magnetic field. We argue further that these instabilities heat predominantly the electrons. We conclude that the efficiency of conversion of accretion energy to radiation can be much smaller than unity only for the unlikely condition that the Ohmic heating of the electrons is negligible.

Motivated by observations of disk galaxies with counter-rotating stars, we have run two-dimensional, collisionless N-body simulations of disk galaxies with significant counter-rotating components. For all our simulations the initial value of Toomre's stability parameter was Q = 1.1. The percentage of counter-rotating particles ranges from 25% to 50%. A stationary one-arm spiral wave is observed to form in each run, persisting from a few to five rotation periods, measured at the half-mass radius. In one run, the spiral wave was initially a leading arm which subsequently transformed into a trailing arm. We also observed a change in spiral direction in the run initially containing equal numbers of particles orbiting in both directions. The results of our simulations support an interpretation of the one armed waves as due to the two stream instability.

The present study of the two-stream instability in stellar disks with counter-rotating components of stars and/or gas is stimulated by recently discovered counter-rotating spiral and S0 galaxies. Strong linear two-stream instability of tightly-wrapped spiral waves is found for one and two-armed waves with the pattern angular speed of the unstable waves always intermediate between the angular speed of the co-rotating matter ($+\Omega$) and that of the counter-rotating matter ($-\Omega$). The instability arises from the interaction of positive and negative energy modes in the co- and counter-rotating components. The unstable waves are in general convective - they move in radius and radial wavenumber space - with the result that amplification of the advected wave is more important than the local growth rate. For a galaxy of co-rotating stars and counter-rotating stars of mass-fraction $\xi_* < {1\over 2}$, or of counter-rotating gas of mass-fraction $\xi_g < {1\over 2}$, the largest amplification is usually for the one-armed leading waves (with respect to the co-rotating stars). For the case of both counter-rotating stars and gas, the largest amplifications are for $\xi_*+\xi_g \approx {1\over 2}$, also for one-armed leading waves. The two-armed trailing waves usually have smaller amplifications. The growth rates and amplifications all decrease as the velocity spreads of the stars and/or gas increase. It is suggested that the spiral waves can provide an effective viscosity for the gas causing its accretion.

We consider accretion disks consisting of counter-rotating gaseous components with an intervening shear layer. Configurations of this type may arise from the accretion of newly supplied counter-rotating gas onto an existing co-rotating gas disk. For simplicity we consider the case where the gas well above the disk midplane is rotating with angular rate $+\Omega$ and that well below has the same properties but is rotating with rate $-\Omega$. Using the Shakura-Sunyaev alpha turbulence model, we find self-similar solutions where a thin (relative to the full disk thickness) equatorial layer accretes very rapidly, essentially at free-fall speed. As a result the accretion speed is much larger than it would be for an alpha disk rotating in one direction. Counter-rotating accretion disks may be a transient stage in the formation of counter-rotating galaxies and in the accretion of matter onto compact objects.

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.