Bi-Maxwellian Plasma Dynamics

• Bi-Maxwellian plasma is a kinetic system with a Maxwellian velocity distribution along the magnetic field and a different perpendicular temperature, leading to pressure anisotropy.
• It modifies electromagnetic wave propagation by altering phase velocities and resonance conditions, thereby supporting instabilities such as firehose, mirror, and whistler modes.
• Multi-moment fluid closures and advanced numerical techniques are used to model bi-Maxwellian effects, enhancing our understanding of solar wind, laboratory plasma, and astrophysical applications.

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}$ (Huang et al., 2019, Bashir et al., 2012).

2. Linear Modes, Wave Propagation, and Dielectric Response

The anisotropy inherent to bi-Maxwellian plasmas modifies the electromagnetic dielectric tensor. The relevant permittivity tensor $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$0, as derived from the linearized Vlasov–Maxwell equations, contains explicit dependency on the $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$1 ratio through Bessel-function sums and plasma dispersion functions. Key effects include:

• R/L (Right/Left) waves (parallel propagation): Anisotropy enters as a correction $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$2 to the phase velocities and modifies resonance conditions.
• Whistler (low-frequency R) branch: For $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$3 the whistler phase speed is increased; for $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$4 it is decreased.
• O (ordinary) mode: For perpendicular propagation, anisotropy corrections arise at finite-Larmor-radius orders, especially via terms like $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$5.
• X and Bernstein modes: To leading order, these remain unaffected by anisotropy, as the explicit combination terms cancel (Bashir et al., 2012).

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 $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$6. In dimensionless terms, instability requires $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$7 or, at the kinetic level, $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$8. This produces negative effective magnetic tension and destabilizes parallel modes (Huang et al., 2019, Bashir et al., 2012).
• Mirror instability: Occurs for $\mathbf{b} = \mathbf{B}/|\mathbf{B}|$9. It preferentially drives nonpropagating oblique modes (Huang et al., 2019).
• Whistler (Weibel-like) instability: For $$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),$$0, non-resonant and resonant branches both support instabilities; field-free cases reduce to the classic Weibel instability ($$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),$$1). Ambient magnetic field and relativistic corrections generally suppress growth. The instability is quenched entirely in the ultra-relativistic regime ($$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),$$2) (Bashir et al., 2013).
• 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 (Bashir et al., 2012).

A summary of instability criteria is provided below:

Instability	Threshold Condition	Physical Driver
Firehose	$$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),$$3	$$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),$$4
Mirror	$$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),$$5	$$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),$$6
Weibel/Whistler	$$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),$$7	$$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),$$8

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 $$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),$$9 comprises:

• Mass continuity:

$v_\parallel = \mathbf{v} \cdot \mathbf{b}$0

• Momentum (with full pressure tensor):

$v_\parallel = \mathbf{v} \cdot \mathbf{b}$1

• Separate parallel and perpendicular pressure evolution equations with coupling to velocity shear and expansion/compression (Huang et al., 2019).

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 (Huang et al., 2019).

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 (Taran et al., 2019).

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

$v_\parallel = \mathbf{v} \cdot \mathbf{b}$2

reducing to the non-relativistic bi-Maxwellian in the $v_\parallel = \mathbf{v} \cdot \mathbf{b}$3 limit. Relativistic corrections generally suppress growth rates and instability thresholds, quenching anisotropy-driven instabilities as $v_\parallel = \mathbf{v} \cdot \mathbf{b}$4 approaches $v_\parallel = \mathbf{v} \cdot \mathbf{b}$5 (Bashir et al., 2013).

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 $v_\parallel = \mathbf{v} \cdot \mathbf{b}$6 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 (Taran et al., 2019).

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 (Huang et al., 2019).

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 (Taran et al., 2019, Huang et al., 2019).

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 (Bashir et al., 2013, Bashir et al., 2012).

---

For detailed functional forms, instability thresholds, and calculation techniques in specific anisotropy-driven regimes, see (Bashir et al., 2012, Bashir et al., 2013, Huang et al., 2019), and (Taran et al., 2019).