Weibel Instability: A Simple Example of Plasma Filamentation
The Weibel instability is a fundamental plasma phenomenon where an anisotropic distribution of particle velocities leads to the spontaneous growth of magnetic fields. This process is critical in astrophysics and laboratory plasma physics, as it explains how magnetic fields can be generated from unmagnetized plasma. To understand the mechanics of this instability, we can examine a simplified model involving counter-streaming electron beams.
The Physical Model
Consider a system consisting of an electron beam with density nb0 moving with an initial velocity v0z. This beam propagates through a plasma of equal density (np0 = nb0) moving in the opposite direction with velocity -v0z. For the sake of simplicity, this analysis assumes a non-relativistic plasma where the initial background electric and magnetic fields are zero (B0 = E0 = 0).
To analyze the stability of this system, we introduce an electromagnetic perturbation in the form of a plane wave propagating along the x-axis (k = kx̂). The electric field perturbation is defined as E1 = Aei(kx-ωt)ẑ. Using Faraday's Law, we can derive the corresponding perturbation magnetic field:
B1 = ŷ(k/ω)E1
ไม่มีภาพประกอบ
Mathematical Derivation of the Instability
Beam Dynamics and Current Density
By linearizing the velocity (vb = vb0 + vb1) and density (nb = nb0 + nb1) of the electron beam, we can determine the perturbation current density Jb1. Using the fluid momentum equation and the fluid continuity equation, the non-zero components of the beam velocity perturbation are found to be:
- vb1z = eE1 / (miω)
- vb1x = (eE1 / miω) * (kvb0 / ω)
When combining the contributions from both the beam and the counter-streaming plasma, the x-components of the net current density vanish. However, the z-components add together, resulting in a net current density perturbation:
J1 = -2nb0e2E1 / (imω) * (1 + k2vb02 / ω2)ẑ
The Dispersion Relation
By applying Maxwell's Equations and defining the effective plasma frequency as ωp2 = 2nb0e2 / (ε0m), we arrive at a bi-quadratic dispersion relation:
ω4 - ω2(ωp2 + k2c2) - ωp2k2v02 = 0
Solving for ω2, we look for modes where the imaginary part of ω is non-zero (ℑ(ω) ≠ 0), which indicates an instability. Under the non-relativistic assumption (v0 ≪ c), the dispersion relation simplifies to:
ω2 = - (ωp2k2v02) / (ωp2 + k2c2)
Since ω2 is negative, ω is purely imaginary (ω = iγ). The growth rate γ is positive, confirming that the perturbation will grow exponentially over time.
Key Facts
- Nature of Perturbation: The instability is primarily magnetic, as the ratio of the magnetic field to the electric field (|B1|/|E1|) is proportional to c/v0, which is much greater than 1.
- Phase Relationship: The electric and magnetic fields are 90° out of phase.
- Physical Result: The growth of the magnetic field leads to the characteristic filamentation structure of the Weibel instability.
- Saturation Point: Saturation occurs when the growth rate γ reaches the order of the electron cyclotron frequency (ωc).
Summary of System Parameters
| Parameter | Symbol/Formula | Description |
|---|---|---|
| Plasma Frequency | ωp2 = 2nb0e2 / (ε0m) | Effective frequency of plasma oscillations |
| Growth Rate | γ = ωp(v0/c) / (1 + ωp2/k2c2)1/2 | The rate at which the instability grows |
| Saturation Field | B ∼ (m/e)ωp(v0/c) | Magnetic field strength at saturation |
| Velocity Condition | v0 ≪ c | Non-relativistic assumption |
Frequently Asked Questions
What causes the Weibel instability?
The Weibel instability is caused by an anisotropy in the particle velocity distribution, such as counter-streaming beams of plasma, which triggers the growth of transverse electromagnetic perturbations.
Why is it called a "filamentation" instability?
It is called filamentation because the growing magnetic fields deflect the particles into concentrated streams or "filaments," creating a structured pattern in the plasma.
Is the resulting perturbation purely magnetic?
While there is a non-zero electric perturbation, the instability is primarily magnetic because the magnetic field magnitude is significantly larger than the electric field magnitude in the non-relativistic limit.
When does the growth of the instability stop?
The instability reaches saturation when the growth rate γ becomes comparable to the electron cyclotron frequency, at which point the magnetic field strength stabilizes.
What is the role of the dispersion relation in this analysis?
The dispersion relation relates the angular frequency (ω) to the wavenumber (k). By finding solutions where ω has a positive imaginary component, physicists can prove that a perturbation will grow exponentially rather than oscillate.