Diffusion in solids is the thermally activated migration of atoms through a crystal lattice, primarily via vacancy exchange or interstitial hopping. Fick's first law relates steady-state flux to a concentration gradient; Fick's second law describes time-dependent concentration profiles. The diffusivity D follows an Arrhenius relationship D = D₀ exp(−Qd/RT), where Qd is the activation energy for diffusion. Engineering processes such as carburization (adding carbon to steel surfaces) and dopant diffusion in semiconductors are directly governed by these principles.
Solve Fick's second law for the semi-infinite solid boundary condition (using the complementary error function solution) applied to carburization problems. Plot concentration vs. depth at different times to build intuition.
Diffusion in solids is superficially similar to diffusion in liquids or gases, but the rigid crystal lattice changes everything. Atoms in a solid are not free to wander — they are trapped in potential wells at lattice sites. For a substitutional atom (one occupying a regular lattice site) to move, it must jump into an adjacent vacancy, and vacancies are rare. This is why substitutional diffusion is slow: the atom must wait for both a thermally activated jump and a neighboring empty site. Interstitial atoms — like carbon squeezed into the gaps of an iron lattice — face a different situation: the interstitial sites are always "available," so the only barrier is the activation energy to squeeze through the lattice. Interstitial diffusion is therefore much faster than substitutional diffusion, even in the same material.
The temperature dependence of diffusivity is captured by the Arrhenius equation D = D₀ exp(−Qd / RT), where Qd is the activation energy, R is the gas constant, and T is the absolute temperature. This is the same form you encountered in chemical kinetics, and for the same reason: both processes require thermal energy to surmount an energy barrier. The exponential sensitivity to temperature means that small changes in T translate to large changes in D — a 50°C increase can change diffusivity by an order of magnitude. In practice, this is why heat-treatment temperatures are tightly controlled.
Fick's first law J = −D(dC/dx) describes the steady-state flux of atoms down a concentration gradient. But most engineering problems involve time-dependent concentration profiles, which requires Fick's second law: ∂C/∂t = D ∂²C/∂x². For the standard carburization setup — a semi-infinite steel bar with a fixed surface carbon concentration Cs exposed at t = 0 — the solution is C(x,t) = Cs − (Cs − C₀)·erf(x / 2√(Dt)), where erf is the error function and C₀ is the initial uniform carbon content. This solution encodes the idea that the "diffusion front" propagates inward as √(Dt): doubling time moves carbon not twice as deep, but only √2 times as deep.
To use this solution, you identify x (depth below surface), t (exposure time), D (diffusivity at the treatment temperature, calculated from the Arrhenius formula), and the boundary/initial conditions. The practical goal in carburization is to achieve a target carbon concentration at a target depth — for example, 0.4 wt% C at 1 mm depth — and you solve for the required time or temperature. This links your abstract understanding of Fick's law back to the hardness profile of a manufactured gear tooth.