Spotlight #29 – Plasma Physics: Debye Shielding, Aurora, Fusion Confinement and Magnetohydrodynamics

A plasma is an ionised gas: a collection of free electrons and positive ions in which long-range electromagnetic forces play the dominant role, rather than short-range collisions. Plasma exists wherever thermal or electromagnetic energy is sufficient to ionise atoms — in stellar interiors at 10 million degrees, in fluorescent lamps at 20 000 K, in lightning channels at 30 000 K, and in fusion tokamaks at 150 million degrees.

The key property that distinguishes a plasma from a random mix of ions and electrons is its collective behaviour: charged particles respond to each other’s fields over distances much larger than their separation, creating waves, instabilities, and a capacity to confine currents and magnetic fields. Understanding plasma behaviour requires both single-particle physics (cyclotron motion, drifts) and fluid theory (magnetohydrodynamics, MHD).

1. Debye Shielding — The Plasma’s Self-Screening Length

Place a positive test charge in a plasma. Nearby electrons are attracted and ions repelled; the resulting charge cloud partially screens the bare Coulomb field. The characteristic distance over which the screening is effective is the Debye length λ_D. Beyond this distance the plasma appears electrically neutral; within it, individual particles experience significant electric fields. The Debye length is the most fundamental length scale of any plasma.

Electric potential around a test charge q in a plasma:
  φ(r) = (q / 4πε_0 r) · exp(−r/λ_D)   (Yukawa / screened Coulomb potential)
  λ_D = Debye length

Debye length (electron screening dominates):
  λ_D = √(ε_0 k_B T_e / n_e e²)
  T_e = electron temperature, n_e = electron number density

Numerical values:
  Solar corona:  T_e ~ 10^6 K,  n_e ~ 10^12 m^−3  → λ_D ~ 0.07 m
  Ionosphere:    T_e ~ 1500 K,  n_e ~ 10^11 m^−3  → λ_D ~ 0.003 m
  Tokamak:       T_e ~ 10^8 K,  n_e ~ 10^20 m^−3  → λ_D ~ 7 × 10^−5 m
  Fluorescent lamp: T_e ~ 2 eV, n_e ~ 10^18 m^−3  → λ_D ~ 10^−5 m

Plasma criteria (conditions for collective behaviour):
  1. Quasi-neutrality: system size L ≫ λ_D   (bulk is neutral; only sheaths deviate)
  2. Many particles per Debye sphere: N_D = (4/3)πλ_D³ n ≫ 1
     N_D ≫ 1 ensures the continuous distribution assumption is valid.
  3. Plasma frequency ω_p ≫ collision frequency ν_collision
     ω_p = √(n_e e² / ε_0 m_e)   (plasma oscillation frequency)
     At f < f_p = ω_p/2π, electromagnetic waves cannot propagate (reflection).

Plasma sheath:
  At a boundary the quasi-neutrality breaks down over a distance ~ few λ_D.
  The sheath is a thin positive-space-charge layer that ensures no net current flows.
  Sheath physics governs plasma-material interactions in chip fabrication and fusion tokamaks alike.

2. Charged Particle Motion — Spirals, Drifts and Mirrors

A charged particle in a uniform magnetic field moves in a circle: the Lorentz force qv×B acts perpendicular to the velocity, providing centripetal acceleration. Combined with any velocity component along the field, the motion is a helix. Non-uniform fields generate systematic drifts perpendicular to both the field and its gradient. These single-particle drifts underlie the large-scale plasma flows in Earth’s magnetosphere and fusion devices.

Cyclotron (gyro) motion:
  Equation: m dv/dt = qv × B
  Circular motion perpendicular to B with:
    Cyclotron frequency: ω_c = |q|B/m   (angular frequency)
    Larmor radius:       r_L = mv_⊥ / |q|B = v_⊥ / ω_c
  Combined with v_∥ along B: helical trajectory

Larmor radii (in 1 T field):
  Electron  (m = 9.1×10^−31 kg, T_e = 100 eV → v_⊥ = 1.87×10^8 m/s):  r_L ~ 1 mm
  Proton    (m = 1.67×10^−27 kg, T_i = 100 eV → v_⊥ = 1.38×10^5 m/s):  r_L ~ 1.4 mm
  Alpha ion (m = 6.64×10^−27 kg, T = 3.5 MeV → v_⊥ = 1.3×10^7 m/s):   r_L ~ 20 cm

E×B drift (in combined electric and magnetic fields):
  v_{E×B} = (E×B) / B²   (independent of particle charge and mass!)
  This is responsible for the overall bulk plasma flow; both ions and electrons drift together.

Gradient drift (in non-uniform B):
  v_{∇B} = (1/2) v_⊥ r_L (∇B × B) / B²   (charge-dependent)
  Direction: opposite for ions and electrons → creates a net current.

Curvature drift:
  In a curved magnetic field (radius of curvature R_c):
  v_c = (mv_∥² r_L) / (qB R_c²) (&hat;R_c × B)
  Curvature + gradient drifts drive the ring current in Earth’s Van Allen belts.

Magnetic mirror effect:
  A charged particle moving from weak to strong B is reflected if pitch angle α < α_mirror:
  sin²α_mirror = B_min / B_max   (mirror ratio)
  The Van Allen radiation belts are magnetic bottles: particles bounce between mirror points.

3. The Aurora Borealis — Charged Particles Painting the Sky

The aurora is a direct visual consequence of single-particle plasma physics. Solar wind electrons and protons are guided along Earth’s dipole magnetic field lines, spiralling inward to the polar regions. Where the field lines plunge into the upper atmosphere, the energetic particles collide with oxygen and nitrogen molecules, exciting them to higher energy states. The emission of light as these states decay creates the aurora’s signature colours.

Solar wind parameters (typical):
  v_sw ~ 400–700 km/s (speed)
  n_sw ~ 5–15 cm^−3 (proton density at 1 AU)
  B_IMF ~ 5–10 nT (interplanetary magnetic field)
  Proton energy: ½m_p v² ~ 1–3 keV

Chapman–Ferraro current layer:
  The solar wind is deflected by Earth’s dipole field at the magnetopause (~10 R_E sunward).
  Magnetic reconnection at the dayside magnetopause injects particles into the magnetosphere.
  Reconnection also occurs in the magnetotail, accelerating particles toward Earth.

Particle precipitation altitudes and colours:
  O (atomic oxygen) green line 557.7 nm:   90–150 km altitude  (most common)
  O (atomic oxygen) red line  630.0 nm:   200–300 km altitude  (softer electrons, <1 keV)
  N&sub2;+ blue/violet 391.4 nm:             60–100 km altitude  (hard electrons, >50 keV)
  Emission altitude correlates with particle energy: harder particles penetrate deeper.

Kp index and auroral oval:
  Kp 0–2: aurora visible only above 65° geomagnetic latitude
  Kp 5–6: visible to ~55° (Scotland, Canada, Germany)
  Kp 9 (storm): visible to <40° (southern US, central Europe, Japan)

Substorm cycle:
  1. Growth phase (~1 h): magnetic flux loaded onto magnetotail
  2. Onset: explosive reconnection in the near-Earth tail
  3. Expansion: bright surge propagates westward along auroral oval
  4. Recovery: gradual dissipation; quiet arcs remain
  Energy input during major storm: ~ 10^17 J (comparable to thousands of nuclear weapons)

4. Tokamak Fusion — Confining a Star in a Bottle

Nuclear fusion — the power source of stars — requires plasma temperatures of 100–200 million degrees Celsius. At these temperatures, deuterium and tritium nuclei overcome their electrostatic repulsion and fuse, releasing helium plus an energetic neutron. The challenge is holding such extreme plasma long and hot enough for the energy output to exceed the energy input. The tokamak, a toroidal magnetic confinement device, is currently the leading approach.

D–T fusion reaction:
  ₁²H + ₁³H → ₂&sup4;He (3.52 MeV) + n (14.06 MeV)
  Total: 17.59 MeV per reaction (vs ~200 MeV for fission of U-235)
  Cross-section peak: ~5 barn at E_cm ~ 100 keV (T ~ 10^8 K)

Lawson criterion (ignition / break-even condition):
  n τ_E T ≥ n τ_E T|_ignition ~ 3 × 10^21 m^−3 s keV  (for D–T)
  n = plasma density (m^−3)
  τ_E = energy confinement time (s)
  T = temperature (keV, where 1 keV = 11.6 million K)
  Fusion power: P_fus ~ n² ⟨σv⟩ E_fus / 4  (for 50–50 D–T mix)

JET world record (2022):
  59 MJ of fusion energy in 5 s → 11.8 MW sustained
  n τ_E T ~ 1.5 × 10^20 m^−3 s keV  (factor ~20 below ignition)

ITER (under construction):
  Major radius R = 6.2 m, minor radius a = 2 m, plasma volume ~840 m³
  Target: Q = P_fus / P_heat ≥ 10  (500 MW out, 50 MW in)
  Expected n τ_E T ~ 6 × 10^21 (close to ignition)

MHD equilibrium (tokamak):
  Grad–Shafranov equation (axisymmetric equilibrium):
    Δ*ψ = −μ_0 R² dp/dψ − f df/dψ
    ψ = poloidal magnetic flux function
    p = pressure, f = RB_φ (toroidal field function)
  Solution gives nested flux surfaces; plasma is confined on tori.
  Safety factor: q = r B_φ / (R B_θ)  (ratio of toroidal to poloidal field line turns)
  Stability requires q > 1 to avoid kink instability; typical q ~ 3–5 at plasma edge.

Key instabilities:
  Kink instability (q < 1): field lines kink, leading to disruption
  Ballooning mode: plasma bulges outward at high β (plasma / magnetic pressure ratio)
  ELMs (Edge-Localised Modes): periodic bursts at plasma edge; must be controlled to protect the wall.

5. Magnetohydrodynamics — Plasma as a Conducting Fluid

When the length and time scales of interest are much larger than the Debye length and cyclotron period, a plasma can be treated as a single conducting fluid — magnetohydrodynamics (MHD). MHD describes solar wind propagation, geomagnetic field generation, astrophysical jet formation and the large-scale stability of fusion plasmas.

Ideal MHD equations (single-fluid, perfect conductor):
  Continuity:     ∂ρ/∂t + ∇·(ρv) = 0
  Momentum:       ρ Dv/Dt = −∇P + J×B
  Energy:         D/Dt (P/ρ²) = 0   (adiabatic)
  Induction:      ∂B/∂t = ∇×(v×B)   (frozen-in flux theorem: B moves with fluid)
  Ohm’s law:     E + v×B = ηJ      (η = resistivity, η = 0 for ideal MHD)
  Closure:        ∇×B = μ_0 J,   ∇·B = 0

Alfvén wave (transverse wave along B):
  Phase speed: v_A = B / √(μ_0 ρ)   (Alfvén speed)
  In the solar corona: B ~ 100 G, ρ ~ 10^−12 kg/m³ → v_A ~ 3 × 10^6 m/s ~ 0.01 c
  Alfvén waves transport energy along field lines and are thought to heat the corona.

Frozen-in flux theorem (ideal MHD):
  Magnetic field lines are “frozen” to the plasma: as plasma moves, field lines move with it.
  Consequence: magnetic topology is preserved in ideal MHD.

Magnetic reconnection (resistive MHD):
  Where field lines of opposite polarity meet, resistivity η ≠ 0 allows field lines to
  break and reconnect in a new topology, releasing stored magnetic energy explosively.
  Sweet–Parker rate: v_rec ~ v_A / √(S)   (S = Lundquist number = μ_0 v_A L / η ~ 10^12 in solar corona)
  Observed reconnection rate ~ 10^−2 v_A  (much faster: “fast reconnection” problem)
  Drives solar flares (up to 10^25 J), coronal mass ejections (CMEs) and geomagnetic storms.

Magnetic pressure and tension:
  Magnetic pressure: P_B = B² / (2μ_0)  (acts like gas pressure, isotropic)
  Magnetic tension:  T_B = B² / μ_0  (acts like a tensioned string along field lines)
  Plasma beta: β = P_plasma / P_B = n k_B T / (B²/2μ_0)
  Tokamak core: β ~ 0.03–0.1;  Solar corona: β ~ 0.01–0.1;  Solar wind: β ~ 1

Magnetic reconnection is one of the central unsolved problems in plasma physics. In the resistive MHD framework, the predicted reconnection rate is far too slow to explain observed solar flares, which can release their energy in minutes. Modern theories — including Hall MHD, electron-scale physics and turbulence-driven reconnection — can match the observed fast rates but remain subject to active experimental study. The Parker Solar Probe mission, launched in 2018 and now flying within 9 solar radii, is collecting in-situ measurements of reconnection regions that will test these theories directly.

Plasma: The Universe’s Dominant State

More than 99% of all visible matter in the universe is in the plasma state. Every star, every active galactic nucleus, every ionised nebula and every solar wind particle is plasma. Yet on Earth, plasma is a laboratory curiosity, an industrial tool (chip etching, welding, neon signs) or a transient phenomenon (lightning, flames). The simulations in this wave trace the bridge from Earth-scale phenomena — fluorescent tubes and auroral curtains — to the astrophysical plasma processes that govern stellar evolution, magnetic field amplification and the energy budget of galaxy clusters. All from the same equations: Maxwell’s laws and Newton’s second law for charged particles.