Point Defects in Metals — Vacancies, Interstitials and Substitutional Atoms

A perfect crystal — every atom occupying its designated lattice site with no impurity, no missing atom, no misplaced atom — exists only at absolute zero and only in thermodynamic theory. In every real metal at every temperature above 0 K, thermodynamic equilibrium demands the presence of point defects: zero-dimensional disruptions to the periodic lattice involving a single atomic site or a small cluster of sites. Far from being mere imperfections, point defects control diffusion rates, solid-solution strengthening, precipitation kinetics, radiation damage, creep behaviour, and electrical resistivity — making their quantitative understanding indispensable across every domain of physical metallurgy and materials engineering.

1. Classification of Point Defects

Point defects are classified by the nature of the departure from perfect crystal periodicity at a single lattice site or its immediate neighbourhood:

1.1 Vacancies

A vacancy is a lattice site from which the host atom is absent. It is the simplest and most important point defect in pure metals. When an atom leaves its lattice site to migrate to the crystal surface, it creates a Schottky defect — the standard vacancy encountered in thermally equilibrated metals. The surrounding atoms relax inward by approximately 1–5% of the nearest-neighbour distance, creating a localised compressive zone around the vacancy. Vacancies are the primary vehicle for solid-state diffusion and therefore govern every diffusion-limited process: precipitation, carburising, sintering, and creep.

1.2 Self-Interstitials and Frenkel Defects

A self-interstitial is a host atom displaced from its lattice site into the interstice between regular atomic positions, producing a large compressive strain field extending over several atomic diameters. The simultaneous creation of a vacancy and a self-interstitial is a Frenkel defect. Self-interstitials have very high formation energies in close-packed metals (typically 2–4 eV compared to 0.5–1.5 eV for vacancies), so thermal Frenkel defects are negligible at equilibrium. They become the dominant defect type under radiation damage, where energetic neutron or ion collisions knock atoms off lattice sites into interstitial positions, creating vacancy-interstitial Frenkel pairs in displacement cascades.

1.3 Substitutional Solute Atoms

When a foreign (solute) atom of similar size replaces a host atom on a lattice site, it forms a substitutional solid solution. The size mismatch between solute and solvent creates a spherically symmetric strain field — compressive if the solute is larger than the host, tensile if smaller. This strain field is the principal mechanism of substitutional solid-solution strengthening. Examples include: Ni substituting for Fe in austenitic stainless steels; Cu substituting for Al in age-hardenable 2xxx alloys; Cr, Mo, and Mn in carbon steel matrices.

1.4 Interstitial Solute Atoms

Small atoms (C, N, H, B, O) with atomic radii less than approximately 60% of the host atom radius can occupy interstitial positions — the gaps between regular lattice sites — forming an interstitial solid solution. In FCC austenite, carbon occupies the large octahedral hole at body-centre and edge-centre positions. In BCC ferrite, carbon is forced into the much smaller tetrahedral hole, creating a strongly asymmetric (tetragonal) distortion rather than a spherical one. This tetragonal distortion interacts far more strongly with dislocations than a spherical substitutional strain field, accounting for the exceptional strengthening per unit carbon in martensite, where quenching traps carbon in supersaturation in the BCT (body-centred tetragonal) lattice.

2. Thermodynamics of Vacancy Formation

The equilibrium vacancy concentration is not zero at any temperature above absolute zero. The thermodynamic argument is precise: introducing n_v vacancies into a crystal of N sites increases the internal energy by n_vH_f (H_f = vacancy formation enthalpy), but simultaneously increases the configurational entropy by k_B ln(W), where W = N!/(N−n_v)!n_v! is the number of ways of distributing the vacancies over available sites. The Gibbs free energy of the crystal is minimised at a specific equilibrium vacancy fraction:

Equilibrium vacancy mole fraction:

  n_v / N = exp(−Q_v / RT)           [simplified, ignoring vibrational entropy]

More precisely:
  n_v / N = exp(ΔS_v / R) × exp(−ΔH_v / RT)

where:
  n_v    = number of vacancies
  N      = total number of lattice sites
  Q_v    = ΔH_v = vacancy formation enthalpy (J/mol or eV/atom)
  ΔS_v   = vacancy formation entropy (~1–3 k_B — usually small)
  R      = 8.314 J/(mol·K)
  T      = absolute temperature (K)

Typical ΔH_v values:
  Aluminium   :  0.68 eV (65.6 kJ/mol)
  Copper      :  1.07 eV (103  kJ/mol)
  Iron (BCC)  :  1.60 eV (154  kJ/mol)
  Nickel      :  1.74 eV (168  kJ/mol)
  Gold        :  0.94 eV (90.7 kJ/mol)

Example — copper at 1000°C (1273 K):
  n_v/N = exp(−103,000 / (8.314 × 1273)) ≈ 1 × 10⁻⁴
  → approximately 1 vacancy per 10,000 atoms near the melting point

At 20°C (293 K):
  n_v/N = exp(−103,000 / (8.314 × 293)) ≈ 10⁻¹⁸
  → vacancies negligible under equilibrium at room temperature

The practical consequence of this strong temperature dependence is quench-retained vacancies: rapid cooling from near the melting point freezes in the high-temperature vacancy population before it can equilibrate to the negligible room-temperature value. In aluminium alloys (2xxx, 6xxx, 7xxx), quenching from solution treatment temperature (~500–535°C) retains ~10⁻⁵ excess vacancies per site, which greatly accelerate solute diffusion during subsequent ageing — by reducing the jump activation energy for vacancy-solute atom exchanges. This is why precipitation hardening in Al alloys proceeds far faster after a water quench than after air cooling. The link between vacancy concentration, diffusivity, and precipitation kinetics is covered in more detail in the annealing and normalising article.

2.1 Vacancy Migration and Diffusion

A vacancy moves through the crystal when an adjacent atom jumps into the vacant site — the vacancy mechanism of diffusion. The atom requires enough thermal energy to surmount the saddle-point configuration between its initial site and the vacant site. The migration enthalpy Q_m is typically 0.5–1.0 eV, comparable in magnitude to the formation enthalpy Q_f. The overall activation energy for self-diffusion Q_D is therefore the sum of both:

Self-diffusion coefficient:
  D = D₀ × exp(−Q_D / RT)

  where  Q_D = Q_f + Q_m   (formation + migration enthalpies)

For substitutional diffusion (vacancy mechanism):
  Q_D = Q_f + Q_m + binding enthalpy (solute–vacancy interaction)

Example — carbon diffusion in FCC austenite (interstitial mechanism, no vacancy needed):
  D₀ = 2.3 × 10⁻⁵ m²/s
  Q_D = 148 kJ/mol
  D at 1000°C ≈ 2.5 × 10⁻¹¹ m²/s

Carbon diffusion in BCC ferrite at 700°C:
  D ≈ 3 × 10⁻¹³ m²/s
  → ~100× faster in austenite than ferrite at the same temperature
     (despite lower T): larger holes reduce migration barrier

3. Interstitial Solid Solutions: The C-Fe Case

The iron-carbon system provides the most consequential example of interstitial solid solution behaviour in engineering metals. The dramatic difference in carbon solubility between FCC austenite and BCC ferrite — 2.14 wt% vs. 0.022 wt% at their respective eutectic/eutectoid temperatures — is entirely explained by the geometry of interstitial holes in the two crystal structures.

3.1 Octahedral vs. Tetrahedral Holes

Every close-packed (FCC) crystal has two types of interstitial site: octahedral holes (at body-centre and edge-centre positions, coordinated by 6 host atoms) and tetrahedral holes (coordinated by 4 host atoms). BCC crystals have tetrahedral and octahedral holes, but the geometry is different — the BCC octahedral hole is smaller than the BCC tetrahedral hole in terms of the effective radius it can accommodate, and both are smaller than the FCC octahedral hole.

Carbon atoms (r_C = 0.077 nm) are slightly too large even for FCC octahedral holes (ideal r = 0.052 nm for close-packed iron), so they always produce some lattice expansion. However, the mismatch in FCC is far less severe than in BCC, allowing far greater solubility before the elastic strain energy destabilises the solid solution. In BCC ferrite, the interstitial carbon produces a tetragonal distortion — the two nearest-neighbour iron atoms along the ⟨100⟩ direction in which the carbon sits are pushed further apart than the four atoms in the perpendicular plane. This asymmetric distortion interacts strongly with dislocation stress fields, providing a large strengthening increment and underpinning the Cottrell atmosphere locking that produces the upper and lower yield point in low-carbon steel. For the consequences of this on the Fe-C phase diagram, see the iron-carbon phase diagram article and the eutectoid reaction article.

4. Substitutional Solid Solutions and the Hume-Rothery Rules

William Hume-Rothery established four empirical rules that predict whether a solute metal will have extensive substitutional solid solubility (>~10 at%) in a given solvent metal. All four must be satisfied simultaneously for a complete or extensive solid solution:

4.1 The Four Rules

Hume-Rothery rules are necessary but not sufficient conditions — they identify candidates for solid solution but do not predict exact solubility limits. First-principles calculations (DFT) and CALPHAD thermodynamic databases are used for quantitative prediction in alloy design. The rules are most useful as rapid screening tools: if any rule is clearly violated, extensive solid solubility can be ruled out. The Cu-Ni system (all four rules satisfied) exhibits complete FCC solid solubility across the entire composition range; the Cu-Ag system (15.2% size difference, same structure, similar electronegativity) shows only terminal solid solutions due to the size factor violation.

4.2 Solid Solution Strengthening Mechanisms

Every solute atom introduces a strain field that interacts with moving dislocations. The two principal interaction types are size misfit (parelastic) interaction — due to the volume misfit between solute and host — and modulus interaction (dielastic) — due to the difference in elastic modulus between solute and host regions. For substitutional solutes, the size misfit dominates. The increase in critical resolved shear stress (CRSS) scales approximately as:

Fleischer model (parelastic + dielastic):
  Δτ = A × G × ε_s^(3/2) × c^(1/2)

where:
  A   = numerical constant (~1/700 for FCC)
  G   = shear modulus of matrix
  c   = solute mole fraction
  ε_s = effective misfit parameter = |ε_a − (η/16)|

  ε_a = atomic size misfit = (1/a)(da/dc)      [lattice parameter change per unit c]
  η   = modulus misfit     = (1/G)(dG/dc)

Interstitial solutes (C, N in BCC iron) — Cochardt model:
  Δτ = const × G × ε_t × c
  (stronger: linear in c, not c^(1/2), due to tetragonal strain → stronger dislocation interaction)

Practical order of strengthening in steel (per wt%):
  C (interstitial)  : ~5000 MPa/wt% in martensite
  N (interstitial)  : ~4500 MPa/wt% in solution
  Si                :   ~83 MPa/wt%
  Mn                :   ~32 MPa/wt%
  Mo                :   ~11 MPa/wt%
  Cr                :    ~8 MPa/wt%

The enormous strengthening coefficient of interstitial carbon in martensite — approximately 5000 MPa per wt%C — relative to substitutional elements arises from the tetragonal distortion of the BCT lattice, which couples strongly to both edge and screw dislocation stress fields. At 0.4 wt% C, the CRSS increment alone accounts for approximately 2000 MPa of the martensite hardness, explaining why untempered martensite in medium-carbon steel is harder than any substitutionally strengthened alloy of comparable composition. The relationship between carbon content, martensite tetragonality (c/a ratio), and hardness is examined in depth in the martensite formation article.

5. Radiation Damage and Point Defects

Nuclear reactor environments produce point defects at rates orders of magnitude above thermal equilibrium through displacement cascades. When a fast neutron (energy ~1 MeV) collides with a lattice atom, it imparts kinetic energy sufficient to displace it from its site — the primary knock-on atom (PKA). The PKA then collides with its neighbours, creating a cascade that displaces thousands of atoms within a ~10 nm volume in picoseconds. Most displaced atoms recombine with vacancies within the cooling cascade, but a fraction (~1%) survive as isolated Frenkel pairs that migrate through the lattice and produce measurable property changes.

Displacement damage quantification:
  Displacements per atom (dpa) = Φ × σ_d × t

  where:
    Φ   = neutron flux (neutrons/m²s)
    σ_d  = displacement cross-section (m²)
    t    = irradiation time (s)

  Typical PWR reactor pressure vessel:
    Φ ~ 10¹⁵ n/cm²/s × 40 years ~ 10²⁴ n/cm² total fluence
    Δdpa ~ 0.01–0.1 dpa at vessel wall

Consequences in reactor steels:
  Radiation hardening: Δσ_y ~ A × √dpa   (dislocation loop obstacle strengthening)
  DBTT shift: ΔT_DBTT ∝ √(fluence)   (Charpy 41 J transition temperature rises)
  Swelling: voids from vacancy cluster coalescence → density decrease at T > 0.4 T_m

Radiation-induced segregation (RIS) is an additional consequence of point defect fluxes in irradiated alloys. Solute atoms that bind preferentially to vacancies are carried with the vacancy flux toward sinks (grain boundaries, dislocations), enriching boundaries in elements such as Ni, Si, and P in austenitic stainless steels. This boundary enrichment can trigger sensitisation and intergranular stress corrosion cracking — a failure mode of significance in boiling water reactor (BWR) internals and light water reactor coolant piping. For more on grain boundary effects, see the grain boundaries article.

6. Point Defects in Engineering Practice

6.1 Age Hardening in Aluminium Alloys

The precipitation hardening sequence in 6xxx Al alloys (Al-Mg-Si) proceeds through a series of metastable phases — GP zones, β″″ needles, β″ rods, equilibrium β (Mg₂Si) — whose nucleation and growth kinetics are governed by vacancy-assisted solute diffusion. Solution treatment at ~540°C dissolves all precipitates into a supersaturated solid solution and simultaneously creates a high vacancy concentration. Water quenching retains both solute supersaturation and excess vacancies. During natural ageing (room temperature) or artificial ageing (160–200°C), vacancies enhance Mg and Si diffusivity by orders of magnitude relative to what would occur at those temperatures in the absence of excess vacancies, enabling GP zone formation and β″″ precipitation within hours rather than years. Alloy designers exploit this by specifying quench rates, pre-ageing treatments, and ageing temperatures to optimise the precipitate size distribution for peak hardness or fracture toughness. See the connection to quenching concepts for analogous treatment in steel systems.

6.2 Carburising and Nitriding

Thermochemical surface hardening by carburising and nitriding exploits the interstitial solid solution capacity of austenite. In gas carburising, the steel component is held at 900–950°C in a carbon-rich atmosphere; carbon dissolves interstitially into the surface austenite and diffuses inward down the concentration gradient according to Fick’s second law. Case depth x is approximately:

Case depth estimation (semi-infinite solid, constant surface concentration):

  C(x,t) = C_s − (C_s − C_0) × erf(x / (2√(Dt)))

  where:
    C_s = surface carbon concentration (wt%)
    C_0 = initial bulk carbon (wt%)
    D   = carbon diffusivity in austenite at process T (m²/s)
    t   = carburising time (s)
    x   = depth from surface (m)

  Approximate case depth (to C = 0.3 wt%):
    x ≈ 2 √(Dt)

  At 925°C: D(C in austenite) ≈ 1.0 × 10⁻¹¹ m²/s
  For 4h (14400 s): x ≈ 2√(1.0×10⁻¹¹ × 14400) ≈ 0.76 mm

After carburising, quenching transforms the carbon-enriched surface layer to martensite (high hardness, ~700–900 HV) while the lower-carbon core transforms to bainite or tempered martensite (tougher). The entire strategy depends on the interstitial solubility of carbon in austenite — if it were as limited as in ferrite, carburising would be thermochemically impossible. For the microstructural outcome at the carburised surface and its relation to bainite formation, see the bainite microstructure article.

6.3 Hydrogen Embrittlement: Interstitial Hydrogen

Hydrogen (r_H = 0.053 nm) is the smallest interstitial atom and dissolves readily in both FCC and BCC metals. In steels, atomic hydrogen at the surface (from corrosion, cathodic protection, electroplating, or welding) enters as H and diffuses interstitially at high rates due to its small size (D_H in α-Fe at 25°C ≈ 10⁻⁹ m²/s — far exceeding carbon or nitrogen diffusivities). Hydrogen accumulates at stress concentrations, grain boundaries, inclusions, and dislocation cores — the trap sites that reduce hydrogen chemical potential. Above a critical local hydrogen concentration, hydrogen embrittlement mechanisms operate: hydrogen-enhanced decohesion (HEDE) at grain boundaries and hydrogen-enhanced localised plasticity (HELP) by dislocation mobility enhancement. The result is brittle fracture at stresses well below the yield strength in sustained-load conditions. Hydrogen damage is a primary concern in high-strength steels (σ_y > 1000 MPa), sour service pipelines (H₂S environments), and welded joints — and is examined in detail in the hydrogen-induced cracking article.

7. Comparative Summary of Point Defect Types

Recommended Reference Books

Disclosure: MetallurgyZone participates in the Amazon Associates programme. If you purchase through these links, we may earn a small commission at no extra cost to you. This helps support free technical content on this site.