Atomic Diffusion in a Crystal Lattice

Atoms don't sit still, they hop

A crystal lattice looks static in a textbook diagram, but at any real temperature above absolute zero its atoms are constantly vibrating, and every so often one of them acquires enough thermal energy to jump into a neighbouring empty lattice site — a vacancy. Over many such hops, driven by nothing more than random thermal motion, atoms migrate through the solid. This is solid-state diffusion, and unlike diffusion in a gas or liquid, it depends critically on the crystal having defects: a perfect, defect-free lattice has essentially nowhere for an atom to hop into.

Fick's laws: diffusion as a flow down a gradient

Adolf Fick's first law (1855) says the flux of diffusing atoms is proportional to how steep the concentration gradient is — atoms net-flow from where they are crowded to where they are sparse:

Fick's 1st law:  J = -D · (∂C/∂x)
Fick's 2nd law:  ∂C/∂t = D · (∂²C/∂x²)

J = flux (atoms per area per time)
C = concentration, x = position, t = time
D = diffusivity (diffusion coefficient)

The second law follows from the first by conservation of atoms and describes how a concentration profile evolves over time — it is the same equation that governs heat conduction, with concentration playing the role of temperature. A sharp initial concentration step smooths into an error-function profile whose width grows as √(Dt): diffusion distance scales with the square root of time, not linearly, which is why doubling a diffusion depth takes four times as long, not twice.

The vacancy mechanism

In most metals and many ionic crystals, atoms move by the vacancy mechanism: an atom adjacent to a vacant lattice site swaps places with the vacancy. Since vacancies themselves are relatively rare — their equilibrium concentration follows an Arrhenius law in the vacancy formation energy — diffusivity depends on both how many vacancies exist to hop into and how easily an atom can clear the energy barrier of the hop itself. Other mechanisms exist (interstitial diffusion for small atoms like carbon or hydrogen squeezing between lattice sites, and interstitialcy or ring mechanisms in specific systems) but vacancy diffusion dominates self-diffusion and substitutional-alloy diffusion in most metals.

Arrhenius behaviour: why diffusion is so temperature-sensitive

The diffusivity D itself is not constant — it depends exponentially on temperature through an Arrhenius relation:

D(T) = D₀ · exp( -Q / (R·T) )

D₀ = pre-exponential factor (attempt frequency × geometric factors)
Q  = activation energy for diffusion (vacancy formation + migration energy)
R  = gas constant, T = absolute temperature

Because Q sits in an exponent, diffusivity is ferociously sensitive to temperature — a modest 10% rise in absolute temperature can multiply D several-fold for a typical metallic activation energy. This exponential sensitivity is why heat treatments in metallurgy (carburising, annealing, sintering) are run at carefully controlled elevated temperatures: a small temperature error compounds into a large error in how deep a diffusion process actually reaches within a given time.

Why this matters in metallurgy

Solid-state diffusion via vacancies is the mechanism behind case-hardening steel (carbon diffusing into the surface layer), homogenising cast alloys, sintering powder metallurgy parts into solid components, and the slow microstructural degradation — creep, precipitate coarsening — that limits how hot a jet turbine blade can run. Engineers use the Arrhenius form directly: measure D at two temperatures, extract Q, and predict processing times or in-service degradation rates at any other temperature within the same mechanism's regime.

Frequently asked questions

Try it live

Everything above runs in your browser — open Crystal Diffusion and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.