Diffusion in Solids

Introduction

Atoms and ions in a solid are not permanently frozen. At finite temperature they vibrate and occasionally cross barriers to neighbouring sites. A substitutional atom may move when a neighbouring regular site is vacant; a small interstitial species can hop among spaces between host atoms. Large-scale chemical mixing, oxide growth, sintering and battery operation all depend on these microscopic jumps. Diffusion laws connect those jumps to measured concentration gradients and time-dependent profiles.

Core explanation

For one-dimensional transport down a concentration gradient, Fick's first law is J = −D ∂C/∂x. Here J is flux, for example mol m⁻² s⁻¹, C is concentration in mol m⁻³, x is position and D is diffusion coefficient in m² s⁻¹. The negative sign says particles flow from higher toward lower concentration under the simple ideal concentration-gradient model. In multicomponent nonideal materials, chemical-potential gradients rather than concentration alone are the fundamental thermodynamic force, and cross-coupling may occur. Fick's law is a useful effective description when its assumptions fit.

Mass conservation gives ∂C/∂t = −∂J/∂x. If D is spatially uniform and does not depend on C, substitution yields Fick's second law , ∂C/∂t = D ∂²C/∂x². It predicts how a profile smooths with time. A characteristic diffusion distance scales like √(Dt), with a numerical coefficient depending on the geometry and exactly how distance is defined. The formula is not a claim that every atom travels the same distance; random walks produce a distribution of displacements. For a one-dimensional unbiased walk with uncorrelated steps, the mean-square displacement is ⟨x²⟩ = 2Dt. MIT's diffusion teaching unit develops flux, Fick's laws and diffusion profiles together.

Microscopically, a vacancy mechanism requires an adjacent vacancy and a successful jump over a migration barrier E m. If vacancies are in thermal equilibrium and require formation energy E f, self-diffusion may have an effective activation energy approximately E f + E m in a simple dilute model. If a dopant fixes vacancy concentration over the temperature range, the measured activation can be closer to a migration contribution, with correlation and association effects added. An interstitial mechanism can avoid a vacancy-formation penalty but may have a different migration barrier. Thus the common Arrhenius form D = D₀ exp(−Q/RT) , with molar Q, is empirical over a suitable range; a slope change may signal a different mechanism or defect regime. Equivalently use E a/k BT when E a is per particle. MIT recitation notes contrast vacancy and interstitial paths and write the Arrhenius dependence.

The diffusion coefficient measures random spreading in a specified context. A self-diffusion tracer coefficient, chemical diffusion coefficient and ionic conductivity-derived diffusion coefficient need not be numerically identical because they weight correlations and thermodynamic factors differently. Grain boundaries and dislocations can create pathways with different transport rates from bulk lattice diffusion. A polycrystal's apparent value depends on microstructure and experimental length scale. NIST's diffusivity and mobility data treatment surveys bulk mechanisms and the conditions underlying Fick-like models.

Temperature has two effects in vacancy-controlled diffusion: it can create more vacancies and help atoms jump more frequently. Do not infer an intrinsic activation energy merely by reading a single diffusion coefficient. At high temperature the crystal can undergo phase changes, defect association can shift, or concentration dependence can become strong. Conversely, at very low temperature a metastable defect population may be frozen in, making a simple equilibrium Arrhenius extrapolation unreliable.

Step-by-step reasoning

1. Identify the moving species and whether it uses a vacancy, interstitial or another pathway. 2. Specify flux and concentration units and the sign of the gradient. 3. Use the first law for flux at a known gradient; combine with conservation for time evolution. 4. Estimate a characteristic length from √(Dt) only with geometry understood. 5. Fit ln D versus 1/T over a single mechanism regime and interpret Q cautiously.

Visual explanation

Draw a row of atoms with one empty regular site; an adjacent atom hops into it, making the vacancy appear to move oppositely. Below draw a smaller ion hopping between interstitial pockets. Alongside plot a concentration hill that broadens and lowers with time; arrows point down its sides, and an inset ln D versus 1/T graph slopes downward.

Real-world analogy

People moving through a theatre can shuffle into empty seats, resembling vacancy migration, or small staff members can pass through aisles between occupied seats, resembling interstitial migration. A crowd initially gathered in one area gradually spreads; the average trend obeys a diffusion equation even though individual paths are random and different.

Real-world example

Carbon can diffuse through interstitial positions in iron under appropriate phase and temperature conditions, enabling steel heat treatment. Substitutional alloying elements often move more slowly because they need vacancies. The actual rates depend on the iron phase, composition and temperature; comparing mechanisms clarifies why the atoms' sizes and site types matter.

Why?

Why does diffusion distance grow approximately as the square root of time? Random jumps in opposite directions partly cancel. Mean-square displacement grows in proportion to jump count and time, so the root-mean-square displacement grows as √t, not linearly with t.

Common misconception

"D is the speed of a diffusing atom." D has units of area per time, not length per time. It measures ensemble spreading. Individual atoms jump intermittently and may move rapidly during a jump while having a small long-term net displacement.

Worked example

Question: A species has D = 1.0 × 10⁻¹² m² s⁻¹ in a one-dimensional simple diffusion model. Estimate its root-mean-square displacement after t = 100 s.

Reasoning: For one-dimensional unbiased diffusion, ⟨x²⟩ = 2Dt. Substitution gives 2 × 1.0 × 10⁻¹² × 100 = 2.0 × 10⁻¹⁰ m². Taking the square root gives 1.41 × 10⁻⁵ m, or 14 μm. This is an rms ensemble scale; some particles travel farther and others less, and a different geometry gives a different coefficient.

Answer: Approximately 14 μm root-mean-square displacement in one dimension.

Quick check

1. What is the sign of J if concentration increases toward positive x and D is positive? Answer: Negative; material diffuses toward decreasing x under Fick's first law.

Exam focus

Keep J, C, D and x units consistent and show the minus sign. Derive Fick's second law from mass conservation when D is constant. Distinguish vacancy formation from migration energy and state when Q can include both.

Advanced insight

In a crystal, jumps can be correlated: after an atom moves, it may be more likely to jump back into the vacancy it just left. A correlation factor modifies the tracer diffusion coefficient relative to a naive random-jump count. For charged species, electric fields add drift, giving a flux with diffusion and migration terms; ionic conductivity therefore cannot always be inferred from Fick's concentration-gradient equation alone.

Summary

Solid diffusion occurs through vacancy, interstitial and microstructural pathways. Fick's first law links flux to gradient, and the second follows from mass conservation for constant D. Diffusion coefficients often follow Arrhenius behaviour, but activation energy depends on defect availability and migration as well as the measurement regime.

Practice questions

1. What units does D have in J = −D∂C/∂x? Answer: Square metres per second when J is mol m⁻² s⁻¹ and C is mol m⁻³. 2. Which mechanism requires an empty regular site next to the moving atom? Answer: Vacancy migration. 3. If time quadruples at fixed D, how does a characteristic √(Dt) distance change? Answer: It doubles under the same geometry and definition. 4. Why can an Arrhenius slope change after doping? Answer: Dopants can fix or associate vacancies, changing whether defect formation and migration both contribute to activation.