Every solid is a thermodynamic compromise between the perfect periodic crystal and the maximum-entropy random structure. Real materials are filled with defects — vacancies, dislocations, grain boundaries — that are not imperfections to be eliminated but levers engineers pull to tune hardness, conductivity, and magnetism. Meanwhile, shrinking a material to nanometre dimensions unlocks quantum phenomena that are invisible in bulk: size-tunable optical absorption in quantum dots, chirality-dependent conductivity in carbon nanotubes, and zero-gap Dirac dispersion in graphene.
1. Crystal Defects and Dislocations
A perfect crystal sits at a free-energy minimum only at 0 K; finite temperature guarantees the presence of point defects. Line defects (dislocations) are the primary mechanism by which metals deform plastically: a dislocation allows atoms to slip one at a time rather than requiring an entire plane to move simultaneously, reducing the theoretical shear stress by a factor of ~1000.
Point Defects, Dislocations & Burgers Vector
Point defects:
Vacancy: missing atom; E_f ≈ 0.5–3 eV (Cu: 1.28 eV)
Interstitial: extra atom in non-lattice site (self- or foreign)
Frenkel pair: vacancy + interstitial on same sublattice (ionic crystals)
Schottky pair: cation + anion vacancy pair (electrical neutrality)
Equilibrium vacancy fraction: n_v/N = exp(−E_f / k_B T)
Dislocation types:
Edge: extra half-plane; Burgers vector b ⊥ dislocation line
Screw: helical path around line; b ∥ dislocation line
Mixed: combination of both characters (most real dislocations)
Burgers vector b:
Defined by Burgers circuit enclosing the dislocation core
Closure failure = b (perfect: b = lattice vector)
|b| ≈ a/√2 ≈ 0.25 nm (FCC metals, close-packed ⟨110⟩ direction)
Slip systems (FCC metals: Al, Cu, Ni, Au):
Slip planes: {111} (close-packed planes)
Slip directions: ⟨110⟩ (close-packed rows) → 12 slip systems
Peierls-Nabarro stress: τ_PN ~ G exp(−2πd / b(1−ν))
d = interplanar spacing, ν = Poisson ratio, G = shear modulus
Explains why close-packed planes slip preferentially
Dislocation density ρ (m−²):
Annealed metal: ρ ≈ 10¹&sup0; m−²
Heavily deformed: ρ ≈ 10¹&sup5; m−²
2. Strengthening Mechanisms
All strengthening mechanisms work by impeding dislocation motion. The yield stress contributions from different mechanisms are roughly additive (the square-root addition rule applies when multiple mechanisms are active simultaneously).
Hall-Petch, Taylor Hardening & Other Mechanisms
Work hardening (Taylor, 1934):
σ = σ_0 + M α G b √ρ (Taylor hardening)
M ≈ 3.06 (Taylor factor for FCC polycrystal)
α ≈ 0.2–0.4 (interaction constant)
Strain hardening exponent n in Hollomon: σ = K ε^n
Considere criterion for necking (n = ε_u):
dσ/dε = σ → uniform elongation ε_u = n
Grain boundary strengthening (Hall-Petch, 1951):
σ_y = σ_0 + k_y / √d
d = grain diameter, k_y ≈ 0.7 MPa·m^(1/2) for mild steel
Nanocrystalline metals (d < 20 nm): inverse Hall-Petch possible
Solid solution strengthening:
Δτ ∝ c^(1/2) · ε^(3/2) (Fleischer model)
c = solute fraction, ε = misfit strain from size difference
Precipitation hardening (Orowan, 1948):
Orowan bypass: Δτ_or ≈ G b / (L − 2r) where L = particle spacing, r = radius
Particle shearing (coherent precipitates): Δτ_sh = F_max / (b L)
Critical particle radius r* where bypass becomes easier than shearing
3. Binary Phase Diagrams
A binary phase diagram maps the thermodynamically stable phases as a function of composition and temperature at fixed pressure. It compresses an enormous amount of experimental thermochemical data into a form directly useful for alloy design and heat treatment.
Phase Diagram Features, Lever Rule & Transformations
Gibbs phase rule (for binary alloy at 1 atm):
F = C − P + 1 (C=2 components, 1 pressure fixed)
F = 3 − P → max 2 phases in two-phase region (tie-line), 1 at eutectic point
Lever rule (composition in two-phase α+L region):
x_0 = alloy composition, x_α = α solvus, x_L = liquidus
Fraction of α: f_α = (x_L − x_0) / (x_L − x_α)
Fraction of liquid: f_L = (x_0 − x_α) / (x_L − x_α)
Eutectic reaction (e.g. Pb-Sn at 61.9 wt% Sn, 183°C):
L → α + β (three-phase invariant point: F=0)
Eutectic microstructure: lamellar α-β spacing λ ≈ 1–10 μm
Peritectic reaction:
L + α → β (less common; Fe-C has peritectic at 0.09 wt% C)
Solid-state transformations (TTT diagram of steel):
Pearlite: γ → α + Fe_3C (diffusional; nose at ~550°C)
Bainite: faster quench; finer carbide dispersed in ferrite
Martensite: displacive (no diffusion); below M_s ≈ 200–400°C
Body-centred tetragonal; hardness up to 900 HV; very brittle
Martensite fraction (Koistinen-Marburger):
f_M = 1 − exp[−0.011(M_s − T)] for T < M_s
4. Semiconductor Band Theory
The electronic band structure of a crystal follows from the periodicity of the lattice. Bloch’s theorem guarantees that eigenstates in a periodic potential take the form of a plane wave modulated by a lattice-periodic function, and the dispersion relation E(k) across the first Brillouin zone determines whether a material is a metal, semiconductor, or insulator.
Bloch's Theorem, Effective Mass & p-n Junction
Bloch's theorem:
ψ_k(r) = u_k(r) e^(ik·r)
u_k(r) = lattice-periodic function; k = crystal momentum in first BZ
Effective mass m* (curvature of band):
1/m* = (1/ℏ²) d²E/dk²
Light holes, heavy holes in Si/Ge valence band (anisotropic m*)
GaAs conduction band: m*_e = 0.067 m_e (very light → high mobility)
Gap types:
Direct gap (GaAs, InP): valence max and conduction min at same k
Efficient optical transitions → LEDs, lasers
Indirect gap (Si, Ge): phonon required to conserve momentum
E_g(Si) = 1.12 eV at 300 K; E_g(Ge) = 0.66 eV; E_g(GaAs) = 1.42 eV
p-n junction equilibrium:
Built-in potential: V_bi = (k_BT/e) ln(N_a N_d / n_i²)
Depletion width: W = √(2εε_0 V_bi (1/N_a + 1/N_d) / e)
Shockley diode equation:
I = I_0 (e^(eV/k_BT) − 1)
I_0 = A e n_i² (D_n/(N_a L_n) + D_p/(N_d L_p)) (reverse saturation current)
Ideality factor η ≈ 1 for diffusion-limited; η ≈ 2 for recombination-limited
5. Polymer Physics
Polymers are long-chain molecules whose mechanical and thermodynamic properties emerge from statistical conformations of the backbone. The freely jointed chain model and Flory’s mean-field theory together explain why rubber is elastic, why polymers phase-separate, and how the glass transition depends on chain stiffness.
Chain Statistics, Rubber Elasticity & Flory-Huggins
Freely jointed chain (N segments, length b):
End-to-end distance: 〈R²〉 = Nb²
Persistence length l_p: orientation correlation 〈t(0)·t(s)〉 = exp(−s/l_p)
Worm-like chain: 〈R²〉 = 2 l_p L (1 − l_p/L (1 − e^(−L/l_p)))
Rubber (affine network) elasticity (Kuhn, 1946):
Shear modulus: G = n k_B T (n = crosslink density per unit volume)
Stress σ = G(λ − λ−²) (λ = stretch ratio, incompressible)
Glass transition T_g:
WLF equation: log(a_T/a_T_g) = −C_1(T−T_g)/(C_2+(T−T_g))
C_1 ≈ 17.4, C_2 ≈ 51.6 K (universal WLF constants near T_g)
Free volume theory: T_g increases with stiffer backbone, shorter side chains
Flory-Huggins theory (polymer-solvent mixing):
ΔG_mix / (nRT) = φ ln φ + (1−φ)ln(1−φ)/N + χ φ(1−φ)
χ = Flory-Huggins parameter (enthalpy of mixing)
χ_c = (1/2)(1 + 1/√N)² ≈ 1/2 for large N (critical point)
Spinodal decomposition (phase separation without nucleation) when:
d²ΔG_mix/dφ² < 0
6. Nanomaterials and Quantum Confinement
When a material is reduced to nanometre dimensions, quantum confinement shifts energy levels, increases surface-to-volume ratio dramatically, and changes melting points, magnetic properties, and optical response. Nanostructures have enabled quantum-dot LED displays, single-wall carbon nanotube transistors, and graphene-based sensors.
Quantum Dots, CNTs, Graphene & Superparamagnetism
Quantum dot confinement (particle-in-a-sphere):
Brus equation: E_gap = E_bulk + (ℏ²π² / 2R²)(1/m*_e + 1/m*_h) − 1.8e²/(4πεε_0 R)
R = dot radius; blue-shift as R decreases (size-tunable colour)
CdSe dots: 2 nm → blue, 7 nm → red (covers visible spectrum)
Carbon nanotube (SWCNT) chirality:
Chiral vector: C = n a_1 + m a_2 (n,m integers, a = 0.246 nm lattice const.)
Armchair (n=m): metallic (zero band gap)
Zigzag (n,0) and chiral: metallic if (n−m) mod 3 = 0, else semiconducting
Band gap of semiconducting SWCNT: E_g = 0.9 eV / d(nm)
High carrier mobility: μ ≈ 10&sup5; cm²/(V·s)
Graphene:
Honeycomb lattice; two-atom unit cell (A, B sublattice) → two bands
Dispersion near K,K' Dirac points: E(k) = ± ℏ v_F |k|
Fermi velocity v_F ≈ 10&sup6; m/s (≈ c/300)
Zero-gap semiconductor (semimetal); opens gap by doping, strain, substrate
Carrier mobility (suspended): ≈ 2 × 10&sup5; cm²/(V·s) at room temp
Superparamagnetism:
Single-domain particles below critical radius r_SD ≈ 10–30 nm
Thermal energy randomises magnetisation direction: blocking temperature T_B
τ = τ_0 exp(K V / k_B T) (Néel-Arrhenius law)
K = anisotropy constant, V = particle volume
Above T_B: no remanence (superparamagnetic); below T_B: ferromagnetic
Application: MRI contrast agents (iron oxide nanoparticles, d ≈ 10 nm)
CALPHAD method: Calculation of Phase Diagrams (CALPHAD) uses assessed thermodynamic databases such as SGTE and TCFE to compute multi-component (up to 25-element) phase diagrams that cannot be measured exhaustively by experiment. Gibbs free energies are parameterised as Redlich-Kister polynomials and optimised against experimental data. CALPHAD is now standard in turbine alloy, steel, and high-entropy alloy design; commercial codes include Thermo-Calc and FactSage.
Try These Simulations
Crystal Structures
3D visualiser of BCC, FCC, HCP and other Bravais lattices with Miller plane and direction annotations.
Alloy Phase Diagram
Binary eutectic and peritectic systems: drag a composition/temperature point to read lever-rule fractions live.
Crystal Growth
Kinetic Monte Carlo surface growing model with adjustable temperature and deposition rate, step-flow to island growth crossover.
Lennard-Jones Fluid
MD simulation of noble-gas atoms interacting via the 6-12 potential; phase transition, pair correlation function and pressure display.