We investigate the stability of test-particle equilibrium orbits in axisymmetric, but otherwise arbitrary, gravitational and electromagnetic fields. We extend previous studies of this problem to include a toroidal magnetic field. We find that, even though the toroidal magnetic field does not alter the location of the circular orbits, it enters the problem as a gyroscopic force with the potential to provide gyroscopic stability. This is in essence similar to the situation encountered in the reduced three-body problem where rotation enables stability around the local maxima of the effective potential. Nevertheless, we show that gyroscopic stabilization by a toroidal magnetic field is impossible for axisymmetric force fields in source-free regions because in this case the effective potential does not possess any local maxima. As an example of an axisymmetric force field with sources, we consider the classical problem of a rotating, aligned magnetosphere. By analyzing the dynamics of halo and equatorial particle orbits we conclude that axisymmetric toroidal fields that are antisymmetric about the equator are unable to provide gyroscopic stabilization. On the other hand, a toroidal magnetic field that does not vanish at the equator can provide gyroscopic stabilization for positively charged particles in prograde equatorial orbits.

Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear term, and we find that we do not need an advective gauge as solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited for non-magnetic problems. Benchmarked results of (hydrodynamical) planet-disk interaction and the of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficiently in straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community.

We derive the conductivity tensor for axisymmetric perturbations of a hot, collisionless, and charge-neutral plasma in the shearing sheet approximation. Our results generalize the well-known linear Vlasov theory for uniform plasmas to differentially rotating plasmas and can be used for wide range of kinetic stability calculations. We apply these results to the linear theory of the magneto-rotational instability (MRI) in collisionless plasmas. We show analytically and numerically how the general kinetic theory results derived here reduce in appropriate limits to previous results in the literature, including the low frequency guiding center (or "kinetic MHD") approximation, Hall MHD, and the gyro-viscous approximation. We revisit the cold plasma model of the MRI and show that, contrary to previous results, an initially unmagnetized collisionless plasma is linearly stable to axisymmetric perturbations in the cold plasma approximation. In addition to their application to astrophysical plasmas, our results provide a useful framework for assessing the linear stability of differentially rotating plasmas in laboratory experiments.

We study the linear stability of weakly magnetized differentially rotating plasmas in both collisionless kinetic theory and Braginskii's theory of collisional, magnetized plasmas. We focus on the very weakly magnetized limit that is important for understanding how astrophysical magnetic fields originate and are amplified at high redshift. We show that the single instability of fluid theory - the magnetorotational instability mediated by magnetic tension - is replaced by two distinct instabilities, one associated with ions and one with electrons. Each of these has a different way of tapping into the free energy of differential rotation. The ion instability is driven by viscous transport of momentum across magnetic field lines due to a finite ion cyclotron frequency (gyroviscosity); the fastest growing modes have wavelengths significantly longer than MHD and Hall MHD predictions. The electron instability is a whistler mode driven unstable by the temperature anisotropy generated by differential rotation; the growth time can be orders of magnitude shorter than the rotation period. The electron instability is an example of a broader class of instabilities that tap into the free energy of differential rotation or shear via the temperature anisotropy they generate. We briefly discuss the application of our results to the stability of planar shear flows and show that such flows are linearly overstable in the presence of fluid gyroviscosity. We also briefly describe the implications of our results for magnetic field amplification in the virialized halos of high redshift galaxies.

We develop an analytic theory to describe spiral density waves propagating in a shearing disc in the weakly nonlinear regime. Such waves are generically found to be excited in simulations of turbulent accretion disks, in particular if said turbulence arises from the magneto-rotational instability (MRI). We derive a modified Burgers equation governing their dynamics, which includes the effects of nonlinear steepening, dispersion, and a bulk viscosity to support shocks. We solve this equation approximately to obtain nonlinear sawtooth solutions that are asymptotically valid at late times. In this limit, the presence of shocks is found to cause the wave amplitude to decrease with time as 1/t^2. The validity of the analytic description is confirmed by direct numerical solution of the full nonlinear equations of motion. The asymptotic forms of the wave profiles of the state variables are also found to occur in MRI simulations indicating that dissipation due to shocks plays a significant role apart from any effects arising from direct coupling to the turbulence.

We study the evolution of a protoplanet of a few earth masses embedded in a protoplanetary disc. If we assume that the atmosphere of the protoplanet, i.e. the volume of gas in hydrostatic equilibrium bound to the core, has a surface radius smaller than the Roche lobe radius, we show that it expands as it accretes both planetesimals and gas at a fixed rate from the nebula until it fills in the Roche lobe. The evolution occurs on a timescale shorter than the formation or migration timescales. Therefore, we conclude that protoplanets of a few earth masses have an atmosphere that extends to the Roche lobe surface, where it matches to the nebula. This is true even when the Bondi radius is smaller than the Roche lobe radius. This is in contrast to the commonly used models in which the static atmosphere extends up to the Bondi radius and is surrounded by a cold accretion flow. As a result, any calculation of the tidal torque exerted by the disc onto the protoplanet should exclude the material present in the Roche lobe, since it is bound to the protoplanet.

A rigorous theory for the generation of a large-scale magnetic field by random non-helically forced motions of a conducting fluid combined with a linear shear is presented in the analytically tractable limit of low Rm and weak shear. The dynamo is kinematic and due to fluctuations in the net (volume-averaged) electromotive force. This is a minimal proof-of-concept quasilinear calculation aiming to put the shear dynamo, a new effect recently found in numerical experiments, on a firm theoretical footing. Numerically observed scalings of the wavenumber and growth rate of the fastest growing mode, previously not understood, are derived analytically. The simplicity of the model suggests that shear dynamo action may be a generic property of sheared magnetohydrodynamic turbulence.

We study and elucidate the mechanism of spiral density wave excitation in a differentially rotating turbulent flow. We formulate a set of wave equations with sources that are only non-zero in the presence of turbulent fluctuations. We solve these in a shearing box domain using a WKBJ method. It is found that, for a particular azimuthal wave length, the wave excitation occurs through a sequence of regularly spaced swings during which the wave changes from leading to trailing form. This is a generic process that is expected to occur in shearing discs with turbulence. Trailing waves of equal amplitude propagating in opposite directions are produced, both of which produce an outward angular momentum flux that we give expressions for as functions of the disc parameters and azimuthal wave length. By solving the wave amplitude equations numerically we justify the WKBJ approach for a Keplerian rotation law for all parameter regimes of interest. In order to quantify the wave excitation completely the important wave source terms need to be specified. Assuming conditions of weak nonlinearity, these can be identified and are associated with a quantity related to the potential vorticity, being the only survivors in the linear regime. Under the additional assumption that the source has a flat power spectrum at long azimuthal wave lengths, the optimal azimuthal wave length produced is found to be determined solely by the WKBJ response and is estimated to be 2 pi H, with H being the nominal disc scale height.

We present fully three-dimensional local MRI simulations with the object of studying the excitation of non-axisymmetric spiral density waves that are observed to always be present in such simulations. They are potentially important for affecting protoplanetary migration through the action of associated stochastic gravitational forces and producing residual transport in MHD inactive regions. The simulations we perform are with zero net flux and produce mean activity levels of alpha ~ 0.005. We reveal the nature of the mechanism responsible for the excitation of these waves by determining the time dependent evolution of the Fourier transforms of the participating state variables. The dominant waves are found to have no vertical structure and to be excited during periodically repeating swings in which they change from leading to trailing. The initial phase of the evolution of such a swing is found to be in excellent agreement with that expected from the WKBJ theory developed in a preceding paper by Heinemann & Papaloizou. We demonstrate that the important source terms causing excitation of the waves are related to a quantity that reduces to the potential vorticity for small perturbations from the background state with no vertical dependence. We find that the root mean square density fluctuations associated with the waves are positively correlated with both this quantity and the general level of hydromagnetic turbulence. The mean angular momentum transport associated with spiral density waves generated in our simulations is estimated to be a significant fraction of that associated with the turbulent Reynolds stress.

Numerical simulations of forced turbulence in elongated shearing boxes are carried out to demonstrate that a nonhelical turbulence in conjunction with a linear shear can give rise to a mean-field dynamo. Exponential growth of magnetic field at scales larger than the outer (forcing) scale of the turbulence is found. Over a range of values of the shearing rate S spanning approximately two orders of magnitude, the growth rate of the magnetic field is proportional to the imposed shear, gamma ~ S, while the characteristic spatial scale of the field is l_b ~ S^(-1/2). The effect is quite general: earlier results for the nonrotating case by Yousef et al. 2008 (PRL 100, 184501) are extended to shearing boxes with Keplerian rotation; it is also shown that the shear dynamo mechanism operates both below and above the threshold for the fluctuation dynamo. The apparently generic nature of the shear dynamo effect makes it an attractive object of study for the purpose of understanding the generation of magnetic fields in astrophysical systems.

The feasibility of a mean-field dynamo in nonhelical turbulence with superimposed linear shear is studied numerically in elongated shearing boxes. Exponential growth of magnetic field at scales much larger than the outer scale of the turbulence is found. The charateristic scale of the field is l_B ~ S^{-1/2} and growth rate is gamma ~ S, where S is the shearing rate. This newly discovered shear dynamo effect potentially represents a very generic mechanism for generating large-scale magnetic fields in a broad class of astrophysical systems with spatially coherent mean flows.

In turbulent high-beta astrophysical plasmas (exemplified by the galaxy cluster plasmas), pressure-anisotropy-driven firehose and mirror fluctuations grow nonlinearly to large amplitudes, dB/B ~ 1, on a timescale comparable to the turnover time of the turbulent motions. The principle of their nonlinear evolution is to generate secularly growing small-scale magnetic fluctuations that on average cancel the temporal change in the large-scale magnetic field responsible for the pressure anisotropies. The presence of small-scale magnetic fluctuations may dramatically affect the transport properties and, thereby, the large-scale dynamics of the high-beta astrophysical plasmas.

We discuss a numerical 3D radiation-MHD simulation of penumbral fine structure in a small sunspot. This simulation shows the development of short filamentary structures with horizontal flows, similar to observed Evershed flows, and an inward propagation of these structures at a speed compatible with observations. We conclude that the Evershed flow represents the horizontal flow component of overturning convection in gaps with strongly reduced field strength. The top of the flow is always directed outward--away from the umbra-- because of the broken symmetry due to the inclined magnetic field. Upflows occur in the inner parts of the gaps and most of the gas turns over radially (outwards and sideways), and descends back down again. The ascending, cooling and overturning flow tends to bend magnetic field lines down, forcing a weakening of the field that makes it easier for gas located in an adjacent layer--further in--to initiate a similar sequence of motion, aided by lateral heating, thus causing the inward propagation of the filament.

We study the influence of the choice of transport coefficients (viscosity and resistivity) on MHD turbulence driven by the magnetorotational instability (MRI) in accretion disks. We follow the methodology described in paper I: we adopt an unstratified shearing box model and focus on the case where the net vertical magnetic flux threading the box vanishes. For the most part we use the finite difference code ZEUS, including explicit transport coefficients in the calculations. However, we also compare our results with those obtained using other algorithms (NIRVANA, the PENCIL code and a spectral code) to demonstrate both the convergence of our results and their independence of the numerical scheme. We find that small scale dissipation affects the saturated state of MHD turbulence. In agreement with recent similar numerical simulations done in the presence of a net vertical magnetic flux, we find that turbulent activity (measured by the rate of angular momentum transport) is an increasing function of the magnetic Prandtl number Pm for all values of the Reynolds number Re that we investigated. We also found that turbulence disappears when the Prandtl number falls below a critical value Pm_c that is apparently a decreasing function of Re. For the limited region of parameter space that can be probed with current computational resources, we always obtained Pm_c>1. We conclude that the magnitudes of the transport coefficients are important in determining the properties of MHD turbulence in numerical simulations in the shearing box with zero net flux, at least for Reynolds numbers and magnetic Prandtl numbers that are such that transport is not dominated by numerical effects and thus can be probed using current computational resources.

We present results of numerical 3D MHD simulations with radiative energy transfer of fine structure in a small sunspot of about 4 Mm width. The simulations show the development of filamentary structures and flow patterns that are, except for the lengths of the filaments, very similar to those observed. The filamentary structures consist of gaps with reduced field strength relative to their surroundings. Calculated synthetic images show dark cores like those seen in the observations; the dark cores are the result of a locally elevated $\tau=1$ surface. The magnetic field in these cores is weaker and more horizontal than for adjacent brighter structures, and the core support a systematic outflow. Movies show migration of the dark-cored structures towards the umbra, and fragments of magnetic flux that are carried away from the spot by a large scale `moat flow'. We conclude that the simulations are in qualitative agreement with observed penumbra filamentary structures, Evershed flows and moving magnetic features.

An efficient algorithm for calculating radiative transfer on massively parallel computers using domain decomposition is presented. The integral formulation of the transfer equation is used to divide the problem into a local but compute-intensive part for calculating the intensity and optical depth integrals, and a nonlocal part for communicating the intensity between adjacent processors. The waiting time of idle processors during the nonlocal communication part does not have a severe impact on the scaling. The wall clock time thus scales nearly linearly with the inverse number of processors.