We describe a new theoretical and numerical framework of the magnetohydrodynamic simulation incorporated with an anisotropic pressure tensor, which can play an important role in a collisionless plasma. A classical approach to handle the anisotropy is based on the double adiabatic approximation assuming that a pressure tensor is well described only by the components parallel and perpendicular to the local magnetic field. This gyrotropic assumption, however, fails around a magnetically neutral region, where the cyclotron period may get comparable to or even longer than a dynamical time in a system, and causes a singularity in the mathematical expression. In this paper, we demonstrate that this singularity can be completely removed away by the combination of direct use of the 2nd-moment of the Vlasov equation and an ingenious gyrotropization model. Numerical tests also verify that the present model properly reduces to the standard MHD or the double adiabatic formulation in an asymptotic manner under an appropriate limit.

We study the formation of slow-mode shocks in collisionless magnetic reconnection by using one-and two-dimensional collisionless MHD codes based on the double adiabatic approximation and the Landau closure model. We bridge the gap between the Petschek-type MHD reconnection model accompanied by a pair of slow shocks and the observational evidence of the rare occasion of in-situ slow shock observations.Our results showed that once magnetic reconnection takes place, a firehose-sense (p || > p ⊥ ) pressure anisotropy arises in the downstream region, and the generated slow shocks are quite weak comparing with those in an isotropic MHD. In spite of the weakness of the shocks, however, the resultant reconnection rate is 10 − 30% higher than that in an isotropic case. This result implies that the slow shock does not necessarily play an important role in the energy conversion in the reconnection system, and is consistent with the satellite observation in the Earth's magnetosphere.PACS numbers: 94.05.-a,

This paper presents a series of stratified shearing-box simulations of collisionless accretion disks in the recently developed framework of kinetic magnetohydrodynamics (MHD), which can handle finite non-gyrotropy of a pressure tensor. Although a fully kinetic simulation predicted a more efficient angular-momentum transport in collisionless disks than in the standard MHD regime, the enhanced transport has not been observed in past kinetic MHD approaches to gyrotropic pressure anisotropy.For the purpose of investigating this missing link between the fully kinetic and MHD treatments, this paper pays attention to the role of non-gyrotropic pressure, and makes a first attempt to incorporate certain collisionless effects into disk-scale, stratified disk simulations. When the timescale of gyrotropization was longer than, or comparable to, the disk rotation frequency of the orbit, we found that the finite non-gyrotropy selectively remaining in the vicinity of current sheets contributes to suppressing magnetic reconnection in the shearing-box system. This leads to increases both in the saturated amplitude of the MHD turbulence driven by magnetorotational instabilities and in the resultant efficiency of angular-momentum transport. Our results seem favorable for fast advection of magnetic fields toward the rotation axis of a central object, which is required to launch an

A new type of instability that is expected to drive magnetohydrodynamic (MHD) turbulence from a purely toroidal magnetic field in an accretion disk is presented. It is already known that in a differentially rotating system, the uniform toroidal magnetic field is unstable due to a magnetorotational instability (MRI) under a non-axisymmetric and vertical perturbation, while it is stable under a purely vertical perturbation. Contrary to the previous study, this paper proposes an unstable mode completely confined to the equatorial plane, driven by the expansive nature of the magnetic pressure gradient force under a non-uniform toroidal field. The basic nature of this growing eigenmode, to which we give a name "magneto-gradient driven instability", is studied using linear analysis, and the corresponding nonlinear evolution is then investigated using two-dimensional ideal MHD simulations. Although a single localized magnetic field channel alone cannot provide sufficient Maxwell stress to contribute significantly to the angular momentum transport, we find that the mode coupling between neighboring toroidal fields under multiple localized magnetic field channels drastically generates a highly turbulent state and leads to the enhanced transport of angular momentum, comparable to the efficiency seen in previous studies on MRIs. This horizontally confined mode may play an important role in the saturation of an MRI through complementray growth with the toroidal MHDs and coupling with magnetic reconnection.