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 ultra-relativistic jet from a black-hole accretion system in, for example, a magnetically arrested disk state.

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.

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 observation. Our results showed that once magnetic reconnection takes place, a firehose-sense 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.