| Jellium model | |
|---|---|
| Name | Jellium model |
| Type | Model |
| Field | Condensed matter physics |
| Introduced | 1920s–1930s |
| Inventor | Eugene Wigner (early formulations), J. C. Slater (terminology) |
| Notable examples | Homogeneous electron gas, Wigner crystal |
Jellium model
The Jellium model is a theoretical model in condensed matter physics and quantum mechanics that represents conduction electrons in a metal as a quantum many-body system moving in a uniform positive background charge. It matters because it isolates electron correlation and exchange effects in the simplest translationally invariant setting, providing a benchmark for approximations in density functional theory (DFT), quantum Monte Carlo (QMC) calculations, and the study of collective excitations such as plasmons.
The jellium model idealizes a metal by replacing discrete ionic cores with a continuous positive charge density (the "jelly") that neutralizes the mobile electron charge. This simplification preserves translational symmetry and removes lattice potentials, allowing focus on the interplay of kinetic energy, Coulomb repulsion, and quantum statistics. Historical motivations trace to early quantum treatments of metals by Eugene Wigner and later developments by John C. Slater and others seeking tractable many-electron systems. Jellium underpins conceptual understanding of screening, the Fermi liquid paradigm, and collective modes in electron gases studied at institutions such as Bell Labs and research groups in Cambridge University and Princeton University.
Mathematically, jellium is defined by the Hamiltonian for N electrons in volume V with periodic boundary conditions: H = sum_i (-ħ^2 ∇_i^2 / 2m) + (1/2) sum_{i≠j} e^2 / |r_i - r_j| + background terms. The positive background contributes a constant potential and self-energy terms to ensure overall charge neutrality. Key parameters include the average electron density n and the Wigner–Seitz radius r_s = (3/4πn)^{1/3} in three dimensions. The noninteracting limit corresponds to the Fermi gas, while interactions introduce exchange and correlation energies often expressed via the exchange-correlation functional in DFT. Important exact limits and sum rules derive from canonical quantization, linear response theory, and the random phase approximation (RPA).
The homogeneous electron gas (HEG) is the translationally invariant realization of jellium; variants include two-dimensional (2D) and one-dimensional (1D) electron gases relevant to surface states, quantum wells, and quantum wires. The HEG exhibits phases such as the paramagnetic Fermi liquid and, at low density, the Wigner crystal predicted by Wigner. Spin-polarized (fully or partially) HEG models probe magnetic instabilities and itinerant ferromagnetism studied by theorists including David Bohm and David Pines in collective-excitation theory. Finite-size adaptations impose spherical or slab geometries to model metallic clusters and surfaces.
Jellium serves as a testbed for many computational techniques. The random phase approximation provides analytic screening and plasmon dispersion at high density. Quantum Monte Carlo (QMC), notably variational Monte Carlo and diffusion Monte Carlo, yields benchmark exchange-correlation energies for the HEG; prominent QMC studies were performed by D. M. Ceperley and B. J. Alder. Perturbative approaches include Hartree–Fock for exchange and many-body perturbation theory such as the GW approximation for quasiparticle energies. Density functional approximations—local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGAs—are often parametrized against HEG QMC data to produce practical exchange-correlation functional forms used across computational materials science.
Jellium predicts Fermi surface properties, ground-state energy as function of r_s, and response functions characterizing screening and dielectric behavior. It yields collective excitations (bulk and surface plasmons) and Lindhard-type single-particle response at weak coupling. Correlation energy behavior crosses over from high-density perturbative limits to low-density regimes where correlation dominates and crystallization into a Wigner lattice may occur. Transport coefficients, effective mass renormalization, and quasiparticle lifetimes have been analyzed within jellium-based many-body frameworks and compared with experiments on simple metals and electron-doped semiconductors.
Beyond serving as a conceptual laboratory, jellium underlies practical computational tools: LDA parametrizations in DFT, pseudopotential testing, and surface science models (jellium slabs approximate metal surfaces and work functions). It informs interpretation of electron-energy-loss spectroscopy (EELS), optical properties of simple metals, and behavior of free-electron-like materials such as alkali metals and noble metals under certain conditions. The model guides the design and understanding of low-dimensional electron systems in semiconductor heterostructures and the physics of graphene and oxide interfaces when screening and correlation are prominent.
Jellium neglects lattice periodicity, ionic potentials, and band-structure effects, limiting quantitative applicability to real materials with directional bonding or d- and f-electron physics. It omits electron–phonon coupling and core-level chemistry. Extensions include incorporating pseudopotentials, embedding treatments to reintroduce ionic structure, and multicomponent jellium for alloys or electron–hole plasmas. Advanced many-body extensions combine HEG benchmarks with techniques such as dynamical mean field theory (DMFT) or GW+DMFT to treat strong correlations, while time-dependent DFT (TDDFT) builds on HEG response functions for dynamical phenomena.