Temperature Anisotropy Relaxation in Plasmas

• Temperature anisotropy relaxation is the process by which plasma with unequal parallel and perpendicular temperatures returns to isotropy via collisions, instabilities, and turbulent interactions.
• The mechanisms vary across regimes, with kinetic microinstabilities dominating in collisionless plasmas and Coulomb collisions in strongly coupled or unmagnetized conditions.
• Simulations and theory reveal that suprathermal effects and strong magnetization can alter relaxation times by orders of magnitude, influencing energy partition and stability in both astrophysical and laboratory settings.

Temperature anisotropy relaxation is the process by which a plasma with unequal parallel ($T_\parallel$) and perpendicular ($T_\perp$) kinetic temperatures with respect to the magnetic field returns toward isotropy on microphysical and macroscopic timescales. The relaxation can occur via Coulomb collisions, quasilinear wave–particle interactions, microinstabilities, or anomalous scattering, with the governing mechanisms and timescales strongly dependent on the plasma regime (collisionality, magnetization, particle distributions, and turbulence). Temperature anisotropy relaxation is fundamental in space, laboratory, and astrophysical plasmas, determining kinetic stability, transport, energy partition, and the persistence or dissipation of non-equilibrium states.

1. Fundamental Mechanisms and Regimes

In weakly coupled or collisionless plasmas, anisotropy is primarily relaxed by kinetic microinstabilities—most notably the electromagnetic ion-cyclotron (EMIC), proton firehose (PFH), and mirror instabilities—and by quasilinear interactions with electromagnetic turbulence (Shaaban et al., 2024, Moya et al., 2020, Shaaban et al., 2016, Santos-Lima et al., 2016). In the collisional, unmagnetized limit (one-component plasma, OCP), relaxation is governed by Coulomb scattering and tractable by kinetic theory or molecular dynamics (Baalrud et al., 2017). In strongly magnetized or partially magnetized regimes relevant to laboratory and non-neutral plasmas, the distinction between parallel and perpendicular temperature relaxation rates becomes critical; gyromotion introduces different collisional and wave–particle pathways for each component, often with orders-of-magnitude separation in relaxation timescales (Jose et al., 2024, Welch et al., 12 Nov 2025).

2. Kinetic Descriptions and Instability Thresholds

In space and astrophysical plasmas, temperature anisotropy is typically quantified by $A = T_\perp / T_\parallel$. Plasma microinstabilities can be triggered for sufficient excess or deficit of $A$ away from unity at given plasma beta ($\beta_\parallel = 8\pi n k_B T_\parallel / B^2$). For $A>1$, the parallel EMIC and mirror instabilities grow; for $A<1$, the PFH and aperiodic firehose modes dominate (Shaaban et al., 2016, Shaaban et al., 2024, Santos-Lima et al., 2016):

$$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$

with empirical coefficients $a,b$ dependent on the mode, plasma composition, and suprathermal content (e.g., $a\approx 0.43$, $T_\perp$0 for EMIC (Moya et al., 2020); $T_\perp$1, $T_\perp$2 for PFH (Santos-Lima et al., 2016)).

The stability thresholds can be shifted by non-Maxwellian (kappa) velocity distributions and by electron–proton anisotropy correlations (Shaaban et al., 2024, Shaaban et al., 2016).

3. Relaxation Dynamics: Linear Theory, Quasilinear Evolution, and Simulation

Collisionless and Weakly Collisional Regimes

Linear kinetic theory predicts that growth rates of relevant instabilities ($T_\perp$3) reach a significant fraction of the proton gyrofrequency ($T_\perp$4), corresponding to e-folding relaxation times of $T_\perp$5—vastly faster than collisional or expansion timescales in the solar wind or ICM (Shaaban et al., 2016, Santos-Lima et al., 2016). Quasilinear theory and simulations show that once unstable, these microinstabilities generate growing electromagnetic fluctuations ($T_\perp$6) which scatter particle pitch angles, rapidly driving $T_\perp$7 back toward marginal stability.

Hybrid simulations and quasilinear calculations demonstrate accelerated and deeper anisotropy relaxation in the presence of suprathermal (bi-Kappa) populations, driven by higher $T_\perp$8 and amplified fluctuation power. For example, in an EMIC regime with $T_\perp$9, $A = T_\perp / T_\parallel$0, bi-Kappa distributions ($A = T_\perp / T_\parallel$1) yield a relaxation time $A = T_\perp / T_\parallel$2, a $A = T_\perp / T_\parallel$3 reduction compared to $A = T_\perp / T_\parallel$4 for a bi-Maxwellian, reflecting enhancement factors of 1.5–2 in growth rates and final fluctuation energy (Shaaban et al., 2024).

Collisional and Strongly Coupled OCP

In the OCP, molecular dynamics show exponential relaxation of anisotropy with a rate $A = T_\perp / T_\parallel$5, well described by the Landau formula at weak coupling ($A = T_\perp / T_\parallel$6), and increasingly by effective-potential theory at moderate coupling ($A = T_\perp / T_\parallel$7). Strongly coupled plasmas exhibit oscillatory structure and memory effects. Markovian kinetic theories fail to capture these features at high $A = T_\perp / T_\parallel$8 (Baalrud et al., 2017).

4. Turbulence-Driven Anisotropy and Self-Regulation

Background electromagnetic turbulence can generate temperature anisotropy even in nominally stable, isotropic plasmas. Quasilinear kinetic theory for a turbulent spectrum shows that damping of high-$A = T_\perp / T_\parallel$9 fluctuations (e.g., at $A$0) channels energy into $A$1, promoting $A$2 and reducing $A$3. Sufficiently energetic tails ($A$4) can drive $A$5 above instability thresholds, at which point instabilities re-grow and regulate the distribution near marginality. The final state is a self-regulated equilibrium with $A$6 and reduced $A$7, matching in situ observations in the solar wind (Moya et al., 2020).

In the turbulent intracluster medium (ICM), the instability-driven (quasilinear) relaxation rate dramatically outpaces the anisotropy driving rate from large-scale AMHD turbulence by $A$8 orders of magnitude, confining the plasma to values of $A$9 arbitrarily close to the kinetic instability thresholds (Santos-Lima et al., 2016).

5. Suprathermal Effects and Species Interplay

Suprathermal tails, described by finite-$\beta_\parallel = 8\pi n k_B T_\parallel / B^2$0 bi-Kappa distributions, lower the threshold for instability, enhance growth rates, and broaden unstable wavenumber bands for both EMIC and PFH modes. Hybrid simulations confirm that suprathermals accelerate and deepen anisotropy relaxation—reducing $\beta_\parallel = 8\pi n k_B T_\parallel / B^2$1 by up to 30\%, raising $\beta_\parallel = 8\pi n k_B T_\parallel / B^2$2 by factors of 2–10, and driving $\beta_\parallel = 8\pi n k_B T_\parallel / B^2$3 closer to unity or even overshooting marginal stability (Shaaban et al., 2024).

Electron–proton anisotropy correlations further modify thresholds and relaxation. Direct or anti-correlation (parameterized by $\beta_\parallel = 8\pi n k_B T_\parallel / B^2$4 in $\beta_\parallel = 8\pi n k_B T_\parallel / B^2$5) can either inhibit or enhance instability development, shifting observed bounds in the $\beta_\parallel = 8\pi n k_B T_\parallel / B^2$6 plane (Shaaban et al., 2016).

6. Strong Magnetization: Anisotropic Temperature Relaxation Rates

In strongly magnetized systems ($\beta_\parallel = 8\pi n k_B T_\parallel / B^2$7), Coulomb collisions and wave–particle interactions become highly anisotropic. The parallel ($\beta_\parallel = 8\pi n k_B T_\parallel / B^2$8) and perpendicular ($\beta_\parallel = 8\pi n k_B T_\parallel / B^2$9) temperature relaxation rates for a given species ($A>1$0) can differ by orders of magnitude (Jose et al., 2024, Welch et al., 12 Nov 2025). Characteristic behavior is as follows:

• For electrons with $A>1$1, ions weakly magnetized ($A>1$2):
  • $A>1$3 is suppressed by $A>1$4–$A>1$5 orders of magnitude ($A>1$6), $A>1$7 modestly enhanced ($A>1$8–$A>1$9).
  • Ion relaxation rates both modestly enhanced ($A<1$0–$A<1$1).
• When both electrons and ions are strongly magnetized, ion–electron perpendicular rates collapse ($A<1$2 for $A<1$3), parallel rates remain of order the unmagnetized value.

As a result, T$A<1$4 equilibrates rapidly (on timescales $A<1$5), while T$A<1$6 persists for much longer ($A<1$7), leading to long-lived temperature anisotropies (Jose et al., 2024, Welch et al., 12 Nov 2025).

Relaxation Dynamics and Scaling Example

$A<1$8	$A<1$9	$$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$0
0.5	1.00	1.00
8.7	1.15	0.05
34	1.17	0.013

(Jose et al., 2024)

This anisotropy in temperature relaxation can impose stringent bottlenecks on plasma cooling strategies in traps and astrophysical settings, requiring explicit tracking of T$$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$1 and T$$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$2 instead of assuming isotropy.

7. Implications, Applications, and Self-Regulation in Astrophysical and Laboratory Plasmas

The fast relaxation of temperature anisotropy in collisionless and weakly collisional plasmas via kinetic instabilities ensures that observed distributions in the solar wind and ICM remain tightly bounded near instability thresholds, explaining the apparent proximity to marginal stability observed in situ (Santos-Lima et al., 2016, Shaaban et al., 2016, Moya et al., 2020, Shaaban et al., 2024).

In laboratory contexts with strong magnetization (Penning–Malmberg traps, antihydrogen experiments), the vastly different T$$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$3 and T$$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$4 relaxation timescales necessitate multi-temperature models. The slow relaxation of T$$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$5 at high $$\begin{align*} &\text{Ion-cyclotron (parallel):}\qquad&&A > 1 + a\,\beta_\parallel^{-b}\ &\text{Mirror (oblique):}\qquad&&A > 1 + a'\,\beta_\perp^{-b'}\ &\text{Firehose (parallel):}\qquad&&A < 1 - a''\,\beta_\parallel^{-b''} \end{align*}$$6 can limit the efficiency of sympathetic cooling and the achievement of ultracold conditions (Welch et al., 12 Nov 2025, Jose et al., 2024).

Temperature anisotropy relaxation, whether via Coulomb collisions, microinstabilities, or turbulence-mediated pathways, constitutes a core transport process and a key constraint on plasma non-equilibrium structure in both natural and experimental plasmas.