Plasma and fusion systems
Connect charged-particle kinetics, collective fields, confinement, diagnostics, and reactor-scale energy balance.
Move between single-particle, kinetic, and fluid plasma descriptions; derive characteristic waves and scales; compare magnetic, inertial, and lattice-assisted fusion claims; and calculate diagnostics and system-level gain without merging reaction evidence with power-plant evidence.
Before you begin
- • Electromagnetic fields and potentials
- • Statistical mechanics
- • Experimental methods and error analysis
By the end, you can
- • Calculate plasma frequency, Debye length, gyroradius, and characteristic transport scales.
- • Relate kinetic and fluid models to their validity regimes.
- • Analyze confinement, stability, and fusion reaction rates.
- • Separate nuclear products, plasma gain, engineering gain, and net-electric gain.
Interactive model
Explore before calculating

Live laboratory
Plasma collective-regime map
Compare shielding, collective oscillation, gyro motion, and apparatus scale. A plasma model needs scale ordering—not merely the presence of ionized particles.
Debye length: 2.35e-5 m
Debye population: 5.44e+4
Plasma frequency: 5.64e+10 rad/s
Cyclotron frequency: 1.76e+10 rad/s
Electron gyroradius: 1.07e-4 m
L/λD: 4.25e+4
The system is large compared with its Debye length. Many particles occupy a Debye sphere.
These are necessary teaching checks, not a sufficient classifier; collisions, distributions, boundaries, ion response, and observation timescale still choose the kinetic or fluid model.
Level 4 · Advanced undergraduate teaching kit
Record the investigation. Teach the reasoning.
A learner-facing lab record and a course-specific instructor guide turn the live model into a repeatable classroom investigation.
Learner record
Plasma regime and gain-chain map
Which dimensionless and length-scale comparisons determine the valid plasma model, and where does fusion gain enter the full facility ledger?
Download learner recordInstructor guide
Teach for evidence, not button pushing
Learners classify a plasma using characteristic scales and keep nuclear, plasma, engineering, and net-electric gain definitions distinct.
Download instructor guideLesson 1 of 3
Charged particles, shielding, and collective scales
When does an ionized gas behave as a plasma rather than a collection of independent collisions?
Particles gyrate in magnetic fields and accelerate along electric fields. Collective behavior appears when many particles occupy a Debye sphere and the system size exceeds the shielding length.
Plasma frequency, cyclotron frequency, collision rate, gyroradius, and Debye length define competing time and length scales. A model is credible only inside its ordering assumptions.
Worked example
How does electron Debye length change if temperature quadruples at fixed density?
- 1. Use λ_D∝√(T_e/n_e).
- 2. Hold n_e fixed.
- 3. Take √4.
The Debye length doubles.
Try it
Plasma regime map
Materials: Parameter sets for laboratory, solar-wind, and fusion plasmas
- 1. Compute characteristic frequencies.
- 2. Compute shielding and gyro scales.
- 3. Compare with device size and observation time.
- 4. Choose kinetic, fluid, or collisional descriptions.
Notice: The same equations do not dominate every plasma; scale ordering chooses the model.
Check your understanding: Does quasineutrality mean charge density is exactly zero everywhere?
Answer: No.
Small charge separations remain essential for fields and waves even when bulk positive and negative densities nearly balance.
Lesson 2 of 3
Plasma waves, transport, and instabilities
How do distribution functions and field feedback create waves or rapidly growing departures?
The Vlasov equation evolves a distribution in phase space under self-consistent fields. Fluid moments give density, flow, pressure, and higher-order closures when kinetic details can be compressed.
Waves and instabilities exchange energy with particles through resonances, gradients, currents, and anisotropy. Transport can greatly exceed binary-collision estimates when turbulence develops.
Worked example
What condition makes a perturbation linearly unstable?
- 1. Assume a mode exp[i(kx−ωt)].
- 2. Allow complex ω=ω_r+iγ.
- 3. Inspect amplitude exp(γt).
- 4. Classify the sign of γ.
γ>0 gives exponential growth; γ<0 gives damping.
Try it
Dispersion-relation explorer
Materials: Plotting notebook
- 1. Choose a simple plasma dispersion relation.
- 2. Vary density and magnetic field.
- 3. Plot real frequency and growth/damping rate.
- 4. Mark where model assumptions fail.
Notice: A dispersion curve is both a prediction and a map of the approximation used to derive it.
Check your understanding: Why can a fluid closure miss Landau damping?
Answer: Landau damping depends on resonant structure in velocity space.
Low-order moments may discard the distribution-function detail that carries the resonance.
Lesson 3 of 3
Fusion reactions, confinement, diagnostics, and gain ledgers
Which measured milestone is being claimed: reactions, plasma heating, facility gain, or delivered electricity?
Fusion rate depends on reactant densities and the velocity-averaged cross section. Magnetic and inertial systems pursue different density-time-temperature regimes but both require impurity, radiation, transport, and stability control.
A neutron or charged-product signal establishes reactions under a detector model. Q_plasma, target gain, wall-plug gain, tritium breeding, component survival, and net electricity are separate ledgers.
Worked example
A system releases 3 MJ of fusion energy after 300 MJ of facility input. What does the facility ledger show?
- 1. State fusion output: 3 MJ.
- 2. State facility input: 300 MJ.
- 3. Compute output/input.
- 4. Keep any target-level or plasma-only definition separate.
Facility energy ratio is 0.01; a narrower scientific gain may still be informative but is not net facility energy.
Try it
Fusion claim ledger
Materials: One fusion press release and its technical source
- 1. Extract every input boundary.
- 2. Extract reaction and energy outputs.
- 3. Identify diagnostic calibration.
- 4. Compute each defensible gain ratio separately.
Notice: Most confusion disappears once every numerator and denominator is named.
Check your understanding: Do statistically credible fusion products prove net-electric operation?
Answer: No.
Reaction evidence and plant energy balance are distinct claims with different measurements.
Formula-to-meaning deck
Read the equation in ordinary language.
λ_D=√(ε₀k_BT_e/(n_e e²))
The Debye length estimates how far electrostatic disturbances penetrate before collective shielding.
ω_pe=√(n_e e²/(m_eε₀))
Electron density and inertia set the natural plasma oscillation frequency.
R_f=n_1n_2⟨σv⟩/(1+δ₁₂)
Fusion reaction rate density combines reactant populations with velocity-averaged nuclear reactivity.
Independent practice
Problem set
Work each problem before opening its hint and solution.
1. How does ω_pe change if electron density increases by a factor of nine?
Reveal hint
Use the square-root dependence.
Reveal solution
It triples.
2. For a mode amplitude proportional to e^{γt}, what does γ=−2 s⁻¹ mean?
Reveal hint
Read the sign and timescale.
Reveal solution
The amplitude damps exponentially with an e-folding time of 0.5 s.
3. A detector efficiency is 20% and records 100 signal counts after background subtraction. Estimate emitted particles under the stated geometry.
Reveal hint
Divide by efficiency before geometry corrections.
Reveal solution
500 particles reaching the detector's accepted channel; total source yield still requires solid-angle and transport corrections.
Derivation studio
Build the result, line by line.
Keep the assumptions visible so the mathematics remains auditable.
Starting point
Electron plasma frequency
Displace the electron fluid by x against a fixed ion background
- 1. The displacement creates surface charge and a restoring electric field.
- 2. Use E=−n_e e x/ε₀ for the slab model.
- 3. Apply m_e ẍ=−eE with consistent sign.
- 4. Identify the oscillator coefficient.
ẍ+ω_pe²x=0, ω_pe²=n_e e²/(m_eε₀)
Collective electrostatic oscillation emerges without requiring particle-by-particle collisions.
Starting point
Maxwellian fusion reactivity
Reaction rate is flux times cross section averaged over relative velocities
- 1. Write dR=n_1n_2σ(v)v f_rel(v)d³v.
- 2. Integrate over relative velocity.
- 3. Define ⟨σv⟩.
- 4. Divide by two for identical reactants.
R_f=n_1n_2⟨σv⟩/(1+δ₁₂)
Temperature influences fusion through the full velocity distribution and energy-dependent cross section.
Computational notebook
Turn the model into an experiment.
Fusion parameter and gain laboratory
Which density, temperature, confinement, diagnostic, and efficiency assumptions dominate a proposed fusion system's conclusion?
Inputs
- • Density and temperature histories
- • Confinement and impurity parameters
- • Reactivity table and diagnostic efficiencies
- • Heating and facility-power ledgers
Algorithm
- 1. Compute characteristic plasma scales.
- 2. Integrate expected reaction yield.
- 3. Propagate detector response and background.
- 4. Calculate plasma, target, facility, and electric gain separately.
Evidence to produce
- • Regime and scale table
- • Predicted-versus-observed product yield
- • Boundary-labeled gain ledger and sensitivity chart
Paper-reading studio
Interrogate the source, not its reputation.
Reconstruct the assumptions, reproduce one calculation, and stop at the boundary of the reported evidence.
Fusion-evidence reconstruction
Does the source establish plasma conditions, reaction products, nuclear heat, system gain, or a power-plant pathway?
- 1. Map diagnostics to each claimed quantity.
- 2. Reproduce detector efficiency and background treatment.
- 3. Name every energy-system boundary.
- 4. Separate measured output from modeled extrapolation.
Calculation to reproduce: Recompute one reaction-rate, triple-product, detector-yield, or gain estimate with propagated uncertainty.
Evidence boundary: Nuclear products, excess heat, plasma gain, and net-electric generation are progressively stronger and distinct claims; none should silently substitute for another.
Continue into the evidence
Source-linked next reading
Chapter 9: Plasma and field effects
Collective plasma structures and the boundary between observation and propulsion interpretation.
Lecture 12: LCF and the integrated picture
Reaction signatures, mechanisms, calorimetry, and net-energy claims kept separate.
Chapter 12: Lattice confinement fusion
A measurement-centered nuclear case study across solid-state environments.