Plasma Physics Course

Total: 8 parts covering plasma fundamentals through advanced topics

Part I, Chapter 4

Plasma Parameters

Characterizing plasmas through dimensionless quantities

4.1 Temperature and Density

The plasma state is fundamentally characterized by two primary quantities: temperature T andnumber density n. These parameters span enormous ranges across different plasma environments, from tenuous interstellar gas to ultra-dense stellar interiors, covering more than 30 orders of magnitude in the (n, T) phase space.

4.1.1 Temperature Scales and Units

In plasma physics, temperature is conventionally expressed in electron volts (eV) rather than Kelvin. The conversion factor relates thermal energy to temperature via Boltzmann's constant:

$$k_B T = 1 \text{ eV} \equiv 1.602 \times 10^{-19} \text{ J} \quad \Leftrightarrow \quad T = 11,604.5 \text{ K}$$

This convention arises because the thermal velocity of particles is naturally expressed in terms of eV. For a particle with thermal energy E = kBT, the characteristic thermal velocity is:

$$v_{th} = \sqrt{\frac{2k_B T}{m}} = \sqrt{\frac{2eV}{m}}$$

where V is the temperature in volts. For electrons at T = 1 eV:

$$v_{th,e} = \sqrt{\frac{2 \times 1.602 \times 10^{-19}}{9.109 \times 10^{-31}}} \approx 5.93 \times 10^5 \text{ m/s}$$

4.1.2 Non-Equilibrium Temperature

Unlike neutral gases, plasmas are often far from thermal equilibrium. Electrons and ions can maintain distinct temperatures Te ≠ Ti for extended periods, particularly when:

  • • Energy is preferentially deposited into one species (e.g., RF heating of electrons)
  • • Collision times are long compared to system evolution times
  • • Expansion or compression occurs faster than equilibration

The temperature equilibration rate between electrons and ions is determined by energy exchange in Coulomb collisions. The equilibration timescale is:

$$\tau_{eq} = \frac{3\sqrt{m_i m_e}}{4\sqrt{2\pi} n e^4 \ln\Lambda} \left(\frac{k_B T_e}{m_e}\right)^{3/2}$$

This can be expressed in terms of the electron-ion collision time:

$$\tau_{eq} \sim \left(\frac{m_i}{m_e}\right)^{1/2} \tau_{ei} \approx 43 \tau_{ei} \quad \text{(for hydrogen)}$$

The factor arises because electrons must undergo many collisions to transfer their small mass of momentum to the much heavier ions. For a deuterium plasma at n = 1020 m−3 and Te = 10 keV:

$$\tau_{eq} \approx 0.19 \text{ s}$$

This is comparable to energy confinement times in fusion devices, explaining why Te ≠ Tiis common in tokamak plasmas.

4.1.3 Velocity Distribution Functions

Temperature is fundamentally a statistical concept, defined through the velocity distribution function. For an isotropic Maxwellian plasma, the distribution is:

$$f(\mathbf{v}) = n\left(\frac{m}{2\pi k_B T}\right)^{3/2} \exp\left(-\frac{m v^2}{2k_B T}\right)$$

The temperature is related to the second moment of the distribution:

$$T = \frac{m}{3k_B} \int |\mathbf{v} - \mathbf{u}|^2 f(\mathbf{v}) d^3v$$

where u is the mean flow velocity. Non-Maxwellian distributions require more careful definition of temperature, often distinguished as parallel and perpendicular temperatures in magnetized plasmas:

$$T_\parallel = \frac{m}{k_B} \langle (v_\parallel - u_\parallel)^2 \rangle, \quad T_\perp = \frac{m}{2k_B} \langle (v_\perp - u_\perp)^2 \rangle$$

4.1.4 Density Regimes Across Nature

Plasma densities vary by more than 30 orders of magnitude across different astrophysical and laboratory environments:

Astrophysical Plasmas

  • • Intergalactic medium:
    n ∼ 10−1 – 103 m−3, T ∼ 106 K
  • • Interstellar medium:
    n ∼ 106 m−3, T ∼ 104 K
  • • Solar wind (1 AU):
    n ∼ 5 × 106 m−3, T ∼ 10 eV
  • • Solar corona:
    n ∼ 1014 m−3, T ∼ 100 eV
  • • Solar core:
    n ∼ 1032 m−3, T ∼ 1 keV

Laboratory Plasmas

  • • Ionosphere F-layer:
    n ∼ 1012 m−3, T ∼ 0.1 eV
  • • Glow discharge:
    n ∼ 1016 m−3, Te ∼ 3 eV
  • • Processing plasma:
    n ∼ 1017 m−3, T ∼ 5 eV
  • • Tokamak edge:
    n ∼ 1019 m−3, T ∼ 100 eV
  • • Tokamak core:
    n ∼ 1020 m−3, T ∼ 10 keV
  • • ICF hotspot:
    n ∼ 1031 m−3, T ∼ 10 keV

4.1.5 Quasi-Neutrality Condition

The density ratio of electrons to ions is constrained by the quasi-neutrality condition. For a plasma with multiple ion species:

$$n_e = \sum_i Z_i n_i$$

where Zi is the charge state of ion species i. This condition is satisfied to extraordinary precision (δn/n ∼ 10−6) over scales larger than the Debye length. Deviations create electric fields that rapidly restore neutrality on the plasma oscillation timescale ωpe−1.

4.2 Beta Parameter

The plasma beta (β) is perhaps the most important dimensionless parameter characterizing magnetized plasmas. It quantifies the ratio of thermal (kinetic) pressure to magnetic pressure:

$$\beta = \frac{P_{thermal}}{P_{magnetic}} = \frac{P}{B^2/2\mu_0}$$

where the thermal pressure is given by the ideal gas law summed over all species:

$$P = \sum_s n_s k_B T_s = n_e k_B T_e + \sum_i n_i k_B T_i$$

For a quasineutral plasma with a single ion species and Te = Ti = T:

$$\beta = \frac{2n k_B T}{B^2/2\mu_0} = \frac{4\mu_0 n k_B T}{B^2}$$

4.2.1 Physical Interpretation

Beta determines which force dominates plasma dynamics:

β ≪ 1 (Low-β regime)

Magnetic pressure dominates. Particles are strongly magnetized, following field lines closely. Magnetic field geometry controls plasma behavior. Field lines are nearly rigid.

$$\frac{\nabla P}{P} \ll \frac{|\nabla B|}{B} \quad \Rightarrow \quad \text{Force balance: } \frac{\mathbf{J} \times \mathbf{B}}{c} \approx 0$$

β ∼ 1 (Moderate-β regime)

Thermal and magnetic pressures are comparable. Plasma can significantly distort magnetic field. MHD equilibrium requires careful balance. Most relevant for fusion applications.

$$\nabla P = \mathbf{J} \times \mathbf{B}$$

β ≫ 1 (High-β regime)

Thermal pressure dominates. Magnetic field is passively advected by plasma flow. Field can be expelled from plasma (diamagnetism). Inertial effects dominate dynamics.

$$\rho \frac{D\mathbf{v}}{Dt} \approx -\nabla P$$

4.2.2 Troyon Beta Limit

In toroidal confinement devices, there exists an empirical stability limit on achievable beta. The Troyon limit relates the maximum stable beta to machine parameters:

$$\beta_{max} \leq g \frac{I_p[\text{MA}]}{a[\text{m}] B_0[\text{T}]} \quad (\%)$$

where g ≈ 2.8 is the normalized beta coefficient, Ip is the plasma current,a is the minor radius, and B0 is the toroidal field. This can be rewritten as:

$$\beta_N = \frac{\beta[\%]}{I_p[\text{MA}]/(a[\text{m}] B_0[\text{T}])} \leq g$$

The limit arises from current-driven MHD instabilities (kink modes, ballooning modes) that become unstable when the plasma pressure gradient becomes too steep. Exceeding βN ≈ 3–4 typically triggers disruptions in conventional tokamaks.

4.2.3 Beta Across Plasma Regimes

Astrophysical Plasmas

  • • Solar corona: β ∼ 10−4 – 10−2
    Magnetically dominated, field controls structure
  • • Solar wind (1 AU): β ∼ 0.5 – 2
    Transition from magnetic to kinetic dominance
  • • Magnetosphere: β ∼ 0.01 – 1
    Varies strongly with location
  • • Molecular clouds: β ∼ 0.01 – 0.1
    Sub-Alfvénic, field supports against gravity

Laboratory Plasmas

  • • Tokamak: β ∼ 0.02 – 0.05
    Limited by MHD stability
  • • Spherical tokamak: β ∼ 0.2 – 0.4
    Higher β due to compact geometry
  • • Stellarator: β ∼ 0.03 – 0.05
    External coils provide stability
  • • RFP: β ∼ 0.1 – 0.2
    Relaxation allows higher β
  • • Field-reversed config: β ∼ 1
    Extreme high-β confinement

4.2.4 Beta and Fusion Reactor Performance

For magnetic fusion energy, beta is critical for economic viability. The fusion power density scales as:

$$P_{fusion} \propto n^2 T^2 \langle\sigma v\rangle \propto \beta^2 B^4 T^{-1}$$

Higher beta means higher power density for given magnetic field, reducing the cost and size of magnets. The magnetic energy required for confinement scales as:

$$E_{mag} \sim \frac{B^2}{2\mu_0} V \sim \frac{P V}{\beta}$$

Thus, doubling β halves the required magnetic energy (and cost) for the same thermal energy. This drives the quest for advanced tokamak scenarios with improved β limits.

Numerical Example: ITER

ITER parameters: n = 1020 m−3, T = 10 keV, B0 = 5.3 T

$$\beta = \frac{4\mu_0 \times 10^{20} \times 1.602 \times 10^{-15}}{(5.3)^2} \approx 0.028 = 2.8\%$$

This is close to the Troyon limit for ITER's configuration, representing highly optimized performance.

4.3 Coupling Parameter

The coupling parameter Γ (also called the plasma coupling constant) quantifies the relative importance of Coulomb potential energy to thermal kinetic energy. It determines whether a plasma behaves as an ideal gas or exhibits strong correlations and collective effects.

4.3.1 Definition and Physical Meaning

The coupling parameter for species s is defined as:

$$\Gamma_s = \frac{(Z_s e)^2}{4\pi\epsilon_0 a_s k_B T_s} = \frac{U_{Coulomb}}{U_{thermal}}$$

where as is the Wigner-Seitz radius, defined as the radius of a sphere containing one particle on average:

$$\frac{4\pi a_s^3}{3} = \frac{1}{n_s} \quad \Rightarrow \quad a_s = \left(\frac{3}{4\pi n_s}\right)^{1/3}$$

The numerator represents the typical Coulomb interaction energy between neighboring particles, while the denominator is the thermal energy. Thus:

  • • Γ ≪ 1: Kinetic energy dominates, particles move freely (ideal plasma)
  • • Γ ∼ 1: Potential and kinetic energies comparable (liquid-like)
  • • Γ ≫ 1: Potential energy dominates, particles localized (crystal-like)

4.3.2 Relation to Plasma Parameter

The coupling parameter is inversely related to the plasma parameter Λ = nλD3. Starting from the Debye length:

$$\lambda_D = \sqrt{\frac{\epsilon_0 k_B T}{n e^2}}$$

we can express the plasma parameter as:

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

Comparing with Γ, we find:

$$\Gamma \sim \left(\frac{3}{4\pi}\right)^{1/3} \Lambda^{-2/3} \approx 0.543 \, \Lambda^{-2/3}$$

Therefore, the classical plasma criterion Λ ≫ 1 is equivalent to weak coupling Γ ≪ 1. Most plasmas studied in plasma physics satisfy both conditions.

4.3.3 Regimes of Plasma Coupling

Weakly Coupled Plasma (Γ ≪ 1)

Classical kinetic theory applies. Binary collisions dominate. Debye shielding is effective. The Boltzmann equation or Vlasov equation with collision operator describes the dynamics.

$$\text{Examples: } \begin{cases} \text{Fusion plasmas: } \Gamma \sim 10^{-3} \\ \text{Gaseous discharges: } \Gamma \sim 10^{-2} \\ \text{Solar wind: } \Gamma \sim 10^{-4} \end{cases}$$

Intermediate Coupling (Γ ∼ 1)

Warm Dense Matter regime. Correlations are important but particles still mobile. Neither kinetic theory nor condensed matter theory fully applies. Requires advanced computational methods (quantum molecular dynamics, density functional theory).

$$\text{Examples: } \begin{cases} \text{Shock-compressed matter} \\ \text{Planetary interiors} \\ \text{ICF ablation fronts} \end{cases}$$

Strongly Coupled Plasma (Γ ≫ 1)

Particles form ordered structures. Exhibits liquid or crystalline properties. Radial distribution function shows long-range order. Phase transitions occur at Γ ≈ 175 (BCC lattice).

$$\text{Examples: } \begin{cases} \text{White dwarf cores: } \Gamma \sim 100-1000 \\ \text{Dusty plasma crystals: } \Gamma \sim 10^3 \\ \text{Ultra-cold plasmas: } \Gamma \sim 10-100 \end{cases}$$

4.3.4 One-Component Plasma (OCP) Model

For strongly coupled systems, the One-Component Plasma (OCP) model is commonly used. It consists of identical point charges in a uniform neutralizing background. The OCP equation of state can be written as:

$$\frac{P}{nk_B T} = 1 + \frac{1}{3}\Gamma u_{ex}(\Gamma)$$

where uex(Γ) is the excess internal energy per particle due to interactions. Molecular dynamics simulations show:

$$u_{ex}(\Gamma) \approx \begin{cases} -0.899\Gamma + 0.095\Gamma^{1/4} - 0.18 & \text{(fluid phase, } \Gamma < 175\text{)} \\ -0.896\Gamma + 1.5 & \text{(crystal phase, } \Gamma > 175\text{)} \end{cases}$$

4.3.5 Experimental Realization: Dusty Plasmas

Dusty or complex plasmas provide laboratory access to the strongly coupled regime. Micron-sized dust particles acquire large negative charges (Zd ∼ 103–105e) in low-pressure discharges:

$$\Gamma_d = \frac{(Z_d e)^2}{4\pi\epsilon_0 a_d k_B T_d}$$

With large Zd and low Td (room temperature), Γd can exceed 1000, leading to visible Coulomb crystals that can be directly imaged. These systems allow experimental study of phase transitions, waves, and transport in strongly coupled plasmas.

Numerical Example: Fusion vs. White Dwarf

ITER plasma: n = 1020 m−3, T = 10 keV

$$a = \left(\frac{3}{4\pi \times 10^{20}}\right)^{1/3} \approx 1.34 \times 10^{-7} \text{ m}$$
$$\Gamma = \frac{1.44 \times 10^{-9}}{1.34 \times 10^{-7} \times 1.602 \times 10^{-15}} \approx 6.7 \times 10^{-3} \ll 1$$

White dwarf core: n = 1036 m−3, T = 107 K ≈ 860 eV

$$a = \left(\frac{3}{4\pi \times 10^{36}}\right)^{1/3} \approx 4.3 \times 10^{-13} \text{ m}$$
$$\Gamma = \frac{1.44 \times 10^{-9}}{4.3 \times 10^{-13} \times 1.38 \times 10^{-16}} \approx 243 \gg 1$$

The white dwarf interior forms a crystalline lattice despite the high temperature!

4.4 Magnetization Parameter

In the presence of a magnetic field, charged particles execute helical gyration around field lines. The degree of magnetization determines whether the magnetic field strongly influences particle dynamics or is merely a perturbation. Several related dimensionless parameters characterize magnetization.

4.4.1 Hall Parameter

The most fundamental magnetization parameter is the Hall parameter ωcτ, which compares the cyclotron period to the collision time:

$$\chi = \omega_c \tau = \frac{|q|B}{m} \cdot \frac{1}{\nu}$$

where ωc = |q|B/m is the cyclotron frequency and τ = ν−1is the collision time. The Hall parameter gives the number of gyro-orbits completed between collisions:

χ ≪ 1: Unmagnetized Regime

Collisions interrupt gyro-motion before one orbit completes. Particles diffuse freely across field lines. Classical transport (not modified by magnetic field). Magnetic field is a weak perturbation.

$$D_\perp \approx D_\parallel \approx v_{th}^2 \tau$$

χ ∼ 1: Weakly Magnetized

Gyro-motion partially develops before collision. Significant reduction in cross-field transport. Transition regime between classical and magnetized behavior.

χ ≫ 1: Magnetized Regime

Particles complete many gyro-orbits between collisions. Strongly confined to field lines. Perpendicular transport severely restricted. Gyrokinetic theory applicable.

$$D_\perp \approx \frac{D_\parallel}{\chi^2} = \frac{v_{th}^2 \tau}{\omega_c^2 \tau^2} = \frac{v_{th}^2}{\omega_c^2 \tau}$$

4.4.2 Larmor Radius Parameter

The Larmor (gyro) radius is the radius of the circular motion of a particle with thermal velocity perpendicular to the magnetic field:

$$r_L = \frac{mv_\perp}{|q|B} = \frac{v_{th}}{\omega_c}$$

For a thermal distribution, using vth = √(2kBT/m):

$$r_L = \frac{\sqrt{2mk_B T}}{|q|B} = \frac{1}{|q|B}\sqrt{2mk_B T}$$

The ratio of Larmor radius to system size L is another key magnetization parameter:

$$\delta = \frac{r_L}{L}$$

When δ ≪ 1, particles are tightly bound to field lines and fluid (MHD) description is valid. When δ ≳ 1, kinetic effects dominate and fluid theory breaks down. In tokamaks, typical values are:

$$\delta_e \sim 10^{-4}, \quad \delta_i \sim 10^{-2}$$

4.4.3 Species-Dependent Magnetization

Electrons and ions have vastly different magnetization due to their mass ratio. For the same temperature:

$$\frac{r_{L,i}}{r_{L,e}} = \sqrt{\frac{m_i}{m_e}} \approx 43 \quad \text{(for hydrogen)}$$
$$\frac{\omega_{ci}}{\omega_{ce}} = \frac{m_e}{m_i} \approx \frac{1}{1836}$$

This leads to the common situation where electrons are highly magnetized (χe ≫ 1) while ions are less so (χi may be ∼ 1 or even ≪ 1). Different transport regimes can coexist for different species.

4.4.4 Alfvén Velocity and Beta

Magnetization is also characterized by comparing thermal velocity to the Alfvén velocity:

$$v_A = \frac{B}{\sqrt{\mu_0 \rho}} = \frac{B}{\sqrt{\mu_0 n_i m_i}}$$

The ratio of thermal to Alfvén speed is directly related to beta:

$$\frac{v_{th}^2}{v_A^2} = \frac{2k_B T/m}{B^2/\mu_0\rho} = \frac{2\mu_0 nk_B T}{B^2} = \beta$$

When β ≪ 1, the Alfvén velocity exceeds thermal velocities and magnetic forces dominate plasma dynamics. MHD waves (Alfvén, fast, slow) propagate faster than particle motions.

4.4.5 Quantitative Examples

Tokamak (ITER)

Parameters: B = 5.3 T, n = 1020 m−3, Te = Ti = 10 keV

$$\omega_{ce} = \frac{eB}{m_e} = 9.3 \times 10^{11} \text{ rad/s}$$
$$\nu_{ei} \approx 3 \times 10^{6} \text{ s}^{-1}$$
$$\chi_e = \frac{\omega_{ce}}{\nu_{ei}} \approx 3.1 \times 10^5 \gg 1$$
$$r_{L,e} = 1.2 \text{ mm}, \quad r_{L,i} = 5.1 \text{ cm}$$

Electrons extremely magnetized; ions moderately magnetized

Solar Corona

Parameters: B = 100 G = 0.01 T, n = 1014 m−3, T = 100 eV

$$\omega_{ce} = 1.76 \times 10^{9} \text{ rad/s}$$
$$\nu_{ei} \approx 1 \text{ s}^{-1}$$
$$\chi_e \approx 1.76 \times 10^{9} \gg 1$$
$$r_{L,e} = 3.2 \text{ cm}, \quad r_{L,i} = 1.4 \text{ m}$$

Collisionless, field-aligned plasma; β ∼ 10−3

The extreme magnetization of electrons (χe ∼ 105–109) in most laboratory and astrophysical plasmas is what makes magnetic confinement fusion possible. Without this strong magnetization, cross-field transport would rapidly disperse the plasma.

4.5 Classification of Plasmas

Plasmas exhibit extraordinary diversity, spanning scales from sub-nanometer (quantum plasmas) to megaparsec (intergalactic medium). They can be classified along multiple independent axes, each highlighting different physical regimes and behaviors.

4.5.1 Classification by Temperature

Cold Plasmas (T < 1 eV)

Partial ionization, typically non-equilibrium with Te ≫ Ti ≈ Tgas. Elastic collisions with neutrals dominate energy balance. Electric field required to sustain ionization.

  • • Fluorescent lamps: Te ∼ 1 eV, n ∼ 1017 m−3
  • • Ionosphere D/E layers: T ∼ 0.03–0.1 eV
  • • Semiconductor processing: Te ∼ 2–5 eV
  • • Glow discharge: Te ∼ 2–4 eV, Tgas ∼ 0.03 eV

Thermal Plasmas (1 eV < T < 100 eV)

High ionization fraction, often near local thermodynamic equilibrium (LTE). Te ≈ Ti ≈ Tgas. Self-sustaining through thermal ionization.

  • • Arc discharge: T ∼ 1–2 eV (10,000–20,000 K)
  • • Plasma torches: T ∼ 2–5 eV
  • • Lightning channel: T ∼ 2–3 eV
  • • Solar chromosphere: T ∼ 1 eV

Hot Plasmas (T > 1 keV)

Fully ionized, collisionless or weakly collisional. Kinetic effects important. Often magnetized. Relevant for fusion and astrophysics.

  • • Tokamak core: T ∼ 10–40 keV
  • • Solar corona: T ∼ 100–200 eV
  • • Solar wind: T ∼ 10 eV
  • • Supernova remnant shocks: T ∼ 1–10 keV
  • • Galaxy cluster ICM: T ∼ 2–10 keV

Ultra-Relativistic (T > 511 keV)

Electron rest mass energy exceeded: kBT > mec2. Relativistic corrections essential. Pair production (e+e−) can occur.

  • • Pulsar magnetospheres: T ∼ MeV
  • • Active galactic nucleus jets: T ∼ GeV
  • • Early universe (t < 1 s): T ∼ MeV–GeV
  • • Quark-gluon plasma: T ∼ 200 MeV (RHIC/LHC)

4.5.2 Classification by Density

Density controls collective behavior through the plasma parameter Λ = nλD3:

Ultra-Low Density (n < 1010 m−3)

Collisionless, mean free path ≫ system size. Single-particle dynamics dominate.

  • • Intergalactic medium: n ∼ 10−1–103 m−3
  • • Pulsar magnetospheres: n ∼ 106 m−3

Low Density (1010 < n < 1018 m−3)

Space plasmas, weakly collisional laboratory plasmas. Kinetic theory essential.

  • • Solar wind: n ∼ 106–107 m−3
  • • Ionosphere: n ∼ 1010–1012 m−3
  • • Low-pressure discharges: n ∼ 1015–1017 m−3

Medium Density (1018 < n < 1024 m−3)

Magnetic fusion, atmospheric pressure plasmas. MHD or gyrokinetic descriptions.

  • • Processing plasmas: n ∼ 1017–1018 m−3
  • • Tokamak edge: n ∼ 1018–1019 m−3
  • • Tokamak core: n ∼ 1020 m−3

High Density (n > 1026 m−3)

ICF, stellar interiors, warm dense matter. Strong coupling, quantum effects emerge.

  • • Solid-density plasmas: n ∼ 1028–1029 m−3
  • • ICF implosions: n ∼ 1030–1032 m−3
  • • White dwarf cores: n ∼ 1036 m−3
  • • Neutron star crusts: n ∼ 1044 m−3 (nuclear density)

4.5.3 Classification by Generation Method

DC Discharge

Electric field between electrodes accelerates electrons, causing impact ionization.

  • • Glow discharge (Crookes tube, neon signs)
  • • Arc discharge (welding, lighting)
  • • Corona discharge (high voltage lines)

RF/Microwave Discharge

Time-varying electromagnetic fields heat electrons, avoiding electrode contamination.

  • • Inductively coupled plasma (ICP)
  • • Capacitively coupled plasma (CCP)
  • • Microwave (ECR) plasma
  • • Helicon plasma sources

Magnetic Confinement

Strong magnetic fields confine hot plasma, heated by various mechanisms.

  • • Tokamak (toroidal + poloidal field)
  • • Stellarator (3D twisted field)
  • • Reversed-field pinch (self-organized)
  • • Magnetic mirror (open field lines)
  • • Field-reversed configuration

Inertial Confinement

Rapid compression achieves extreme densities; inertia provides brief confinement.

  • • Laser-driven ICF (NIF, LMJ)
  • • Z-pinch (magnetic compression)
  • • Ion-beam driven
  • • Magnetized liner inertial fusion (MagLIF)

Laser/Beam Produced

High-intensity radiation or particle beams ionize and heat target material.

  • • Laser-produced plasma (ablation)
  • • Electron-beam melting
  • • Ion-beam surface modification
  • • Laser wakefield accelerators

Natural/Astrophysical

Spontaneous formation through UV radiation, shocks, gravitational compression, etc.

  • • Solar/stellar coronae (heating mystery)
  • • Planetary ionospheres (photoionization)
  • • Supernova remnants (shock heating)
  • • Accretion disk plasmas

4.6 The Density-Temperature Phase Space

The complete landscape of plasma physics can be mapped on a logarithmic (n, T) diagram, spanning more than 40 orders of magnitude in density and 10 orders of magnitude in temperature. This phase space is bounded by fundamental physical limits and divided into distinct regimes by important dimensionless parameters.

4.6.1 Classical Plasma Boundary

The fundamental requirement for collective plasma behavior is that the Debye sphere contains many particles:

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

Rearranging this condition gives the classical plasma boundary:

$$n \ll \left(\frac{4\pi\epsilon_0 k_B T}{e^2}\right)^{3/2} \cdot \frac{1}{\Lambda}$$

For Λ = 10 (minimum for plasma behavior), this gives:

$$n[\text{m}^{-3}] \lesssim 1.5 \times 10^{19} \cdot T[\text{eV}]^{3/2}$$

Below this line, the coupling parameter Γ becomes large and the system behaves more like a liquid or solid than a plasma. This boundary separates classical plasmas from strongly coupled Coulomb systems.

4.6.2 Quantum Degeneracy Boundary

At high densities and low temperatures, quantum mechanical effects become important when thethermal de Broglie wavelength becomes comparable to the inter-particle spacing:

$$\Lambda_{th} = \frac{h}{\sqrt{2\pi m k_B T}}$$

Quantum degeneracy occurs when Λth ≳ n−1/3, or equivalently when thedegeneracy parameter exceeds unity:

$$n_Q = \frac{1}{\Lambda_{th}^3} = \frac{(2\pi m k_B T)^{3/2}}{h^3} \lesssim n$$

For electrons, this gives the quantum boundary:

$$n[\text{m}^{-3}] \gtrsim 4.8 \times 10^{27} \cdot T[\text{eV}]^{3/2}$$

Above this line, electrons obey Fermi-Dirac statistics rather than classical Maxwell-Boltzmann statistics. The electron gas becomes degenerate, with Fermi pressure dominating. This regime is crucial in white dwarfs, neutron stars, and warm dense matter physics.

4.6.3 Relativistic Boundary

When the thermal energy approaches the electron rest mass energy, relativistic corrections become essential:

$$k_B T \sim m_e c^2 = 511 \text{ keV} \quad \text{or} \quad T \sim 5.9 \times 10^9 \text{ K}$$

Above this temperature, the momentum distribution becomes relativistic:

$$E = \sqrt{(pc)^2 + (m_e c^2)^2} \approx pc \quad \text{for } p \gg m_e c$$

At even higher temperatures (T ≳ 1 MeV), electron-positron pair production becomes significant:

$$\gamma + \gamma \leftrightarrow e^+ + e^-$$

This regime is relevant for pulsar magnetospheres, gamma-ray bursts, and the early universe.

4.6.4 Ionization Equilibrium

The ionization state of the plasma depends on the balance between ionization and recombination. The Saha equation gives the equilibrium ionization fraction:

$$\frac{n_{i+1} n_e}{n_i} = \frac{2g_{i+1}}{g_i} \left(\frac{2\pi m_e k_B T}{h^2}\right)^{3/2} \exp\left(-\frac{\chi_i}{k_B T}\right)$$

where χi is the ionization energy from state i to i+1, and giare statistical weights. For hydrogen ionization (χ = 13.6 eV):

$$\frac{n_e^2}{n_H} = 2.4 \times 10^{21} \cdot T[\text{K}]^{3/2} \exp\left(-\frac{1.58 \times 10^5}{T[\text{K}]}\right) \quad [\text{m}^{-3}]$$

This determines the ionization fraction as a function of (n, T). Full ionization requires T ≳ χi/kB or extremely low density where recombination is negligible.

4.6.5 Fusion Ignition Contours

For magnetic confinement fusion, the Lawson criterion specifies the minimumnTτE (triple product) required for ignition. For D-T fusion:

$$n T \tau_E \geq 3 \times 10^{21} \text{ m}^{-3} \cdot \text{keV} \cdot \text{s}$$

The fusion power density depends on the fusion reactivity ⟨σv⟩:

$$P_{fusion} = \frac{1}{4} n^2 \langle\sigma v\rangle E_{fusion}$$

For D-T, ⟨σv⟩ peaks at T ≈ 70 keV, but practical operation occurs at T = 10–20 keV where the product of reactivity and energy confinement is optimized.

Inertial confinement fusion operates at much higher density (n ∼ 1031 m−3) but extremely short confinement time (τE ∼ ps–ns), satisfying the same nTτ criterion.

4.6.6 Additional Phase Space Boundaries

Transparency vs. Opacity

When the mean free path for photons λmfp ∼ (nσ)−1 becomes smaller than system size, the plasma becomes optically thick. This boundary depends on photon energy and plasma composition.

$$n \gtrsim \frac{1}{\sigma L} \quad \text{(optical depth } \tau \sim 1\text{)}$$

Collisional vs. Collisionless

The transition from collisional to collisionless behavior occurs when the collision mean free path exceeds the system size L. For Coulomb collisions:

$$\lambda_{mfp} = \frac{v_{th}}{\nu_{ei}} = \frac{12\pi^{3/2} \epsilon_0^2 (k_B T)^2}{n e^4 \ln\Lambda} \gtrsim L$$

Most space plasmas are collisionless; fusion plasmas lie in between.

Equipartition Line

The line Te = Ti divides two-temperature plasmas from equilibrated plasmas. Equilibration time scales as (mi/me)1/2τei, so non-equilibrium is common in rapidly evolving systems.

4.6.7 Plasma Regimes in (n, T) Space

The complete (n, T) diagram reveals several distinct plasma regimes:

Classical Weakly Coupled

Region: Low-moderate n, moderate-high T
Parameters: Λ ≫ 1, Γ ≪ 1
Examples: Fusion plasmas, solar corona

Strongly Coupled

Region: High n, low-moderate T
Parameters: Γ ≳ 1
Examples: Warm dense matter, white dwarf cores

Quantum Degenerate

Region: Very high n, low-moderate T
Parameters: n ≫ nQ
Examples: White dwarfs, neutron star crusts

Relativistic

Region: Any n, very high T (MeV scale)
Parameters: kBT ≳ mec2
Examples: GRBs, AGN jets, early universe

Summary: Key Boundaries in (n, T) Space

Classical plasma:$n \lesssim 1.5 \times 10^{19} T^{3/2}$ (Λ ≳ 10)
Quantum degeneracy:$n \gtrsim 4.8 \times 10^{27} T^{3/2}$ (electrons)
Relativistic regime:T ≳ 511 keV (electrons)
D-T fusion ignition:$nT\tau_E \gtrsim 3 \times 10^{21}$ m$^{-3}$·keV·s
Pair production:T ≳ 1 MeV

Understanding the (n, T) phase space is essential for recognizing which physical processes dominate in different plasma environments and selecting appropriate theoretical descriptions (classical vs. quantum, kinetic vs. fluid, non-relativistic vs. relativistic).

← CollisionsNext: Magnetized Plasmas →