Defect Concentration Calculator

A perfect crystal is never the equilibrium state above absolute zero. Creating a vacancy costs the formation enthalpy Ef, but scattering n vacancies over N sites adds configurational entropy, and minimizing the Gibbs energy G = nEf minus TS with the Stirling approximation yields the simple Boltzmann result n/N = exp(-Ef/kBT). With kB equal to 8.617e-5 eV/K, the exponent is just Ef divided by 0.02585 eV at room temperature, which is why formation energies of even 1 eV give vanishing populations at 300 K yet substantial ones near melting. Measured monovacancy formation energies run from about 0.67 eV in aluminium and 1.28 eV in copper to roughly 3.6 eV in tungsten, and the equilibrium site fraction at the melting point clusters in the narrow range of 1e-4 to 1e-3 for most metals. Ionic crystals conserve charge, so defects appear in pairs and the exponent is halved. Schottky pairs, a cation and anion vacancy together, dominate in alkali halides such as NaCl where the pair formation energy is near 2.3 eV, while Frenkel pairs, a vacancy plus a displaced interstitial, dominate in silver halides where the ion size mismatch is large. Several limitations deserve attention. The formula omits the vibrational entropy prefactor exp(Sf/kB), typically a factor of 1 to 10. Below a crossover temperature the population is extrinsic, fixed by aliovalent impurities rather than thermodynamics. Charged defects in semiconductors depend on Fermi level and chemical potential, requiring Kroger-Vink bookkeeping. Quenched or irradiated samples are far from equilibrium entirely.

Yakov Frenkel proposed in 1926 that atoms in a crystal could leave their normal sites for interstitial positions, introducing the idea of thermally generated lattice disorder and giving a mechanism for ionic conduction in solids that were otherwise thought rigid. Walter Schottky and Carl Wagner developed the statistical thermodynamics of this disorder around 1930 and 1931, treating defect formation as a chemical equilibrium and deriving the mass-action relations that explain non-stoichiometry and ionic transport in compounds such as silver chloride and zinc oxide. Their framework showed why the conductivity of an ionic crystal follows an Arrhenius law with an activation energy combining formation and migration terms, and it earned Wagner recognition as a founder of solid state chemistry. Direct measurement came much later. Robert Simmons and Robert Balluffi published experiments around 1960 to 1963 comparing macroscopic thermal expansion with the X-ray lattice parameter expansion of aluminium, copper, silver and gold, the difference giving the vacancy concentration without any model assumption, and these remain benchmark values. Positron annihilation spectroscopy was developed from the 1960s onward as a sensitive vacancy probe, since positrons localize in the reduced electron density of an open volume. Notation was standardized by Ferdinand Kroger and Hans Vink in 1956, whose symbols for site, species and effective charge are now universal. Since the 1990s density functional theory supercell calculations have supplied formation energies for defects too dilute or too transient to measure.