Skip to content

Plasma Physics: Brief Overview

A plasma screens electric fields over the Debye length:

λD=ε0kBTnee2\lambda_D = \sqrt{\frac{\varepsilon_0 k_B T}{n_e e^2}}

For ne=1018n_e = 10^{18} m3^{-3}, T=104T = 10^4 K: λD=7.4×105\lambda_D = 7.4 \times 10^{-5} m =74μ= 74\,\muM.

The plasma frequency:

ωp=nee2meε0\omega_p = \sqrt{\frac{n_e e^2}{m_e \varepsilon_0}}

For ne=1018n_e = 10^{18} m3^{-3}: ωp=5.64×1010\omega_p = 5.64 \times 10^{10} rad/s, fp=8.98f_p = 8.98 GHz. EM waves with ω<ωp\omega < \omega_p cannot propagate (evanescent).

Small displacements of the electron cloud create restoring forces, leading to Langmuir waves:

ωLangmuir=ωp(1+3kBT2mek2ωp2)1/2\omega_{\text{Langmuir} = \omega_p\left(1 + \frac{3k_BT}{2m_e}\frac{k^2}{\omega_p^2}\right)^{-1/2}}

At long wavelengths (k0k \to 0): ωωp\omega \to \omega_p (undamped). With ion motion: the ion-acoustic wave has ω2=k2cs2/(1+k2λD2)\omega^2 = k^2 c_s^2/(1 + k^2\lambda_D^2) where cs=kBT/mic_s = \sqrt{k_BT/m_i}.

Problem. A uniformly charged sphere of radius RR has total charge QQ. Find EE inside and outside.

Solution. Outside (r>Rr > R): EdA=Q/ε0    E4πr2=Q/ε0    E=Q4πε0r2\oint E \cdot dA = Q/\varepsilon_0 \implies E \cdot 4\pi r^2 = Q/\varepsilon_0 \implies E = \frac{Q}{4\pi\varepsilon_0 r^2}. Inside (r<Rr < R): enclosed charge =Q(r/R)3= Q(r/R)^3. E=Qr4πε0R3E = \frac{Qr}{4\pi\varepsilon_0 R^3}.

\blacksquare

Problem. An EM wave has E0=100V/mE_0 = 100 \mathrm{ V/m} in vacuum. Find the average Poynting vector magnitude.

Solution. S=E022μ0c=10022×4π×107×3×108=1000075413.3W/m2{\langle S \rangle = \frac{E_0^2}{2\mu_0 c} = \frac{100^2}{2 \times 4\pi \times 10^{-7} \times 3 \times 10^8} = \frac{10000}{754} \approx 13.3 \mathrm{ W/m}^2}.

\blacksquare

  • Confusing Gauss’s law applications. Gauss’s law is most useful for systems with high symmetry (spherical, cylindrical, planar). Fix: Choose a Gaussian surface matching the symmetry; the flux through the surface equals the enclosed charge divided by ε0\varepsilon_0.
  • Wrong Maxwell equation sign. Faraday’s law has a negative sign: ×E=Bt\nabla \times \vec{E} = -\frac{\partial \vec{B}}{\partial t}. Fix: The minus sign reflects Lenz’s law — the induced EMF opposes the change in flux.
  • Confusing D\vec{D}and E\vec{E}, H\vec{H}and B\vec{B}. D=ε0E+P\vec{D} = \varepsilon_0\vec{E} + \vec{P}; H=B/μ0M\vec{H} = \vec{B}/\mu_0 - \vec{M}. Fix: In vacuum: D=ε0E\vec{D} = \varepsilon_0\vec{E}, H=B/μ0\vec{H} = \vec{B}/\mu_0.
  • Maxwell’s equations: Gauss’s law, Gauss’s law for magnetism, Faraday’s law, Ampère-Maxwell law.
  • Gauss’s law: EdA=Qenc/ε0\oint \vec{E} \cdot d\vec{A} = Q_{\text{enc}}/\varepsilon_0.
  • EM waves: E0=cB0E_0 = cB_0; c=1/μ0ε0c = 1/\sqrt{\mu_0\varepsilon_0}; Poynting vector S=E×H/μ0\vec{S} = \vec{E} \times \vec{H}/\mu_0.
  • Boundary conditions: tangential EE and normal BB are continuous across interfaces.
TopicSiteLink
[Electromagnetism]A-LevelView
[Electromagnetism]IBView
[Electromagnetism]DSEView
[Electromagnetism]UniversityView
QuantityFormulaPhysical role
Debye lengthλD=ε0kBT/nee2\lambda_D = \sqrt{\varepsilon_0 k_B T / n_e e^2}Distance over which electric fields are screened
Plasma frequencyωp=nee2/meε0\omega_p = \sqrt{n_e e^2 / m_e \varepsilon_0}Natural oscillation frequency of electron gas
Langmuir waveω2=ωp2+3k2vth2\omega^2 = \omega_p^2 + 3k^2 v_{\rm th}^2Electrostatic wave in unmagnetised plasma
Ion-acoustic waveω2=k2cs2/(1+k2λD2)\omega^2 = k^2 c_s^2 / (1 + k^2\lambda_D^2)Low-frequency wave with ion inertia and electron pressure
Electron gyrofrequencyωce=eB/me\omega_{ce} = eB/m_eCyclotron frequency in magnetised plasma
  • Confusing Debye shielding with perfect neutrality. A plasma is quasineutral (ninen_i \approx n_e) on scales large compared to λD\lambda_D, but charge separation exists on Debye-length scales. Fix: Use λD\lambda_D as the scale below which individual charges matter.
  • Assuming all EM waves propagate in a plasma. Waves with ω<ωp\omega < \omega_p are evanescent — they decay exponentially. Fix: The cut-off condition is ω>ωp\omega > \omega_p for propagation; below ωp\omega_p the refractive index becomes imaginary.
  • Forgetting ion motion in low-frequency waves. The ion-acoustic wave requires mobile ions; at frequencies above ωpi\omega_{pi} (ion plasma frequency), ions cannot respond. Fix: Check whether ωωpi\omega \ll \omega_{pi} before using the ion-acoustic dispersion.
  • Treating Coulomb collisions as rare. While high-temperature plasmas are often collisionless, the collision frequency scales as T3/2T^{-3/2}; cold, dense plasmas can be collisional. Fix: Compare the mean free path to system size using νeineTe3/2\nu_{ei} \propto n_e T_e^{-3/2}.
  • Fusion energy (tokamaks): Magnetic confinement of deuterium-tritium plasmas at T108T \sim 10^8 K requires understanding of MHD stability, transport, and wave heating.
  • Space physics: The solar wind (ne107n_e \sim 10^7 m3^{-3}, T105T \sim 10^5 K) is a plasma that interacts with Earth’s magnetosphere, causing aurorae and geomagnetic storms.
  • Semiconductor processing: Low-temperature plasmas (Te104T_e \sim 10^4 K, ne1016n_e \sim 10^{16} m3^{-3}) are used for etching and deposition in microchip fabrication.
  • Radio astronomy: Pulsar signals propagate through the interstellar medium (ISM) plasma; dispersion measurements give the column density DM=nedl{\rm DM} = \int n_e\, dl.

13.6 Worked Example: Debye Length in the Solar Corona

Section titled “13.6 Worked Example: Debye Length in the Solar Corona”

Problem. The solar corona has ne1014n_e \approx 10^{14} m3^{-3} and T106T \approx 10^6 K. Compute the Debye length. How many electrons are in a Debye sphere?

Solution.

λD=ε0kBTnee2=(8.85×1012)(1.38×1023)(106)(1014)(1.6×1019)2\lambda_D = \sqrt{\frac{\varepsilon_0 k_B T}{n_e e^2}} = \sqrt{\frac{(8.85\times10^{-12})(1.38\times10^{-23})(10^6)}{(10^{14})(1.6\times10^{-19})^2}}

λD1.22×10282.56×1024=4.77×1056.9×103  m=6.9  mm\lambda_D \approx \sqrt{\frac{1.22\times10^{-28}}{2.56\times10^{-24}}} = \sqrt{4.77\times10^{-5}} \approx 6.9\times10^{-3}\;\mathrm{m} = 6.9\;\mathrm{mm}

The Debye sphere volume is 43πλD31.38×106\frac{4}{3}\pi\lambda_D^3 \approx 1.38\times10^{-6} m3^3, containing ND=ne43πλD31014×1.38×1061.4×108N_D = n_e \cdot \frac{4}{3}\pi\lambda_D^3 \approx 10^{14} \times 1.38\times10^{-6} \approx 1.4\times10^8 electrons. Since ND1N_D \gg 1, the corona satisfies the plasma criterion for collective behaviour.

\blacksquare

RegimeConditionKey behaviour
Debye shieldingr>λDr > \lambda_DElectric fields screened out
Plasma oscillationsωωp\omega \approx \omega_pCollective electron oscillation (Langmuir)
EM wave propagationω>ωp\omega > \omega_pWave propagates through plasma
EM wave cut-offω<ωp\omega < \omega_pWave is evanescent, reflected
Ion-acoustic wavesTeTiT_e \gg T_iSound-like waves with cs=kBTe/mic_s = \sqrt{k_BT_e/m_i}