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

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'en 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*.

💡 Research Summary

This paper presents the first first‑principles particle‑in‑cell (PIC) investigation of the magnetorotational instability (MRI) in a collisionless plasma, focusing on local two‑dimensional (axisymmetric) shearing‑box simulations. The motivation stems from low‑luminosity accretion flows such as that onto Sgr A*, where Coulomb collision frequencies are far below dynamical rates, rendering the MHD approximation questionable. By treating electrons and ions as individual particles, the authors capture kinetic effects that are inaccessible to fluid models.

The numerical setup employs a 2‑D (x–z) shearing box with a net vertical magnetic field Bz. Several runs explore different initial Alfvén‑to‑light‑speed ratios (vA/c = 0.01, 0.03, 0.1) and ion‑to‑electron mass ratios (mi/me = 25 and 100). The grid resolution and particle number are chosen to resolve thin current sheets and reconnection layers. Periodic radial (x) and shearing‑periodic azimuthal (y) boundary conditions are applied, while the vertical direction is treated as periodic.

Key results are as follows. First, regardless of the initial vA/c, the MRI amplifies the magnetic field continuously until the Alfvén speed approaches the speed of light (vA ≈ c). This mirrors the lack of saturation of channel modes seen in axisymmetric MHD studies. Second, the rapid field growth drives a strong pressure anisotropy with p⊥ > p∥. The anisotropy stays close to the mirror‑instability threshold, and mirror modes are continuously excited, providing a self‑regulating mechanism that limits further anisotropy growth. Third, particle energization proceeds through two complementary channels. Magnetic reconnection in the thin current sheets generated by the MRI produces non‑thermal power‑law tails in both ion and electron energy distributions, indicating efficient acceleration that could power high‑energy emission. Simultaneously, the pressure anisotropy generates a viscous stress that dissipates shear energy into heat, raising the bulk temperature of the plasma. The authors demonstrate that both mechanisms operate robustly across the parameter space examined, and that the overall heating efficiency does not depend sensitively on the chosen mass ratio or box aspect ratio.

The study also discusses limitations. The 2‑D axisymmetric geometry suppresses non‑axisymmetric parasitic modes (e.g., Kelvin‑Helmholtz, Parker) and three‑dimensional turbulence, which are known to affect MRI saturation in fully kinetic simulations. Moreover, realistic astrophysical mass ratios (mi/me ≈ 1836) and larger scale separations remain computationally prohibitive. Nonetheless, the work establishes a baseline kinetic picture: in a collisionless disk, MRI can drive magnetic fields to relativistic strength, generate near‑threshold mirror anisotropy, and simultaneously heat the plasma via reconnection‑driven particle acceleration and anisotropy‑driven viscous dissipation.

Implications for observations are significant. The non‑thermal particle populations produced by reconnection could explain the broad, power‑law spectra observed from radio to gamma‑rays in low‑luminosity accretion flows, while the viscous heating may set the overall radiative efficiency. The authors propose that future three‑dimensional PIC simulations, coupled with radiative transfer calculations, will be essential to make quantitative predictions for sources such as Sgr A*. In summary, this work provides the first kinetic confirmation that the MRI remains a potent driver of magnetic amplification and plasma heating even in the extreme, collisionless regime of astrophysical accretion disks.