Weibel instability

The Weibel instability is a plasma instability present in homogeneous or nearly homogeneous electromagnetic plasmas which possess an anisotropy in momentum (velocity) space. This anisotropy is most generally understood as two temperatures in different directions. Burton Fried showed that this instability can be understood more simply as the superposition of many counter-streaming beams. In this sense, it is like the two-stream instability except that the perturbations are electromagnetic and result in filamentation as opposed to electrostatic perturbations which would result in charge bunching. In the linear limit the instability causes exponential growth of electromagnetic fields in the plasma which help restore momentum space isotropy. In very extreme cases, the Weibel instability is related to one- or two-dimensional stream instabilities.

Consider an electron-ion plasma in which the ions are fixed and the electrons are hotter in the y-direction than in x or z-direction.

To see how magnetic field perturbation would grow, suppose a field B = B cos kx spontaneously arises from noise. The Lorentz force then bends the electron trajectories with the result that upward-moving-ev x B electrons congregate at B and downward-moving ones at A. The resulting current $$j = -en v_e$$ sheets generate magnetic field that enhances the original field and thus perturbation grows.

Weibel instability is also common in astrophysical plasmas, such as collisionless shock formation in supernova remnants and $$\gamma$$-ray bursts.

A Simple Example of Weibel Instability
As a simple example of Weibel instability, consider an electron beam with density $$n_{b0}$$ and initial velocity $$v_0 \mathbf{z}$$ propagating in a plasma of density $$n_{p0} = n_{b0}$$ with velocity $$-v_0 \mathbf{z}$$. The analysis below will show how an electromagnetic perturbation in the form of a plane wave gives rise to a Weibel instability in this simple anisotropic plasma system. We assume a non-relativistic plasma for simplicity.

We assume there is no background electric or magnetic field i.e. $$\mathbf{B_0} = \mathbf{E_0} = 0$$. The perturbation will be taken as an electromagnetic wave propagating along $$ \mathbf{\hat{x}} $$ i.e. $$\mathbf{k} = k \mathbf{\hat{x}}$$. Assume the electric field has the form


 * $$\mathbf{E_1} = A e^{i(kx-\omega t)} \mathbf{z} $$

With the assumed spatial and time dependence, we may use $$ \frac{\partial}{\partial t} \rightarrow -i \omega $$ and $$ \nabla \rightarrow i k \mathbf{\hat{x}} $$. From Faraday's Law, we may obtain the perturbation magnetic field


 * $$ \nabla \times \mathbf{E_1} = - \frac{\partial \mathbf{B_1}}{\partial t} \Rightarrow i \mathbf{k} \times \mathbf{E_1} = i \omega \mathbf{B_1} \Rightarrow \mathbf{B_1} = \mathbf{\hat{y}} \frac{k}{\omega} E_1 $$

Consider the electron beam. We assume small perturbations, and so linearize the velocity $$ \mathbf{v_b} = \mathbf{v_{b0}} + \mathbf{v_{b1}} $$ and density $$ n_b = n_{b0} + n_{b1} $$. The goal is to find the perturbation electron beam current density

$$ \mathbf{J_{b1}} = - e n_b \mathbf{v_b} = - e n_{b0} \mathbf{v_{b1}} - e n_{b1} \mathbf{v_{b0}} $$

where second-order terms have been neglected. To do that, we start with the fluid momentum equation for the electron beam


 * $$ m(\frac{\partial \mathbf{v_b}}{\partial t} + (\mathbf{v_b} \cdot \nabla) \mathbf{v_b}) = -e \mathbf{E} - e \mathbf{v_b} \times \mathbf{B} $$

which can be simplified by noting that $$\frac{\partial \mathbf{v_{b0}}}{\partial t} = \nabla \cdot \mathbf{v_{b0}} = 0 $$ and neglecting second-order terms. With the plane wave assumption for the derivatives, the momentum equation becomes

$$ -i \omega m \mathbf{v_{b1}} = -e \mathbf{E_1} - e \mathbf{v_{b0}} \times \mathbf{B_1} $$

We can decompose the above equations in components, paying attention to the cross product at the far right, and obtain the non-zero components of the beam velocity perturbation:


 * $$ v_{b1z} = \frac{e E_1}{m i \omega } $$


 * $$ v_{b1x} = \frac{e E_1}{m i \omega} \frac{k v_{b0}}{\omega} $$

To find the perturbation density $$ n_{b1} $$, we use the fluid continuity equation for the electron beam

$$ \frac{\partial n_b}{\partial t} + \nabla \cdot (n_b \mathbf{v_b}) = 0 $$

which can again be simplified by noting that $$ \frac{\partial n_{b0}}{\partial t} = \nabla n_{b0} = 0 $$ and neglecting second-order terms. The result is


 * $$ n_{b1} = n_{b0} \frac{k}{\omega} v_{b1x} $$

Using these results, we may use the equation for the beam perturbation current density given above to find


 * $$ J_{b1x} = - n_{b0} e^2 E_1 \frac{k v_{b0}}{i m \omega^2}$$


 * $$ J_{b1z} = - n_{b0} e^2 E_1 \frac{1}{i m \omega}(1+ \frac{k^2 v_{b0}^2}{\omega^2})$$

Analogous expressions can be written for the perturbation current density of the left-moving plasma. By noting that the x-component of the perturbation current density is proportional to $$v_0$$, we see that with our assumptions for the beam and plasma unperturbed densities and velocities the x-component of the net current density will vanish, whereas the z-components, which are proportional to $$v_0^2$$, will add. The net current density perturbation is therefore


 * $$ \mathbf{J_1} = -2 n_{b0} e^2 E_1 \frac{1}{i m \omega}(1+ \frac{k^2 v_{b0}^2}{\omega^2}) \mathbf{\hat{z}} $$

The dispersion relation can now be found from Maxwell's Equations:


 * $$ \nabla \times \mathbf{E_1} = i \omega \mathbf{B_1} $$


 * $$ \nabla \times \mathbf{B_1} = \mu_0 \mathbf{J_1} - i \omega \epsilon_0 \mu_0 \mathbf{E_1} $$


 * $$ \Rightarrow \nabla \times \nabla \times \mathbf{E_1} = -\nabla^2 \mathbf{E_1} + \nabla (\nabla \cdot \mathbf{E_1}) = k^2 \mathbf{E_1} + i \mathbf{k} (i \mathbf{k} \cdot \mathbf{E_1}) = k^2 \mathbf{E_1} = i \omega \nabla \times \mathbf{B_1} = \frac{i \omega}{c^2 \epsilon_0} \mathbf{J_1} + \frac{\omega^2}{c^2} \mathbf{E_1} $$

where $$ c = \frac{1}{\sqrt{\epsilon_0 \mu_0}} $$ is the speed of light in free space. By defining the effective plasma frequency $$ \omega_p^2 = \frac{2 n_{b0} e^2}{\epsilon_0 m} $$, the equation above results in


 * $$ k^2 - \frac{\omega^2}{c^2} = -\frac{\omega_p^2}{c^2}(1+\frac{k^2v_0^2}{\omega^2}) \Rightarrow \omega^4 - \omega^2 (\omega_p^2 + k^2 c^2) - \omega_p^2 k^2 v_0^2 = 0 $$

This bi-quadratic equation may be easily solved to give the dispersion relation

$$ \omega^2 = \frac{1}{2} (\omega_p^2 + k^2 c^2 \pm \sqrt{(\omega_p^2+k^2 c^2)^2 + 4 \omega_p^2 k^2 v_0^2} )$$

In the search for instabilities, we look for $$ Im(\omega) \neq 0 $$ ($$k$$ is assumed real). Therefore, we must take the dispersion relation/mode corresponding to the minus sign in the equation above.

To gain further insight on the instability, it is useful to harness our non-relativistic assumption $$ v_0 << c $$ to simplify the square root term, by noting that


 * $$ \sqrt{(\omega_p^2+k^2 c^2)^2 + 4 \omega_p^2 k^2 v_0^2} = (\omega_p^2 + k^2 c^2)(1+ \frac{4 \omega_p^2k^2v_0^2} {(\omega_p^2+k^2c^2)^2})^{1/2} \approx (\omega_p^2 + k^2 c^2)(1+ \frac{2 \omega_p^2k^2v_0^2}{(\omega_p^2+k^2c^2)^2}) $$

The resulting dispersion relation is then much simpler

$$ \omega^2 = \frac{-\omega_p^2 k^2 v_0^2}{\omega_p^2 + k^2c^2} < 0 $$

$$ \omega $$ is purely imaginary. Writing $$ \omega = i \gamma $$


 * $$ \gamma = \frac{\omega_p k v_0}{(\omega_p^2+k^2 c^2)^{1/2}} = \omega_p \frac{v_0}{c} \frac{1}{(1+\frac{\omega_p^2}{k^2 c^2})^{1/2}} $$

we see that $$ Im(\omega) > 0 $$, indeed corresponding to an instability.

The electromagnetic fields then have the form


 * $$ \mathbf{E_1} = A \mathbf{\hat{z}} e^{\gamma t} e^{i k x} $$


 * $$ \mathbf{B_1} = \mathbf{\hat{y}} \frac{k}{\omega} E_1 = \mathbf{\hat{y}} \frac{k}{i \gamma} A e^{\gamma t} e^{i k x} $$

Therefore, the electric and magnetic fields are $$90^o$$ out of phase, and by noting that

$$ \frac{|B_1|}{|E_1|} = \frac{k}{\gamma} \propto \frac{c}{v_0} >> 1 $$

so we see this is a primarily magnetic perturbation although there is a non-zero electric perturbation. The magnetic field growth results in the characteristic filamentation structure of Weibel instability. Saturation will happen when the growth rate $$ \gamma $$ is on the order of the electron cyclotron frequency

$$ \gamma \sim \omega_p \frac{v_0}{c} \sim \omega_c \Rightarrow B \sim \frac{m}{e} \omega_p \frac{v_0}{c} $$