What is Plasma?

Plasma seems to have the kinds of properties one would like for life. It’s somewhat like liquid water — unpredictable and thus able to behave in an enormously complex fashion. It could probably carry as much information as DNA does. It has at least the potential for organizing itself in interesting ways.

– Freeman Dyson (1986)
Our sun is composed of plasma.

Plasma – colloquially referred to as the `4th state of matter’ – exhibits many interesting phenomena due to the diverse range of physics that govern its behavior, including electromagnetics, fluid/statistical mechanics, and even quantum mechanics. Most are familiar with the states of matter that are common on Earth. The reader is guaranteed to have observed how one type of matter in particular – water – transitions between the three common states of matter (solid ice, liquid water, water vapor), and thus they are familiar with the characteristics of each of those phases. A gas is best understood as a collection of molecules that have a high enough average speed that they overcome the intermolecular forces that make them stick to other alike molecules and fly about, experiencing collisions with one another. Plasma occurs when we inject even more energy into these molecules, breaking up their atomic structure in the process.

Plasma occurs when enough energy is injected into a material to break up its atomic structure.

To proceed with our water example, if we continue to heat up the gas, eventually the covalent bonds that hold together the water molecules will break apart into monatomic particles (the H-O bond energy of water is \sim450 kJ/mol, and O_2 has a bond energy of \sim500 kJ/mol). At this point, the substance is still a gas, albeit a very hot one. As we continue to inject energy into the system, we excite various elevated electronic energy states. and eventually, the electrons in the oxygen and hydrogen molecules will be knocked free of their orbitals (the ionization energies of oxygen and hydrogen are \sim1200 kJ/mol), and we are left with a plasma. Although this state of matter is uncommon on Earth (mostly observed by the naked eye in the form of lightning/electrical arcs and the aurora borealis), it is actually the state of matter that is most common in the universe as it is the phase of matter that stars are composed of. In the figure below, we see an oft-reproduced figure in the world of plasma physics, a plot of temperature vs. density, with photos of examples of plasmas that exist within the regime of density and temperature at that location in the plot. As one can see from this figure, plasma exists in the cold expanse of interstellar space, within comet tails, stellar nebulae, stars, lightning, and in fusion devices like those pictured at the top of the figure (Sandia Z-machine on the left, and an inertial confinement hohlraum on the right). Below what is called the Fermi temperature (a statistical/quantum mechanics-derived concept that arises from the statistical properties of some types of fundamental particles, here denoted by a solid black line), solid-state plasmas exist, which can be found in the cores of stars, and even in the form of metals here on Earth, which themselves are best-modeled as a plasma due to the fact that the electrons within them are more or less free to move.

A plot of electron temperature vs density, with photos of various examples of plasma placed at the location in the domain representative of their temperature and density. The dotted lines correspond to different values of Debye length, and the solid line represents the transition to a degenerate Fermi gas.

Now we can start to get more precise. First, let’s get an idea of how much ionization is present in a given gas at a given temperature. Suppose we want to figure out an approximation for what fraction of air molecules are ionized. Let’s simplify our analysis by assuming that air is roughly equivalent to monatomic nitrogen, and then let’s find our how much of it is singly-ionized. If we assume we are in thermal equilibrium at a temperature that is not close to zero, then we are safe in assuming that the energy distribution of the particles follows the Boltzmann distribution (see the derivation in another learn page) where the number of particles in the jth energy level is N_j^*\sim C_je^{-\epsilon_j/k_BT}/Q, where C_j is the number of available states in the jth energy level, \epsilon_j is the jth energy level relative to the ground state and Q = \sum_jC_je^{-\epsilon_j/k_BT} is the partition function [1]. If we sum over all the energy levels to find the total number of particles in each species (the three species here are neutral atomic nitrogen, singly ionized nitrogen N^+, and electrons) and compute the quotient of these totals, we get an expression reminiscent of the Law of Mass Action:

    \begin{gather*}\frac{N_{\Te{N}^+}^*N_{e^-}^*}{N_{\Te{N}}^*} = \frac{Q_{\Te{N}^+}Q_{e^-}}{Q_{\Te{N}}} e^{-\frac{\epsilon_i}{k_BT}}.\end{gather*}

where \epsilon_i is the ionization energy.

The partition function in general has two contributions, translational and internal; Q=Q_{tr}Q_{int}. For the electrons, the internal contribution is simply Q_{int} = 2 since they are fundamental particles that can either be spin up or spin down. The internal contributions of the atoms and ions are more complicated and have to do with electronic excitation (since we are dealing with monatomic gas, we do not need to concern ourselves with rotation and vibration of molecular bonds), but are mostly second order so we will assume Q_{int}\approx1 for both. The translational contribution to the partition function is, in general,

    \begin{gather*}Q_{tr} = V\left(\frac{2\pi mk_BT}{h^2}\right)^{\frac{3}{2}},\end{gather*}

where h=6.626\times10^{-34} J-s is Planck’s constant. This partition function arises from considering the available translational energy states and assuming we are in the Boltzmann limit (states sparsely populated). Now, we need to note a couple of things. First, m_{\Te{N}<em>+}\approx m</em>{\Te{N}}. Second, we can replace the volume term V with N/n, where n is the number density of the given species. Finally, since each ionization event creates exactly one ion and one electron, we have N_{\Te{N}^+}=N_{e^-}. With these things in mind, we get:

(1)   \begin{gather*}\frac{N_{\Te{N}^+}^*N_{e^-}^*}{N_{\Te{N}}^*} = 2\frac{N_{\Te{N}^+}^*N_{e^-}^*}{N_{\Te{N}}^*}\frac{n_{\Te{N}}^*}{n_{\Te{N}^+}^*n_{e^-}^*} \left(\frac{2\pi m_ek_BT}{h^2}\right)^{\frac{3}{2}} e^{-\frac{\epsilon_i}{k_BT}}\notag\\ \frac{n_{\Te{N}^+}^{2}}{n_{\Te{N}}^*} = \left(\frac{2^{5/3}\pi m_ek_B}{h^2}\right)^{\frac{3}{2}} T^{\frac{3}{2}} e^{-\frac{\epsilon_i}{k_BT}} \approx 4.8\times10^{21} T^{\frac{3}{2}} e^{-\frac{\epsilon_i}{k_BT}}. \end{gather*}

This equation is referred to as the Saha equation, and gives us an estimate for degree of ionization of a gas in thermal equilibrium. For our nitrogen example, if we assume we are at room temperature T\approx300 K where the neutral density is about n_n\approx 3\times 10^{25} m^{-3}, since the ionization energy of nitrogen is \epsilon_i=14.5 eV, the fractional ionization is basically nil; n_i/n_n\sim 10^{-122}. That matches our experience; we do not see much plasma appearing at sea level. But nonetheless, the Saha equation tells us that every gas has some degree of ionization, so where do we draw the line and call it a plasma?

Quantitative and Qualitative Descriptions and Criteria

The definition of plasma given in [2] is that plasma is a quasineutral gas of charges and neutral particles which exhibits collective behavior. Collective behavior here is a bit of a qualitative term; it is referring to the characteristic of plasma that Dyson was describing in the epigraph preceding this section. By collective behavior, we mean that local behavior of the constituents of a plasma is driven by the motion of other constituents, sometimes quite strongly. Plasma gets its name from Irving Langmuir who was reminded of the properties of blood plasma, but it’s a bit of a misnomer since plasma translates in Ancient Greek to moldable `substance’, where plasma is anything but. If we `squeeze’ a plasma in one location, we create some sort of strong local current density and displacement, which then produces the resultant EM fields, which then go on to affect the motion of the constituents of the plasma both near and far from our `squeeze’ location. Because of this, plasma seems to have a mind of its own. Moreover, if the degree of ionization is too low, there is not enough charge motion and separation to create the requisite macroscopic fields to strongly affect particles far away from the `squeeze’.

The `quasi-neutral’ part of the definition has to do with a phenomenon called `shielding’. Imagine we have a plasma background with a decently high density (forgive the imprecise language; I will put a quantitative criterion on this shortly), and then we place a positive test charge in the plasma. Immediately, the electrons in the plasma will rush to meet the test charge, that is, until the ions that they `left behind’ (ions are much more massive than electrons and thus their inertia prevents them from moving on the same time scale) begin to represent an effective cloud of positive charge surrounding the system, pulling the electrons back in the direction from which they came. If one draws a boundary around this area of charge displacement, the view from outside the boundary is as if no test charge is present since the constituents of the plasma have shielded out its potential. If there was no plasma present, an observer form the same distance would see the potential produced by the test charge in vacuum. A diagram of this scenario is shown in the figure below.

Test charge placed in vacuum (left) and within a plasma (right). The electric fields of the test charge in vacuum are felt from far away, while in the plasma case, the redistribution of the charges makes the test charge ‘invisible’ outside of the gray dashed line.

We’re going to derive the length scale over which this shielding occurs. Let’s imagine we have a plane of constant potential \phi_0 at x=0 in an infinite uniform singly ionized plasma, and then let’s see how the potential changes as we move away from the origin. To keep things simple, we ignore the ion motion due to the mismatch in inertia I mentioned before, so there is a uniform background of positive ions in which electrons are free to move about. In this case, Poisson’s equation (Gauss’ law with \Bf{E}=-\nabla\phi) becomes

(2)   \begin{gather*}\ders[\phi]{x} = \frac{e}{\eps_0}(n_e(x)-n_i).\end{gather*}

We will set the density far away from the origin to n_{\infty}=n_i. Assuming that the plasma is in thermal equilibrium, the electron distribution function is of the Boltzmann type;

    \begin{gather*}f_e(u) = Ae^{-(mu^2/2-e\phi)/k_BT_e}.\end{gather*}

If we integrate this distribution over all speeds and note that as \phi\RA0, n_e\RA n_{\infty}, we can obtain the electron density in terms of the potential \phi:

    \begin{gather*}n_e = n_{\infty}e^{\frac{e\phi}{k_BT_e}}.\end{gather*}

We substitute this into Poisson’s equation 2;

    \begin{gather*}\ders[\phi]{x} = \frac{en_{\infty}}{\eps_0}(e^{\frac{e\phi}{k_BT_e}}-1).\end{gather*}

As we discussed earlier, the potential is going to fall quite fast as we move away from the origin, so we can safely assume |e\phi/k_BT_e|\ll1 and expand the exponential in a Taylor series. We will drop the second order term so that

    \begin{gather*}\ders[\phi]{x} = \frac{e^2n_{\infty}}{\eps_0k_BT_e}\phi,\end{gather*}

leading to a solution of exponential decay;

    \begin{gather*}\phi = \phi_0e^{-\sqrt{\frac{e^2n_{\infty}}{\eps_0k_BT_e}}|x|}.\end{gather*}

This defines a length scale which is called the Debye length:

(3)   \begin{gather*}\lambda_D = \sqrt{\frac{\eps_0k_BT_e}{e^2n_{\infty}}}.\end{gather*}

This distance can be understood as a characteristic shielding length in a plasma. If the Debye length is sufficiently small (in other words, if the characteristic length of your experimental scenario is L, \lambda_D\ll L), then any local disruptions of the charge balance in the plasma will be shielded away from the macro perspective; i.e. the plasma can be considered quasi-neutral. Of course, we need to have a plasma density high enough such that the number of charges in the Debye sphere is large, N_D = n_e\frac{4\pi\lambda_D^3}{3} \gg 1.

Another important characteristic quantity of a plasma is the plasma frequency \omega_p. Suppose you have a slab of plasma for which all the electrons in an interval (x_0,x_0+\Delta x) are shifted one \Delta x in the positive direction, resulting in an electric field at x=x_0+\Delta x. If we use the integral form of Gauss’ law with one boundary of the integration surface at x=x_0+\Delta x with surface area L^2, and the other boundary far off at large x where the plasma is neutral, then we can come up with an expression for the magnitude of the electric field as a function of the displacement;

    \begin{gather*}\int_S\Bf{E}\cdot\Te{d}\Bf{S} = \frac{Q_{encl}}{\eps_0}\\-L^2E = \frac{-L^2\Delta xn_e e}{\eps_0}\\E = \frac{\Delta xn_e e}{\eps_0}.\end{gather*}

We know that the equation of motion for a single electron is \ders[x]{t} = -eE/m_e, so then the equation that describes the dynamics of the displaced electrons in the plasma is

    \begin{gather*}\ders[\Delta x]{t} = -\frac{n_e e^2}{\eps_0m_e}\Delta x,\end{gather*}

which denotes an oscillation with frequency

(4)   \begin{gather*}\omega_p = \sqrt{\frac{e^2n_e}{\eps_0 m_e}},\end{gather*}

known as the plasma frequency. Another key criterion for plasma is that the frequency of collisions between charged particles and neutral atoms \nu is small compared to the plasma frequency, \omega_p \gg \nu. This ensures that the dynamics of the plasma do not end up most closely resembling plain-old hydrodynamics.

Bibliography

This article is an excerpt from Prof. Rodríguez’s PhD thesis. If you would like to cite it for your own work, please cite:

– J. A. Rodríguez, Electromagnetic Wave Manipulation with Plasma Metamaterials, Ph.D. thesis, Stanford University (2023).

[1]: Walter G. Vincenti and Jr. Charles H. Kruger. Introduction to Physical Gas Dynamics. Krieger Publishing, Malabar, FL, 1986.

[2]: Francis F. Chen. Introduction to Plasma Physics and Controlled Fusion, Third Edition. Springer Publishing, Switzerland, 2015.