Fick's Laws of Diffusion

Fick's laws of diffusion are a pair of fundamental constitutive equations formulated in 1855 by the German physiologist and physicist Adolf Eugen Fick that describe macroscopic mass transport by diffusion in gases, liquids, solids, and biological tissues. Recognizing the mathematical analogy between the random thermal migration of dissolved salt ions in water and Joseph Fourier's 1822 law of heat conduction (as well as Georg Ohm's 1827 law of electrical current), Fick established a quantitative bridge between local spatial concentration gradients and net molar or mass flux.

Fick's First Law applies to steady-state systems (or instantaneous flux profiles), stating that the diffusive flux vector J is directly proportional to the negative spatial gradient of concentration (−∇*c*), scaled by the diffusion coefficient D. Fick's Second Law combines this constitutive flux relationship with the conservation of mass (the continuity equation) to yield a parabolic partial differential equation—mathematically identical to the heat equation—that governs how a non-steady-state concentration field c(r, t) evolves over space and time. Together, these equations underpin quantitative models of respiratory gas exchange, cell membrane permeability, synapse neurotransmitter clearance, semiconductor dopant profiles, and Alan Turing's reaction–diffusion morphogenesis.

Mathematical Formulation in One and Three Dimensions

Fick's First Law postulates that in an isotropic medium without bulk fluid motion, the net diffusive flux J (the amount of substance passing perpendicularly through a unit cross-sectional area per unit time, with SI units of mol·m⁻²·s⁻¹ or particles·m⁻²·s⁻¹) is proportional to the spatial derivative of molar concentration c (mol·m⁻³):

`` One Dimension: J_x = -D (∂c / ∂x) Three Dimensions: J = -D ∇c = -D (x̂ ∂c/∂x + ŷ ∂c/∂y + ẑ ∂c/∂z) ``

The negative sign indicates that net transport proceeds down the concentration gradient, from regions of high concentration to regions of low concentration. In non-ideal solutions where activity coefficients deviate from unity, the thermodynamic driving force is the spatial gradient of the chemical potential (μ = μ° + RT ln a), yielding the generalized Onsager–Fick relation:

`` J = - (D · c / RT) ∇μ ``

When the diffusion coefficient D is independent of concentration and a steady-state gradient is maintained across a planar slab or membrane of thickness Δ*x* bounded by constant concentrations c₁ and c₂, integrating Fick's First Law yields a linear concentration profile with constant flux J = D(c₁ − c₂) / Δ*x*.

Membrane Permeability and Physiological Gas Exchange

When a solute diffuses across a biological cell membrane or respiratory epithelium of surface area A and thickness Δ*x*, the intramembrane concentration equals the aqueous concentration multiplied by the dimensionless lipid–water partition coefficient K. Integrating Fick's First Law across the barrier gives the total molar transport rate (dn/dt = *J*·*A*):

`` dn / dt = - P_m · A · Δc = - (D · K · A / Δx) (c_in - c_out) ``

where P_m = DK/Δ*x* is the membrane permeability coefficient (m·s⁻¹). For dissolved respiratory gases such as O₂ and CO₂, Henry's Law relates molar concentration to partial pressure (c = *S*·*P*, where S is gas solubility). Substituting Henry's Law yields the physiological form of Fick's First Law for pulmonary or gill gas exchange:

`` V̇_gas = (A · D · S / Δx) · (P₁ - P₂) = D_L · ΔP ``

where V̇_gas is the volumetric gas uptake rate and D_L = ADS/Δ*x* is the diffusing capacity of the lung. This relationship explains why vertebrate lungs and teleost fish gills maximize surface area A (~70–100 m² in adult human alveoli), minimize barrier thickness Δ*x* (~0.2–0.5 μm across type I pneumocytes and capillary endothelium), and employ countercurrent exchange in fish gills to sustain a positive partial pressure difference Δ*P* along the entire capillary length.

Derivation from the Continuity Equation

When solute concentrations vary with time, mass conservation within an infinitesimal control volume requires that the local rate of concentration change (∂*c*/∂*t*) equals the negative divergence of the flux vector (−∇·J), provided no chemical reactions generate or consume the species:

`` Continuity Equation: ∂c / ∂t + ∇ · J = 0 ``

Substituting Fick's First Law (J = −*D*∇*c*) into the continuity equation yields Fick's Second Law in its general heterogeneous form:

`` ∂c / ∂t = ∇ · (D ∇c) ``

In an isotropic, homogeneous medium where the diffusion coefficient D is spatially invariant, D factors out of the differential operator to give the linear diffusion equation:

`` One Dimension: ∂c / ∂t = D (∂²c / ∂x²) Three Dimensions: ∂c / ∂t = D ∇²c = D (∂²c/∂x² + ∂²c/∂y² + ∂²c/∂z²) ``

Physically, the second spatial derivative (∂²c/∂*x*², the curvature or concavity of the concentration profile) determines whether a local region gains or loses solute: where the concentration profile is concave upward (∂²c/∂*x*² > 0), more solute enters than leaves and ∂*c*/∂*t* > 0; where the profile is linear (∂²c/∂*x*² = 0), steady state prevails and ∂*c*/∂*t* = 0.

Fundamental Gaussian Propagator and Error-Function Solutions

Two analytical solutions of Fick's Second Law serve as foundational models in biophysics and analytical chemistry:

1. Instantaneous Point-Source Solution (Green's Function): If N moles (or particles) of solute are deposited at t = 0 as a Dirac delta pulse at the origin x = 0 in an unbounded one-dimensional domain (c(x, 0) = N δ(x)), solving Fick's Second Law via Fourier transformation yields a spreading Gaussian normal distribution:

`` c(x, t) = [ N / √(4 π D t) ] · exp( -x² / (4 D t) ) `` In n spatial dimensions (n = 1, 2, 3), the pre-exponential factor generalizes to N / (4π*Dt*)^n/2. The variance of this Gaussian distribution is σ² = 2*Dt*, recovering Einstein's mean-squared displacement relation ⟨*x*²⟩ = 2*Dt*.

1. Semi-Infinite Constant-Surface Boundary Solution: If a semi-infinite medium (x ≥ 0) initially at concentration c₀ has its boundary at x = 0 stepped to a constant concentration c_s for all t > 0, applying the Boltzmann similarity variable η = x / (2√(Dt)) transforms the partial differential equation into an ordinary differential equation whose solution involves the complementary error function (erfc):

`` [ c(x, t) - c₀ ] / [ c_s - c₀ ] = erfc( x / [2 √(D t)] ) = 1 - erf( x / [2 √(D t)] ) ``

Coupled Chemical Kinetics and Turing Patterns

In living tissues, diffusing signaling molecules simultaneously undergo enzymatic synthesis, receptor capture, or proteolytic degradation. Incorporating a chemical reaction rate term R(c) into Fick's Second Law produces the reaction–diffusion equation:

`` ∂c / ∂t = D ∇²c + R(c) ``

For a developmental morphogen (such as the transcription factor Bicoid in the Drosophila syncytial embryo) secreted at x = 0 with constant flux J₀ while undergoing first-order linear degradation at rate constant k (R(c) = −*k*·*c*), setting ∂*c*/∂*t* = 0 yields a steady-state profile that decays with an exponential distribution over space:

`` c(x) = (J₀ / √(D · k)) · exp( -x / λ ), where λ = √(D / k) ``

Here, the characteristic decay length λ = √(D/k) determines the spatial range over which target genes are activated at distinct concentration thresholds. In 1952, Alan Turing demonstrated that when two morphogens—a short-range autocatalytic activator (u, with small diffusion coefficient D_u) and a long-range inhibitor (v, with D_v ≫ D_u)—interact via coupled nonlinear reaction–diffusion equations, a homogeneous steady state can undergo diffusion-driven instability, spontaneously generating periodic stationary stripes, spots, and anatomical patterns.

Why does Fick's First Law contain a negative sign?

The spatial gradient vector ∇*c* (or derivative ∂*c*/∂*x* in one dimension) points in the direction of increasing concentration. Because random thermal motion causes more molecules to leave a densely populated region than a sparsely populated one, the net macroscopic flux J always points in the opposite direction—from high concentration toward low concentration. The minus sign ensures that when ∂*c*/∂*x* is negative, the flux J_x in the positive x-direction is a positive quantity for a positive diffusion coefficient D.

How do clinical diseases like pulmonary edema and fibrosis impair gas exchange via Fick's Law?

According to the physiological form of Fick's First Law (V̇_gas = (*A*·*D*·*S* / Δ*x*)·Δ*P*), the rate of oxygen transfer into pulmonary capillaries is directly proportional to alveolar surface area A and inversely proportional to barrier thickness Δ*x*. Pulmonary edema (fluid accumulation in the interstitial space) and pulmonary fibrosis (collagen deposition) increase Δ*x* from ~0.3 μm to several micrometers, whereas emphysema destroys alveolar septa and reduces A. Both pathologies lower the lung diffusing capacity D_L, causing arterial hypoxemia during exercise.

When do Fick's laws break down in biological or polymeric systems?

Classical Fickian diffusion assumes a dilute solute moving through a homogeneous, memoryless fluid where the mean-squared displacement scales linearly with time (⟨*r*²⟩ ∝ t¹). Non-Fickian (anomalous) diffusion arises in glassy polymers where solvent penetration is coupled to slow viscoelastic polymer chain relaxation (Case II transport), and inside crowded intracellular matrices or chromatin networks where transient binding traps and cytoskeletal obstacles produce subdiffusive power-law scaling (⟨*r*²⟩ ∝ t^α with 0.5 < α < 0.9) governed by fractional diffusion equations.

1. Fick, A. (1855). "Ueber Diffusion." Annalen der Physik und Chemie, 170(1), 59–86. https://doi.org/10.1002/andp.18551700105
2. Crank, J. (1975). The Mathematics of Diffusion (2nd ed.). Oxford University Press. https://doi.org/10.1021/ja01565a079
3. Turing, A. M. (1952). "The chemical basis of morphogenesis." Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 237(641), 37–72. https://doi.org/10.1098/rstb.1952.0012
4. Weibel, E. R. (1984). The Pathway for Oxygen: Structure and Function in the Mammalian Respiratory System. Harvard University Press. https://doi.org/10.1016/0034-5687(73)90030-390030-3)