Fick's laws of diffusion describe diffusion and were derived by Adolf Fick in 1855.

They can be used

Contents

[hide]

 1Fick's first law

 2Fick's second law

 3Derivation of Fick's laws

o 3.1Fick's first law

o 3.2Fick's second law

 4Derivation

o 4.1Example solution in one dimension: diffusion length

o 4.2Generalizations

 5Applications

o 5.1Biological perspective

o 5.2Fick's flow in liquids

o 5.3Semiconductor fabrication applications

 6History

 7See also

 8Notes

 9References

 10External links

Fick's first law[edit]

Fick's first law relates the diffusive flux to the concentration under the assumption of steady state. It
postulates that the flux goes from regions of high concentration to regions of low concentration,
with a magnitude that is proportional to the concentration gradient (spatial derivative), or in
simplistic terms the concept that a solute will move from a region of high concentration to a region
of low concentration across a concentration gradient. In one (spatial) dimension, the law is:

where

 J is the "diffusion flux," of which the dimension is amount of substance per unit area per unit
time, so it is expressed in such units as mol m−2 s−1. J measures the amount of substance that
will flow through a unit area during a unit time interval.
  D is the diffusion coefficient or diffusivity. Its dimension is area per unit time, so typical
units for expressing it would be m2/s.

 φ (for ideal mixtures) is the concentration, of which the dimension is amount of substance
per unit volume. It might be expressed in units of mol/m3.

 x is position, the dimension of which is length. It might thus be expressed in the unit m.

D is proportional to the squared 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.6 × 10−9 to 2 × 10−9 m2/s. For biological
molecules the diffusion coefficients normally range from 10−11 to 10−10 m2/s.

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

where J denotes the diffusion flux vector.

The driving force for the one-dimensional diffusion is the quantity −∂φ/∂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:

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

If the primary variable is mass fraction (yi, given, for example, in kg/kg), then the equation changes
to:

where ρ is the fluid density (for example, in kg/m3). Note that the density is outside
the gradient operator.

Fick's second law[edit]

Fick's second law predicts how diffusion causes the concentration to change with time. It is a partial
differential equation which in one dimension reads:

where

 φ is the concentration in dimensions of [(amount of substance) length−3], example

 t is time [s]
  D is the diffusion coefficient in dimensions of [length2 time−1], example m2/s

 x is the position [length], example m

In two or more dimensions we must use the Laplacian Δ = ∇2, which generalises the second
derivative, obtaining the equation

Derivation of Fick's laws[edit]

Fick's first law[edit]

In one dimension, the following derivation is based on a similar argument made in Berg 1977
(see references).

Consider a collection of particles performing a random walk in one dimension with length scale
Δx and time scale Δt. Let N(x,t) be the number of particles at position x at time t.

At a given time step, half of the particles would move left and half would move right. Since half of
the particles at point x move right and half of the particles at point x + Δx move left, the net
movement to the right is:

The flux, J, is this net movement of particles across some area element of area a, normal to the
random walk during a time interval Δt. Hence we may write:

Multiplying the top and bottom of the righthand side by (Δx)2 and rewriting, one obtains:

Concentration is defined as particles per unit volume, and hence

In addition, (Δx)2/2Δt is the definition of the diffusion constant in one dimension, D. Thus our
expression simplifies to:

In the limit where Δx is infinitesimal, the righthand side becomes a space derivative:

Fick's second law[edit]

Fick's second law can be derived from Fick's first law and the mass conservation in absence of any
chemical reactions:
Assuming the diffusion coefficient D to be a constant, one can exchange the orders of the
differentiation and multiply by the constant:

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

For the case of diffusion in two or more dimensions Fick's second law becomes

which is analogous to the heat equation.

If the diffusion coefficient is not a constant, but depends upon the coordinate or concentration,
Fick's second law yields

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

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

Derivation[edit]

Fick's second law is a special case of the convection–diffusion equation in which there is
no advective flux and no net volumetric source. It can be derived from the continuity equation:

where is the total flux and R is a net volumetric source for . The only source of flux in this
situation is assumed to be diffusive flux:

Plugging the definition of diffusive flux to the continuity equation and assuming there is no source
(R = 0), we arrive at Fick's second law:

If flux were the result of both diffusive flux and advective flux, the convection–diffusion equation is
the result.
Example solution in one dimension: diffusion length[edit]

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

where erfc is the complementary error function. This is the case when corrosive gases diffuse
through the oxidative layer towards the metal surface (if we assume that concentration of gases in
the environment is constant and the diffusion space (i. e., corrosion product layer) is semi-infinite –
starting at 0 at the surface and spreading infinitely deep in the material). If, in its turn, the diffusion
space is infinite (lasting both through the layer with n(x,0) = 0, x > 0 and that with n(x,0) = n0, x ≤ 0),
then the solution is amended only with coefficient  1⁄2 in front of n0 (this might seem obvious, as the
diffusion now occurs in both directions). This case is valid when some solution with
concentration n0 is put in contact with a layer of pure solvent. (Bokstein, 2005) The length 2√Dt is
called the diffusion length and provides a measure of how far the concentration has propagated in
the x-direction by diffusion in time t (Bird, 1976).

As a quick approximation of the error function, the first 2 terms of the Taylor series can be used:

If D is time-dependent, the diffusion length becomes

This idea is useful for estimating a diffusion length over a heating and cooling cycle, where D varies
with temperature.

Generalizations[edit]

1. In non-homogeneous media, the diffusion coefficient varies in space, D = D(x). This dependence
does not affect Fick's first law but the second law changes:

2. In anisotropic media, the diffusion coefficient depends on the direction. It is a

it is the product of a tensor and a vector:

For the diffusion equation this formula gives

3. For inhomogeneous anisotropic media these two forms of the diffusion equation should be
combined in

4. The approach based on Einstein's mobility and Teorell formula gives the following generalization
of Fick's equation for the multicomponent diffusion of the perfect components:

where φi are concentrations of the components and Dij is the matrix of coefficients. Here,
indices i, j are related to the various components and not to the space coordinates.

The Chapman–Enskog formulae for diffusion in gases include exactly the same terms. These physical
models of diffusion are different from the test models ∂tφi = ΣjDij Δφj which are valid for very small
deviations from the uniform equilibrium. Earlier, such terms were introduced in the Maxwell–Stefan
diffusion equation.

For anisotropic multicomponent diffusion coefficients one needs a rank-four tensor, for example
Dij,αβ, where i, j refer to the components and α, β = 1, 2, 3 correspond to the space coordinates.

Applications[edit]

Equations based on Fick's law have been commonly used to model transport processes in
foods, neurons, biopolymers, pharmaceuticals, poroussoils, population dynamics, nuclear
materials, plasma physics, and semiconductor doping processes. Theory of all voltammetric methods
is based on solutions of Fick's equation. Much 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". It can be shown that the Fick's law can be obtained from the Maxwell–
Stefan equations[1] of multi-component mass transfer. The Fick's law is limiting case of the Maxwell–
Stefan equations, when the mixture is extremely dilute and every chemical species is interacting only
with the bulk mixture and not with other species. To account for the presence of multiple species in
a non-dilute mixture, several variations of the Maxwell–Stefan equations are used. See also non-
diagonal coupled transport processes (Onsager relationship).

Biological perspective[edit]

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

in which,

 P is the permeability, an experimentally determined membrane "conductance" for a given

 c2 − c1 is the difference in concentration of the gas across the membrane for the direction of
flow (from c1 to c2).
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.

Fick's flow in liquids[edit]

When two miscible liquids are brought into contact, and diffusion takes place, the macroscopic (or
average) concentration evolves following Fick's law. On a mesoscopic scale, that is, between the
macroscopic scale described by Fick's law and molecular scale, where molecular random walkstake
place, fluctuations cannot be neglected. Such situations can be successfully modeled with Landau-
Lifshitz fluctuating hydrodynamics. In this theoretical framework, diffusion is due to fluctuations
whose dimensions range from the molecular scale to the macroscopic scale.[3]

In particular, fluctuating hydrodynamic equations include a Fick's flow term, with a given diffusion
coefficient, along with hydrodynamics equations and stochastic terms describing fluctuations. When
calculating the fluctuations with a perturbative approach, the zero order approximation is Fick's law.
The first order gives the fluctuations, and it comes out that fluctuations contribute to diffusion. This
represents somehow a tautology, since the phenomena described by a lower order approximation is
the result of a higher approximation: this problem is solved only by renormalizing the fluctuating
hydrodynamics equations.

Semiconductor fabrication applications[edit]

Integrated circuit fabrication technologies, model processes like CVD, thermal oxidation, wet
oxidation, doping, etc. use diffusion equations obtained from Fick's 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)