PIC Simulations of Velocity-Space Instabilities in a Decreasing Magnetic Field: Viscosity and Thermal Conduction

Abstract

We use particle-in-cell (PIC) simulations of a collisionless, electron-ion plasma with a decreasing background magnetic field, $B$, to study the effect of velocity-space instabilities on the viscous heating and thermal conduction of the plasma. If $B$ decreases, the adiabatic invariance of the magnetic moment gives rise to pressure anisotropies with $p_{||,j} > p_{\perp,j}$ ($p_{||,j}$ and $p_{\perp,j}$ represent the pressure of species $j$ ($=i$ or $e$) parallel and perpendicular to the magnetic field). Linear theory indicates that, for sufficiently large anisotropies, different velocity-space instabilities can be triggered. These instabilities, which grow on scales comparable to the electron and ion Larmor radii, in principle have the ability to pitch-angle scatter the particles, limiting the growth of the anisotropies. Our PIC simulations focus on the nonlinear, saturated regime of the instabilities. This is done through the permanent decrease of the magnetic field by an imposed shear in the plasma. Our results show that, in the regime $2 \lesssim β_j \lesssim 20$ ($β_j \equiv 8πp_j/B^2$), the saturated ion and electron pressure anisotropies are controlled by the combined effect of the oblique ion firehose (OIF) and the fast magnetosonic/whistler (FM/W) instabilities. These instabilities grow preferentially on the ion Larmor radius scale, and make the ion and electron pressure anisotropies nearly equal: $Δp_e/p_{||,e} \approx Δp_i/p_{||,i}$ (where $Δp_j=p_{\perp,j} - p_{||,j}$). We also quantify the thermal conduction of the plasma by directly calculating the mean free path of electrons along the mean magnetic field, which we find strongly depends on whether $B$ decreases or increases. Our results can be applied in studies of low collisionality plasmas such as the solar wind, the intracluster medium, and some accretion disks around black holes.

Published in the ApJ