Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

TL;DR

This paper investigates electron heat transport regulation in high-beta plasmas via whistler turbulence, confirming inverse-beta scaling and developing an effective collision operator to model turbulence effects.

Contribution

It provides a comprehensive set of PIC simulations across various parameters and introduces methods to construct an effective collision operator for whistler turbulence.

Findings

Heat flux scales inversely with electron beta.

Simulation results support the robustness of the inverse-beta scaling.

Developed methods for effective collision operator construction.

Abstract

The regulation of electron heat transport in high-beta, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as the inverse of electron beta. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of electron beta and temperature-gradient length scales under two different physical setups. The saturated heat flux in all of our simulations follows the expected inverse-electron-beta scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90-degree scattering problem in the cosmic ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.