Fully kinetic shearing-box simulations of magnetorotational turbulence in 2D and 3D. I. Pair plasmas

Abstract

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e. kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic Particle-in-Cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e. MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron-ion simulations, on which we will focus in future work.

Submitted to ApJ

References in corpus (27)

• First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring
• First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• MHD simulations of the magnetorotational instability in a shearing box with zero net flux. II. The effect of transport coefficients
• MHD simulations of the magnetorotational instability in a shearing box with zero net flux. I. The issue of convergence
• Improved Constraints on the Acceleration History of the Universe and the Properties of the Dark Energy
• Kinetic turbulence in relativistic plasma: from thermal bath to non-thermal continuum
• Magnetorotational Turbulence and Dynamo in a Collisionless Plasma
• Hybrid-Kinetic Simulations of Ion Heating in Alfvénic Turbulence
• Aspect ratio dependence in magnetorotational instability shearing box simulations
• The Radiative Efficiency and Spectra of Slowly Accreting Black Holes from Two-Temperature GRRMHD Simulations
• Sgr A* X-ray flares from non-thermal particle acceleration in a magnetically arrested disc
• Angular Momentum Transport and Particle Acceleration during Magnetorotational Instability in a Kinetic Accretion Disk
• On the viability of the shearing box approximation for numerical studies of MHD turbulence in accretion disks
• Radiation GRMHD Simulations of the Hard State of Black Hole X-ray Binaries and the Collapse of a Hot Accretion Flow
• On stochastic heating and its phase-space signatures in low-$\beta$ kinetic turbulence
• Reconnection and particle acceleration in 3D current sheet evolution in moderately-magnetized astrophysical pair plasma
• On the convergence of Magnetorotational turbulence in stratified isothermal shearing boxes
• Axisymmetric Magnetorotational Instability in Viscous Accretion Disks
• Shearingbox-implementation for the central-upwind, constraint-transport MHD-code NIRVANA
• On the iterated Crank-Nicolson for hyperbolic and parabolic equations in numerical relativity
• Nonthermal emission from the plunging region: a model for the high-energy tail of black hole X-ray binary soft states
• Energy partition between Alfvénic and compressive fluctuations in magnetorotational turbulence with near-azimuthal mean magnetic field
• The Super-Alfvénic Rotational Instability in accretion disks about black holes
• Stabilization of Magnetic Reconnection in Relativistic Current Sheet
• Stratified Simulations of Collisionless Accretion Disks
• Triggering tearing in a forming current sheet with the mirror instability