Thermionic emission

Release of Electrons from a heated material; one source of Electrons in an X-ray tube.

Heating a metal gives some Electrons enough energy to escape its surface. A heated diode cathode can therefore supply an electric current across a vacuum gap to a collecting diode anode. This makes temperature a control for the supply of Electrons. Conversely, measuring that supply at several temperatures lets us estimate the surface work function: the energy barrier measured from the Fermi level to the nearby vacuum level. A Richardson plot turns the strong exponential function dependence into a slope.

Three questions must be kept separate: how many Electrons cross the surface, how many traverse the gap, and how many the collecting diode anode records. The Richardson-Dushman law answers the first under specified surface assumptions. Space charge can obstruct the second, while incomplete collection can alter the third.

The Richardson-Dushman law models the magnitude of emitted electric current density asHere $T$ is absolute thermodynamic temperature, $W$ is the work function, $k_B$ is the Boltzmann constant, and $A$ is the Thermionic Richardson constant. The signed electric current density points opposite to the outgoing Electron motion; $J$ here is positive. For uniform emission over area $S$, the emitted electric current magnitude is $I=SJ$.

The Richardson-Dushman law assumes an emitter close to thermal equilibrium, a surface with a well-defined and approximately constant work function over the fitted temperature range, and predominantly thermal escape over the barrier. Interpreting a measured collector electric current with this law additionally requires an emission-limited thermionic current, known emitting area, and negligible return emission from the collector. Strong electric fields can lower the surface barrier or enable quantum tunnelling; then this zero-field expression is insufficient.

For a deliberately chosen teaching model, take $W=2.0\mathrm{eV}$ and $A=120\mathrm{A}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{K}}^{-2}$, expressing the barrier in electronvolts. These are illustrative parameters, not measured properties of a named material. Using $k_B=8.617333262\times 10^{-5}\mathrm{eV}{\mathrm{K}}^{-1}$, at $1000\mathrm{K}$ the Boltzmann factor is about $8.33\times 10^{-11}$, yetA 20 percent increase in temperature raises the predicted electric current density by about 69 times. For a chosen emitting area $S=0.010\mathrm{c}{\mathrm{m}}^2$, these correspond to approximately $0.100\mathrm{mA}$ and $6.89\mathrm{mA}$ if collection is complete. The small escaping fraction is replenished by the external electrical circuit.

The sensitivity is especially clear from the logarithmic derivative at fixed $A,W$:In the chosen model at $1000\mathrm{K}$, a $1\mathrm{K}$ rise changes $J$ by about 2.5 percent to first order. Accurate temperature measurement matters even when electric current measurement is easy.

The coefficient $A$ has dimensions of electric current density divided by squared temperature. For ideal nonrelativistic particles with the free electron mass, two spin states, a planar uniform surface, and unit escape transmission, the thermionic tail-flux derivation giveswhere $e$ is the positive elementary charge and $h$ is the Planck constant. The two area units differ by a factor of $10^4$. A fitted $A$ need not equal $A_0$: reflection at the surface, the electronic band structure, and uncertain emitting area affect the inferred coefficient. A constant escape transmission $\tau$ would give $A=\tau A_0$ in this simple model; an energy-dependent transmission generally requires an integral instead.

The name joins Owen Willans Richardson with Saul Dushman. Dushman's 1923 paper develops a thermodynamic treatment of Electron emission analogous to evaporation and evaluates the prefactor using quantum ideas.^[ref]

Choose a fixed positive reference $A_{\mathrm{ref}}$ with the same units as $A$. Taking the natural logarithm of a dimensionless ratio givesThus the Richardson plot is a straight line when $A$ and $W$ are constant, with slope $-W/k_B$ and intercept $\ln (A/A_{\mathrm{ref}})$. The familiar notation $\ln (J/T^2)$ silently chooses units; the reference makes that choice explicit. Changing $A_{\mathrm{ref}}$ shifts the intercept and leaves the inferred work function unchanged. If the horizontal axis instead uses $1000\mathrm{K}/T$, its slope is $-W/(1000\mathrm{K}{k}_B)$.

For the chosen $2.0\mathrm{eV}$ example, choose $A_{\mathrm{ref}}=120\mathrm{A}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{K}}^{-2}$. Then $y(1000\mathrm{K})=-23.2090$ and $y(1200\mathrm{K})=-19.3409$. Their slope is approximately $-23209\mathrm{K}$, giving $W=2.0\mathrm{eV}$. Two temperatures determine a line algebraically; several temperatures allow curvature and drift to be detected. These values are calculated, rather than experimental observations.

A straight trace alone does not establish a unique surface work function. Adsorption, changing surface coverage, or a mixture of regions with different barriers can change the inferred slope. For independent emitting patches with area fractions $f_i$, constant $A_i$ and barriers $W_i$,Lower barriers contribute disproportionately. Even with fixed patches, the slope of the Richardson plot gives a temperature-dependent weighted mean of the $W_i$, rather than their simple area average. If $A$ or $W$ themselves vary with temperature, the apparent slope-derived barrier becomesCheck collection, temperature calibration, and surface stability before assigning curvature to a new emission mechanism.

Consider a planar surface with outward normal along $z$. For an auxiliary equilibrium electron gas just outside it, let $n_v$ be the total number density, counting both signs of $v_z$. Its Maxwell-Boltzmann velocity distribution has normal-component densityAn Electron with $v_z>0$ crosses more area per unit time when it moves faster. The outgoing number flux therefore uses $v_z{f}_z$, rather than just the probability of $v_z>0$:This integral contributes only $T^{1/2}$. Holding a classical bulk number density fixed and imposing a normal kinetic energy threshold $W$ would instead give $\Gamma =n\sqrt{k_{B}T/(2\pi m_e)}{e}^{-W/(k_{B}T)}$; that argument alone does not produce $T^2$.

For the ideal surface model, measure external Electron energies from the metal's chemical potential $\mu$. A freely propagating external state has $E-\mu =W+p^2/(2m_e)$. When $W\gg k_{B}T$, its Fermi-Dirac distribution reduces to the Maxwell–Boltzmann statistics tail. Counting two spin states per phase space cell of volume $h^3$ givesEach of the three Gaussian integrals supplies $T^{1/2}$. Combining this $T^{3/2}$ state count with the $T^{1/2}$ crossing factor givesThe auxiliary atmosphere supplies an equilibrium state-counting argument; an operating collector removes the outgoing Electrons. It need not maintain a real equilibrium gas throughout the gap. Unit transmission and the ideal free Electron state count are substantive assumptions, so this calculation does not determine an arbitrary real surface's measured $A$.

In particular, Electrons inside a metal generally form a degenerate Fermi gas. A classical distribution for the escaping tail does not imply a classical distribution for the whole interior. Richardson himself emphasized this distinction in his lecture.^[ref]

An emission-limited thermionic current is set by the rate at which the heated surface supplies Electrons, with nearly complete collection. At fixed temperature, increasing the attracting electric potential difference then produces an approximate collection plateau. That plateau is the useful regime for fitting the Richardson-Dushman law, provided the electric field has not significantly lowered the surface barrier.

The emitting surface still matters in this regime. Adsorption changes the work function; different crystal faces can have different barriers. Heating may desorb contaminants or change a coating, so repeated measurements at the same temperature need not agree unless the surface preparation is controlled. For fixed $A,T$, a barrier change $\Delta W$ multiplies $J$ by $e^{-\Delta W/(k_{B}T)}$. In the chosen model at $1000\mathrm{K}$, raising the barrier by $0.10\mathrm{eV}$ reduces $J$ to about 0.313 of its previous value.

A space-charge-limited thermionic current occurs when the negative space charge of emitted Electrons creates an electric field that obstructs further transport across the gap. Raising the emitter temperature may then increase the surface supply without producing the corresponding collector electric current. The gap geometry and collecting electric potential difference matter even when the surface is stable.

To distinguish this from an emission-limited thermionic current, vary the collecting electric potential difference at fixed emitter temperature, then repeat at another temperature. A voltage-dependent transport region followed by a collection plateau separates a gap constraint from a supply constraint. A weak collector electric current alone cannot establish a high work function. Nor does simply choosing a very large collecting electric potential difference guarantee a valid zero-field measurement: the surface barrier may then change.

In his Nobel lecture, delivered on 12 December 1929 for the 1928 prize, Owen Willans Richardson recalled heating wires for weeks while hand-pumping away gas released by the apparatus. With platinum, he found a finite emission per unit area that rose rapidly with temperature. His initial 1901 formula used $T^{1/2}$; he preferred the $T^2$ form in 1911.^[ref]

The early apparatus surrounded a heated platinum wire with a metal cylinder connected to earth through a galvanometer. Wire electrical resistance supplied the temperature readout. After checking collection saturation, Richardson measured at a wire potential of $-120\mathrm{V}$, over roughly $1300$–$1600\mathrm{K}$; reversing polarity suppressed the negative discharge.^[ref]

The inference rests on several controls. A temperature sweep probes thermal escape; reversing the collecting electric field tests the sign of the mobile electric charge; varying collection voltage tests whether transport is limiting. Prolonged pumping and heating test whether a transient gas release is responsible. These controls address different alternatives, and none by itself establishes a clean, uniform surface.

For a modern repetition, stabilize each temperature, measure the collector electric current against collecting electric potential difference, and record the pressure and surface preparation. Repeat selected temperatures after cooling to reveal drift. Use the collection plateau and the actual emitting area to construct a Richardson plot. The resulting slope estimates the work function only within the stated surface and transport assumptions; it does not turn an uncontrolled wire into a reference material.