---
title: Bi-Maxwellian Plasma Dynamics
url: https://www.emergentmind.com/topics/bi-maxwellian-plasma
type: topic
---

# Bi-Maxwellian Plasma Dynamics

A bi-Maxwellian plasma is a kinetic system in which the particle velocity distribution is Maxwellian along the direction of the ambient magnetic field but possesses a generally different (and independent) temperature in the perpendicular direction. This results in pressure anisotropy and fundamentally alters the wave modes, closure relations, stability thresholds, and transport processes, compared to isotropic Maxwellian plasmas. The bi-Maxwellian description underpins several multi-moment closures and serves as the starting point for the kinetic analysis of anisotropic wave propagation, plasma heating, and temperature-driven instabilities.

## 1. The Bi-Maxwellian Distribution and Pressure Anisotropy

The velocity-space distribution for a plasma species $s$ in a uniform magnetic field direction $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$ is
\[
f_s(\mathbf{v}) = \frac{n_s}{\pi^{3/2} v_{s\perp}^2 v_{s\parallel}} \exp\left(-\frac{v_\perp^2}{v_{s\perp}^2} - \frac{v_\parallel^2}{v_{s\parallel}^2}\right),
\]
where $v_\parallel = \mathbf{v} \cdot \mathbf{b}$ and $v_\perp^2 = |\mathbf{v}|^2 - v_\parallel^2$ define the velocities parallel and perpendicular to the field, and $v_{s\parallel}, v_{s\perp}$ are the thermal speeds. The perpendicular and parallel kinetic pressures, defined as moments of the distribution, are
\[
p_{s\parallel} = m_s \int v_\parallel^2 f_s \,\mathrm{d}^3v, \qquad
p_{s\perp} = \frac{m_s}{2} \int v_\perp^2 f_s \,\mathrm{d}^3v,
\]
yielding the pressure tensor:
\[
\mathbf{P}_s = p_{s\perp}\,\mathbf{I} + (p_{s\parallel} - p_{s\perp})\,\mathbf{b}\mathbf{b}.
\]
The defining feature is the distinction $p_{s\parallel} \neq p_{s\perp}$, leading to pressure anisotropy, encapsulated by the anisotropy parameter $A_s = T_{s\perp} / T_{s\parallel}$ [1902.06816, 1205.6896].

## 2. Linear Modes, Wave Propagation, and Dielectric Response

The anisotropy inherent to bi-Maxwellian plasmas modifies the electromagnetic dielectric tensor. The relevant permittivity tensor $\varepsilon_{ij}$, as derived from the linearized Vlasov–Maxwell equations, contains explicit dependency on the $T_{\perp} / T_{\parallel}$ ratio through Bessel-function sums and plasma dispersion functions. Key effects include:

- **R/L (Right/Left) waves (parallel propagation):** Anisotropy enters as a correction $\propto (A_e - 1) \beta_{\parallel e}$ to the phase velocities and modifies resonance conditions.
- **Whistler (low-frequency R) branch:** For $A_e > 1$ the whistler phase speed is increased; for $A_e < 1$ it is decreased.
- **O (ordinary) mode:** For perpendicular propagation, anisotropy corrections arise at finite-Larmor-radius orders, especially via terms like $(A_\alpha^{-1} - 1) k_\perp^2 v_{th\perp\alpha}^2/(2\Omega_{0\alpha}^2)$.
- **X and Bernstein modes:** To leading order, these remain unaffected by anisotropy, as the explicit combination terms cancel [1205.6896].

Modes propagating obliquely to the magnetic field, such as kinetic Alfvén waves or the generalized fast mode, exhibit two separate acoustic corrections: one due to parallel anisotropy ("parallel acoustic effect") and one due to FLR corrections.

## 3. Instabilities: Firehose, Mirror, Whistler/Weibel

Pressure anisotropy in bi-Maxwellian plasmas enables a spectrum of distinct microinstabilities:

- **Firehose instability:** Triggered when $p_{\parallel} > p_{\perp} + B^2/\mu_0$. In dimensionless terms, instability requires $\beta_{\parallel} - \beta_{\perp} > 2$ or, at the kinetic level, $A < 1 - 2/\beta_{\parallel}$. This produces negative effective magnetic tension and destabilizes parallel modes [1902.06816, 1205.6896].

- **Mirror instability:** Occurs for $p_{\perp}/p_{\parallel} > 1 + 1/\beta_{\parallel}$. It preferentially drives nonpropagating oblique modes [1902.06816].

- **Whistler (Weibel-like) instability:** For $A = T_\perp / T_\parallel > 1$, non-resonant and resonant branches both support instabilities; field-free cases reduce to the classic Weibel instability ($\gamma \propto k_\parallel v_{th\parallel} \sqrt{A-1}$). Ambient magnetic field and relativistic corrections generally suppress growth. The instability is quenched entirely in the ultra-relativistic regime ($T_\perp \rightarrow m c^2 \Rightarrow \gamma \rightarrow 0$) [1304.2661].

- **Fast-mode instability:** For highly oblique propagation, the magnetosonic term can suppress firehose-type growth, and the instability thresholds become strongly dependent on anisotropy and propagation angle [1205.6896].

A summary of instability criteria is provided below:

| Instability      | Threshold Condition           | Physical Driver                     |
|------------------|-----------------------------|-------------------------------------|
| Firehose         | $p_\parallel - p_\perp > B^2/\mu_0$     | $T_\parallel > T_\perp$          |
| Mirror           | $p_\perp/p_\parallel > 1 + 1/\beta_\parallel$ |  $T_\perp > T_\parallel$           |
| Weibel/Whistler  | $A \equiv T_\perp/T_\parallel > 1$   |  $T_\perp > T_\parallel$           |

## 4. Multi-Moment Fluid Closures: Six-Moment and Higher

To achieve tractable descriptions while retaining key kinetic anisotropy effects, the six-moment closure assumes a bi-Maxwellian phase space, advancing equations for parallel and perpendicular pressures independently. The six-moment system for each species $s$ comprises:

- Mass continuity:  
  $\partial_t \rho_s + \nabla\cdot(\rho_s \mathbf{u}_s) = 0$
- Momentum (with full pressure tensor):  
  $\partial_t (\rho_s \mathbf{u}_s) + \nabla\cdot\left[ \rho_s \mathbf{u}_s\mathbf{u}_s + p_{s\perp} \mathbf{I} + (p_{s\parallel} - p_{s\perp}) \mathbf{b}\mathbf{b} \right] = \frac{q_s}{m_s}\rho_s (\mathbf{E} + \mathbf{u}_s\times\mathbf{B})$
- Separate parallel and perpendicular pressure evolution equations with coupling to velocity shear and expansion/compression [1902.06816].

This model is computationally less intensive than a full ten-moment treatment (which advects all six unique components of the pressure tensor), while providing significantly greater fidelity than isotropic five-moment MHD [1902.06816].

The moment-hierarchy formalism is employed in two-fluid models of the solar wind, capturing both bi-Maxwellian protons and kappa-Maxwellian electrons; the proton moments are closed at the level of pressure anisotropy and associated heat flux components [1908.09198].

## 5. Kinetic Effects in Relativistic and Non-Relativistic Regimes

In relativistic contexts (e.g., gamma-ray burst sources, relativistic jets), the semi-relativistic generalization of the bi-Maxwellian distribution becomes essential. The equilibrium is
\[
f_0(p_\perp, p_\parallel) = \frac{1}{(2\pi m T_\parallel)^{1/2}(2\pi m T_\perp)}(1 + p_\perp^2/m^2 c^2)^{-1/2} \exp\left[-\sqrt{m^2 c^2 + p_\perp^2}/T_\perp - p_\parallel^2/(2m T_\parallel)\right],
\]
reducing to the non-relativistic bi-Maxwellian in the $T_\perp \ll m c^2$ limit. Relativistic corrections generally suppress growth rates and instability thresholds, quenching anisotropy-driven instabilities as $T_\perp$ approaches $m c^2$ [1304.2661].

## 6. Applications and Observational Relevance

Bi-Maxwellian models are standard in solar wind, magnetospheric, and laboratory plasmas where field-aligned and perpendicular heating or acceleration are decoupled. In fast solar wind modeling, the bi-Maxwellian proton component yields evolution equations matching observed power-law profiles for density, temperature, and heat flux over $0.3-1$ AU, provided the kinetic-closure terms and instability-modified collision rates are retained. The parallel and perpendicular temperature components cross over as a function of heliocentric distance, compatible with in situ measurements; instability thresholds regulate the attainable anisotropy and shape turbulence dissipation and heating [1908.09198].

The six-moment closure is robust across a broad parameter space, quantitatively capturing wave propagation, linear instability thresholds, and anisotropic shock structure, barring regions of strong non-gyrotropic pressure, such as the electron diffusion region in magnetic reconnection [1902.06816].

## 7. Numerical and Theoretical Approaches

Numerical implementations exploit implicit schemes for stiffness control and explicit second-order fluxes for advection. Maxwell’s equations with divergence cleaning (via hyperbolic-parabolic schemes and auxiliary damped wave variables) enforce solenoidal and charge conservation constraints. In solar wind applications, an iterated Crank-Nicolson method handles multi-moment evolution along flux tubes [1908.09198, 1902.06816].

Theoretical analyses employ either closure-moment expansions or direct kinetic solution of the Vlasov-Maxwell system, with the latter needed for threshold and growth-rate calculations in the presence of strong anisotropy or relativistic corrections [1304.2661, 1205.6896].

---

For detailed functional forms, instability thresholds, and calculation techniques in specific anisotropy-driven regimes, see [1205.6896], [1304.2661], [1902.06816], and [1908.09198].

Source: https://www.emergentmind.com/topics/bi-maxwellian-plasma