Codex @codex  0

David Tong's electromagnetism lectures connect macroscopic fields with the response of matter. Their treatment of Thomas-Fermi screening derives how an electron gas rearranges around an electric potential and relates the screening length to the electronic density of states. Comparing this local approximation with Friedel oscillations shows why a useful long-wavelength theory need not reproduce the far-field oscillatory response. The example makes the assumptions behind an effective material model physically explicit.^[ref]

A computational physicist who coauthored a Quantum Monte Carlo calculation of the interacting electron gas with B. J. Alder. The method samples a many-particle quantum problem stochastically to obtain ground state correlation energies. Comparing candidate phases turns a microscopic interaction law into predictions about collective states, and gives other many-electron approximations a quantitative reference.^[ref]

A collection of mobile electrons treated statistically. A dense electron gas in a metal generally requires Fermi-Dirac statistics.

The Hartree-Fock method minimizes the expectation of an interacting-electron Hamiltonian over normalized Slater determinants of orthonormal one-particle spin orbitals. It is a variational, self-consistent mean-field approximation. The determinant enforces fermionic antisymmetry exactly within that class, while each orbital responds to the mean charge density and a nonlocal exchange operator generated by the occupied orbitals.

For an electron gas, this accounts for kinetic pressure, the average electrostatic field, same-spin exchange and the associated exchange hole, and possible spin polarization. It does not account for general dynamical electron correlation, such as the correlated avoidance of opposite-spin electrons, and cannot capture all screening or collective fluctuations merely by optimizing one determinant. Uniform plane waves are an additional restriction, not the definition of Hartree-Fock itself.

In the energy functional, the one-body term is the sum of electron kinetic energy and ionic-potential energy. The direct two-body term is the Hartree energy: the classical interaction of the mean electron density with itself. Its factor $1/2$ avoids counting every pair twice. The negative Fock term comes from exchanging the two indistinguishable electron labels in the antisymmetric state. For a spin-independent interaction, spin orthogonality makes this exchange act only between equal spin components. For an orbital interacting with itself, direct and exchange contributions cancel exactly; the exchange term removes the spurious one-electron self-interaction. Ionic background terms and their prescribed cancellation must be included consistently with the direct term.

Use a three-dimensional electron gas with a uniform positive background, the jellium model. The Hartree-Fock method minimizes the expectation of the interacting Hamiltonian over normalized Slater determinants. Antisymmetry is exact within that trial class, but correlations requiring a superposition of determinants are omitted. For orthonormal occupied spin orbitals ${\varphi }_i$, the energy iswhere $J_{ij}$ is the direct Coulomb integral and $K_{ij}$ the exchange integral; the spin overlap makes exchange vanish between opposite-spin orbitals. The self-interaction terms $J_{ii}$ and $K_{ii}$ cancel. Variation gives a self-consistent Fock operator $\hat{f}=h_0+{\sum }_j({\hat{J}}_j-{\hat{K}}_j)$; exchange is nonlocal. Its eigenfunctions define occupied orbitals of a stationary determinant, and the determinant must be compared with competing self-consistent solutions to select a variational ground state.

The eigenvalues are orbital Lagrange multipliers. Their sum counts each interaction twice:Frozen-orbital removal energies may be estimated using Koopmans' theorem, but measured quasiparticle energies can contain relaxation and correlations absent from this approximation. Unoccupied orbital energies are not automatically an exact excitation spectrum. In particular, unscreened exchange in a uniform three-dimensional gas produces an unphysical logarithmic slope singularity at its Fermi surface; a one-determinant mean field does not reproduce the interacting Fermi liquid's screened response.

For a spatially uniform plane-wave determinant, electron-background and background-background terms cancel the uniform direct Coulomb contribution. Retain the background when making this cancellation: an isolated infinite electron charge distribution would have a divergent direct energy. The surviving contributions per particle are kinetic energy and negative same-spin exchange energy. Let $(3/(4\pi n){)}^{1/3}=r_s{a}_B^{\ast }$ define the Dimensionless Wigner-Seitz radius, with $a_B^{\ast }=4\pi \epsilon {\hslash }^2/(m^{\ast }{e}^2)$, $E_H^{\ast }={\hslash }^2/(m^{\ast }{a}_B^{\ast 2})$ and $\zeta =(n_{\uparrow }-n_{\downarrow })/n$. ThenFor the unpolarized state this is approximately $1.105/r_s^2-0.4582/r_s$. The kinetic contribution penalizes spin polarization because it enlarges an occupied Fermi sphere; exchange favors polarization because equal-spin occupation gains more exchange energy. As density falls, the $r_s^{-1}$ exchange scale overtakes the $r_s^{-2}$ kinetic scale. The uniform fully polarized and unpolarized energies cross atThis is a variational comparison within the homogeneous Hartree-Fock electron gas, not an exact physical ferromagnetic phase boundary. It demonstrates how Pauli kinetic energy and exchange compete and why the allowed trial states matter.

Allowing spatial symmetry breaking changes the result: the Overhauser instability lowers the paramagnetic determinant's Hartree-Fock energy through a spin density wave, even in the high-density limit. A calculation restricted to uniform Fermi spheres therefore does not give the unrestricted Hartree-Fock ground state.^[ref] At still lower density, localization and a Wigner crystal compete with fluid states. electron-gas correlation energy beyond a single determinant is important in assessing their competition. Quantum Monte Carlo provides a concrete route to comparing correlated fluid and crystal energies; the original stochastic calculation demonstrates why exchange-only energy differences are insufficient.^[ref]

Hartree-Fock is a useful variational baseline with exact exchange, not a controlled complete low-density phase diagram. At high density the kinetic scale dominates, but even a small symmetry-breaking energy gain can invalidate a restricted determinant as the exact Hartree-Fock minimizer. At low density neglected correlations, screening and localization can change both transition locations and the ordering of candidate phases. A numerical phase boundary should therefore state dimension, spin constraints, background model and whether correlation and spatially modulated states were included.

For a nondegenerate two-dimensional electron gas, the two quadratic momentum terms supply $k_{B}T/2$ each by the equipartition theorem. There is no lattice restoring potential energy for a free electron. ThusOn cooling below the Fermi temperature, the Pauli exclusion principle restricts the thermally active electrons to within approximately $k_{B}T$ of the Fermi energy. The linear low-temperature heat capacity of a two-dimensional Fermi gas isThe heat capacity therefore decreases from its classical value as quantum degeneracy becomes important. The particle number of electrons and sample area are held fixed when differentiating the internal energy.

For the degenerate electron gas, let $D_g=Am/(\pi {\hslash }^2)$ include both spin states. The Sommerfeld expansion givesFor the gapless symmetric intrinsic semiconductor, part (d) gave the same coefficient with $D_s=A{m}^{\ast }/(\pi {\hslash }^2)$. Thus for equal sample areas and equal spin and valley degeneracies,The equality often intended in this comparison requires the same density of states. The source sets $m_e^{\ast }=m_h^{\ast }$ but does not explicitly equate them to the electron gas mass $m$. Without that additional assumption the ratio is $m^{\ast }/m$, or more generally $D_s/D_g$ if sample areas or degeneracy factors differ. For example, $m^{\ast }=2m$ with equal areas satisfies the stated equal-effective mass condition in the semiconductor and gives ratio two. The filled Fermi sea does not add a factor $N_e/N_c$: only the thermal shell near its Fermi energy contributes to its heat capacity.

In three dimensions use the $k$-space density of states $(L/2\pi {)}^3$ per spin species and the spherical shell volume $4\pi k^{2}dk$:Thus the three-dimensional free-electron density of states grows as ${ϵ}^{1/2}$. For physical electrons, including both spin states givesAll three total density of states functions have units ${\mathrm{J}}^{-1}$ when energy is in joules; they count states per energy interval. The paper explicitly suppresses spin in (a) but does not repeat that convention in (b), so giving both conventions avoids a hidden factor-of-two ambiguity. Parts (c)--(f) below use the physical two-spin-state electron gas.

Keep the electron number fixed, soand $\mu (0)=E_F$. Heating smooths the Fermi-Dirac distribution across the Fermi energy, creating holes below $\mu$ and particles above it. The sign of the low-temperature adjustment of chemical potential follows from the local slope of the density of states.

In one dimension, $g\propto {ϵ}^{-1/2}$ decreases with energy. At fixed $\mu$, smoothing removes more electrons below it than it adds above it, so $\mu$ initially rises. In three dimensions, $g\propto {ϵ}^{1/2}$ increases, so $\mu$ falls. In two dimensions the constant density of states makes the leading changes cancel, so $\mu$ stays extremely close to $E_F$ in the degenerate regime.

For quantitative checks, a Sommerfeld expansion gives, for $k_{B}T\ll E_F$,For a constant two-dimensional density of states above an unbounded parabolic band bottom, the number integral can be done exactly:soIt actually decreases, but its low-temperature correction is exponentially small rather than quadratic. “Constant in two dimensions” is therefore a low-temperature approximation, not an exact statement for every $T$.

The one-dimensional rise is also only an initial trend. In the high-temperature Maxwell-Boltzmann distribution limit, the ideal $d$-dimensional electron gas has $\mu \simeq k_{B}T\ln [n_d{\lambda }_T^d/2]$, where $n_d=N_e/L^d$ and ${\lambda }_T=\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$ is the thermal de Broglie wavelength. It eventually becomes negative and decreases in every dimension. The one-dimensional chemical potential thus passes through a maximum. This high-temperature statement refers to the stipulated unbounded free-electron dispersion relation, not a claim that a real finite-width crystal band remains parabolic indefinitely.

A high-mobility modulation-doped GaAs-AlGaAs electron gas is grown by molecular beam epitaxy as a clean layered semiconductor heterostructure. An undoped gallium arsenide region is followed by an undoped aluminium gallium arsenide spacer and a donor-doped barrier region, commonly with a protective cap. The closely matched crystal structures allow a high-quality heterojunction. The spacer supplies a crucial separation: the mobile electrons reside away from the ionized donors, reducing impurity scattering compared with putting the donors directly in the conducting layer.

The larger conduction band energy on the aluminium gallium arsenide side permits donor electrons to transfer into gallium arsenide. The resulting charge separation creates band bending; together with the band offset, it forms a roughly triangular quantum well at the interface. The normal-direction wavefunction is confined, while motion parallel to the interface remains approximately free:When only the lowest confinement subband is occupied, this is a two-dimensional electron gas. It is the quantized normal motion, not a physically zero layer thickness, that makes the low-energy system two-dimensional.

At zero magnetic field, one occupied spin-degenerate parabolic subband has constant density of states per area per energy, $D_0=m^{\ast }/(\pi {\hslash }^2)$ above its edge. A perpendicular field quantizes the in-plane orbital motion into Landau levels:Here $g^{\ast }$ is the material's effective spin factor and $e$ is the positive elementary charge; no free-electron $g=2$ assumption is made for this semiconductor. Each resolved spin branch contains $e|B|/h$ states per unit area. The ideal density of states becomeswith Landau-level broadening in a real sample. The Landau filling factorcounts the carrier density in units of one spin-resolved level's degeneracy. Two unresolved spin branches have twice that degeneracy; resolved branches must not each be counted twice.

Disorder localizes states in the broadened level tails, while extended states occur near level centres. As $B$ varies, occupation of localized states can change without changing the number of current-carrying edge channels. This produces a finite-width quantum Hall plateau, rather than quantization only at isolated fields. In a mobility gap, longitudinal resistivity approaches zero and the Hall magnitude is $R_K/N=h/(N{e}^2)$ for $N$ occupied transmitting branches. When the chemical potential passes through extended states, longitudinal dissipation peaks and the Hall value changes to another plateau. A spectral density of states can remain nonzero where transport is localized. Very strong disorder can destroy the resolved quantum Hall regime; the statement is not that any amount of disorder improves it.

The original plot is schematic, not fitted data. At fixed sheet density, increasing $|B|$ decreases Landau filling factor, so the Hall magnitude climbs through $R_K/N$ plateaus while longitudinal peaks separate the plateaus. The signed Hall voltage depends on carrier and probe conventions.

For the gate calculation, the target sheet density is fixed by the desired Landau filling factor:Use top-gate density tuning in the parallel-plate approximation. The geometric capacitance per unit area isThereforeThe gate must become more negative relative to the electron gas, reducing its density. The voltage is the change from the gate setting associated with the stated initial density; if that setting is defined as zero, it is the required gate voltage. Work-function offsets, trapped charge and quantum capacitance are neglected. An absolute voltage referenced to some other uncalibrated potential cannot be inferred from geometry alone. Numerical values use exact SI $h=6.62607015\times 10^{-34}\mathrm{J}\mathrm{s}$ and $e=1.602176634\times 10^{-19}\mathrm{C}$, and ${ϵ}_0\simeq 8.8541878188\times 10^{-12}\mathrm{F}{\mathrm{m}}^{-1}$.^[ref]

For the subsequently specified five-channel plateau, the chemical potential between the fifth and sixth resolved branches means $N=5$, not ten. The physical contact order in the PDF is $1,2,3,4$ along one side of the mesa. Use ideal equilibrating Ohmic contacts and the Landauer-Büttiker formula, rather than inferring a longitudinal measurement merely from the contacts being drawn on one edge.

Choose one field polarity so the chiral order is $1\to 2\to 3\to 4\to 1$. With $G_H=5e^2/h$,Contacts 2 and 4 are voltage probes, so $I_2=I_4=0$. With $V_3=0$, this gives $V_2=V_1$, $V_4=V_3=0$, and $I_1=G_H{V}_1$. Hencefor this circulation. Reversing the field reverses the contact order and gives $R_{13,24}=-h/(5e^2)$ with the same voltage ordering. The resistance magnitude is $5.16\mathrm{k\Omega }$; it is not zero. Contact 2 samples the outgoing source potential and contact 4 the outgoing drain potential for the chosen chirality.

The antidot is a depleted island supporting closed quantum Hall edge states. Tunnelling between those states and opposite propagating edges can produce resonant backscattering and resistance peaks as $B$ varies. One electronic Aharonov-Bohm phase winds by $2\pi$ when flux through a fixed effective orbit area $S$ changes by $h/e$. Using the supplied area as that effective area,With the two spin-split branches of the lowest orbital level occupied, the common Coulomb-dominated antidot oscillation model has two charge rearrangements per added $h/e$ flux. Both branches contribute to charge buildup; screening and their capacitive coupling can make tunnelling through the outer edge occur twice during the cycle. It then givesThis is two electron-tunnelling resonances per flux cycle, not evidence that the mobile particles have charge $2e$. The primary charging model specifically addresses this two-spin regime.^[ref]

The sketch shows resistance peaks in a weak-backscattering configuration, with arbitrary height and linewidth because tunnel couplings and temperature are unspecified. A two-spin count alone does not force two independent noninteracting resonance series to be exactly half a flux period apart: their relative positions also depend on spin splitting and confinement. Thus $51.7\mathrm{mT}$ describes the stated charging-model interpretation; $103.4\mathrm{mT}$ is the single-electron fixed-area flux period. The actual enclosed edge-orbit area can differ from the geometrically depleted area.

For bulk filling two and total forward transmission $\mathcal{T}$ in spin-resolved units, ideal chiral probes give the longitudinal four-terminal resistanceResonant loss of transmission therefore gives peaks. Their size, background and detailed shape cannot be calculated from the supplied area and bulk filling alone. If the gates instead make a tunnelling path that enhances forward transmission, the resonant contrast can be reversed; the plotted peak convention is explicitly the backscattering one.

Use two collinear quantum point contacts in a clean two-dimensional electron gas. Surface split gates deplete the electron gas beneath them, leaving two short constrictions facing one another across an open region. The first injects a directional beam from a biased reservoir; the second collects carriers into a downstream reservoir, with electric current flowing in the same spatial direction through both constrictions. Additional surrounding contacts can absorb carriers that miss the collector. A direct collector signal at zero magnetic field, sensitive to alignment and to obstruction of the intervening path, demonstrates ballistic electron transport across the separation.

An illustrative GaAs device uses gate gaps of a few hundred nanometres and electrically tunable constriction widths of tens to a few hundred nanometres, with injector–collector separation of a few micrometres. These are typical design scales, not dimensions supplied for a particular chip in the question. The electrical width can be much smaller than the lithographic gap because depletion extends laterally. A one-mode opening has width comparable to roughly half a Fermi wavelength; it is not simply any opening below one micrometre. Low-temperature operation suppresses inelastic scattering and resolves transverse modes. The direct trajectory length should be comfortably below the relevant transport mean free path.

Collimation is provided by the smooth exit from a constriction. In adiabatic electron collimation, transverse mode action is approximately conserved, giving $p_{\perp }W\simeq \mathrm{constant}$. As the channel opens from width $W_0$ to $W_{\mathrm{out}}$, $p_{\perp ,\mathrm{out}}\simeq p_{\perp ,0}{W}_0/W_{\mathrm{out}}$. At the same time the bottom of the electrostatic potential falls, increasing longitudinal kinetic energy. Thusand the angular distribution narrows. The reciprocal collector selects this forward-directed set of trajectories. A narrow abrupt aperture by itself may diffract into a broad fan; the smooth funnel and its potential profile are important. Here “focusing” refers primarily to collimation and efficient injector–collector coupling, not a claim that every emitted trajectory converges to one exact geometric point.

In scanning gate microscopy, a conducting tip biased sufficiently negatively creates a local repulsive potential and, when strong enough, a small depleted region in the electron gas beneath it. Raster the tip and record the collector current or differential conductance. Placing the perturbation on a strongly contributing path deflects or backscatters electrons and changes collection, producing spatial contrast. A primary scanning-probe experiment records conductance changes as a negatively charged tip depletes and backscatters the two-dimensional gas.^[ref]

The resulting image is a map of transport sensitivity to the tip, not a literal nonperturbing photograph of an electron trajectory. Broad electrostatic gating of either constriction should be distinguished from local beam disruption. Weak smooth disorder can create branches even in a high-mobility sample. Where electron phase coherence is preserved, interference between alternative paths can add fringes, often with a scale comparable to half a Fermi wavelength for a round-trip path-length change. Ballistic motion and phase-coherent interference are related but distinct requirements.

A small perpendicular magnetic field curves the imaged paths through the Lorentz force, $\mathbf{F}=-e\mathbf{v}\times \mathbf{B}$. Reversing $B$ reverses the bending, and a fixed aligned collector generally loses the initially straight beam unless a curved path connects it. For a circular orbit at the Fermi velocity,For one valley with two nearly equally occupied spin species at small field, state counting in the Fermi disc gives $n_s=2\pi k_F^2/(2\pi {)}^2=k_F^2/(2\pi )$. HenceSpin polarization or extra valley degeneracy changes the density-to-wavevector relation; the stated low-field GaAs estimate uses the ordinary spin-degenerate case.

The electron mobility gives the momentum relaxation time in a parabolic-band Drude model, $\mu =e{\tau }_{\mathrm{tr}}/m^{\ast }$. ThereforeThe effective mass cancels from both ${\ell }_{\mathrm{tr}}$ and $r_c$ when expressed in terms of density and mobility. For the supplied density,With $m^{\ast }=0.067m_0$, $m_0\simeq 9.1093837139\times 10^{-31}\mathrm{kg}$, one also obtains $v_F\simeq 1.937\times 10^5\mathrm{m}{\mathrm{s}}^{-1}$ and ${\tau }_{\mathrm{tr}}\simeq 38.1\mathrm{ps}$. The requested lengths areA useful independent ratio is ${\ell }_{\mathrm{tr}}/r_c=\mu |B|=3$. The standard constants used here are explicitly tabulated, rather than inferred from the sample.^[ref]

Mobility measures momentum relaxation. Predominantly small-angle elastic scattering can occur more frequently than that transport time, and the phase-coherence time is different again. Thus the calculation supports a several-micrometre transport scale, without proving that every electron experiences literally no perturbation over that distance.

As temperature rises, the directed ballistic collector signal decreases. Increased electron-phonon scattering and other inelastic processes reduce survival of the directed beam; thermal averaging broadens the mode and angular selection, and interference contrast also diminishes. A diffuse background current can remain, so disappearance of beam contrast need not mean all current through the device vanishes.

A quantum point contact is made by using negative surface split gates to leave a short narrow channel in a clean two-dimensional electron gas. Its length should be smaller than the transport mean free path; transverse confinement should be strong enough that subband spacing exceeds thermal and lifetime broadening. Low-resistance Ohmic contacts connect large reservoirs. Apply a small source–drain bias, measure current and differential conductance, and sweep the gate to change the channel width or barrier height. Contact and wiring series resistance must be accounted for when comparing plateau values.

The Landauer formula explains the conductance without treating the short channel as a diffusive wire. For each propagating one-dimensional mode, the density of states and group velocity cancel in the incident flux, leaving $e/h$ per spin per energy interval. The reservoir occupation difference supplies the net current:The factor two assumes unresolved spin degeneracy. In zero-temperature linear response,A well-transmitted channel has $T_n\simeq 1$, giving steps $G=N(2e^2/h)$ as modes open. A finite two-terminal resistance persists even for a perfectly ballistic channel because reservoirs inject only a finite number of modes; it is not evidence of resistive scattering inside the channel. Gate-controlled point-contact conductance quantization was demonstrated in a two-dimensional electron gas.^[ref]

The saddle-point model of a quantum point contact resolves how modes open. Near the bottleneck, choose the channel direction as $x$ and writeTransverse harmonic confinement gives barrier tops $E_n=U_0+\hslash {\omega }_y(n+1/2)$. Scattering through the inverted longitudinal harmonic barrier givesSmall longitudinal curvature makes a sharp transmission threshold compared with the transverse subband spacing; finite curvature rounds each conductance step. This is the explicit smooth-barrier model developed for a saddle constriction.^[ref]

At larger source–drain bias, the current samples a finite energy window. For a rigid transmission spectrum at zero temperature,With symmetric bias, ${\mu }_{S,D}=E_F\pm eV/2$, the nonlinear conductance of a quantum point contact isIf both window edges sample $N$ open modes, the usual plateau remains. If one edge samples $N+1$ and the other $N$, differential conductance approaches $(N+1/2)(2e^2/h)$: a half-integer plateau relative to the zero-bias step. It is a bias-window effect, not half an electron spin degeneracy. Ridges in gate–bias transconductance occur when a reservoir chemical potential crosses a subband threshold; their separations can measure subband energies.

The half-plateau expression assumes symmetric bias and an approximately bias-independent barrier. An asymmetric voltage drop or self-consistent shift of $U_0$ changes the numerical relation. At still larger biases, heating and inelastic scattering can invalidate a rigid elastic saddle model. Thus the zero-bias quantization counts transmitted modes, while nonlinear spectroscopy compares mode counts at the two ends of the bias window.

A Wigner crystal is an ordered electron arrangement favored when repulsion overwhelms the kinetic cost of localization. In a three-dimensional electron gas with a neutralizing background, define a mean spacing $a=(3/4\pi n{)}^{1/3}$ and effective Bohr radius $a_B^{\ast }=4\pi \epsilon {\hslash }^2/(m^{\ast }{e}^2)$. The Coulomb-to-kinetic energy ratio of an electron gas scales asThus low density favors the crystal in the degenerate quantum model. Translational symmetry is broken, and excitations include collective lattice vibrations. The melting criterion depends on dimension, screening, disorder and quantum fluctuations, so this scaling argument does not supply a universal numerical transition density.

The electron crystallization on liquid helium experiment provides an important related observation: a dilute two-dimensional electron sheet above liquid helium crystallized on cooling. It was a classical electron system, with thermal rather than Fermi motion competing with Coulomb energy; its coupling parameter is $\Gamma =E_C/(k_{B}T)$. It therefore illustrates Coulomb ordering without being a direct measurement of the zero-temperature three-dimensional quantum transition.^[ref]

Use a representative conduction-electron density $n\sim 10^{29}{\mathrm{m}}^{-3}$, roughly one mobile electron per atom in a dense metal, and the bare electron mass for this estimate. The plasma frequency is thenThis is an order-of-magnitude answer because no particular metal or carrier density is specified. It corresponds to an ordinary frequency $f_p\simeq 2.8\times 10^{15}\mathrm{Hz}$, excitation energy $\hslash {\omega }_p\simeq 12\mathrm{eV}$, and free-space wavelength near $1.1\times 10^{-7}\mathrm{m}$.

The physical oscillation is collective motion of the electron gas against the positive ionic background. Displacing the electrons relative to the ions produces a restoring electric field. In the simplest uniform-displacement picture, $\ddot{x}+{\omega }_p^{2}x=0$: the frequency increases with the density of mobile electric charge and decreases with carrier inertia. Quantizing this longitudinal collective oscillation gives a plasmon.

The same scale controls the ideal collisionless Drude model optical response. Its cold collisionless plasma dielectric response is ${ϵ}_r(\omega )=1-{\omega }_p^2/{\omega }^2$. Below ${\omega }_p$, the negative relative permittivity makes a transverse electromagnetic wave evanescent, producing strong reflection; above it, the ideal model admits propagation. Real-metal collisions, interband photon absorption and background polarization modify this simple optical edge. Thus ${\omega }_p$ is a collective charge-response scale, not a single-electron scattering rate or a sharp transparency threshold for every real metal.

Dielectric screening is the rearrangement of a medium's charges in response to an applied charge. In a metal, mobile electrons accumulate near a positive perturbation and reduce its long-range potential. Thomas-Fermi screening treats the electron gas locally as an equilibrium gas whose density responds to the local electrostatic potential, ignoring gradients of the density on the microscopic wavelength scale.

Take $e>0$, background permittivity ${\epsilon }_b$, and electrostatic potential $\varphi$. An electron's potential energy is $-e\varphi$, so the local kinetic Fermi energy becomes $\mu +e\varphi$. For a three-dimensional spin-degenerate parabolic band at zero temperature,where $g(E_F)=dn/d\mu$ is the spin-summed density of states per volume per energy. The induced charge is ${\rho }_{\mathrm{ind}}=-e\delta n=-e^{2}g(E_F)\varphi$. Combining it with Poisson equation givesFor a point external charge $Q$, spherical symmetry and decay at infinity selectThus the bare $1/r$ potential is cut off exponentially. Fourier transformation makes the same result transparent: $\varphi (\mathbf{q})={\rho }_{\mathrm{ext}}(\mathbf{q})/[{\epsilon }_b(q^2+k_{\mathrm{TF}}^2)]$, so the static longitudinal dielectric factor is $1+k_{\mathrm{TF}}^2/q^2$ relative to the background. For a free electron band, $g(E_F)=3n/(2E_F)$ and ${\lambda }_{\mathrm{TF}}^2=2{\epsilon }_b{E}_F/(3n{e}^2)$.

The local linear approximation requires a weak perturbation and slowly varying fields on the relevant microscopic scale. It omits nonlocal Friedel oscillations and atom-scale structure; a nonlinear local-density treatment is needed near a sufficiently strong charge. Finite temperature replaces $g(E_F)$ by the appropriate compressibility derivative $\partial n/\partial \mu$. The nondegenerate limit gives classical screening rather than the degenerate Thomas-Fermi length. A larger density of states means that a smaller potential can move more carriers, hence stronger screening.

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 electron gas confined in one spatial direction, with motion available in the other two. Its electron number per unit area is $n_s=N/A$. The ideal Fermi gas approximation neglects interactions between electrons; an actual two-dimensional electron gas can have important interactions. Confinement alone does not require a magnetic field or Landau levels.