Fick's law of diffusion

From Wikipedia, the free encyclopedia

For the technique of measuring cardiac output, see Fick principle.

Fick's laws of diffusion describe diffusion and can be used to solve for the diffusion coefficient D. They were derived by Adolf Fick in the year 1855.

1 First law
2 Second law
3 Example solution in one dimension: diffusion length
4 Applicability
5 Temperature dependence of the diffusion coefficient
6 Biological perspective
7 Semiconductor fabrication applications
8 See also
9 References
- 9.1 Notations
- 9.2 Footnotes
10 External links

[edit] First law

Fick's first law is used in steady-state diffusion, i.e., when the concentration within the diffusion volume does not change with respect to time $(\, J_\mathrm{in} = J_\mathrm{out})$ . In one (spatial) dimension, this is

$J = - D \frac{\partial \phi}{\partial x}$

where

$J$ is the diffusion flux in dimensions of [(amount of substance) length⁻² time^-1], example $\bigg(\frac{\mathrm{mol}}{ m^2\cdot s}\bigg)$

$\, D$ is the diffusion coefficient or diffusivity in dimensions of [length² time⁻¹], example $\bigg(\frac{m^2}{s}\bigg)$

$\, \phi$ (for ideal mixtures) is the concentration in dimensions of [(amount of substance) length⁻³], example $\bigg(\frac\mathrm{mol}{m^3}\bigg)$

$\, x$ is the position [length], example $\,m$

$\, D$ is proportional to the velocity of the diffusing particles, which depends on the temperature, viscosity of the fluid and the size of the particles according to the Stokes-Einstein relation. In dilute aqueous solutions the diffusion coefficients of most ions are similar and have values that at room temperature are in the range of 0.6x10^-9 to 2x10^-9 m²/s. For biological molecules the diffusion coefficients normally range from 10^-11 to 10^-10 m²/s.

In two or more dimensions we must use $\nabla$ , the del or gradient operator, which generalises the first derivative, obtaining

$J=- D\nabla \phi$ .

The driving force for the one-dimensional diffusion is the quantity $- \frac{\partial \phi}{\partial x}$

which for ideal mixtures is the concentration gradient. In chemical systems other than ideal solutions or mixtures, the driving force for diffusion of each species is the gradient of chemical potential of this species. Then Fick's first law (one-dimensional case) can be written as:

$J_i = - \frac{D c_i}{RT} \frac{\partial \mu_i}{\partial x}$

where the index i denotes the ith species, c is the concentration (mol/m³), R is the universal gas constant (J/(K mol)), T is the absolute temperature (K), and μ is the chemical potential (J/mol).

[edit] Second law

Fick's second law is used in non-steady or continually changing state diffusion, i.e., when the concentration within the diffusion volume changes with respect to time.

$\frac{\partial \phi}{\partial t} = D\,\frac{\partial^2 \phi}{\partial x^2}\,\!$

Where

$\,\phi$ is the concentration in dimensions of [(amount of substance) length^-3], [mol m^-3]
$\, t$ is time [s]
$\, D$ is the diffusion coefficient in dimensions of [length² time^-1], [m² s^-1]
$\, x$ is the position [length], [m]

It can be derived from the Fick's First law and the mass balance:

$\frac{\partial \phi}{\partial t} =-\,\frac{\partial}{\partial x}\,J = \frac{\partial}{\partial x}\bigg(\,D\,\frac{\partial}{\partial x}\phi\,\bigg)\,\!$

Assuming the diffusion coefficient D to be a constant we can exchange the orders of the differentiating and multiplying by the constant:

$\frac{\partial}{\partial x}\bigg(\,D\,\frac{\partial}{\partial x} \phi\,\bigg) = D\,\frac{\partial}{\partial x} \frac{\partial}{\partial x} \,\phi = D\,\frac{\partial^2\phi}{\partial x^2}$

and, thus, receive the form of the Fick's equations as was stated above.

For the case of diffusion in two or more dimensions the Second Fick's Law is:

$\frac{\partial \phi}{\partial t} = D\,\nabla^2\,\phi\,\!$ ,

which is analogous to the heat equation.

If the diffusion coefficient is not a constant, but depends upon the coordinate and/or concentration, the Second Fick's Law becomes:

$\frac{\partial \phi}{\partial t} = \nabla \cdot (\,D\,\nabla\,\phi\,)\,\!$

An important example is the case where $φ$ is at a steady state, i.e. the concentration does not change by time, so that the left part of the above equation is identically zero. In one dimension with constant $\, D$ , the solution for the concentration will be a linear change of concentrations along $\, x$ . In two or more dimensions we obtain

$\nabla^2\,\phi =0\!$

which is Laplace's equation, the solutions to which are called harmonic functions by mathematicians.

[edit] Example solution in one dimension: diffusion length

A simple case of diffusion with time t in one dimension (taken as the x-axis) of a density $n (x, t)$ from a boundary located at position $x = 0$ where the density is maintained at a value $n (0)$ is

$n \left(x,t \right)=n(0)$ $\mathrm{erfc} \left(\begin{matrix}\frac{x}{\sqrt{4Dt}}\end{matrix}\right)$ .

where erfc is the complementary error function. The length $\sqrt{(4Dt)}$ is called the diffusion length and provides a measure of how far the density has propagated in the x-direction by diffusion in time t.

For more detail on diffusion length, see these examples.

[edit] Applicability

Equations based on Fick's law have been commonly used to model transport processes in foods, neurons, biopolymers, pharmaceuticals, porous soils, population dynamics, semiconductor doping process, etc. A large amount of experimental research in polymer science and food science has shown that a more general approach is required to describe transport of components in materials undergoing glass transition. In the vicinity of glass transition the flow behavior becomes "non-Fickian". See also non-diagonal coupled transport processes (Onsager relationship).

[edit] Temperature dependence of the diffusion coefficient

The diffusion coefficient at different temperatures is often found to be well predicted by

$D = D_0\cdot e^{-\frac{E_{A}}{R\cdot T}},$

where

$\, D$ is the diffusion coefficient
$\, D_0$ is the maximum diffusion coefficient (at infinite temperature)
$\, E_A$ is the activation energy for diffusion in dimensions of [energy (amount of substance)⁻¹]
$\, T$ is the temperature in units of [absolute temperature] (kelvins or degrees Rankine)
$\, R$ is the gas constant in dimensions of [energy temperature⁻¹ (amount of substance)⁻¹]

An equation of this form is known as the Arrhenius equation.

Typically, a compound's diffusion coefficient is ~10,000x greater in air than in water. Carbon dioxide in air has a diffusion coefficient of 16 mm²/s, and in water, its coefficient is 0.0016 mm²/s [1].

An approximate dependence of the diffusion coefficient on temperature can often be found using Stokes-Einstein equation, which predicts that:

$\frac {D_{T1}} {D_{T2}} = \frac {T_1} {T_2} \frac {\mu_{T2}} {\mu_{T1}}$

where:

T1 and T2 denote temperatures 1 and 2, respectively

D is the diffusion coefficient (m²/s)

T is the absolute temperature (K),

μ is the dynamic viscosity of the solvent (Pa·s)

[edit] Biological perspective

The first law gives rise to the following formula:^[1]

$Flux = {-P\cdot A\cdot (c_2 - c_1)}\,\!$

in which,

$\, P$ is the permeability, an experimentally determined membrane "conductance" for a given gas at a given temperature.
$\, A$ is the surface area over which diffusion is taking place.
$\, c_2 - c_1$ is the difference in concentration of the gas across the membrane for the direction of flow (from $c 1$ to $c 2$ ).

Fick's first law is also important in radiation transfer equations. However, in this context it becomes inaccurate when the diffusion constant is low and the radiation becomes limited by the speed of light rather than by the resistance of the material the radiation is flowing through. In this situation, one can use a flux limiter.

The exchange rate of a gas across a fluid membrane can be determined by using this law together with Graham's law.

[edit] Semiconductor fabrication applications

IC Fabrication technologies, model processes like CVD, Thermal Oxidation, and Wet Oxidation, Doping etc using Diffusion equations obtained from Ficks law.

In certain cases, the solutions are obtained for boundary conditions such as constant source concentration diffusion, limited source concentration, or moving boundary diffusion (where junction depth keeps moving into the substrate).

[edit] See also

[edit] References

[edit] Notations

A. Fick, Phil. Mag. (1855), 10, 30.
A. Fick, Poggendorff's Annel. Physik. (1855), 94, 59.
W.F. Smith, Foundations of Materials Science and Engineering 3^rd ed., McGraw-Hill (2004)
H.C. Berg, Random Walks in Biology, Princeton (1977)

[edit] Footnotes

^ Physiology at MCG 3/3ch9/s3ch9_2

[edit] External links

Diffusion fundamentals
Diffusion coefficient
Ficks Law Calculator and a host of other science tools - Rex Njoku & Dr.Anthony Steyermark

Categories: Diffusion | Statistical mechanics | Physical chemistry