Illustrated Physics
physics, drawn
Series

Plasma oscillation — the series

12 essays on one idea, from the one that introduces it to the one that assumes the rest.
  1. One line, 21 decades of density, and every mirror on it. Plasma frequency against electron density, both logarithmic, with the horizontal rules at frequencies a reader already has a feel for. A wave is reflected by everything to the right of where its own rule meets the line and passes through everything to the left. a fluorescent tube at 1.0e+17 m⁻³ cuts off at 2.84 GHz; ionosphere, D layer at night at 1.0e+8 m⁻³ cuts off at 90 kHz; ionosphere, E layer at 1.0e+11 m⁻³ cuts off at 2.84 MHz; ionosphere, F2 layer at 1.0e+12 m⁻³ cuts off at 8.98 MHz; aluminium's conduction electrons at 1.8e+29 m⁻³ cuts off at 3819.9 THz. So the F2 layer turns back a 1 MHz broadcast and lets a 100 MHz one straight out, which is why one of them is heard across an ocean at night and the other stops at the horizon; and aluminium's cutoff sits in the far ultraviolet, which is why it is a mirror for everything visible and a window above 78 nm. The dependence is a square root, so the line has slope one half: a hundredfold denser plasma reflects only ten times the frequency. Arriving at an angle helps by a factor of sec θ — 1.00 at 0°, 1.15 at 30°, 2.00 at 60° — because only the component of the motion along the density gradient has to be turned round.

    The frequency below which nothing gets in

    Free charges give a medium a permittivity that is negative, and a negative permittivity is not an absorbing medium — it is one in which no wave exists at all. Below that frequency the reflection is total, exactly rather than nearly, because there is no transmitted wave and nothing to absorb. The same expression puts the number at 9 MHz for the ionosphere and 3.8 PHz for aluminium.

    part 1 · astrophysics
  2. A click that arrives sorted by pitch. Arrival time against frequency for a broadband pulse travelling 40 megametres along a magnetic field line through a plasma of 1e+8 electrons per cubic metre in a field of 300 nanotesla — the magnetosphere, roughly. The gyrofrequency is 8.4 kilohertz and the plasma frequency 90 kilohertz, so this is the regime where a right-hand circularly polarised wave travels happily below both of them. It does not travel at one speed. The group velocity is 2c√ω(ω_c − ω)^(3/2)/(ω_p ω_c), which is zero at zero frequency and zero again at the gyrofrequency and peaks in between — at 2.10 kilohertz here, which is a quarter of the gyrofrequency, found by searching the drawn function rather than quoted. Everything either side of that peak is slower, so the curve has a nose: 3.21 s at 500 Hz, 2.20 s at 2000 Hz, 3.59 s at 5000 Hz. The branch below the nose is the classic whistler — a lightning stroke in one hemisphere reaching a receiver in the other as a note gliding downward over about a second, which is what gave the phenomenon its name — and the branch above it arrives as a rising tone at the same time. Both are observed, and a recording that shows the two joined at the nose is how the gyrofrequency along the path is read off directly. Read the other way it is an instrument: the product of the delay and the square root of the frequency is 79 s·Hz^½ here, nearly constant across the low end of the band, and it measures the electron content of a path through the magnetosphere that nothing else could reach.

    The whistle that arrives sorted

    A lightning stroke in one hemisphere reaches a receiver in the other as a note gliding downward over about a second. Nothing dispersed the sound; there was no sound. A radio pulse travelled forty megametres along a magnetic field line through a plasma whose group velocity depends on frequency, and arrived with its frequencies separated by up to three seconds.

    part 2 · astrophysics
  3. The charge a plasma hides. The potential around a charge in a plasma, as a fraction of the bare Coulomb potential at the same distance, against distance measured in screening lengths. The electrons crowd toward the charge and the ions move away until the rearrangement cancels the field, and what is left falls as exp(−r/λ_D) on top of the ordinary 1/r. At one screening length the potential is already down to 37 per cent of the bare value, at three to 5 per cent, and at ten to 4 × 10⁻⁵ — so a charge in a plasma is invisible beyond a few λ_D, and the long range of the Coulomb force, which is what makes electrostatics awkward everywhere else, is simply gone. The screening length at 10¹¹ m⁻³ is 3.8 mm at 300 K, 6.9 mm at 1000 K, 21.8 mm at 10000 K, rising as the square root of the temperature because a hotter electron is harder to hold in place. Drawn this way the three curves coincide exactly: the shape is universal and the only thing a plasma's density and temperature decide is the length written on the axis.

    The long-range force that does not reach

    The Coulomb force falls off as slowly as gravity does, which is what makes electrostatics awkward — every charge is in principle in contact with every other. Put the same charge into a plasma and it becomes invisible beyond a few millimetres, because the mobile charges around it rearrange until its field is cancelled. What is left is a screened potential with a definite range, and that range is what makes a plasma a plasma.

    part 3 · astrophysics
  4. A wave that dies with nothing to rub against. The electric field of a plasma wave at kλ = 0.5, against time in plasma periods, on a logarithmic scale, obtained by integrating the collisionless kinetic equation as an initial-value problem. There are no collisions in the equation, no viscosity and no resistance; the only operator acting on the distribution is a rotation of phase whose rate depends on the particle's speed. The field nevertheless falls exponentially, at 0.1534 per plasma period, against the published root of the kinetic dispersion relation at this wavenumber, 0.1534, and Landau's asymptotic formula's 0.1514. Meanwhile the free energy of the perturbation — the weighted norm of the distribution plus the field energy, which the equation conserves exactly — moves by 2.4e-10. So nothing has been dissipated: every joule the field loses is still in the distribution, and the accounting closes to a part in ten thousand million. The energy has gone into the particles' ordered motion, and the information about the wave is wound into structure at finer and finer scales in velocity.

    The wave that dies with nothing to rub against

    Every damping in this collection so far removes energy from a wave and puts it somewhere warmer. This one removes it and produces no heat at all: there are no collisions in the equation, the entropy is unchanged, the whole thing runs backwards perfectly, and the wave still dies exponentially. What it dies into is structure in velocity too fine for a field to see.

    part 4 · astrophysics
  5. The one place in a plasma where the charges do not balance. The potential and the two densities through a sheath, in Debye lengths from the sheath edge, for a plasma of hydrogen ions. The potential is measured downward in units of the electron temperature and reaches 2.84 at the wall — the value at which the two fluxes balance. The layer is 15.1 Debye lengths thick — for 3 eV electrons at 1e+16 per cubic metre, a Debye length of 129 µm and a sheath of 1.9 mm. Inside it the electron density falls as the Boltzmann factor while the ion density falls only as the ions speed up, so the two separate: at the wall the ion density exceeds the electron density by 85 per cent of itself. That is the whole of the difference between a sheath and the plasma it borders — quasineutrality is not an assumption that can be made here, and Poisson's equation has to be solved instead of replaced.

    The wall a plasma builds against itself

    A plasma is quasineutral everywhere except where it touches something. Electrons are faster than ions by the square root of the mass ratio, so any surface is struck by far more of them, charges negative, and goes on charging until the two arrivals balance — leaving a layer a few Debye lengths thick in which the charges do not cancel at all.

    part 5 · astrophysics
  6. The electrode that walks itself negative. An electrode driven through a blocking capacitor by a symmetric waveform of amplitude 20 kTe, in an argon-like plasma whose electrons are collected 108 times as readily as its ions at saturation, starting uncharged. Potentials are relative to the plasma, in electron temperatures. On the first positive half-cycle the electrode draws a flood of electrons and nothing on the negative half-cycle can return the charge, because the ions arrive at a fixed, small rate; the electrode walks down, and within 6 cycles its mean potential is within half an electron temperature of the self-bias. That bias, found by integrating the charging, is −22.27 kTe, and the charge balance with a Bessel function gives −22.27. The floating potential of the same wall undriven is −4.68. At 3 eV, a 60 V drive makes a 67 V bias with no direct current supplied anywhere.

    The bias no battery supplies

    Drive an electrode in a plasma through a capacitor with a perfectly symmetric waveform and it charges to a steady negative voltage almost as large as the waveform's own amplitude. No direct current flows anywhere and no battery is connected. The plasma's sheath lets electrons in far more easily than it lets them out, and every semiconductor wafer etched in the last forty years was held at a voltage made that way.

    part 6 · astrophysics
  7. The self-bias set by a phase alone. The dc self-bias against the phase between the fundamental and its second harmonic, each of amplitude A, for three chambers: symmetry parameter ε = 0.5, 1, 2, where ε = 1 is a chamber with equal sheaths on both sides. In the symmetric chamber the bias runs from −0.438 A at 0° through zero at 45° to 0.438 at 90°, and back — −43.8 V to 43.8 V at A = 100 V — staying within 0.011 A of a straight line between 0° and 90°. An asymmetric chamber adds a fixed offset of the kind unequal electrodes produce and leaves the phase control in place, so the geometry and the waveform add rather than compete.

    The bias a waveform can choose

    A plasma between two identical electrodes, driven by a sine wave, has no reason to prefer either side, and it does not: no direct voltage appears. Add the second harmonic and it does, although the chamber is still symmetric, no direct current flows and the power spectrum does not depend on the phase between the two. The phase alone sets a bias of up to 0.44 of each harmonic's amplitude, of either sign, and moves the energy of the ions from one electrode to the other while their number stays put.

    part 7 · astrophysics
  8. Rays at a frequency the ionosphere cannot stop head-on. Ray paths at 14 MHz launched at angles from 10° to 75° from the vertical into a Chapman layer peaking at 300 km with a critical frequency of 8 MHz, over a flat Earth. A wave sent straight up at 14 MHz passes through: it is above the layer's plasma frequency everywhere. A ray at φ₀ from the vertical turns back where the plasma frequency equals f cos φ₀ — the secant law — so only rays more oblique than arccos(8/14) = 55.2° return. Those that do land beyond about 1133 km; nothing comes down closer than that: the skip zone, a ring of silence round a transmitter that is audible much further away. Vertical scale exaggerated about six times.

    The frequency the sky returns only at a slant

    The ionosphere turns back any radio wave below its plasma frequency, about eight megahertz on a summer afternoon, and lets anything above it out into space — if the wave is sent straight up. Send a fourteen-megahertz wave up at a slant and it comes back down, a thousand kilometres away. Only the part of the wave's motion across the layer has to be stopped, and that part is the frequency times the cosine of the angle. The rule turns one frequency into a range of them, sets the highest frequency that can reach any given distance, and leaves a ring of silence round every short-wave transmitter, inside which it cannot be heard at all.

    part 8 · astrophysics
  9. A sound that survives only in a hot-electron plasma. The damping of an ion-acoustic wave at kλD = 0.3, as the fraction of its amplitude lost per radian of oscillation, −γ/ω, against the ratio of electron to ion temperature, for hydrogen and argon plasmas, on logarithmic axes: from the kinetic dispersion relation (solid) and from the weak-damping formula (dashed), which is valid only where the damping is small. With equal temperatures the wave loses 0.45 of its amplitude per radian in hydrogen — it is gone within a cycle. At Tₑ/Tᵢ = 10 the loss is 0.031; at 100 it is 0.0126, set by the electrons alone. Argon's floor is lower, 0.0020, because its ions are heavier and its electrons a smaller fraction of a sound wave's speed.

    The sound a plasma carries on its electrons' heat

    Sound in a gas needs two things: something heavy to be pushed, and a pressure to push it. In an ordinary gas the same molecules supply both. In a plasma they can be split — the ions carry the inertia and the electrons, a thousand times lighter and often a hundred times hotter, carry the pressure — and the result is a sound that travels at a speed set by one particle's temperature and the other's mass. It travels only while the two temperatures stay far apart. Warm the ions to match the electrons and the sound dies within a cycle, with nothing to rub against.

    part 9 · astrophysics
  10. Two beams that each stay straight, and wind each other into whirlpools. Phase space of a one-dimensional particle simulation: every dot is an electron, plotted at its position along a periodic box one fastest-growing wavelength long (10.26 v₀/ωₚ) and its velocity, in units of the beam speed. One beam starts at +v₀, the other at −v₀, with a 0.0001 ripple in position. At t = 18/ωₚthe beams are visibly rippled; by 26 the ripple has grown into overturning, and by 40 the two beams have wound round each other into a single vortex, filamented into strands. When growth first stops, at ωₚt = 25.2, the wave's field is 0.56 m v₀ωₚ/e; electrons trapped in it bounce at 0.59 ωₚ, faster than the 0.354 ωₚ at which it was growing, and it traps any electron within ±1.92 v₀ of the wave's speed — zero, by symmetry — so both beams are inside its separatrix; total energy was conserved to 0.31%.

    The two beams that unravel each other

    A cold beam of electrons, alone, is perfectly stable: disturb it and it rings at its plasma frequency and carries the ringing downstream. Send a second beam through it the other way and the pair is not stable at all. A ripple of the right length grows by a factor of e in less than half a plasma period, takes a third of the beams' drift energy into an electric field, and then winds both beams into a single whirlpool in which their drift is gone. Neither beam can do it. The growth exists only in the pair, and a beam a millionth as dense as the plasma it crosses is enough to start it.

    part 10 · astrophysics
  11. The same grains as a crystal and as a liquid. 160 grains repelling one another by a screened Coulomb force, held together by a bowl-shaped confinement as a levitated layer is, and jostled by gas collisions at two temperatures; each picture is the end of a simulation of their motion. Grains with exactly six close neighbours are coloured. Left: strongly coupled — thermal energy about a thousandth of the interaction at the spacing — the grains settle into a hexagonal lattice, 79 per cent of those away from the edge with six neighbours. Right: the same grains with their thermal energy raised seventy-fivefold, the lattice melted, 9 per cent. In the laboratory the grains are a hundred micrometres apart and big enough to film with an ordinary camera, so the melting of a crystal can be watched particle by particle.

    The dust a plasma freezes into a crystal

    Sprinkle microscopic plastic spheres into a glowing discharge and each one charges itself to tens of thousands of electrons, because the plasma's light, fast electrons reach it far more often than its heavy ions. The charge holds the grains up against gravity in the sheath above the electrode, and it makes their mutual repulsion so strong compared with their jiggling that, at room temperature, they freeze into a hexagonal crystal with a lattice spacing a camera can resolve. A plasma, the hottest and most disordered state of matter, is where the most strongly ordered one was found to form most easily.

    part 11 · astrophysics
  12. The product a fusion plasma has to beat. The density, temperature and confinement time multiplied together at which a deuterium–tritium plasma heats itself, against temperature. Multiplying by T turns the previous curve into a pressure times a time, which is what a magnet can actually hold, since the field limits the pressure nT and not the density alone. The curve has a clear minimum of 3.0 × 10²¹ keV s/m³ near 14 keV, and that is why magnetic machines are designed to run hot but not hotter. ITER's design point, a density of 10²⁰ per cubic metre at 8.9 keV and 3.7 s, gives 3.3 × 10²¹: a little under this line, because ITER is meant to be driven by external heating with a tenfold gain and not to ignite outright.

    The product a fusion plasma has to beat

    A star makes less power per litre than a compost heap, and it has burned for four and a half thousand million years. A fusion plasma makes thousands of times more and has to be held for seconds, and it has taken seventy years to get close. Both are obeying one inequality, written down in a classified report in 1955, in which density and time appear only as their product.

    part 12 · astrophysics

All series