Plasma Physics Course

Total: 8 parts covering plasma fundamentals through advanced topics

Part I, Chapter 1

Basic Properties & Definition

What is plasma and when does matter behave collectively?

1.1 Definition of Plasma

A plasma is a quasi-neutral gas of charged and neutral particles that exhibits collective behavior. This definition contains three key criteria:

Three Criteria for Plasma

  1. Sufficient ionization: Enough charged particles for collective effects
  2. Quasi-neutrality: Equal positive and negative charge densities on macroscopic scales
  3. Collective behavior: Long-range Coulomb interactions dominate over binary collisions

1.1.1 Ionization Mechanisms

Plasmas are created when sufficient energy is supplied to ionize neutral atoms. The ionization energy for hydrogen is Ei = 13.6 eV. Common ionization mechanisms include:

• Thermal Ionization

At high temperatures, thermal collisions provide sufficient energy. The Saha equationrelates ionization fraction to temperature:

$$\frac{n_i n_e}{n_n} = 2.4 \times 10^{21} \frac{T^{3/2}}{g_i/g_n} e^{-E_i/k_B T} \quad \text{[m}^{-3}\text{]}$$

where gi and gn are statistical weights. For T ≈ Ei/kB ≈ 105 K, appreciable ionization occurs.

Derivation of the Saha Equation

Consider the ionization equilibrium reaction A ⇌ A⁺ + e⁻. In thermal equilibrium, the chemical potentials balance:

$$\mu_A = \mu_{A^+} + \mu_e$$

Step 1. For an ideal gas species α, the chemical potential is related to its single-particle partition function Zα via:

$$\mu_\alpha = -k_B T \ln\!\left(\frac{Z_\alpha}{N_\alpha}\right), \quad Z_\alpha = Z_\alpha^{\text{trans}} \cdot Z_\alpha^{\text{int}}$$

Step 2. The translational partition function for a free particle of mass mαin volume V is:

$$Z_\alpha^{\text{trans}} = V\left(\frac{2\pi m_\alpha k_B T}{h^2}\right)^{3/2} = \frac{V}{\Lambda_\alpha^3}$$

where Λα = h/√(2πmαkBT) is the thermal de Broglie wavelength. The internal partition function reduces to the statistical weight of the ground state at moderate temperatures: Zintα ≈ gα.

Step 3. Substituting chemical potentials into the equilibrium condition and using nα = Nα/V, the law of mass action gives:

$$\frac{n_i n_e}{n_n} = \frac{g_i}{g_n} \cdot \frac{Z_e^{\text{int}}}{\Lambda_e^3} \cdot e^{-E_i / k_B T}$$

The electron has spin degeneracy Zinte = ge = 2, and the ion-to-atom mass ratio is ≈ 1 so their translational contributions cancel. Inserting:

$$\frac{n_i n_e}{n_n} = \frac{2 g_i}{g_n}\left(\frac{2\pi m_e k_B T}{h^2}\right)^{3/2} e^{-E_i/k_B T}$$

Step 4. Evaluating the prefactor numerically with T in Kelvin:

$$\left(\frac{2\pi m_e k_B}{h^2}\right)^{3/2} = 2.415 \times 10^{21} \text{ m}^{-3}\text{K}^{-3/2}$$

This yields the Saha equation. The exponential factor e−Ei/kBT shows ionization is activated — it rises steeply once kBT reaches a fraction of Ei. Note that significant ionization occurs already at T ∼ Ei/20 because of the large T3/2prefactor from phase-space (translational degrees of freedom of the freed electron).

• Electron Impact Ionization

Energetic electrons ionize neutrals via collisions. The ionization rate coefficient is:

$$\langle \sigma v \rangle \approx 10^{-13} \left(\frac{T_e}{E_i}\right)^{1/2} e^{-E_i/T_e} \quad \text{[m}^3\text{/s]}$$

• Photoionization

Photons with energy hν > Ei ionize atoms directly. Important in stellar atmospheres and the ionosphere. Photoionization cross-section scales as:

$$\sigma_{\text{photo}} \sim 10^{-22} \left(\frac{\lambda}{100 \text{ nm}}\right)^3 \quad \text{[m}^2\text{]}$$

1.1.2 Ionization Fraction

The degree of ionization α is the fraction of atoms ionized:

$$\alpha = \frac{n_i}{n_i + n_n}$$

Plasmas range from weakly ionized (α ≪ 1, dominated by neutral collisions) tofully ionized (α ≈ 1, no neutrals present):

Classification by Ionization

  • • Weakly ionized (α < 10−4): Gas discharges, ionosphere F-region
  • • Partially ionized (10−4 < α < 0.9): Tokamak edge, solar chromosphere
  • • Fully ionized (α > 0.9): Tokamak core, solar corona, stellar interiors

1.1.3 Collective Behavior

Unlike neutral gases where molecules interact via short-range collisions, plasmas exhibitcollective behavior because the Coulomb force is long-ranged. A single charged particle influences many neighbors simultaneously.

The Coulomb potential energy between two charges separated by distance r is:

$$U(r) = \frac{e^2}{4\pi\epsilon_0 r} \approx \frac{1.44 \text{ eV·nm}}{r}$$

The mean interparticle spacing is n−1/3. For a typical laboratory plasma with ne = 1019 m−3 (n−1/3 ≈ 5 nm) and Te = 10 eV, the potential energy at nearest-neighbor distance is comparable to the kinetic energy:

$$\frac{U(n^{-1/3})}{k_B T} = \frac{e^2}{4\pi\epsilon_0 k_B T n^{-1/3}} \sim 0.3$$

This shows that Coulomb interactions significantly affect particle dynamics. However, because the plasma is quasi-neutral, the potential is shielded at distances beyond the Debye length, making the effective interaction range finite.

The Fourth State of Matter

Plasmas comprise ~99% of the visible universe. From neon signs (1016 m−3) to stellar cores (1032 m−3), plasma physics spans 16 orders of magnitude in density and 8 orders in temperature!

1.2 Debye Shielding

Consider a test charge q placed in a plasma. Mobile charges will rearrange to screen the electric field. For an electron-ion plasma with equal temperatures Te = Ti = T, the potential around the test charge is:

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

where λD is the Debye length, the characteristic shielding distance:

$$\boxed{\lambda_D = \sqrt{\frac{\epsilon_0 k_B T}{n_e e^2}}} = \sqrt{\frac{\epsilon_0 k_B T_e}{n_e e^2} + \frac{\epsilon_0 k_B T_i}{n_i Z^2 e^2}}$$

1.2.1 Detailed Derivation from Poisson-Boltzmann

Starting with Poisson's equation:

$$\nabla^2 \phi = -\frac{\rho}{\epsilon_0} = -\frac{e}{\epsilon_0}(Z n_i - n_e)$$

Assuming the plasma responds in thermal equilibrium, particles follow Boltzmann distributions:

$$n_e(r) = n_0 \exp\left(\frac{e\phi}{k_B T_e}\right), \quad n_i(r) = n_0 \exp\left(\frac{-Ze\phi}{k_B T_i}\right)$$

For |eφ| ≪ kBT (weak potential), expand to first order:

$$n_e \approx n_0\left(1 + \frac{e\phi}{k_B T_e}\right), \quad n_i \approx n_0\left(1 - \frac{Ze\phi}{k_B T_i}\right)$$

Substituting into Poisson's equation:

$$\nabla^2 \phi = \frac{e n_0}{\epsilon_0}\left[\left(1 + \frac{e\phi}{k_B T_e}\right) - Z\left(1 - \frac{Ze\phi}{k_B T_i}\right)\right]$$

Using quasi-neutrality (n0 = Zni0), the constant terms cancel:

$$\nabla^2 \phi = \frac{e n_0}{\epsilon_0}\left(\frac{e}{k_B T_e} + \frac{Z^2 e}{k_B T_i}\right)\phi = \frac{\phi}{\lambda_D^2}$$

where the generalized Debye length for multiple species is:

$$\frac{1}{\lambda_D^2} = \sum_\alpha \frac{n_\alpha q_\alpha^2}{\epsilon_0 k_B T_\alpha} = \frac{n_e e^2}{\epsilon_0 k_B T_e} + \frac{n_i Z^2 e^2}{\epsilon_0 k_B T_i}$$

For spherical symmetry with boundary conditions φ(r→0) = q/(4πε0r) and φ(r→∞) = 0:

$$\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{d\phi}{dr}\right) = \frac{\phi}{\lambda_D^2}$$

Solution:

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

1.2.2 Numerical Values and Scaling

For a hydrogen plasma with Te = Ti = T:

$$\lambda_D \approx 7.43 \times 10^3 \sqrt{\frac{T[\text{eV}]}{n_e[\text{m}^{-3}]}} \quad \text{meters}$$

Example: Debye Lengths in Nature

  • • Tokamak core: ne = 1020 m−3, Te = 10 keV → λD ≈ 70 μm
  • • Gas discharge: ne = 1016 m−3, Te = 2 eV → λD ≈ 100 μm
  • • Ionosphere: ne = 1012 m−3, Te = 0.1 eV → λD ≈ 7 mm
  • • Solar corona: ne = 1015 m−3, Te = 100 eV → λD ≈ 7 cm
  • • Interstellar medium: ne = 106 m−3, Te = 1 eV → λD ≈ 7 m

1.2.3 Total Shielding Charge

The total charge in the shielding cloud can be calculated by integrating the charge density:

$$Q_{\text{shield}} = \int_0^\infty 4\pi r^2 \rho(r)\, dr = -q$$

This confirms that the shielding cloud exactly neutralizes the test charge. The exponential falloff ensures convergence of the integral, unlike the bare Coulomb potential.

1.2.4 Validity of Linear Approximation

The linearization |eφ| ≪ kBT is valid for:

$$r \gg r_c = \frac{e^2}{4\pi\epsilon_0 k_B T}$$

where rc is the distance where potential energy equals thermal energy. For T = 10 eV, rc ≈ 1.4 nm. Since λD ≫ rc in most plasmas, linear theory is valid except very close to the test charge.

Physical Interpretation

The Debye length represents the fundamental length scale in plasma physics. For r ≪ λD, individual particle interactions dominate. For r ≫ λD, collective effects dominate and the plasma is quasi-neutral. The exponential decay e−r/λDeffectively converts the long-range Coulomb interaction into a short-range Yukawa potential.

1.2.5 Time-Dependent Shielding

The above derivation assumes instantaneous response. For time-varying fields with frequency ω, shielding is incomplete when ω ≳ ωpe. The dynamic Debye length is:

$$\lambda_D(\omega) = \lambda_D \sqrt{1 + \frac{\omega^2}{\omega_{pe}^2}}$$

For ω ≫ ωpe, shielding breaks down and electromagnetic waves can propagate through the plasma.

1.3 Plasma Frequency

If charge neutrality is perturbed, the plasma responds by oscillating at the plasma frequency. This fundamental timescale characterizes how quickly a plasma can respond to disturbances.

1.3.1 Simple Derivation: Cold Plasma Oscillation

Consider a one-dimensional plasma slab where electrons are displaced by δx from a uniform ion background. This creates a charge separation and electric field:

$$\sigma = n_e e \delta x \quad \Rightarrow \quad E = \frac{\sigma}{\epsilon_0} = \frac{n_e e \delta x}{\epsilon_0}$$

The restoring force on electrons (ignoring thermal pressure) gives Newton's equation:

$$m_e \frac{d^2 \delta x}{dt^2} = -e E = -\frac{n_e e^2}{\epsilon_0} \delta x$$

This is simple harmonic motion with the electron plasma frequency:

$$\boxed{\omega_{pe} = \sqrt{\frac{n_e e^2}{\epsilon_0 m_e}} \approx 56.4 \sqrt{n_e[\text{m}^{-3}]} \text{ rad/s}}$$

Or in frequency units:

$$f_{pe} = \frac{\omega_{pe}}{2\pi} \approx 8.98 \sqrt{n_e[\text{m}^{-3}]} \text{ Hz}$$

1.3.2 Ion Plasma Frequency

Similarly, ions oscillate at the ion plasma frequency:

$$\omega_{pi} = \sqrt{\frac{n_i Z^2 e^2}{\epsilon_0 m_i}} = \sqrt{\frac{m_e}{m_i}} \omega_{pe}$$

For hydrogen, mi/me ≈ 1836, so:

$$\omega_{pi} \approx \frac{\omega_{pe}}{43} \ll \omega_{pe}$$

This mass ratio creates timescale separation: electrons respond much faster than ions, leading to distinct high-frequency (electron) and low-frequency (ion) phenomena.

1.3.3 Dispersion Relation and EM Wave Propagation

The plasma frequency determines whether electromagnetic waves can propagate. We derive the dispersion relation from Maxwell's equations coupled to the cold electron fluid.

Derivation of EM Dispersion Relation

Step 1. Start with Maxwell's equations (no free charges, transverse wave):

$$\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}, \qquad \nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t}$$

Step 2. The cold electron fluid equation of motion (linearized, ions immobile):

$$m_e \frac{\partial \mathbf{v}_e}{\partial t} = -e\mathbf{E} \quad \Rightarrow \quad \mathbf{J} = -n_e e \mathbf{v}_e$$

Step 3. Take the curl of Faraday's law and substitute Ampère's law:

$$\nabla \times (\nabla \times \mathbf{E}) = -\frac{\partial}{\partial t}(\nabla \times \mathbf{B}) = -\mu_0 \frac{\partial \mathbf{J}}{\partial t} - \frac{1}{c^2}\frac{\partial^2 \mathbf{E}}{\partial t^2}$$

Using ∇ × (∇ × E) = −∇²E (for transverse waves with ∇·E = 0):

$$\nabla^2 \mathbf{E} - \frac{1}{c^2}\frac{\partial^2 \mathbf{E}}{\partial t^2} = \mu_0 \frac{\partial \mathbf{J}}{\partial t}$$

Step 4. From the equation of motion, ∂J/∂t = (nee²/me)E. Substituting:

$$\nabla^2 \mathbf{E} - \frac{1}{c^2}\frac{\partial^2 \mathbf{E}}{\partial t^2} = \frac{\mu_0 n_e e^2}{m_e}\mathbf{E} = \frac{\omega_{pe}^2}{c^2}\mathbf{E}$$

Step 5. Assume plane-wave solution E = E0 ei(k·r − ωt), so ∇² → −k² and ∂²/∂t² → −ω²:

$$-k^2 + \frac{\omega^2}{c^2} = \frac{\omega_{pe}^2}{c^2}$$

Rearranging yields the EM wave dispersion relation in an unmagnetized plasma:

$$\boxed{\omega^2 = \omega_{pe}^2 + c^2 k^2}$$

For propagating waves (real k), we require:

$$\boxed{\omega > \omega_{pe}}$$

EM waves with ω < ωpe are evanescent—they cannot propagate and decay exponentially. The plasma acts as a high-pass filter. The wave becomes evanescent with spatial decay rate κ = √(ωpe² − ω²)/c, giving a skin depth δ = 1/κ.

Application: Ionospheric Radio Reflection

The ionosphere has ne ∼ 1012 m−3, giving fpe ∼ 9 MHz. AM radio waves (f < 1.6 MHz) can penetrate and reflect off higher layers, enabling long-distance communication. FM radio (f ∼ 100 MHz > fpe) passes through, limiting range to line-of-sight.

1.3.4 Refractive Index

The phase velocity vp = ω/k defines the refractive index. We derive it directly from the dispersion relation.

Derivation of Refractive Index and Critical Density

Step 1. From the dispersion relation ω² = ωpe² + c²k², solve for k:

$$k = \frac{1}{c}\sqrt{\omega^2 - \omega_{pe}^2}$$

Step 2. The phase velocity is vp = ω/k, and the refractive index is n = c/vp = ck/ω:

$$n = \frac{ck}{\omega} = \frac{c}{\omega} \cdot \frac{\sqrt{\omega^2 - \omega_{pe}^2}}{c} = \sqrt{1 - \frac{\omega_{pe}^2}{\omega^2}}$$

Step 3. The cutoff condition n = 0 (wave cannot propagate) occurs when ω = ωpe. Substituting ωpe² = nee²/(ε₀me), the critical density is the electron density at which ωpe = ω:

$$n_c = \frac{\epsilon_0 m_e \omega^2}{e^2} = \frac{4\pi^2 \epsilon_0 m_e f^2}{e^2}$$

Step 4. Note that vp = c/n > c (superluminal phase velocity), while the group velocity carries energy and information:

$$v_g = \frac{d\omega}{dk} = \frac{c^2 k}{\omega} = c\sqrt{1 - \frac{\omega_{pe}^2}{\omega^2}} < c$$

The product vp · vg = c² satisfies causality. The group velocity vanishes at the cutoff ω = ωpe, meaning the wave energy is reflected.

$$\boxed{n = \frac{c}{v_p} = \sqrt{1 - \frac{\omega_{pe}^2}{\omega^2}} = \sqrt{1 - \frac{n_e}{n_c}}}$$

where the critical density is:

$$n_c = \frac{\epsilon_0 m_e \omega^2}{e^2} \approx 1.1 \times 10^{10} \left[\frac{f[\text{Hz}]}{1\text{ Hz}}\right]^2 \text{ m}^{-3}$$

For laser fusion (λ = 351 nm, f = 8.5 × 1014 Hz):

$$n_c \approx 8 \times 10^{27} \text{ m}^{-3}$$

The laser cannot penetrate plasma regions where ne > nc, depositing energy at the critical density surface.

1.3.5 Plasma Dielectric Function

The plasma response to EM waves is described by the dielectric function, which we derive from the electron equation of motion.

Derivation of Dielectric Function

Step 1. Apply an oscillating electric field E = E0e−iωtto the cold electron fluid. Newton's law gives:

$$m_e \frac{d\mathbf{v}_e}{dt} = -e\mathbf{E} \quad \Rightarrow \quad -i\omega m_e \mathbf{v}_e = -e\mathbf{E} \quad \Rightarrow \quad \mathbf{v}_e = \frac{e\mathbf{E}}{i\omega m_e}$$

Step 2. The current density is J = −neeve, defining the AC conductivity σ(ω):

$$\mathbf{J} = -n_e e \mathbf{v}_e = \frac{n_e e^2}{i\omega m_e}\mathbf{E} = \sigma(\omega)\mathbf{E}, \quad \sigma(\omega) = \frac{i n_e e^2}{\omega m_e}$$

Step 3. In a dielectric medium, the displacement field is D = ε₀ε(ω)E. Ampère's law in the frequency domain reads:

$$\nabla \times \mathbf{B} = \mu_0\mathbf{J} + \mu_0\epsilon_0 \frac{\partial \mathbf{E}}{\partial t} = \mu_0\epsilon_0\left[\frac{\sigma}{-i\omega\epsilon_0} + 1\right](-i\omega)\mathbf{E}$$

Identifying the dielectric function as ε(ω) = 1 + iσ/(ωε₀):

$$\epsilon(\omega) = 1 + \frac{i\sigma(\omega)}{\omega \epsilon_0} = 1 + \frac{i}{\omega\epsilon_0}\cdot\frac{in_e e^2}{\omega m_e} = 1 - \frac{n_e e^2}{\epsilon_0 m_e \omega^2}$$

Recognizing ωpe² = nee²/(ε₀me), we obtain the cold-plasma dielectric function. Note the purely imaginary conductivity — the plasma is a reactive (lossless) medium; adding collisions via a friction term −meνve gives ε(ω) = 1 − ωpe²/[ω(ω + iν)], introducing absorption.

$$\boxed{\epsilon(\omega) = 1 - \frac{\omega_{pe}^2}{\omega^2}}$$

When ε < 0 (ω < ωpe), waves are evanescent. When ε > 0 (ω > ωpe), waves propagate. The zero crossing ε = 0 at ω = ωpe defines the cutoff.

Typical Plasma Frequencies

Plasma Typene [m−3]fpe
Interstellar medium10690 kHz
Solar wind (1 AU)107280 kHz
Ionosphere10129 MHz
Solar corona101490 MHz
Gas discharge1016900 MHz
Lab plasma10189 GHz
Tokamak core102090 GHz
Laser target10289 × 106 GHz

1.3.6 Plasma Period and Timescale Ordering

The plasma period τpe = 2π/ωpe sets the fundamental timescale:

$$\tau_{pe} = \frac{1}{f_{pe}} \approx \frac{0.11}{\sqrt{n_e[\text{m}^{-3}]}} \text{ seconds}$$

Plasma phenomena separate into distinct timescale regimes:

$$\omega_{ce} \gtrsim \omega_{pe} \gg \omega_{ci} \sim \omega_{pi} \gg \nu_c$$

where ωce,ci are cyclotron frequencies and νc is the collision frequency. This hierarchy enables reduced descriptions on different timescales.

1.4 Plasma Parameter

A crucial dimensionless quantity characterizing plasma behavior is the number of particles in a Debye sphere (volume VD = (4π/3)λD3):

$$N_D = n_e \frac{4\pi}{3} \lambda_D^3 = \frac{4\pi}{3} n_e \left(\frac{\epsilon_0 k_B T_e}{n_e e^2}\right)^{3/2}$$

Simplifying:

$$\boxed{N_D = \frac{4\pi}{3} \frac{(\epsilon_0 k_B T_e)^{3/2}}{e^3 \sqrt{n_e}} = 1.72 \times 10^9 \frac{T_e^{3/2}[\text{eV}]}{n_e^{1/2}[\text{m}^{-3}]}}$$

Note the scaling: ND ∝ T3/2/√n. Higher temperature or lower density favors collective behavior.

1.4.1 Physical Interpretation

ND measures the relative importance of collective vs. binary interactions:

• ND ≫ 1 (Weakly Coupled Plasma)

Many particles shield each test charge → collective electromagnetic fields dominate → plasma behavior. Individual Coulomb collisions are rare events. Kinetic theory and fluid descriptions apply.

• ND ∼ 1 (Strongly Coupled Plasma)

Few particles available for shielding → particle correlations important → liquid-like behavior. Examples: ultracold plasmas, white dwarf interiors, warm dense matter.

• ND ≪ 1 (Not Really a Plasma)

Shielding breaks down → gas-like behavior with binary collisions dominating. The system doesn't exhibit true plasma collective phenomena.

1.4.2 Plasma Coupling Parameter

The inverse quantity is the plasma coupling parameter Γ:

$$\Gamma = \frac{1}{N_D} = \frac{U_{\text{pot}}}{U_{\text{kin}}} = \frac{(e^2/4\pi\epsilon_0 a)}{k_B T}$$

where a = n−1/3 is the mean interparticle spacing. More precisely, using the Wigner-Seitz radius aws = (3/4πn)1/3:

$$\Gamma = \frac{e^2}{4\pi\epsilon_0 a_{\text{ws}} k_B T} = \frac{e^2}{4\pi\epsilon_0 k_B T}\left(\frac{4\pi n}{3}\right)^{1/3}$$

Numerically:

$$\Gamma \approx 2.3 \times 10^{-3} \frac{n_e^{1/3}[\text{m}^{-3}]}{T[\text{eV}]}$$

1.4.3 Criterion for Plasma Behavior

For ideal plasma (weakly coupled) behavior, we require:

$$\boxed{\Gamma \ll 1 \quad \Leftrightarrow \quad N_D \gg 1}$$

This ensures collective effects dominate over two-body correlations. Most fusion and space plasmas satisfy Γ ≲ 0.01, justifying kinetic and fluid theories.

1.4.4 Strongly Coupled Plasmas

When Γ ≳ 1, the plasma enters the strongly coupled regime where particle correlations are essential:

Phase Transitions in Strongly Coupled Plasmas

  • • Γ ≈ 2: Liquid-like correlations develop
  • • Γ ≈ 170: Crystallization into Coulomb crystal (BCC lattice)
  • • Γ > 170: Solid phase with phonon modes

Strongly coupled plasmas occur in:

  • Dusty plasmas: Micron-sized charged grains form visible Coulomb crystals
  • White dwarf interiors: Γ ∼ 10−100, electron fluid exhibits liquid/solid behavior
  • Warm dense matter: Matter at solid density but T ∼ 1 eV (inertial fusion, planetary cores)
  • Ultracold neutral plasmas: T ∼ 1 K, created by photoionizing laser-cooled atoms

1.4.5 Quantum Effects

At high density or low temperature, quantum effects become important. The quantum parametercompares thermal de Broglie wavelength to interparticle spacing:

$$\Theta = \frac{\lambda_{\text{dB}}}{a} = \frac{\hbar}{a\sqrt{2m k_B T}} = \frac{n^{1/3}\hbar}{\sqrt{2m k_B T}}$$

When Θ ≳ 1, quantum effects (Pauli exclusion, zero-point motion) are significant:

• Classical Regime (Θ ≪ 1)

Thermal energy dominates quantum effects. Most fusion and astrophysical plasmas.

• Degenerate Regime (Θ ≫ 1)

Quantum statistics dominate. Fermi-Dirac distribution applies. Examples: white dwarf electrons (ne ∼ 1036 m−3, T ∼ 107 K, Θe ∼ 100).

The Fermi temperature TF = EF/kB is the temperature scale below which quantum degeneracy dominates.

Derivation of Fermi Energy and Temperature

Step 1. For a free electron gas, the density of states per unit volume is:

$$g(\varepsilon) = \frac{1}{2\pi^2}\left(\frac{2m_e}{\hbar^2}\right)^{3/2}\sqrt{\varepsilon}$$

(including spin degeneracy factor 2). This counts the number of quantum states in a thin energy shell [ε, ε + dε] per unit volume.

Step 2. At T = 0 (complete degeneracy), all states up to the Fermi energy EFare filled. The electron density is:

$$n_e = \int_0^{E_F} g(\varepsilon)\, d\varepsilon = \frac{1}{2\pi^2}\left(\frac{2m_e}{\hbar^2}\right)^{3/2} \int_0^{E_F}\sqrt{\varepsilon}\, d\varepsilon$$

Step 3. Evaluating the integral ∫₀EF √ε dε = (2/3) EF3/2:

$$n_e = \frac{1}{3\pi^2}\left(\frac{2m_e E_F}{\hbar^2}\right)^{3/2}$$

Step 4. Solving for EF:

$$E_F = \frac{\hbar^2}{2m_e}(3\pi^2 n_e)^{2/3}$$

Defining TF = EF/kB and inserting constants gives:

$$\boxed{T_F = \frac{\hbar^2}{2m_e k_B}(3\pi^2 n_e)^{2/3} \approx 5.9 \times 10^{-11}\, n_e^{2/3}[\text{m}^{-3}] \quad \text{[K]}}$$

For T ≪ TF, the plasma is degenerate: electrons fill states up to EFwith occupation number ≈ 1, and the electron pressure is the quantum degeneracy pressure Pe = (2/5)neEF, independent of temperature.

Example: Phase Diagram of Matter

The (Γ, Θ) plane provides a unified phase diagram for all matter:

  • • Γ ≪ 1, Θ ≪ 1: Classical ideal plasma (tokamaks, solar corona)
  • • Γ ≪ 1, Θ ≫ 1: Quantum ideal plasma (stellar interiors)
  • • Γ ≫ 1, Θ ≪ 1: Classical strongly coupled (dusty plasmas, ultracold plasmas)
  • • Γ ≫ 1, Θ ≫ 1: Quantum strongly coupled (white dwarfs, neutron star crusts)

1.5 Quasi-Neutrality

On length scales L ≫ λD, a plasma is quasi-neutral:

$$\boxed{n_e \approx Z n_i \quad \text{for } L \gg \lambda_D}$$

This doesn't mean the plasma is perfectly neutral—small deviations create the electric fields that govern plasma dynamics. The charge imbalance is typically:

$$\left|\frac{n_e - Z n_i}{n_e}\right| \sim \frac{\lambda_D}{L}$$

1.5.1 Physical Mechanism

Quasi-neutrality is enforced dynamically by the strong restoring forces that arise from charge separation. If a region develops net charge density ρ ≠ 0:

$$\nabla \cdot \mathbf{E} = \frac{\rho}{\epsilon_0} \quad \Rightarrow \quad E \sim \frac{\rho L}{\epsilon_0}$$

This electric field accelerates electrons with force F = −eE, creating a restoring current that neutralizes the charge on timescale ∼ τpe. Any attempted charge separation excitesplasma oscillations at ωpe, which rapidly restore quasi-neutrality.

1.5.2 Quantitative Estimate

Consider a charge imbalance parameter:

$$\delta = \frac{|n_e - Z n_i|}{n_e}$$

The electric field from this imbalance (over length L) is:

$$E \sim \frac{e n_e \delta L}{\epsilon_0}$$

The corresponding electrostatic potential energy is:

$$U_E \sim e E L = \frac{e^2 n_e \delta L^2}{\epsilon_0}$$

Comparing to thermal energy kBT, quasi-neutrality (δ ≪ 1) holds when:

$$\frac{e^2 n_e L^2}{\epsilon_0 k_B T} \gg 1 \quad \Rightarrow \quad L \gg \lambda_D$$

In this limit, the energetic cost of charge separation is so large that:

$$\delta \sim \left(\frac{\lambda_D}{L}\right)^2$$

1.5.3 Numerical Example: Laboratory Plasma

Consider a plasma with ne = 1018 m−3, Te = 10 eV:

$$\lambda_D \approx 7.4 \times 10^3 \sqrt{\frac{10}{10^{18}}} \approx 74 \,\mu\text{m}$$

For a system size L = 1 m:

$$\delta \sim \frac{\lambda_D}{L} = \frac{74 \times 10^{-6}}{1} \approx 10^{-4}$$

This means |ne − ni|/ne ∼ 0.01%, yet this tiny imbalance produces electric fields of order:

$$E \sim \frac{e n_e \delta L}{\epsilon_0} \sim 10 \text{ V/m}$$

These "small" fields drive significant plasma flows and currents!

1.5.4 Implications for Plasma Theory

The quasi-neutrality approximation simplifies Poisson's equation:

$$\nabla^2 \phi = -\frac{e}{\epsilon_0}(Z n_i - n_e) \approx 0 \quad \Rightarrow \quad n_e \approx Z n_i$$

This reduces the number of independent variables in plasma fluid equations. However, quasi-neutrality is an approximation, not an exact constraint. At boundaries (sheaths, double layers) and on scales ∼ λD, charge separation is important.

1.5.5 Breakdown of Quasi-Neutrality

Quasi-neutrality breaks down in several important situations:

• Plasma Sheaths

Near walls, charge separation over ∼few λD creates potential drops of ∼kBTe/e that confine electrons and accelerate ions. Essential for wall-plasma interactions.

• Double Layers

Localized potential structures with strong charge separation. Observed in auroral acceleration regions and laboratory plasmas.

• Electrostatic Waves

Plasma waves involve oscillating charge density perturbations. For wavelength k−1 ∼ λD, charge separation is significant (Langmuir waves, ion acoustic waves).

• Relativistic Plasmas

In ultra-high-intensity laser-plasma interactions or pulsar magnetospheres, strong EM fields can drive significant charge separation.

Ambipolar Diffusion

Quasi-neutrality leads to ambipolar diffusion: electrons and ions diffuse together to maintain charge neutrality. We derive the ambipolar diffusion coefficient from the coupled fluid equations.

Derivation

Step 1. The particle flux for each species in steady state with friction:

$$\Gamma_e = -D_e \nabla n - n\mu_e \mathbf{E}, \qquad \Gamma_i = -D_i \nabla n + n\mu_i \mathbf{E}$$

where Dα = kBTα/(mανα) is the diffusion coefficient and μα = |qα|/(mανα) is the mobility (Einstein relation: Dα = μαkBTα/|qα|).

Step 2. Quasi-neutrality requires equal fluxes: Γe = Γi = Γ. Setting them equal:

$$-D_e \nabla n - n\mu_e E = -D_i \nabla n + n\mu_i E$$

Step 3. Solve for the ambipolar electric field:

$$E = \frac{D_i - D_e}{\mu_i + \mu_e}\frac{\nabla n}{n}$$

Since De ≫ Di, the field points inward (retards electrons, accelerates ions).

Step 4. Substitute E back into either flux equation:

$$\Gamma = -\frac{D_i \mu_e + D_e \mu_i}{\mu_e + \mu_i}\nabla n = -D_a \nabla n$$

Step 5. Using μe ≫ μi (since me ≪ mi) and the Einstein relation De/μe = kBTe/e:

$$D_a \approx \frac{D_i \mu_e}{\mu_e} + \frac{D_e \mu_i}{\mu_e} \approx D_i + D_i\frac{T_e}{T_i} = D_i\left(1 + \frac{T_e}{T_i}\right)$$
$$\boxed{D_a = \frac{D_i \mu_e + D_e \mu_i}{\mu_e + \mu_i} \approx D_i \left(1 + \frac{T_e}{T_i}\right)}$$

The ambipolar rate is controlled by the slower species (ions) but enhanced by the electron temperature ratio. For Te = Ti, Da ≈ 2Di. For Te ≫ Ti (common in gas discharges), Da ≈ DiTe/Ti ≫ Di, significantly enhancing ion transport.

1.6 Summary: Three Plasma Criteria

For a gas to behave as a plasma, it must satisfy three criteria simultaneously:

1. Many particles in Debye sphere

$$N_D = n_e \lambda_D^3 \gg 1 \quad \Leftrightarrow \quad \lambda_D \gg n_e^{-1/3}$$

Ensures collective effects dominate over binary collisions. The Debye sphere must contain enough particles for statistical shielding. Equivalently, the Debye length must exceed the mean interparticle spacing.

2. System size exceeds Debye length

$$L \gg \lambda_D$$

Ensures quasi-neutrality on macroscopic scales. The plasma can shield out internal electric fields and support collective oscillations. Without this criterion, boundary effects dominate and collective behavior cannot develop.

3. Plasma period shorter than observation time

$$\omega_{pe} \tau_{\text{obs}} \gg 1 \quad \Leftrightarrow \quad \tau_{\text{obs}} \gg \omega_{pe}^{-1}$$

Ensures the plasma can respond collectively on relevant timescales. If observations occur faster than τpe, the system behaves like a collection of individual particles rather than a plasma. Also requires collision time τcoll ≫ τpe for kinetic theory to apply.

1.6.1 Combined Plasma Criteria

These three criteria can be combined into a single inequality. Using λD ∝ √(T/n) and ωpe ∝ √n:

$$N_D = \frac{4\pi}{3}n\lambda_D^3 = \frac{4\pi}{3}\frac{(\epsilon_0 k_B T)^{3/2}}{e^3 n^{1/2}} \gg 1$$

This single condition encapsulates the essence of plasma behavior: high temperature and low density favor collective effects.

1.6.2 Plasma Parameter Space

Natural plasmas span an enormous range in the (n, T) plane:

Plasma Typene [m−3]Te [eV]λDNDΓ
Interstellar medium10617 m10910−7
Solar wind (1 AU)1071070 m101110−8
Ionosphere10120.17 mm10610−5
Solar corona10151007 cm101010−8
Gas discharge10162100 μm10510−4
Tokamak edge10191007 μm10710−6
Tokamak core102010470 μm101010−8
ICF target10301031 nm1030.1
White dwarf interior10361060.01 nm1010

Note how ND ranges from 103 to 1011, yet all satisfy ND ≫ 1. The coupling parameter Γ ranges from 10−8 (ideal plasma) to 10 (strongly coupled).

1.6.3 Key Length and Time Scales

Plasma physics is characterized by a hierarchy of scales:

Length Scales

$$r_L = \frac{v_\perp}{\omega_c} \quad \text{(Larmor radius)}$$
$$\lambda_D = \sqrt{\frac{\epsilon_0 k_B T}{n e^2}} \quad \text{(Debye length)}$$
$$c/\omega_{pe} \quad \text{(Collisionless skin depth)}$$
$$\lambda_{mfp} = v_{\text{th}}/\nu \quad \text{(Mean free path)}$$
$$L \quad \text{(System size)}$$

Time Scales

$$\omega_{ce}^{-1} = \frac{m_e}{eB} \quad \text{(Electron cyclotron)}$$
$$\omega_{pe}^{-1} = \sqrt{\frac{\epsilon_0 m_e}{n e^2}} \quad \text{(Plasma period)}$$
$$\omega_{ci}^{-1} = \frac{m_i}{ZeB} \quad \text{(Ion cyclotron)}$$
$$\nu^{-1} \quad \text{(Collision time)}$$
$$\tau_{\text{conf}} \quad \text{(Confinement time)}$$

Typical ordering: ωce−1 ≪ ωpe−1 ≪ ωci−1 ≪ ν−1 ≪ τconf

The Fourth State of Matter

These three criteria define plasma behavior across 20 orders of magnitude in density (106−1030 m−3) and 8 orders of magnitude in temperature (0.1−106 eV)—from neon signs to stellar cores, from the interstellar medium to inertial fusion targets. Plasmas comprise ~99% of the visible universe, making plasma physics truly the study of cosmic matter!

1.6.4 Plasma Definition Revisited

We can now give a more precise definition:

Formal Definition of a Plasma

A plasma is an ionized gas satisfying:

(1) Weak coupling: Γ = (Upot/Ukin) ≪ 1, or equivalently ND ≫ 1

(2) Quasi-neutrality: L ≫ λD, ensuring |ne − Zni|/ne ≪ 1

(3) Collective behavior: ωpeτobs ≫ 1, and typically ωpe ≫ νcoll

These conditions ensure the gas exhibits collective electromagnetic phenomena—waves, instabilities, self-consistent fields—that distinguish plasma from ordinary neutral gas behavior.

← Part I OverviewNext: Single Particle Motion →