PhysSandbox
Classical MechanicsWaves & SoundElectricity & MagnetismOptics & LightGravity & OrbitsLabs
🌙Astronomy & The Sky🌡️Thermodynamics🌍Biophysics, Fluids & Geoscience📐Math Visualization🔧Engineering🧪Chemistry

More from Electricity & Magnetism

Other simulators in this category — or see all 62.

View category →
NewUniversity / research

Langmuir Plasma Oscillations

Collective electron oscillations at ω_p = √(n_e e²/ε₀m_e): animated δn wave, Bohm–Gross dispersion ω(k), and ω_p vs density.

Launch Simulator
NewUniversity / research

Alfvén Wave on a Magnetic Field Line

Transverse MHD Alfvén wave: v_A = B/√(μ₀ρ), animated field-line displacement, standing modes and pulse reflections at fixed boundaries.

Launch Simulator
NewUniversity / research

Magnetic Reconnection (Sweet–Parker)

Toy Sweet–Parker X-point: opposite B-regions, current sheet, inflow v_in ≈ v_A/√S and outflow, Lundquist number S and sheet thickness δ/L.

Launch Simulator
NewUniversity / research

Parker Spiral Solar Wind

Rotating Sun plus radial solar wind: interplanetary B spirals with tan ψ = Ωr/v_sw, animated field lines (top view), ψ(r) and |B| falloff.

Launch Simulator
NewUniversity / research

Radiation Belts Drift Shells

Van Allen belts: gradient and curvature drift on dipole L-shells, mirror latitude, equatorial loss cone α_lc(L); electrons west, protons east.

Launch Simulator
NewSchool

Eddy Current Tube

Magnet fall: dv/dt = g in air vs g − k v in copper pipe (toy drag model).

Launch Simulator
PhysSandbox

Interactive physics, chemistry, and engineering simulators for students, teachers, and curious minds.

Physics

  • Classical Mechanics
  • Waves & Sound
  • Electricity & Magnetism

Science

  • Optics & Light
  • Gravity & Orbits
  • Astronomy & The Sky

More

  • Thermodynamics
  • Biophysics, Fluids & Geoscience
  • Math Visualization
  • Engineering
  • Chemistry

© 2026 PhysSandbox. Free interactive science simulators.

PrivacyTermsContact
Home/Electricity & Magnetism/Debye Shielding in a Plasma

Debye Shielding in a Plasma

In a plasma, mobile electrons respond to a localized electric potential and rearrange to partially cancel external fields — Debye shielding. For a test point charge Q immersed in a quasineutral electron plasma with density n_e and temperature T_e, linearized Poisson–Boltzmann theory gives the Yukawa potential φ(r) = (Q/4πε₀r) exp(−r/λ_D), where the Debye length λ_D = √(ε₀ k_B T_e / (n_e e²)) sets the exponential decay scale. At r ≫ λ_D the potential is strongly suppressed compared with the bare Coulomb 1/r law; at r ≪ λ_D the field resembles vacuum Coulomb behavior. The screening parameter κ = 1/λ_D grows with density and falls with temperature. This simulator plots screened versus bare φ(r), the radial electric field magnitude, a two-dimensional color map of φ around the charge with a circle at r = λ_D, and curves of λ_D versus T (fixed n) and versus n (fixed T). It uses SI units with sliders for T_e, log₁₀ n_e, and test charge Q. Ion dynamics, magnetic fields, and nonlinear shielding at very large potentials are omitted.

Who it's for: Undergraduate plasma physics, space physics, or electromagnetism students after Coulomb's law and before kinetic theory or Langmuir waves.

Key terms

  • Debye shielding
  • Debye length
  • Yukawa potential
  • Plasma quasineutrality
  • Screening
  • Poisson–Boltzmann
  • Electron temperature

Debye shielding

10000K
18
1nC

Yukawa screened potential φ(r) = (Q/4πε₀r) exp(−r/λ_D) with λ_D = √(ε₀k_BT/n_e e²). Mobile electrons rearrange to cancel the test field beyond a few Debye lengths — the basis of quasineutrality in plasma.

Measured values

Debye length λ_D6.90 µm
κ = 1/λ_D1.45e+5m⁻¹
n_e1.00×10¹⁸ m⁻³
T_e10000K

How it works

Debye shielding of a test charge in a plasma: Yukawa potential φ ∝ e^{-r/λ_D}/r, λ_D from T and n_e, 2D map and λ_D scaling curves.

Frequently asked questions

What is the physical meaning of λ_D?
It is the characteristic distance over which mobile charges can rearrange to cancel an imposed electric field. Beyond a few λ_D, the plasma is nearly neutral on average.
Why does φ(r) look like (1/r) exp(−r/λ_D)?
The Laplacian of a spherically symmetric Yukawa potential solves Poisson's equation with an exponential screening term from linearized charge response, giving exponential decay on top of the Coulomb geometry.
How do T and n affect shielding?
Higher n_e means more electrons available to screen → smaller λ_D. Higher T_e means less willingness to stay in the potential well → larger λ_D.
When does this linear model break down?
If eφ ≳ k_B T_e the response is nonlinear; very low density or collisionless kinetic effects also modify the profile. This page is the standard introductory Yukawa picture.