Understanding Diffusion in Gases and Liquids: A Comprehensive Overview

1. Lecture 17: Diffusion PHYS 430/603 material Laszlo Takacs UMBC Department of Physics

2. Diffusion in gases and liquids

3. Types of diffusion • Interstitial diffusion - C in Fe, H in Ti, etc. • Vacancy diffusion • Self diffusion, tracer diffusion - healing of defects, time-dependent mechanical phenomena • Chemical diffusion - diffusion of one element through another, solid state chemical reactions • Bulk diffusion • Diffusion in grain boundaries, lattice defects • Surface diffusion • Diffusion is a statistical phenomenon, driven by the random motion of atoms.

4. Phenomenological description Fick’s law: Continuity: Combined: If D = const. jD in cm2/s; c in particles/cm3 To be solved for c(x,t) or c(r,t) Typical solution depends on x/(Dt)1/2. One can define a diffusion front that travels as x2 ~ 6Dt If defined that way, about 99% of the diffusing element is still behind the diffusion front.

5. Evaluating a diffusion couple Assume constant D and symmetrical solution to x = 0: Introduce new variables:

6. We can guess the solution Adjust G and H to the boundary conditions

7. Diffusion is a thermally activated process, resulting in an exponential temperature dependence: D = D0 exp(-Q/kT) It can also depend on orientation relative to the crystal axes. Interstitials diffuse much faster than substitutional impurities and matrix atoms.

8. Rough correlations • D0 for self diffusion in metals ~ 1 cm2/s (0.1 to 10) • ∆HD / Tm ~ 1.5*10-3 eV/K • ∆HD / Lm ~ 15 These relationships are useful if no measured data are available. Typically estimates of diffusion are quite rough, there is a huge variation, high accuracy is not needed. Also, defects and impurities influence diffusion very much, it is difficult to control it.

9. Examples Temperature is uniformly T = 912°C, the alpha-gamma transition of iron. Diffusion time t = 1 h =3600 s, distance x2 = 6Dt Self diffusion in alpha Fe: D0 = 2.0 cm2/s, ∆HD = 2.5 eV D = 2.0 exp(-2.5/8.6*10-5*11850) = 4.4*10-11 cm2/s, x ~ 10 µm Self diffusion in gamma Fe: D0 = 0.4 cm2/s, ∆HD = 2.8 eV D = 4.6*10-13 cm2/s, x ~ 1 µm 60Co in alpha Fe: D0 = 0.2 cm2/s, ∆HD = 2.4 eV D = 1.2*10-11 cm2/s, x ~ 5 µm C in alpha Fe: D0 = 0.004 cm2/s, ∆HD = 0.83 eV D = 1.2*10-6 cm2/s, x ~ 1.6 mm C in gamma Fe: D0 = 0.67 cm2/s, ∆HD = 1.6 eV D = 1.0*10-7 cm2/s, x ~ 0.5 mm Notice the very large differences.

10. Vacancy diffusion in a solid anddislocation climb due to capturing vacancies

11. A simple atomistic view of diffusion Let’s consider two neighboring lattice planes of area A (left and right) at a distance d, perpendicular to the x direction. Nl = Ad c and Nr = Ad [c + (dc/dx)*d] If  is the jump frequency and jump in 6 mutually perpendicular directions has the same probability: jr = /6 (Adc) / A and jl = /6 [Ad [c + (dc/dx)d] / A j = jr - jl = - d2/6 dc/dx = - D dc/dx, thus D = d2/6

12. Problems • Geometry is more complicated, even for the simplest case of interstitial diffusion. 1/6 assumption is good for simple cubic, bcc gives the same result. An anisotropic hcp lattice is certainly different. • At higher concentration, the jump probability depends on the occupation of the potential target sites. • For substitutional atoms (self diffusion included) jump can only occur, if the target site is vacant. Estimate for D must be multiplied by the vacancy concentration, cv = exp( - Gf/kT ). This factor is not needed for vacancy diffusion, it is much faster. • Probability of jump depends on previous jump - correlation effects. • Estimating jump frequency is very difficult - but needed for quantitative estimates:  =  exp( - Gm /kT ) Gm migration free energy,  Debye frequency. • Diffusion current should be related to the Gibbs free energy, rather than concentration.

13. More than one diffusing species Diffusion is usually asymmetric, the boundary between regions is shifting, often difficult to define. Very importantly, solid state reactions always take place at interfaces via diffusion. The products form intermediate layers that often act as diffusion barriers. Decreasing the particle size - e.g. by ball milling - accelerates the reactions. Diffusion relates to volume changes, formation of voids, etc. According to the experts, diffusion is much more complicated than it looks for a superficial first look. Rather than focusing on the difficult delicacies of the theoretical description, let’s consider a few examples.

14. The Kirkendall experiment Zn diffuses faster than Cu in this system. As a result, the Cu - brass interface shifts, the Mo wires marking the interface drift toward each other.

15. After 15 h at 500°C layered structure developed that exhibits every compound in the phase diagram in order of composition. The layer thickness depends on the diffusion coefficients. Concentration gradient exists within the layer, if the phase diagram allows. Each layer is single-phase, two phase regions do not form.

17. Surface and interface diffusion • More open, diffusion can take place faster. In fact, dislocations provide routes for fast diffusion also. Diffusion in grain boundaries makes predicting the behavior of polycrystals difficult. • Activation energy is lowest for surface diffusion, larger for grain boundary diffusion, and the highest for bulk diffusion. Therefore, bulk diffusion dominates at high temperature (larger number of possible sites) but grain boundary and surface diffusion takes over at lower temperature. • Technological importance: Thin film deposition, weakening of grain boundaries by impurities, shape changes of particles, etc. • Diffusion can be accelerated by introducing grain boundaries and other lattice defects. This is an important benefit from MA (mechanical alloying).