Search

We study the evolution of eccentricity and inclination of massive planets in low-density cavities of protoplanetary discs using three-dimensional (3D) simulations. When the planet's orbit is aligned with the equatorial plane of the disc, the eccentricity increases to high values of 0.7-0.9 due to the resonant interaction with the inner parts of the disc. For planets on inclined orbits, the eccentricity increases due to the Kozai-Lidov mechanism, where the disc acts as an external massive body that perturbs the planet's orbit. At small inclination angles, < 30 degrees, the resonant interaction with the inner disc strongly contributes to the eccentricity growth, while at larger angles, eccentricity growth is mainly due to the Kozai-Lidov mechanism. We conclude that planets inside low-density cavities tend to acquire high eccentricity if favorable conditions give sufficient time for growth. The final value of the planet's eccentricity, after the disc dispersal depends on the planet's mass and properties of the cavity and protoplanetary disc., Comment: 19 pages, 17 figures, submitted to MNRAS

We present results of global three-dimensional (3D) magnetohydrodynamic (MHD) simulations of accretion onto magnetized stars where both the magnetic and rotational axes of the star are tilted about the rotational axis of the disc. We observed that initially the inner parts of the disc are warped, tilted, and process due to the magnetic interaction between the magnetosphere and the disc. Later, larger tilted discs form with the size increasing with the magnetic moment of the star. The normal vector to the discs are tilted at different angles, from 5-10 degrees up to 30-40 degrees. Small tilts may result from the winding of the magnetic field lines about the rotational axis of the star and the action of the magnetic force which tends to align the disc. Another possible explanation is the magnetic Bardeen-Petterson effect in which the disc settles in the equatorial plane of the star due to precessional and viscous torques in the disc. Tilted discs slowly precess with the time scale of the order of 50 Keplerian periods at the reference radius (approx. 3 stellar radii). Our results can be applied to different types of stars where signs of tilted discs and/or slow precession have been observed., Comment: 14 pages, 13 figures, MNRAS, accepted. The new version of the paper is more complete compared with the earlier version

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., Comment: 17 pages, 14 figures, accepted by MNRAS

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., Comment: 11 pages, 4 figures, accepted by MNRAS

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., Comment: 11 pages, 12 figures, submitted to MNRAS

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., Comment: 16 pages, 12 figures

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., Comment: 23 pages, 17 figures, submitted to New Astronomy

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., Comment: 17 pages, 14 figures, published in MNRAS

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., Comment: 11 pages, 4 figures

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., Comment: 11 pages, 8 figures, conference proceedings: "Physics at the Magnetospheric Boundary", Geneva, Switzerland, 25-28 June, 2013

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., Comment: 14 pages, 11 figures, submitted to MNRAS

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., Comment: 4 pages, 4 figures, conference proceedings: "Physics at the Magnetospheric Boundary", Geneva, Switzerland, 25-28 June, 2013

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., Comment: 14 pages, 8 figures. arXiv admin note: substantial text overlap with arXiv:0907.3394

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., Comment: 3 pages

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., Comment: Accepted for publication in MNRAS

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)$., Comment: 5 pages, 5 figures

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., Comment: 8 pages, 6 figures

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., Comment: 14 pages, 17 figures

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., Comment: 26 pages, 25 figures

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., Comment: 6 pages, 4 figures

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)., Comment: 6 pages, 1 figure. To appear in Proc. of Xth International Conference on Gravitation, Astrophysics, and Cosmology, Qui-Nhon, Vietnam, Dec. 17-22, 2011

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., Comment: 15 pages, 13 figures, accepted by MNRAS

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., Comment: 9 pages, 9 figures, accepted for publication in MNRAS

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$., Comment: 5 pages, 3 figures

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., Comment: 7 pages, 5 figures, proceedings of the 25th Texas Symposium on Relativistic Astrophysics, held in Heidelberg in December 6-10, 2010, submitted to Proceedings of Science (POS)

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., Comment: 6 pages, 4 figures

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., Comment: 10 pages, 11 figures

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., Comment: 7 pages, 7 figures

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., Comment: 8 pages, 7 figures

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., Comment: 29 pages and 29 figures. This version has a major expansion after comments by a referee. The 1-st version is correct but mainly describes the conical wind. This version describes in greater detail both the conical winds and the propeller regime. Accepted to the MNRAS

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., Comment: 7 pages, 5 figures

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$., Comment: 10 pages, 7 figures

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., Comment: 6 pages, 5 figures. Published in: Triggering Relativistic Jets (Eds. William H. Lee & Enrico Ramirez-Ruiz) Revista Mexicana de Astronomia y Astrofisica (Serie de Conferencias) Vol. 27, Contents of Supplementary CD, p.240; 2007RMxAC..27..240L

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., Comment: 7 pages, 9 figures; Published in: Triggering Relativistic Jets (Eds. William H. Lee & Enrico Ramirez-Ruiz) Revista Mexicana de Astronomia y Astrofisica (Serie de Conferencias) Vol. 27, Contents of Supplementary CD, p.245

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., Comment: 14 pages, 14 figures, submitted to the MNRAS. See animations at http://www.astro.cornell.edu/~romanova/conical.htm ; http://www.astro.cornell.edu/~romanova/propeller.htm ; http://www.astro.cornell.edu/us-rus/disk_prop.htm

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., Comment: Invited review, conference proceedings of the meeting "Protostellar Jets in Context", 7-12 July 2008, island of Rhodes, Greece; editors: profs. Tom Ray and Kanaris Tsinganos; 10 pages, 10 figures, see animations at http://www.astro.cornell.edu/~romanova/conical.htm and http://www.astro.cornell.edu/~romanova/propeller.htm

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., Comment: 8 pages, 6 figures

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., Comment: 12 pages, 10 figures, Star-Disk Interaction in Young Stars, Proceedings of the International Astronomical Union, IAU Symposium, Volume 243. See animations at http://astro.cornell.edu/~romanova/projects.htm and at http://astro.cornell.edu/us-rus/

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., Comment: 13 pages, 11 figures, submitted to the MNRAS after revision, see animations at http://www.astro.cornell.edu/us-rus/shock.htm

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., Comment: 13 pages, 21 figures, accepted for publication in MNRAS

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., Comment: 6 pages, 2 figures

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., Comment: 5 pages, 2 figures

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., Comment: 24 pages, 8 figures

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., Comment: 11 pages, 10 figures. Accepted for publication in MNRAS

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., Comment: 6 figures, accepted for publication in Astrophysics and Space Science

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., Comment: Published in Proceedings of the 13th Young Scientists' Conference on Astronomy and Space Physics, held in Kyiv, Ukraine, April 25-29, 2006, Eds.: Golovin, A.; Ivashchenko, G.; Simon, A

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., Comment: 10 pages, 4 figures, accepted by MNRAS. See version with better resolution figures and animation at http://astrosun2.astro.cornell.edu/us-rus/propeller.htm