A nonlinear theory of the parallel firehose and gyrothermal instabilities in a weakly collisional plasma

Weakly collisional plasmas dynamically develop pressure anisotropies with respect to the magnetic field. These anisotropies trigger plasma instabilities at scales just above the ion Larmor radius \rho_i and much below the mean free path \lambda_{mfp}. They have growth rates of a fraction of the ion cyclotron frequency - much faster than either the global dynamics or local turbulence. The instabilities dramatically modify the transport properties and, therefore, the macroscopic dynamics of the plasma. Their nonlinear evolution drives pressure anisotropies towards marginal stability, controlled by the plasma beta \beta_i. Here this nonlinear evolution is worked out for the simplest analytically tractable example - the parallel firehose instability. In the nonlinear regime, both analytical theory and the numerical solution predict secular growth of magnetic fluctuations. They develop a k^{-3} spectrum, extending from scales somewhat larger than \rho_i to the maximum scale that grows secularly with time (~t^{1/2}); the relative pressure anisotropy (\pperp-\ppar)/\ppar tends to the marginal value -2/\beta_i. The marginal state is achieved via changes in the magnetic field, not particle scattering. When a parallel ion heat flux is present, the firehose mutates into the new gyrothermal instability (GTI), which continues to exist up to firehose-stable values of pressure anisotropy, which can be positive and are limited by the heat flux. The nonlinear evolution of the GTI also features secular growth of magnetic fluctuations, but the spectrum is eventually dominated by modes around the scale ~\rho_i l_T/\lambda_{mfp}, where l_T is the scale of the parallel temperature variation. Implications for momentum and heat transport are speculated about. This study is motivated by the dynamics of galaxy cluster plasmas.

💡 Research Summary

The paper addresses how weakly collisional, high‑beta plasmas develop pressure anisotropies (p⊥≠p∥) that quickly excite micro‑scale instabilities at wavelengths just above the ion Larmor radius (ρi) but far below the mean free path (λmfp). These instabilities grow at a fraction of the ion cyclotron frequency (Ωi), far faster than the bulk dynamics or turbulent motions, and they dramatically alter the plasma’s transport properties.
The authors focus on the simplest analytically tractable case: the parallel firehose instability. Using a kinetic‑MHD framework that includes a weak collision operator, they first recover the linear threshold Δp/p∥ ≤ ‑2/βi and the fastest‑growing mode with k∥ρi ≈ 1, k⊥≈0. In the nonlinear stage they discover that the instability does not primarily scatter particles; instead, the growing magnetic fluctuations δB remodel the background field B0. This “magnetic‑field‑adjustment” mechanism drives the pressure anisotropy toward the marginal value –2/βi. The fluctuation spectrum evolves into a secular k‑3 power law extending from scales slightly larger than ρi up to a maximum scale that grows as Lmax ∝ t½. Consequently, magnetic energy is transferred to ever larger scales while the anisotropy remains pinned at the marginal firehose limit.
When a parallel ion heat flux q∥ is present, the firehose morphs into a new gyrothermal instability (GTI). The GTI can persist even when the pressure anisotropy is firehose‑stable (i.e., positive), because the heat flux supplies a free‑energy source. Its most unstable wavelength is set by the competition between the ion Larmor radius and the temperature‑gradient scale lT, yielding k∥ρi ≈ (λmfp/lT)‑½, or equivalently a physical scale ∼ ρi lT/λmfp. In the nonlinear regime the GTI also exhibits secular growth of δB, but the spectrum eventually becomes dominated by modes around this characteristic scale rather than by a broad k‑3 cascade. The final anisotropy settles at a value that can be less negative (or even positive) than –2/βi, limited by the magnitude of the heat flux.
Numerical 1‑D kinetic simulations confirm the analytical predictions: magnetic fluctuations grow secularly, the pressure anisotropy relaxes to the marginal state without significant particle scattering, and the GTI spectrum peaks at the predicted scale when q∥ is non‑zero.
The authors discuss the implications for macroscopic transport. Because the instabilities regulate Δp through magnetic‑field changes, the effective viscosity and thermal conductivity are reduced far below the classical Braginskii values. This provides a natural explanation for the suppressed heat conduction and anomalous velocity shear observed in the intracluster medium of galaxy clusters. Moreover, the predicted k‑3 magnetic spectrum could be observable as small‑scale magnetic turbulence in high‑resolution Faraday‑rotation or X‑ray spectroscopic data.
In conclusion, the paper delivers a complete nonlinear theory for the parallel firehose and its heat‑flux‑driven extension, the gyrothermal instability, in weakly collisional plasmas. It highlights how micro‑instabilities self‑regulate pressure anisotropy, generate secular magnetic fluctuations, and reshape transport coefficients, offering a robust framework for interpreting observations of high‑beta astrophysical plasmas such as those in galaxy clusters. Future work should extend the analysis to fully three‑dimensional turbulence, incorporate realistic collision operators, and compare directly with observational diagnostics.