Abstract:The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the rate of Coulomb collisions between particles is very small, making the disk evolve essentially as a collisionless plasma. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell (PIC) plasma simulations. In this initial study we focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field until the Alfvén velocity, v_A, is comparable to the speed of light, c (independent of the initial value of v_A/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the perpendicular pressure being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission -- from the radio to the gamma-rays -- of systems such as Sgr A*.

High-Resolution Particle-In-Cell Simulations of Two-Dimensional Bernstein-Greene-Kruskal Modes

High-resolution particle-in-cell simulations of two-dimensional Bernstein–Greene–Kruskal modes

Three-Dimensional Particle-In-Cell Simulations of Two-Dimensional Bernstein-Greene-Kruskal Modes

Parallel Ion Compressibility Effects on Kinetic Ballooning Mode for Different Magnetic Shears

A "Galactic Disk"-Model for Three-Dimensional Bernstein-Greene-Kruskal Modes in a Finite Magnetic Field

Sensitivity of Kinetic Ballooning Mode Instability to Tokamak Equilibrium Implementations

Field Theory and a Structure-Preserving Geometric Particle-in-cell Algorithm for Drift Wave Instability and Turbulence

Particle-in-cell Simulation of Plasma Wave Excitation in Two-Dimensional Electron Gas

Particle-in-Cell simulation of two-dimensional electron velocity shear driven instability in relativistic domain

The role of noise in PIC and Vlasov simulations of the Buneman instability

3D cylindrical BGK model of electron phase-space holes with finite velocity and polarization drift

Geometric Electrostatic Particle-In-Cell Algorithm on Unstructured Meshes

High Frequency Gyrokinetic Particle-In-Cell Simulation: Application to Heating of Magnetically Confined Plasmas

Canonical Symplectic Particle-In-Cell Method For Long-Term Large-Scale Simulations Of The Vlasov-Maxwell Equations

A Particle Simulation of Current Sheet Instabilities under Finite Guide Field

3D Electrostatic Gyrokinetic Electron and Fully Kinetic Ion Simulation of Lower-Hybrid Drift Instability of Harris Current Sheet

Gyrokinetic Electron and Fully Kinetic Ion Simulations of Fast Magnetosonic Waves in the Magnetosphere

Gyrokinetic Particle-in-Cell simulation of the effect of compressional magnetic perturbations on the microscopic instabilities in Tokamak

Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code

Variational symplectic particle-in-cell simulation of nonlinear mode conversion from extraordinary waves to Bernstein waves

Local 2D Particle-in-cell simulations of the collisionless MRI