Theoretical astrophysics and plasma physics at RPC

Abstract:

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 3D 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.

Abstract:

Plasmas whose Coulomb-collision rates are very small may relax on shorter timescales to non-Maxwellian quasi-equilibria, which, nevertheless, have a universal form, with dependence on initial conditions retained only via an infinite set of Casimir invariants enforcing phase-volume conservation. These are distributions derived by Lynden-Bell (Mon. Not. R. Astron. Soc., vol. 136, 1967, p. 101) via a statistical-mechanical entropy-maximisation procedure, assuming perfect mixing of phase-space elements. To show that these equilibria are reached dynamically, one must derive an effective 'collisionless collision integral' for which they are fixed points - unique and inevitable provided the integral has an appropriate H-theorem. We describe how such collision integrals are derived and what assumptions are required for them to have a closed form, how to prove the H-theorems for them, and why, for a system carrying sufficiently large electric-fluctuation energy, collisionless relaxation should be fast. It is suggested that collisionless dynamics may favour maximising entropy locally in phase space before converging to global maximum-entropy states. Relaxation due to interspecies interaction is examined, leading, inter alia, to spontaneous transient generation of electron currents. The formalism also allows efficient recovery of 'true' collision integrals for both classical and quantum plasmas.