Electron-Scale Weibel Filaments in Plasmas

• Electron-scale Weibel filaments are quasi-magnetostatic, self-organized current structures emerging from electron anisotropy in plasmas.
• They involve rapid exponential growth, nonlinear saturation, and merging processes that are modeled through PIC simulations and kinetic theories.
• These filaments, scaling with the electron skin depth, are key to magnetic field generation and energy conversion in HED experiments and astrophysical shocks.

Electron-scale Weibel filaments are quasi-magnetostatic, self-organized current structures generated via the Weibel instability in plasmas exhibiting electron temperature or momentum anisotropy. These filaments are fundamental to the rapid generation and self-organization of small-scale magnetic fields in a wide range of environments, including high-energy-density (HED) laser plasmas, relativistic beam–plasma systems, and astrophysical shocks. Their characteristic scale is set by the electron skin depth, $\delta_e = c/\omega_{pe}$, leading to the formation of magnetic and current structures on sub-micron to micron scales in laboratory plasmas, and up to macro-kilometer scales in astrophysical regimes.

1. Physical Basis: Linear Weibel Instability and Scaling Laws

The electron Weibel instability is rooted in the aperiodic growth of transverse magnetic modes in an anisotropic electron distribution, typically quantified by temperature anisotropy $A_e = T_\parallel/T_\perp - 1$ (in thermal systems) or by counterstreaming velocities (in beam-driven systems). For a bi-Maxwellian plasma, kinetic and cold-fluid models yield the general dispersion relation,

$$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$

where the instability operates for $A_e > (k c/\omega_{pe})^2$ and peaks at $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$ (Zhang et al., 2022, Stockem et al., 2010, Garasev et al., 2021, Kocharovsky et al., 2023). The fastest growing mode has wavelength

$$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$

with growth rate $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$ (cold limit), but reduced to $\gamma_{max}/\omega_{pe} \sim 0.1$–$0.3$ for realistic laboratory anisotropies or kinetic corrections (Zhang et al., 2022, Shukla et al., 2017). Similar scaling holds for relativistic counterstreaming-beam systems, where the anisotropy is set by the drift velocity $u_0$ and Lorentz factor $A_e = T_\parallel/T_\perp - 1$0, giving $A_e = T_\parallel/T_\perp - 1$1 and filaments with radius $A_e = T_\parallel/T_\perp - 1$2 (Huynh et al., 2022, Parsons et al., 2023).

2. Structure and Temporal Evolution of Electron-Scale Filaments

Electron-scale filaments manifest as self-pinched current channels (for $A_e = T_\parallel/T_\perp - 1$3 or $A_e = T_\parallel/T_\perp - 1$4) encircled by transverse, quasi-static magnetic fields ($A_e = T_\parallel/T_\perp - 1$5, $A_e = T_\parallel/T_\perp - 1$6), with spatial cross-sections typically $A_e = T_\parallel/T_\perp - 1$7–$A_e = T_\parallel/T_\perp - 1$8 $A_e = T_\parallel/T_\perp - 1$9 in radius. Two-dimensional (2D) and three-dimensional (3D) PIC simulations and experimental diagnostics consistently observe the following metrics:

Filament Parameter	Typical Value (HED lab plasma)	Reference
Filament radius $$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$0	$$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$1–$$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$2 (micron scale)	(Kocharovsky et al., 2023, Garasev et al., 2021, Parsons et al., 2023)
Spacing $$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$3	$$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$4–$$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$5 (few µm–100 µm)	(Zhang et al., 2022, Garasev et al., 2021, Kocharovsky et al., 2023)
Peak $$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$6-field	$$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$7–$$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$8 T (lab)	(Kocharovsky et al., 2023, Ngirmang et al., 2019, Garasev et al., 2021)
Current density $$\gamma^2(k) = \omega_{pe}^2 \left[ A_e - \left( \frac{k c}{\omega_{pe}} \right)^2 \right]$$9	$A_e > (k c/\omega_{pe})^2$0–$A_e > (k c/\omega_{pe})^2$1 A/m$A_e > (k c/\omega_{pe})^2$2	(Kocharovsky et al., 2023)
Energy fraction $A_e > (k c/\omega_{pe})^2$3	few percent	(Zhang et al., 2022, Kocharovsky et al., 2023)

Filament growth proceeds through distinct phases:

• Linear phase: Exponential amplification of B$A_e > (k c/\omega_{pe})^2$4 at rates determined by $A_e > (k c/\omega_{pe})^2$5, with initial fluctuations seeded by anisotropy or counterstreaming.
• Nonlinear saturation: Self-generated B-fields confine the electron trajectories such that their Larmor radius $A_e > (k c/\omega_{pe})^2$6, halting linear growth. The magnetic energy density $A_e > (k c/\omega_{pe})^2$7 1–5% of the electron thermal energy (Zhang et al., 2022, Kocharovsky et al., 2023).
• Merging/coalescence phase: Like-currents attract and merge into larger structures, shifting power to lower $A_e > (k c/\omega_{pe})^2$8 and increasing filament size. This occurs on the ion-acoustic or Alfvénic timescale in laser plasmas (Kocharovsky et al., 2023, Shukla et al., 2017).
• Late-time decay and turbulence: Filaments persist, decay via merging, and can seed quasi-turbulent magnetic fields relevant for downstream plasma dynamics (Kocharovsky et al., 2023).

3. Mechanisms of Energy Conversion and Field Structure

The Weibel process efficiently converts a fraction of the plasma thermal energy into magnetic field energy. In the UCLA experiment, up to $A_e > (k c/\omega_{pe})^2$9 of the initial electron thermal energy is converted into quasi-static magnetic fields near saturation (Zhang et al., 2022). This is consistent with broader PIC and experimental results in both thermal (Kocharovsky et al., 2023, Zhang et al., 2022, Ngirmang et al., 2019) and beam-driven (Parsons et al., 2023, Kumar et al., 2015) systems.

Electron heating is primarily driven by the longitudinal electric field $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$0 within filaments, which comprises both inductive (Faraday) and electrostatic components. Quantitative decomposition yields inductive contributions of 60% and electrostatic 40% in relativistic counterstreaming scenarios (Kumar et al., 2015). The net work $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$1 transfers energy predominantly to electrons, raising their mean kinetic energy to 25–30% of the initial ion energy in relativistic cases (Kumar et al., 2015).

The force balance at saturation involves both magnetic-pressure-gradient ($k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$2) and magnetic-tension ($k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$3) terms and, in multidimensional systems, the tension term is comparably important, facilitating the growth of nonlinear electric fields with force amplitudes on par with those from $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$4 (Stockem et al., 2010). Further, fully 2D/3D simulations reveal the emergence of “eddy” (transverse) current structures, which in turn generate out-of-plane B-field components, enabling the transition to fully 3D magnetic turbulence (Stockem et al., 2010).

4. Nonlinear Dynamics: Merging, Reconnection, and Turbulence

Post-saturation evolution is dominated by filament coalescence and magnetic reconnection (Shukla et al., 2017, Kocharovsky et al., 2023). The merging of like-current filaments increases the transverse scale $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$5 as $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$6, with current per filament and associated magnetic energy scaling as $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$7, $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$8. Once the filament current approaches the Alfvén threshold $k_{max} c/\omega_{pe} \simeq \sqrt{A_e}$9, further increases in current are suppressed, and subsequent mergers are mediated by collisionless reconnection events (Shukla et al., 2017).

These reconnection events release energetic jets of electrons in the plane transverse to the currents, producing rapid spikes in transverse kinetic energy and stepwise drops in magnetic energy. This mechanism is responsible for the rapid transverse heating observed experimentally and in simulation (Shukla et al., 2017, Kocharovsky et al., 2023).

Late-time dynamics include a cascade of magnetic energy from the initial filament scale to longer wavelengths, gradual decay of anisotropy, and slow energy transfer to larger-scale turbulence (Kocharovsky et al., 2023, Garasev et al., 2021). The resulting small-scale fields persist long after the initial instability ceases, providing a seed for further kinetic processes.

5. Laboratory and Astrophysical Manifestations

Electron-scale Weibel filaments have been observed directly in expanding laser-produced plasmas via ultrafast shadowgraphy and proton radiography, with spatial scales (diameters) of 1–2 µm and inter-filament spacings of 3–100 µm depending on plasma parameters (Ngirmang et al., 2019, Sutcliffe et al., 2022, Kocharovsky et al., 2023, Zhang et al., 2022). Self-generated magnetic fields reach tens to hundreds of Tesla, matching synthetic probe diagnostics from PIC simulations (Ngirmang et al., 2019). In high-energy-density (HED) experiments, observed filament wavelengths ($$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$0–$$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$1 µm) and growth rates ($$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$2–$$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$3 ns$$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$4) are consistent with electron-driven Weibel instability for modest anisotropy $$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$5, exceeding predictions from ion-driven or magnetothermal mechanisms (Sutcliffe et al., 2022). The magnetic power spectra may match analytic gyrokinetic predictions, with $$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$6 at small scales (Sutcliffe et al., 2022).

In collisionless relativistic shocks, Weibel filaments determine the microphysics of particle reflection and injection into diffusive shock acceleration (DSA). In PIC simulations, the radii $$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$7–$$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$8 $$\lambda_{max} = \frac{2\pi}{k_{max}} \simeq 2\pi \frac{c}{\omega_{pe} \sqrt{A_e}} = 2\pi \delta_e / \sqrt{A_e}$$9 and spacing $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$0–$\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$1 $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$2 determine reflection probabilities for electrons/positrons. Reflection and filament-hopping survival set the efficiency for particle injection into nonthermal tails at $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$3%–5% (Parsons et al., 2023).

Astrophysically, electron-scale filaments are relevant for the pre-magnetization and thermalization of collisionless shock fronts found in gamma-ray bursts, pulsar wind nebulae, and planetary magnetosheath current sheets (Kocharovsky et al., 2023, Garasev et al., 2021, Grassi et al., 2016).

6. Influence of External Magnetic Fields and Plasma Inhomogeneities

The presence of a pre-existing, flow-aligned magnetic field modifies the linear and nonlinear evolution of electron Weibel filaments (Grassi et al., 2016, Garasev et al., 2021, Kocharovsky et al., 2023). The linear growth rate and scale are suppressed for long wavelengths, with the critical field $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$4 above which all Weibel modes are quenched. However, at higher electron temperatures or short-wavelength modes, the nonlinear saturation amplitude and final filament structure are only weakly affected (Grassi et al., 2016). Guide fields aligned perpendicular to the anisotropy axis preferentially suppress the formation and persistence of z-pinch filaments, more strongly inhibiting small-scale turbulence and driving the system toward larger-scale current sheets or sheath currents (Garasev et al., 2021, Kocharovsky et al., 2023).

Inhomogeneous plasma density profiles also modulate filamentation properties: filament size scales as $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$5, and steep density gradients lead to localized filament formation (Kocharovsky et al., 2023). For expanding HED plasmas, filament characteristics are tied directly to the spatial and temporal evolution of electron anisotropy and local plasma density.

7. Spectral Signatures and Scaling Relations

The spectral evolution of magnetic fields generated by the Weibel instability exhibits characteristic behaviors:

• During the linear and early nonlinear regime, the wavenumber spectrum of B-fields is broad, but rapidly narrows as merging and coalescence select a dominant mode (Zhang et al., 2022).
• In the nonlinear coalescence regime, the spectral peak shifts toward lower $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$6 with a narrowing half-width, indicative of inverse cascading (Zhang et al., 2022, Shukla et al., 2017, Kocharovsky et al., 2023).
• Observationally, magnetic power spectra may follow a power law $\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$7 at scales below the electron Larmor radius, matching analytic gyrokinetic predictions (Sutcliffe et al., 2022).

Scaling relations extracted from experiment and simulation summarize the hierarchy of length, time, and field strength in both laboratory and astrophysical environments (Kocharovsky et al., 2023, Garasev et al., 2021):

Quantity	Formula
Skin depth	$\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$8
Filament radius	$\gamma_{max} \simeq \omega_{pe} \sqrt{A_e}$9–$\gamma_{max}/\omega_{pe} \sim 0.1$0
Growth rate (cold/rel.)	$\gamma_{max}/\omega_{pe} \sim 0.1$1; or $\gamma_{max}/\omega_{pe} \sim 0.1$2
B-field at saturation	$\gamma_{max}/\omega_{pe} \sim 0.1$3
Merging time	$\gamma_{max}/\omega_{pe} \sim 0.1$4–$\gamma_{max}/\omega_{pe} \sim 0.1$5
Energy fraction	$\gamma_{max}/\omega_{pe} \sim 0.1$6–$\gamma_{max}/\omega_{pe} \sim 0.1$7

Electron-scale Weibel filaments are thus a universal feature of anisotropic, collisionless plasmas, providing a robust mechanism for small-scale magnetization and energy transformation relevant to both laboratory HED experiments and astrophysical shock environments. Their detailed structure and evolution depend sensitively on anisotropy, plasma density, external fields, and system geometry (Zhang et al., 2022, Kocharovsky et al., 2023, Garasev et al., 2021, Ngirmang et al., 2019, Grassi et al., 2016, Parsons et al., 2023, Huynh et al., 2022, Stockem et al., 2010, Shukla et al., 2017, Kumar et al., 2015).