Debye Shielding

A plasma is a state in which a large number of charged particles move about independently yet, as a whole, behave like a single collective. The best entry point to that “collective behavior” is the topic of this page: Debye shielding. Once you understand Debye shielding, you gain a length scale called the Debye length, and from there you can follow a single continuous path through plasma oscillations, quasineutrality, and finally the rigorous definition of “what a plasma even is.”

Start with Intuition (High School)

Imagine placing a single tiny sphere carrying a positive charge inside a plasma. Around this sphere there are already many electrons (negative charge) and ions (positive charge) that were flying around within the plasma.

The positive sphere attracts electrons and repels ions. As a result, a “cloud of electrons,” where negative electrons gather in slight excess, forms around the sphere. This cloud makes the positive charge of the sphere harder to see from the outside. It is just like covering a dazzling light bulb with a black mesh, so that its brightness seen from far away is softened. This is Debye shielding.

Because electrons are much lighter than ions, they move quickly and form the cloud first. As a result, seen from another particle a little way off, the influence of the sphere’s charge almost disappears. The measure of “how far the electric influence reaches” is the length called the Debye length.

The important point here is why the electric influence does not vanish completely, but remains out to some distance. If the particles were cold (low temperature) and sat still, the electrons would stick tightly to the sphere and completely cancel its charge. But real particles carry heat and are constantly moving around. This thermal motion blurs the cloud and makes the shielding incomplete. In other words, the Debye length is the length set by the balance between the electric force that pulls particles in and the thermal agitation that tries to scatter them apart. The higher the temperature, the blurrier the cloud and the longer the Debye length; the higher the particle density, the denser the cloud and the shorter the Debye length.

Understand the Physics (Undergraduate)

Let us express Debye shielding with equations. In vacuum, the potential around a charge $q$ follows the Coulomb potential

$$\phi_{\text{vac}}(r) = \frac{q}{4\pi\varepsilon_0 r}$$

and decays slowly, inversely proportional to the distance $r$. Here $\varepsilon_0$ is the permittivity of free space. In a plasma, however, the surrounding charged particles rearrange to cancel the potential, so it changes into a Yukawa-type (screened Coulomb) potential

$$\phi(r) = \frac{q}{4\pi\varepsilon_0 r}\, e^{-r/\lambda_D}$$

The exponential factor $e^{-r/\lambda_D}$ takes effect, and once the distance exceeds $\lambda_D$ the potential rapidly approaches zero. This $\lambda_D$ is the Debye length.

The Debye length can be written using the electron temperature $T_e$ and electron density $n_e$ as

$$\lambda_D = \sqrt{\frac{\varepsilon_0 k_B T_e}{n_e e^2}}$$

where $k_B$ is the Boltzmann constant and $e$ is the elementary charge. The fact that the temperature $T_e$ appears in the numerator and the density $n_e$ in the denominator matches the earlier intuition (longer at higher temperature, shorter at higher density).

The line of derivation is as follows. Within a potential $\phi$, the electron density follows a Boltzmann distribution $n_e \exp(e\phi/k_B T_e)$. Assuming the potential is small ($e\phi \ll k_B T_e$), we expand the exponential to first order and substitute it into Poisson’s equation $\nabla^2 \phi = -\rho/\varepsilon_0$, obtaining an equation of the form $\nabla^2 \phi = \phi/\lambda_D^2$. Its solution is the Yukawa-type potential above, and $\lambda_D$ naturally appears as the length scale of the shielding.

In a fusion plasma (temperature around 10 keV, density around $10^{20}\ \text{m}^{-3}$), the Debye length is about 70 micrometers. This length is far smaller than the device size (meters) and, conversely, far larger than the mean interparticle distance (about 0.2 micrometers at this density), satisfying these two inequalities simultaneously. This relation of being “sufficiently smaller than the device and sufficiently larger than the interparticle spacing” is the heart of the definition of a plasma discussed later.

Another important quantity is plasma oscillation. When the electron population in a uniform plasma is displaced slightly as a whole by $x$, charge separation arises with respect to the positive ions, and an electric field (restoring force) acts to bring them back. The electrons begin simple harmonic motion under this restoring force, and its angular frequency is given by the electron plasma frequency

$$\omega_{pe} = \sqrt{\frac{n_e e^2}{\varepsilon_0 m_e}}$$

where $m_e$ is the electron mass. In a fusion plasma the frequency $f_{pe} = \omega_{pe}/2\pi$ is about 90 GHz. Because ions are much more massive, their response is slower, and the ion plasma frequency is smaller than the electron case by a factor of $\sqrt{m_e/m_i}$. For hydrogen ions it is about 1/43 that of electrons.

The Debye length and the plasma frequency are linked through the electron thermal velocity $v_{th,e} = \sqrt{k_B T_e/m_e}$ by the relation $\lambda_D \sim v_{th,e}/\omega_{pe}$. This means that “during one plasma oscillation an electron travels roughly one Debye length by thermal motion,” showing that shielding and oscillation arise from the same physics.

Deepen the Theory (Graduate)

The quantity that expresses how many particles are contained within a Debye sphere (a sphere of radius $\lambda_D$) is the plasma parameter. The number of particles inside the Debye sphere is defined as

$$N_D = \frac{4}{3}\pi \lambda_D^3\, n_e$$

and $N_D \gg 1$ is the condition for collective behavior. In a fusion plasma $N_D$ reaches about $10^8$. The larger $N_D$ is, the more particles cooperate to carry the shielding and oscillation, so collective effects dominate over individual two-body Coulomb collisions. Conversely, a system with $N_D \sim 1$ is called a “strongly coupled plasma,” where the statistical treatment breaks down and a different theory is required. Fusion plasmas are classified as “weakly coupled, ideal plasmas” with $N_D \gg 1$.

With this in mind, the three rigorous conditions for an ionized gas to behave as a plasma can be organized as follows.

First, the Debye length must be sufficiently smaller than the system size $L$ ($\lambda_D \ll L$). This ensures that charge is shielded throughout most of the system and quasineutrality holds.

Second, the number of particles inside the Debye sphere must be large ($N_D \gg 1$). This gives the shielding a statistical meaning, so that collective effects exceed two-body collisions.

Third, the electron plasma frequency must be sufficiently higher than the collision frequency $\nu$ between electrons and neutral particles ($\omega_{pe}\tau > 1$, $\tau = 1/\nu$). This is the condition that electrons are not disrupted by collisions with neutrals before completing a plasma oscillation, so that electromagnetic collective motion prevails over gas-like collisional motion.

These are not independent of one another; each expresses a single idea, that “electromagnetic collective effects prevail over individual particle motion and thermal disturbances,” from the three aspects of length, particle number, and time.

Quasineutrality is the approximation that, on scales larger than the Debye length, the electron density and the ion density are nearly equal ($n_e \approx Z n_i$) and the net charge density can be regarded as nearly zero. Many plasma theories, starting with MHD (magnetohydrodynamics), take this quasineutrality as their starting point. Quasineutrality is not universal, however: it breaks down on small scales of order the Debye length, in high-frequency electric fields, and at the boundary between a plasma and a wall.

A representative example where quasineutrality breaks down is the sheath. When a plasma contacts a solid wall, the fast-moving electrons flow into the wall first and charge it negatively, and a positive space-charge layer with a thickness of a few times the Debye length forms in front of the wall. This is the sheath, and quasineutrality does not hold within it. The sheath has an important condition: the Bohm criterion, which states that a stable sheath cannot form unless the ions flowing into the wall are accelerated to at least the ion sound speed $c_s = \sqrt{k_B T_e/m_i}$. The sheath governs the flow of particles and energy from the plasma to the wall, and determines how the motion of charged particles terminates near the wall. The motion of particles within the bulk plasma itself is treated in Charged Particle Motion.

Research Frontier (PhD)

The physics of Debye shielding and sheaths, though a classical theory, remains directly connected to the frontier of fusion research today.

One of the most active subjects of study is plasma-wall interaction. In a tokamak divertor, the energy and angle of the ions that reach the wall through the sheath govern the sputtering, erosion, and heat load of the wall material. The energy distribution of ions passing through the sheath, and the structure of the magnetized sheath, where the magnetic field lines strike the wall at a shallow angle, are being studied through both experiment and kinetic simulation. For the engineering design of and challenges facing divertors, see Divertor.

In sheath modeling, phenomena that cannot be captured by a fluid treatment alone become important. For secondary electron emission at the wall surface, non-Maxwellian velocity distributions, and transient phenomena (the bursty heat loads brought about by edge instabilities called ELMs), kinetic simulations by the PIC (particle-in-cell) method play an important role. How to generalize the Bohm criterion to non-Maxwellian distributions and magnetized situations is also a continuing research topic.

The Debye length is also deeply involved in the principles of plasma diagnostics. In a Langmuir probe, the theory of the sheath that forms on the probe surface determines the measurement accuracy. In Thomson scattering, the relative magnitude of the observed wavelength scale (the inverse of the wavenumber) and the Debye length distinguishes whether the scattering is “incoherent scattering” from individual electrons or “collective scattering” that reflects the collective motion of electrons. The quantity that expresses this boundary is the Salpeter parameter $\alpha = 1/(k\lambda_D)$, which is the foundation for deriving density and temperature from measurements.

English keywords that appear frequently when reading papers include Debye length, plasma parameter, quasineutrality, sheath, presheath, Bohm criterion, plasma-wall interaction, and magnetized sheath. All of these lie on the extension of the basic concepts of shielding, oscillation, and quasineutrality built up on this page.

Check Your Understanding

Q1. Around a charge placed in a plasma, why does the electric influence reach only out to a certain distance?

A Because the particles constantly move about through thermal motion, the shielding cloud blurs and the shielding becomes incomplete B Because the particles come completely to rest and stick tightly to the charge, canceling the electricity C Because the charge itself disappears with time D Because there are no charged particles in a plasma

Q2. Looking at the equation for the Debye length, how does it change when the temperature rises and when the density rises?

A It becomes shorter when the temperature rises and longer when the density rises B It becomes longer when the temperature rises and shorter when the density rises C It changes with neither temperature nor density D It always becomes longer with either temperature or density

Q3. By what physical relation are the Debye length and the plasma frequency linked?

A The two are entirely unrelated and arise from separate physics B The Debye length grows in proportion to the plasma frequency C They are linked through the electron thermal velocity, and during one plasma oscillation an electron travels roughly one Debye length by thermal motion D The higher the plasma frequency, the longer the Debye length must always be

Q4. Which combination is correct as the three conditions for an ionized gas to behave as a plasma?

A The Debye length is larger than the system size, the number of particles in the Debye sphere is small, and the plasma frequency is lower than the collision frequency B The Debye length is sufficiently smaller than the system size, the number of particles in the Debye sphere is large, and the electron plasma frequency is higher than the collision frequency with neutral particles C The temperature is near absolute zero, the density is near zero, and all charge is neutralized D The Debye length is smaller than the interparticle spacing, all particles are at rest, and the magnetic field is zero

Q5. What is the sheath, a representative example where quasineutrality breaks down?

A A layer of concentrated negative charge that forms in the center of the plasma B A region where charge is completely neutralized on scales much larger than the Debye length C A positive space-charge layer with a thickness of a few times the Debye length that forms at the boundary where the plasma contacts the wall D A steady state in which electrons and ions move at the same velocity

Related Topics

• Plasma Physics Overview - Introduction to plasma parameters
• Charged Particle Motion - Motion of particles in fields
• Divertor - Engineering of wall heat and particle handling involving the sheath
• MHD - Fluid description of plasmas
• Transport Phenomena - Particle and heat transport