advanced32 min studyLesson 291 in Physics

Debye shielding and the plasma parameter

Derive the screening length and separate weak coupling, quasineutrality, and response-time assumptions.

Before you begin

On this page

The question for today

What you’ll understand

Linearize the density response of specified mobile species.
Compute the number of electrons in an electron Debye sphere.
Identify where a quasineutral approximation fails.

A gas containing charged particles does not automatically behave like a smooth, neutral conducting fluid on every scale. Screening, statistical smoothness, and collective response impose different conditions. A useful plasma description states which species can rearrange, which observation scale is being resolved, and whether individual interactions are weak. These distinctions prevent a screening formula from becoming an unrestricted definition of all plasmas.

Screening is a density response

Consider electrons at temperature $T_e$ in kelvin, with initially uniform density $n_0$ and stationary singly charged ions. A weak positive potential attracts extra electrons. Their equilibrium Boltzmann density is $n_e=n_0\exp(e\phi/k_BT_e)$ , where $e$ is the positive elementary charge. Expanding the exponential requires $|e\phi|\ll k_BT_e$ ; this is a local condition on the potential energy, not merely a statement that the test charge is small.

Substitution into Poisson's equation gives a screened rather than a bare Coulomb potential. Define the electron Debye length using electrons alone. If ions also have time to achieve a Boltzmann response, their inverse squared screening length adds to the electron contribution. Equal electron and ion temperatures therefore shorten the combined length by $\sqrt{2}$ ; replacing the electron length silently would change the physical assumption about mobile species.

\lambda_{De}=\sqrt{\frac{\epsilon_0 k_BT_e}{n_0e^2}},\qquad (\nabla^2-\lambda_{De}^{-2})\phi=-\frac{\rho_{\rm ext}}{\epsilon_0},\qquad \phi(r)=\frac{Qe^{-r/\lambda_{De}}}{4\pi\epsilon_0r}

For classical isothermal electrons and immobile ions, the linearized potential outside a localized test source decays exponentially. Length is in metres, density in m⁻³, external charge density in C/m³, and potential in volts. The point-source expression applies only where linearization is valid, not arbitrarily close to the singular source.↔ Scroll wide equations horizontally.

Many particles and a large system are separate tests

Define $N_D$ as the mean electron count inside a sphere of radius $\lambda_{De}$ . A large value supports a smooth collective field with relatively small number fluctuations and weak pair correlations in the classical plasma being considered. Some texts call its reciprocal the plasma parameter, so state the convention. It is not enough to count particles in the entire apparatus: a large total population can still have few particles within its local screening volume.

Bulk quasineutrality also needs a macroscopic length $L$ much larger than the relevant screening length. Poisson's equation then allows a small fractional charge imbalance to sustain a finite electric field. Under a thermal-scale potential estimate, that imbalance scales as $(\lambda_{De}/L)^2$ . Quasineutrality means a controlled small imbalance, not exact cancellation imposed before deriving the electric field. A material sheath deliberately resolves the scales on which this approximation stops working.

N_D=\frac{4\pi}{3}n_0\lambda_{De}^3\gg1,\qquad \frac{|n_i-n_e|}{n_0}\sim\left(\frac{\lambda_{De}}L\right)^2\frac{|e\phi|}{k_BT_e}

The second relation is a scale estimate from Poisson's equation, not an exact identity or a bound. L is a potential-variation length in metres. Large N_D identifies the weakly coupled continuum regime used here; strongly coupled plasmas require other descriptions.↔ Scroll wide equations horizontally.

Screening has a clock and a domain

A static Boltzmann calculation describes a relaxed response. Electrons cannot rearrange instantaneously following an arbitrarily rapid perturbation. Their collective response rate is set by the electron plasma frequency, developed next. To observe weakly damped electron oscillations, electron collision frequencies should be small compared with that frequency. This dynamical requirement is different from the equilibrium statement that a Debye sphere contains many particles. Neither condition alone fixes the behaviour at a boundary.

The screened potential is particularly useful as a diagnostic of scale dependence. Far outside a screening cloud the test charge is strongly hidden; within it the Coulomb factor remains important. The mathematical point-charge solution cannot justify linear response near its divergent origin. Likewise, finding a Debye-scale electric field does not establish that the material has ceased being a plasma. It establishes that a quasineutral bulk model is insufficient for the particular structure being resolved.

Work it through

Check a small laboratory column

Take $n_0=10^{16}\,\mathrm{m^{-3}}$ , electron thermal energy $k_BT_e=1.60\times10^{-18}\,\mathrm J$ (10 eV), and column radius $L=1.0\,\mathrm{cm}$ . Use $\epsilon_0=8.85\times10^{-12}\,\mathrm{F/m}$ and $e=1.60\times10^{-19}\,\mathrm C$ . Find the electron Debye length and electron Debye-sphere population.

The squared length is $(8.85\times10^{-12})(1.60\times10^{-18})/[10^{16}(1.60\times10^{-19})^2]=5.53\times10^{-8}\,\mathrm{m^2}$ . Keep the density in SI units before taking a square root.
Thus $\lambda_{De}=2.35\times10^{-4}\,\mathrm m=0.235\,\mathrm{mm}$ and $\lambda_{De}/L=0.0235$ . The radius spans about 43 electron screening lengths, supporting a bulk-scale approximation away from the wall.
Using extra digits of the length, the sphere count is $(4\pi/3)10^{16}(2.352\times10^{-4})^3\simeq5.45\times10^5$ . This is comfortably above unity. It checks weak coupling in the stated classical model but does not measure collision rates or prove a specific equilibrium profile.

The result

The electron screening length is about 0.235 mm and its sphere contains about 545,000 electrons. Both spatial tests favour a smooth quasineutral interior; a wall layer still requires its own treatment.

The living notebook

Explore the model

Interactive lab

Interactive lab · sampling

A statistical analogy for a smooth density

Reset and set both controls to p = 0.50 and n = 20; predict the mean's variance. Reset and set both controls to p = 0.50 and n = 80; predict a fourfold smaller variance, then inspect the readouts. These are exact independent Bernoulli means, not counts in a correlated Debye cloud, and playback highlights outcomes rather than collecting measurements.

Exact binomial massNormal interval massHighlighted outcome / mean

Mean sample proportion: 0.5
Sample-mean variance: 0.0125
Uniform CDF error bound: ≤ 0.1062

Success probability p0.5Independent samples n20

t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forLarger independent populations reduce relative statistical variation; the lab does not calculate screening or Coulomb correlations.

Pause & reason

Check your understanding

Interactive checks need JavaScript. You can still read the explanation.

Read the checkpoint explanation

Each mobile species contributes positively to the inverse squared screening length. Equal temperatures and densities give two equal contributions, shortening the combined length by √2. This is a change of response model, not a change in the electron-only definition.

Your turn

Try it yourself

A second classical plasma has four times the electron density and the same electron temperature as the worked example. Its radius remains 1 cm. Find its electron Debye length, Debye-sphere electron count, and the direction of the change in scale separation.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The length scales as the inverse square root of density; then recompute density times length cubed.

Compare with the solution

The length halves to 0.1175 mm. The sphere population is multiplied by $4(1/2)^3=1/2$ , giving approximately $2.72\times10^5$ . The radius now spans about 85 electron Debye lengths, so bulk spatial separation improves even though the local sphere population falls. Both sphere populations remain large. Increasing density at fixed temperature therefore improves this geometric test while moving the classical coupling test in the opposite direction; the two conditions are logically independent.

Keep the idea

Take it with you

Specify mobile species before assigning a Debye screening length.
Large Debye-sphere population and small Debye-to-system length ratio test different approximations.
Quasineutrality retains small charge imbalances and excludes resolved sheaths.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.