Materials Science & Engineering —
From Crystal Lattices to Smart Alloys

Materials science sits at the intersection of physics, chemistry, and engineering: it asks why iron rusts, why aluminium is light and strong, why Nitinol springs back after enormous deformation, and how superconductors expel magnetic fields. This spotlight tours the interactive simulations on mysimulator.uk that bring these phenomena to life — from the geometry of crystal packing to the quantum physics of band gaps and Bose-Einstein condensation in Cooper pairs.

Crystal Structures & Packing

Almost all solid metals, ceramics, and many minerals are crystalline: their atoms sit on a periodic lattice. The three most common metallic structures are:

The diamond cubic structure (C, Si, Ge) places atoms on an FCC lattice plus 4 interior atoms at ¼,¼,¼ tetrahedral positions. Packing is only 34.0%, but the strong directional covalent bonds make diamond the hardest natural material. The NaCl (rock salt) structure is a paradigmatic ionic crystal: alternating Na+ and Cl- ions on an FCC lattice, each coordinated by 6 ions of opposite charge.

Packing efficiency = (n × Vatom) / Vunit cell
FCC: n = 4   (8×⅛ corners + 6×½ faces)
BCC: n = 2   (8×⅛ + 1 centre)
SC:  n = 1   (8×⅛ corners only)
🔮

Crystal Structures 3D Viewer renders SC, BCC, FCC, HCP, Diamond, and NaCl unit cells using the painter’s algorithm. Drag to rotate, scroll to zoom. Toggle the 2×2×2 supercell to see periodic boundary conditions in action. Mouseover atoms to see coordination numbers.

Alloy Phase Diagrams

Pure metals are rarely used in engineering. Alloys — mixing two or more metallic elements — offer dramatically improved mechanical properties. The equilibrium behaviour of a binary alloy is captured in a phase diagram: a map of which phases (liquid, solid solution, intermetallic compounds) are stable as a function of temperature and composition.

The most important features of a binary phase diagram are:

Lever rule for a 2-phase region (phases α, L):

Weight fraction liquid: WL = (C0 − Cα) / (CL − Cα)
Weight fraction solid: Wα = (CL − C0) / (CL − Cα)

where C0 = overall composition, CL = liquidus composition, Cα = solidus composition at temperature T
🧪

Alloy Phase Diagram Simulator lets you explore binary eutectic, isomorphous, and peritectic systems. Click any point on the phase diagram to read off phases, compositions, and lever-rule fractions. Animate slow cooling of an alloy to see how microstructure evolves.

The Stress-Strain Curve

Applying a uniaxial load to a metal specimen produces a characteristic stress-strain curve that encodes nearly everything important about the material’s mechanical behaviour:

  1. Elastic region (Hooke’s law, σ = Eε): stress is proportional to strain; deformation is reversible. The slope is Young’s modulus E. For steel: E ≈ 200 GPa.
  2. Yield pointy): permanent plastic deformation begins. Dislocations start to move. Below σy, parts spring back; above it, they stay bent.
  3. Plastic hardening: as dislocations pile up and tangle, stress must increase to continue deformation. Described by the power law σ = σy + Kεpn (n = strain hardening exponent, typically 0.1–0.5).
  4. Ultimate Tensile Strength (UTS): peak engineering stress. Beyond UTS, necking localisation occurs and the cross-sectional area drops faster than hardening can compensate.
  5. Fracture: the specimen breaks. Ductile materials show significant plastic strain before fracture; brittle materials (ceramics, cast iron) show almost none.
Elastic: σ = Eε   (E = Young’s modulus)
Plastic: σ = σy + Kεpn
Toughness (area under curve) = ∫0εf σ dε     [J/m³]
Resilience (elastic area) = σy²/(2E)     [J/m³]
🟰

Stress-Strain Curve Simulator animates the tension test from zero load to fracture. An animated specimen shows elastic stretching, yielding, necking, and rupture. Compare five materials: steel, aluminium, rubber, bone, and polymer. The shaded toughness area updates in real time.

Shape Memory Alloys & Superelasticity

Shape memory alloys (SMAs) are a class of smart materials capable of recovering large strains (>6%) that would permanently deform any conventional metal. They exploit a reversible solid-state phase transformation between two crystallographic phases:

The transformation is reversible and thermoelastic: on heating, the “shape memory effect” drives martensite back to austenite, recovering the parent shape. On cooling below Af without applied stress, the material transforms back to martensite.

Superelastic (pseudoelastic) behaviour

At temperatures above Af, applying sufficient stress mechanically induces the martensitic transformation. Releasing the stress reverses the transformation and the material springs back completely — despite having undergone 6–8% strain. This is superelasticity, familiar from eyeglass frames that can be bent without breaking.

The key thermodynamic relationship coupling temperature and stress to the transformation is the Clausius-Clapeyron equation:

dσ / dT = −ΔH / (ε0 · T0)

ΔH = latent heat of transformation     ε0 = transformation strain (6–8% for Nitinol)
T0 = equilibrium transformation temperature

Typical slope: 5–8 MPa/K for Nitinol NiTi

The stress-strain hysteresis loop encloses an area equal to the energy dissipated per cycle as latent heat of transformation. This makes SMAs excellent vibration dampers as well as actuators.

🔧

Shape Memory Alloy Simulator shows the austenite-martensite lattice transformation alongside the superelastic stress-strain hysteresis loop. Adjust temperature and stress to move between phases. Switch modes to see shape memory vs. superelastic cycles. Presets for Nitinol, Cu-Al-Ni, Fe-Mn-Si, and Cu-Zn-Al.

Material presets

Key SMA systems

Nitinol NiTi: Ms ≈ −20..+80°C (tunable by Ni:Ti ratio), recovery strain 8%, biocompatible. Medical: stents, guidewires, orthodontic archwires.
Cu-Al-Ni: Ms up to 200°C, cheaper than Ni, brittler, used in high-temp actuators.
Fe-Mn-Si: ferromagnetic, large hysteresis, pipeline coupling applications.
Cu-Zn-Al: low cost, Ms −100..+120°C, used in thermal actuators.

Electronic Band Structure

The electrical, optical, and magnetic properties of a solid are largely determined by its electronic band structure: the allowed energy levels for electrons as a function of crystal momentum k. Quantum mechanical treatment (Bloch’s theorem) shows that in a periodic potential, electron wave functions take the form of a plane wave modulated by the lattice periodicity:

ψn,k(r) = un,k(r) · eik·r     (Bloch theorem)

un,k has the same periodicity as the crystal lattice
n = band index     k = crystal momentum in first Brillouin zone

Solving the Schrödinger equation with a periodic potential produces energy bands (allowed energy ranges) separated by band gaps (forbidden ranges). The electrical character of a material is determined by how these bands are filled:

Band Structure Simulator visualises electron energy E(k) curves along the first Brillouin zone for free electrons, the nearly-free electron model (weak periodic potential), and tight-binding. Adjust the lattice potential strength to watch the gap open at zone boundaries.

Superconductivity

Below a critical temperature Tc, certain materials exhibit two hallmark quantum phenomena: (1) zero electrical resistance and (2) expulsion of magnetic flux (Meissner effect). Classical physics cannot explain either. BCS (Bardeen-Cooper-Schrieffer) theory (1957, Nobel Prize 1972) provides the microscopic mechanism:

In a metal, positively charged ions respond to the passage of an electron by being slightly displaced. The resulting lattice polarisation can attract a second electron, forming a weakly bound Cooper pair with momentum +k and −k and opposite spins. Cooper pairs have integer spin (bosonic) and can condense into a single macroscopic quantum state described by an order parameter ψ = |ψ|e.

BCS energy gap: Δ(0) ≈ 1.764 kBTc

Coherence length: ξ = ℏvF/(πΔ)
London penetration depth: λL = √(m*/(nsμ0e²))

Ginzburg-Landau parameter: κ = λ/ξ
Type I: κ < 1/√2     Type II: κ > 1/√2

High-temperature superconductors (Tc > 30 K) like YBCO (93 K) involve a more complex pairing mechanism still under active research, with d-wave symmetry in the copper-oxide planes.

🧲

Superconductivity Simulator visualises BCS Cooper pair formation in a 2D crystal lattice, the Meissner effect with field lines expelled from the bulk, and Type II Abrikosov vortex lattice. Switch between Hg, Pb, Nb, Nb3Sn, and YBCO presets to explore how Tc varies across material families.

Adsorption on Surfaces

Surface science is a crucial branch of materials science: heterogeneous catalysis, corrosion, sensor design, and thin-film deposition all depend on how molecules interact with surfaces. The simplest quantitative model of adsorption is the Langmuir isotherm:

Langmuir: θ = KP / (1 + KP)

θ = fractional surface coverage    K = adsorption equilibrium constant    P = partial pressure

BET (multi-layer): q/qm = Cx / [(1−x)(1−x+Cx)]    x = P/P0

Freundlich: q = K · P1/n     (empirical, heterogeneous surfaces)

The BET (Brunauer-Emmett-Teller) isotherm extends Langmuir to multi-layer physisorption and is the standard method for measuring specific surface area of porous materials (N2 adsorption at 77 K). A BET surface area of 1000 m²/g is common for activated carbons used in water purification and energy storage.

🧪

Adsorption Isotherm Simulator animates gas molecules landing on and desorbing from a 28-site surface and plots all three isotherm models (Langmuir, BET, Freundlich) simultaneously so you can see exactly where they agree and diverge. Adjust affinity K, BET constant C, and the Freundlich 1/n to map out the parameter space.

Materials Science in the Platform

The materials science category on mysimulator.uk is growing rapidly. Here is a summary of what is currently available:

Further Exploration

Materials science connects naturally to neighbouring fields covered elsewhere on the platform:

← Spotlight #47: Galaxies & Dark Matter