Non-neutral plasma

In physics a non-neutral plasma is a plasma whose net charge creates an electric field large enough to play an important or even dominant role in the plasma dynamics. The simplest non-neutral plasmas are plasmas consisting of a single charge species. Examples of single species non-neutral plasmas that have been created in laboratory experiments are plasmas consisting entirely of electrons, pure ion plasmas, positron plasmas, and antiproton plasmas.

Non-neutral plasmas are used for research into basic plasma phenomena such as cross-magnetic field transport, nonlinear vortex interactions, and plasma waves and instabilities. They have also been used to create cold neutral antimatter, by carefully mixing and recombining cryogenic pure positron and pure antiproton plasmas. Positron plasmas are also used in atomic physics experiments that study the interaction of antimatter with neutral atoms and molecules. Cryogenic pure ion plasmas have been used in studies of strongly coupled plasmas and quantum entanglement. More prosaically, pure electron plasmas are used to produce the microwaves in microwave ovens, via the magnetron instability.

Neutral plasmas in contact with a solid surface (that is, most laboratory plasmas) are typically non-neutral in their edge regions. Due to unequal loss rates to the surface for electrons and ions, an electric field (the "ambipolar field") builds up, acting to hold back the more mobile species until the loss rates are the same. The electrostatic potential (as measured in electron-volts) required to produce this electric field depends on many variables but is often on the order of the electron temperature.

Non-neutral plasmas for which all species have the same sign of charge have exceptional confinement properties compared to neutral plasmas. They can be confined in a thermal equilibrium state using only static electric and magnetic fields, in a Penning trap configuration (see Fig. 1). Confinement times of up to several hours have been achieved. Using the "rotating wall" method, the plasma confinement time can be increased arbitrarily.

Such non-neutral plasmas can also access novel states of matter. For instance, they can be cooled to cryogenic temperatures without recombination (since there is no oppositely charged species with which to recombine). If the temperature is sufficiently low (typically on the order of 10 mK), the plasma can become a non-neutral liquid or a crystal. The body-centered-cubic structure of these plasma crystals has been observed by Bragg scattering in experiments on laser-cooled pure beryllium plasmas.

Equilibrium of a single species non-neutral plasma
Non-neutral plasmas with a single sign of charge can be confined for long periods of time using only static electric and magnetic fields. One such configuration is called a Penning trap, after the inventor F. M. Penning. The cylindrical version of the trap is also sometimes referred to as a Penning-Malmberg trap, after Prof. John Malmberg. The trap consists of several cylindrically symmetric electrodes and a uniform magnetic field applied along the axis of the trap (Fig 1). Plasmas are confined in the axial direction by biasing the end electrodes so as to create an axial potential well that will trap charges of a given sign (the sign is assumed to be positive in the figure). In the radial direction, confinement is provided by the v × B Lorentz force due to rotation of the plasma about the trap axis. Plasma rotation causes an inward directed Lorentz force that just balances the outward directed forces caused by the unneutralized plasma as well as the centrifugal force. Mathematically, radial force balance implies a balance between electric, magnetic and centrifugal forces:

where m is particle mass, q is particle charge, r is radial distance from the trap axis and Er is the radial component of the electric field. This quadratic equation can be solved for the rotational velocity $$ v_\theta$$, leading to two solutions, a slow-rotation and a fast-rotation solution. The rate of rotation $$\omega = -v_\theta/r$$ for these two solutions can be written as


 * $$ \omega = \frac{\Omega_c}{2} \pm\sqrt{{\Omega_c}^2/4-q E_r/{m r}} $$,

where $$\Omega_c = q B/m$$ is the cyclotron frequency. Depending on the radial electric field, the solutions for the rotation rate fall in the range $$ 0 \le \omega/\Omega_c \le 1 $$. The slow and fast rotation modes meet when the electric field is such that $$ q E_r/{m r} = {\Omega_c}^2/4 $$. This is called the Brillouin limit; it is an equation for the maximum possible radial electric field that allows plasma confinement.

This radial electric field can be related to the plasma density n through the Poisson equation,


 * $$ \frac{1}{r} \frac{\partial}{\partial r}(r E_r)= qn/\epsilon_0, $$

and this equation can be used to obtain a relation between the density and the plasma rotation rate. If we assume that the rotation rate is uniform in radius (i.e. the plasma rotates as a rigid body), then Eq. (1) implies that the radial electric field is proportional to radius r. Solving for Er from this equation in terms of $$\omega$$ and substituting the result into Poisson's equation yields

This equation implies that the maximum possible density occurs at the Brillouin limit, and has the value


 * $$ n_B= \frac{\epsilon_0 m {\Omega_c}^2}{2 q^2}=\frac{B^2/(2 \mu_0)}{m c^2}$$

where $$ c=1/\sqrt{\mu_0 \epsilon_0}$$ is the speed of light. Thus, the rest energy density of the plasma, n·m·c2, is less than or equal to the magnetic energy density $$ B^2/(2\mu_0)$$ of the magnetic field. This is a fairly stringent requirement on the density. For a magnetic field of 10 tesla, the Brillouin density for electrons is only nB = $$.

The density predicted by Eq.(2), scaled by the Brillouin density, is shown as a function of rotation rate in Fig. (2). Two rotation rates yield the same density, corresponding to the slow and fast rotation solutions.

Plasma loss processes; the rotating wall method
In experiments on single species plasmas, plasma rotation rates in the tens of kHz range are not uncommon, even in the slow rotation mode. This rapid rotation is necessary to provide the confining radial Lorentz force for the plasma. However, if there is neutral gas in the trap, collisions between the plasma and the gas cause the plasma rotation to slow, leading to radial expansion of the plasma until it comes in contact with the surrounding electrodes and is lost. This loss process can be alleviated by operating the trap in an ultra high vacuum. However, even under such conditions the plasma rotation can still be slowed through the interaction of the plasma with "errors" in the external confinement fields. If these fields are not perfectly cylindrically symmetric, the asymmetries can torque on the plasma, reducing the rotation rate. Such field errors are unavoidable in any actual experiment, and limit the plasma confinement time.

It is possible to overcome this plasma loss mechanism by applying a rotating field error to the plasma. If the error rotates faster than the plasma, it acts to spin up the plasma (similar to how the spinning blade of a blender causes the food to spin), counteracting the effect of field errors that are stationary in the frame of the laboratory. This rotating field error is referred to as a "rotating wall", after the theory idea that one could reverse the effect of a trap asymmetry by simply rotating the entire trap at the plasma rotation frequency. Since this is impractical, one instead rotates the trap electric field rather than the entire trap, by applying suitably phased voltages to a set of electrodes surrounding the plasma.

Cryogenic non-neutral plasmas: correlated states
When a non-neutral plasma is cooled to cryogenic temperatures, it does not recombine to a neutral gas as would a neutral plasma, because there are no oppositely charged particles with which to recombine. As a result, the system can access novel strongly coupled non-neutral states of matter, including plasma crystals consisting solely of a single charge species. These strongly coupled non-neutral plasmas are parametrized by the coupling parameter &Gamma;, defined as


 * $$\Gamma = \frac{q^2}{4 \pi\epsilon_0 k_B T a},$$

where $$ T$$ is the temperature and $$a$$ is the Wigner–Seitz radius (or mean inter-particle spacing), given in terms of the density $$n$$ by the expression $$ 4\pi a^3 n/3 =1$$. The coupling parameter can be thought of as the ratio of the mean interaction energy between nearest-neighbor pairs, $$ q^2/(4\pi\epsilon_0 a)$$, and the mean kinetic energy of order $$ k_b T$$. When this ratio is small, interactions are weak and the plasma is nearly an ideal gas of charges moving in the mean-field produced by the other charges. However, when $$\Gamma > 1$$ interactions between particles are important and the plasma behaves more like a liquid, or even a crystal if $$ \Gamma$$ is sufficiently large. In fact, computer simulations and theory have predicted that for an infinite homogeneous plasma the system exhibits a gradual onset of short-range order consistent with a liquid-like state for $$\Gamma \approx 2 $$, and there is predicted to be a first-order phase transition to a body-centered-cubic crystal for $$\Gamma\simeq 175$$.

Experiments have observed this crystalline state in a pure beryllium ion plasma that was laser-cooled to the millikelvin temperature range. The mean inter-particle spacing in this pure ion crystal was on the order of 10-20 μm, much larger than in neutral crystalline matter. This spacing corresponds to a density on the order of 108-109 cm−3, somewhat less than the Brillouin limit for beryllium in the 4.5 tesla magnetic field of the experiment. Cryogenic temperatures were then required in order to obtain a $$ \Gamma$$ value in the strongly coupled regime. The experiments measured the crystal structure by the Bragg-scattering technique, wherein a collimated laser beam was scattered off of the crystal, displaying Bragg peaks at the expected scattering angles for a bcc lattice (See Fig. 3).

When small numbers of ions are laser-cooled, they form crystalline "Coulomb clusters". The symmetry of the cluster depends on the form of the external confinement fields. An interactive 3D view of some of the clusters can be found here.