Point Defects, Defect Thermodynamics & Non-Stoichiometry

At any finite temperature ($T > 0\text{ K}$), an ideal, defect-free crystal lattice is thermodynamically unstable. The spontaneous formation of point defects increases the internal enthalpy ($\Delta H_d > 0$), but this energetic penalty is overcome by the enormous increase in configurational entropy ($\Delta S_{\text{conf}} > 0$).

Free Energy Formulation

Consider creating $n$ isolated vacancies in a crystal containing $N$ regular atomic lattice sites. The Gibbs free energy change of defect formation is:

where $\Delta H_v$ is the enthalpy of vacancy formation and $\Delta S_{\text{vib}}$ is the local vibrational entropy change due to altered phonon frequencies around the vacancy. The number of distinct ways $W$ to distribute $n$ vacancies among $N+n$ available sites is:

Applying Boltzmann's entropy formula $S = k_B \ln W$ and Stirling's approximation ($\ln x! \approx x\ln x - x$):

Minimization and Equilibrium Concentration

Thermodynamic equilibrium at constant temperature and pressure requires:

Differentiating $\Delta G$:

Since $n \ll N$, $N+n \approx N$:

Because $\left( \frac{\partial^2 \Delta G}{\partial n^2} \right) = \frac{k_B T}{n} > 0$, this represents a true global minimum. Point defects are equilibrium thermodynamic entities, unlike dislocations or grain boundaries.

In ionic crystals, point defects must preserve macroscopic electrical neutrality. Walter Schottky (1935) identified the stoichiometric stoichiometric defect pair known as the Schottky defect.

Definition and Mechanism

A Schottky defect consists of a stoichiometric pair of vacancies: one cation vacancy and one anion vacancy (for a $1:1$ binary crystal $MX$).

• To maintain crystal density and surface integrity, the displaced cation and anion migrate to the exterior crystal surface, creating new regular lattice sites:

• Both sublattices lose ions in equal proportions, preserving charge neutrality and stoichiometric composition.
• Crystal density decreases because the crystal volume increases while mass remains constant.

Equilibrium Mass-Action Law

Let $N$ be the number of cation sites (equal to anion sites in $MX$). If $n_S$ Schottky pairs are formed, the number of ways to arrange $n_S$ cation vacancies on $N$ cation sites is $W_c = \frac{N!}{(N-n_S)! n_S!}$, and similarly for anion vacancies $W_a = \frac{N!}{(N-n_S)! n_S!}$. Total configurational entropy:

Minimizing $\Delta G = n_S \Delta H_S - T \Delta S_{\text{conf}}$:

Assuming $n_S \ll N$:

Notice the crucial factor of $1/2$ in the exponent: because the defect is a binary pair, the effective activation energy is half the total pair formation enthalpy $\Delta H_S/2$.

• Typical formation enthalpies in alkali halides are $\Delta H_S \approx 2.0 - 2.5\text{ eV}$, yielding $n_S/N \sim 10^{-5}$ near the melting point and $n_S/N \sim 10^{-16}$ at room temperature.

Yakov Frenkel (1926) formulated an alternative stoichiometric defect mechanism: the Frenkel defect.

Microscopic Mechanism

A Frenkel defect consists of an ion that leaves its normal crystallographic lattice site and squeezes into an unoccupied interstitial void, creating a vacancy-interstitial pair:

• The stoichiometry and total mass remain unaltered.
• The overall macroscopic crystal volume remains virtually unchanged, so crystal density does not decrease appreciably.
• Frenkel defects typically occur on the cation sublattice (cation Frenkel defect) because cations are substantially smaller than anions ($r^+ \ll r^-$) and fit much more readily into tight interstitial cavities.
• Anion Frenkel defects ($\text{V}_X^\bullet + X_i'$) are rare due to severe steric overlap, but occur in fluorite-type crystals ($\text{CaF}_2, \text{UO}_2$) where empty interstitial cubes accommodate fluoride or oxide ions.

Thermodynamic Formulation

Let $N$ be the number of normal lattice sites and $N_i$ be the number of available interstitial sites. The number of ways to distribute $n_F$ vacancies among $N$ normal sites is $\frac{N!}{(N-n_F)! n_F!}$, and $n_F$ interstitials among $N_i$ interstitial sites is $\frac{N_i!}{(N_i-n_F)! n_F!}$. Minimizing $\Delta G$:

Again, the exponent contains the factor $1/2$.

• Frenkel Archetype: Silver halides ($\text{AgCl}, \text{AgBr}$). The polarizable $4d^{10}$ shell of $\text{Ag}^+$ enables easy deformation and stabilization in tetrahedral interstitial sites of the FCC bromide lattice ($\Delta H_F \approx 1.2\text{ eV}$), providing the photographic and fast ionic conduction mechanism of silver salts.

F. A. Kröger and H. J. Vink (1956) developed a systematic chemical notation for writing point defect equilibria in crystalline solids.

The Kröger-Vink Notation Syntax: $A_S^C$

A defect species is represented by three components:

1. Main Symbol ($A$): Indicates the entity currently occupying the site.

• Element symbol (e.g., $\text{Na}, \text{Cl}, \text{Fe}$).
• $\text{V}$ represents a vacancy (empty site).
• $e$ or $e'$ represents a free conduction electron.
• $h$ or $h^\bullet$ represents an electron hole in the valence band.

2. Subscript ($S$): Indicates the crystallographic site type.

• Element symbol whose site is occupied (e.g., $M, X, \text{Na}, \text{Cl}$).
• $i$ represents an interstitial site.

3. Superscript ($C$): Indicates the effective (virtual) charge relative to the perfect, unperturbed crystal lattice:

• $\times$ represents neutral effective charge ($q_{\text{eff}} = 0$).
• $\bullet$ represents positive effective charge ($q_{\text{eff}} = +1$).
• $'$ (prime) represents negative effective charge ($q_{\text{eff}} = -1$).

Rules for Balancing Defect Equations

1. Conservation of Mass: Chemical atoms must balance identically on both sides (vacancies $\text{V}$ have zero mass).

2. Conservation of Effective Charge: The algebraic sum of effective charges on the left must equal the sum on the right:

3. Conservation of Site Ratio: The ratio of cation sites to anion sites in the crystal host must remain constant according to the host stoichiometry. Sites cannot be created or destroyed independently on one sublattice alone.

Examples of Standard Defect Reactions

• Schottky defect in $\text{NaCl}$:

• Frenkel defect in $\text{AgBr}$:

• Electron-hole thermal excitation across the band gap:

Many transition metal oxides, sulfides, and halides deviate measurably from ideal integer stoichiometry depending on temperature and ambient oxygen partial pressure ($P_{\text{O}_2}$). These are classified into two broad regimes:

1. Metal-Deficient Oxides ($M_{1-x}O$, e.g., $\text{Fe}_{1-x}\text{O}, \text{Ni}_{1-x}\text{O}, \text{CoO}$)

• Structural Nature: Cation vacancies are created when the crystal is exposed to oxidizing atmospheres (high $P_{\text{O}_2}$).
• Defect Reaction: Oxygen gas incorporates into the lattice, creating regular oxygen sites and an equivalent number of cation vacancies:

• To preserve charge neutrality, two metal cations are oxidized to higher valence states (e.g., $2\text{Fe}^{2+} \to 2\text{Fe}^{3+}$, represented as two electron holes $2h^\bullet$).
• By the law of mass action:

Under electroneutrality $[h^\bullet] = 2[\text{V}_M'']$:

The vacancy concentration and hole electrical conductivity ($\sigma = p e \mu_h$) scale as $P_{\text{O}_2}^{1/6}$ ($p$-type semiconductor).

2. Anion-Deficient Oxides ($MO_{2-x}$ or $M_{1+x}O$, e.g., $\text{ZnO}, \text{TiO}_{2-x}, \text{CeO}_{2-x}$)

• Structural Nature: Oxygen leaves the lattice into the gas phase under reducing conditions (low $P_{\text{O}_2}$), leaving behind oxygen vacancies and conduction electrons:

• Mass action relation:

Under electroneutrality $[e'] = 2[\text{V}_O^{\bullet\bullet}]$:

Conductivity scales inversely with oxygen pressure as $P_{\text{O}_2}^{-1/6}$ ($n$-type semiconductor).

When an alkali halide crystal (e.g., colorless $\text{NaCl}$) is heated in an excess alkali metal vapor (e.g., sodium vapor) and rapidly quenched, it develops intense coloration (yellow for $\text{NaCl}$, violet for $\text{KCl}$, blue for $\text{KBr}$). This is caused by the formation of an F-center (Farbzentrum, German for color center).

Microscopic Architecture of an F-Center

1. Gaseous sodium atoms deposit on the crystal surface:

1. The sodium atom ionizes to $\text{Na}^+$, incorporating into the cation lattice and creating an empty chloride vacancy ($\text{V}_{\text{Cl}}^\bullet$).
2. The liberated electron diffuses rapidly into the crystal bulk and is captured by the electrostatic potential well of the positive anion vacancy:

An F-center is an electron trapped in an anion vacancy, surrounded symmetrically by six regular nearest-neighbor cations ($\text{Na}^+$).

Particle-in-a-Spherical-Box Model & The Mollwo-Ivey Law

The trapped electron behaves quantum mechanically as a particle in a three-dimensional potential well of diameter approximately equal to the lattice parameter $a$. The ground state is an $s$-like orbital ($1s$), and optical absorption promotes the electron to an excited $p$-like state ($2p$):

Empirically, E. Mollwo (1931) and H. Ivey (1947) demonstrated that the optical absorption maximum wavelength $\lambda_{\max}$ across all alkali halides follows the power law:

where $a$ is the cubic lattice parameter and the empirical exponent $n \approx 1.8 - 2.0$. For example, $\text{NaCl}$ ($a = 5.64\text{ Å}$) absorbs at $\lambda = 465\text{ nm}$ (blue light absorption, reflecting complementary yellow).

A Brouwer diagram (or Kröger-Vink diagram) is a double-logarithmic plot displaying the concentrations of all point defects ($\log [D_j]$) as a continuous function of thermodynamic activity (most commonly $\log P_{\text{O}_2}$) at constant temperature.

Brouwer Approximation Method

The full electroneutrality equation for a metal oxide with vacancies, electrons, and holes is:

Because solving this high-order polynomial analytically is cumbersome, G. Brouwer (1954) proposed dividing the diagram into distinct pressure regimes where only the two largest dominant defect species are balanced:

1. Regime I: High Oxygen Partial Pressure (Oxidizing / p-type):

• Dominant neutrality condition: $2[\text{V}_M''] \approx [h^\bullet]$.
• Defect slopes:

• Minority defects: $[e'] \propto P_{\text{O}_2}^{-1/6}$, $[\text{V}_O^{\bullet\bullet}] \propto P_{\text{O}_2}^{-1/3}$.

2. Regime II: Intermediate Oxygen Partial Pressure (Stoichiometric / Intrinsic):

• Thermal electronic excitation dominates: $[e'] \approx [h^\bullet] = \sqrt{K_i} = \text{constant}$ (independent of $P_{\text{O}_2}$).
• Or Schottky equilibrium dominates: $[\text{V}_M''] \approx [\text{V}_O^{\bullet\bullet}] = \sqrt{K_S} = \text{constant}$.

3. Regime III: Low Oxygen Partial Pressure (Reducing / n-type):

• Dominant neutrality condition: $[e'] \approx 2[\text{V}_O^{\bullet\bullet}]$.
• Defect slopes:

• Minority defects: $[h^\bullet] \propto P_{\text{O}_2}^{+1/6}$, $[\text{V}_M''] \propto P_{\text{O}_2}^{-1/3}$.

In each regime, defect concentrations plot as straight lines with rational fractional slopes ($0, \pm 1/6, \pm 1/4, \pm 1/2$).

When the departure from stoichiometry becomes large ($x > 0.01$), isolated point defects interact electrostatically and elastically, condensing into complex defect clusters and crystallographic shear structures.

Defect Clusters in Wüstite (Fe_{1-x}O)

In $\text{Fe}_{1-x}\text{O}$ ($0.05 \le x \le 0.15$), cation vacancies do not remain randomly dispersed.

• Roth Clusters (4:1 Cluster): Four octahedral $\text{Fe}^{2+}$ vacancies group together around a central tetrahedral interstitial site occupied by an oxidized $\text{Fe}^{3+}$ cation.
• Koch-Cohen Clusters (13:4 Cluster): Aggregations of 13 octahedral vacancies and 4 tetrahedral $\text{Fe}^{3+}$ ions surrounded by relaxed neighboring ions, forming a localized inverse-spinel sub-microdomain (pre-nucleation of magnetite $\text{Fe}_3\text{O}_4$).

Crystallographic Shear (CS) & Magnéli Phases

In early transition metal oxides with high formal oxidation states ($\text{WO}_3, \text{TiO}_2, \text{MoO}_3$), large non-stoichiometry is accommodated without creating point vacancies. Instead, the crystal eliminates entire planes of oxygen atoms through crystallographic shear (CS):

1. Corner-sharing octahedra collapse along a specific crystallographic plane (e.g., $(121)$ or $(132)$ in rutile).
2. The collapsed region slips, transforming corner-sharing octahedra into denser edge-sharing or face-sharing octahedra.
3. Periodic recurrence of these shear planes generates discrete homologous series of fully ordered, stoichiometric compounds known as Magnéli phases:

• Titanium oxides: $\text{Ti}_n\text{O}_{2n-1}$ ($4 \le n \le 10$, e.g., $\text{Ti}_4\text{O}_7$).
• Tungsten oxides: $\text{W}_n\text{O}_{3n-1}$ and $\text{W}_n\text{O}_{3n-2}$ (e.g., $\text{W}_{20}\text{O}_{58}$).

Magnéli phases exhibit high electrical conductivities and chemical inertness, functioning as corrosion-resistant electrocatalyst supports.