Fick's laws of diffusion: Difference between revisions

Latest revision as of 21:47, 13 November 2025

Template:Short description Script error: No such module "For". Template:Use dmy dates

Molecular diffusion from a microscopic and macroscopic point of view. Initially, there are solute molecules on the left side of a barrier (purple line) and none on the right. The barrier is removed, and the solute diffuses to fill the whole container. Top: A single molecule moves around randomly. Middle: With more molecules, there is a clear trend where the solute fills the container more and more uniformly. Bottom: With an enormous number of solute molecules, randomness becomes undetectable: The solute appears to move smoothly and systematically from high-concentration areas to low-concentration areas. This smooth flow is described by Fick's laws.

Fick's laws of diffusion describe diffusion and were first posited by Adolf Fick in 1855 on the basis of largely experimental results. They can be used to solve for the diffusion coefficient, Template:Mvar. Fick's first law can be used to derive his second law which in turn is identical to the diffusion equation.

Fick's first law: Movement of particles from high to low concentration (diffusive flux) is directly proportional to the particle's concentration gradient.^[1]

Fick's second law: Prediction of change in concentration gradient with time due to diffusion.

A diffusion process that obeys Fick's laws is called normal or Fickian diffusion; otherwise, it is called anomalous diffusion or non-Fickian diffusion.

History

In 1855, physiologist Adolf Fick first reported^[2] his now well-known laws governing the transport of mass through diffusive means. Fick's work was inspired by the earlier experiments of Thomas Graham, which fell short of proposing the fundamental laws for which Fick would become famous. Fick's law is analogous to the relationships discovered at the same epoch by other eminent scientists: Darcy's law (hydraulic flow), Ohm's law (charge transport), and Fourier's law (heat transport).

Fick's experiments (modeled on Graham's) dealt with measuring the concentrations and fluxes of salt, diffusing between two reservoirs through tubes of water. It is notable that Fick's work primarily concerned diffusion in fluids, because at the time, diffusion in solids was not considered generally possible.^[3] Today, Fick's laws form the core of our understanding of diffusion in solids, liquids, and gases (in the absence of bulk fluid motion in the latter two cases). When a diffusion process does not follow Fick's laws (which happens in cases of diffusion through porous media and diffusion of swelling penetrants, among others),^[4]^[5] it is referred to as non-Fickian.

Fick's first law

Fick's first law relates the diffusive flux to the Script error: No such module "anchor".gradient of the concentration. 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 can be written in various forms, where the most common form (see^[6]^[7]) is in a molar basis: $J = - D \frac{d φ}{d x},$ where

Template:Mvar is the diffusion flux, of which the dimension is the amount of substance per unit area per unit time. Template:Mvar measures the amount of substance that will flow through a unit area during a unit time interval,
Template:Mvar is the diffusion coefficient or diffusivity. Its dimension is area per unit time,
$\frac{d φ}{d x}$ is the concentration gradient,
Template:Mvar (for ideal mixtures) is the concentration, with a dimension of amount of substance per unit volume,
Template:Mvar is position, the dimension of which is length.

Template:Mvar 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. The modeling and prediction of Fick's diffusion coefficients is difficult. They can be estimated using the empirical Vignes correlation model^[8] or the physically-motivated entropy scaling.^[9] In dilute aqueous solutions the diffusion coefficients of most ions are similar and have values that at room temperature are in the range of Template:Val. For biological molecules the diffusion coefficients normally range from 10⁻¹⁰ to 10⁻¹¹ m²/s.

In two or more dimensions we must use Template:Math, the del or gradient operator, which generalises the first derivative, obtaining $𝐉 = - D \nabla φ,$ where Template:Math denotes the diffusion flux.

The driving force for the one-dimensional diffusion is the quantity Template:Math, which for ideal mixtures is the concentration gradient.

Variations of the first law

Another form for the first law is to write it with the primary variable as mass fraction (Template:Mvar, given for example in kg/kg), then the equation changes to $𝐉_{i} = - \frac{ρ D}{M_{i}} \nabla y_{i},$ where

the index Template:Mvar denotes the Template:Mvar^th species,
Template:Math is the diffusion flux of the Template:Mvar^th species (for example in mol/m²/s),
Template:Math is the molar mass of the Template:Mvar^th species,
Template:Mvar is the mixture density (for example in kg/m³).

The $ρ$ is outside the gradient operator. This is because $y_{i} = \frac{ρ_{s i}}{ρ},$ where Template:Mvar is the partial density of the Template:Mvarth species.

Beyond this, in chemical systems other than ideal solutions or mixtures, the driving force for the diffusion of each species is the gradient of chemical potential of this species. Then Fick's first law (one-dimensional case) can be written $J_{i} = - \frac{D c_{i}}{R T} \frac{\partial μ_{i}}{\partial x},$ where

the index Template:Mvar denotes the Template:Mvar^th species,
Template:Mvar is the concentration (mol/m³),
Template:Mvar is the universal gas constant (J/K/mol),
Template:Mvar is the absolute temperature (K),
Template:Mvar is the chemical potential (J/mol).

The driving force of Fick's law can be expressed as a fugacity difference: $J_{i} = - \frac{D}{R T} \frac{\partial f_{i}}{\partial x},$

where $f_{i}$ is the fugacity in Pa. $f_{i}$ is a partial pressure of component Template:Math in a vapor $f_{i}^{G}$ or liquid $f_{i}^{L}$ phase. At vapor liquid equilibrium the evaporation flux is zero because $f_{i}^{G} = f_{i}^{L}$ .

Derivation of Fick's first law for gases

Four versions of Fick's law for binary gas mixtures are given below. These assume: thermal diffusion is negligible; the body force per unit mass is the same on both species; and either pressure is constant or both species have the same molar mass. Under these conditions, Ref.^[10] shows in detail how the diffusion equation from the kinetic theory of gases reduces to this version of Fick's law: $V_{𝐢} = - D \nabla \ln y_{i},$ where Template:Math is the diffusion velocity of species Template:Mvar. In terms of species flux this is $J_{𝐢} = - \frac{ρ D}{M_{i}} \nabla y_{i} .$

If, additionally, $\nabla ρ = 0$ , this reduces to the most common form of Fick's law, $J_{𝐢} = - D \nabla φ .$

If (instead of or in addition to $\nabla ρ = 0$ ) both species have the same molar mass, Fick's law becomes $J_{𝐢} = - \frac{ρ D}{M_{i}} \nabla x_{i},$ where $x_{i}$ is the mole fraction of species Template:Mvar.

Fick's second law

Fick's second law predicts how diffusion causes the concentration to change with respect to time. It is a partial differential equation which in one dimension reads $\frac{\partial φ}{\partial t} = D \frac{\partial^{2} φ}{\partial x^{2}},$ where

Template:Mvar is the concentration in dimensions of $[N L^{- 3}]$ , example mol/m³; Template:Math is a function that depends on location Template:Mvar and time Template:Mvar,
Template:Mvar is time, example s,
Template:Mvar is the diffusion coefficient in dimensions of $[L^{2} T^{- 1}]$ , example m²/s,
Template:Mvar is the position, example m.

In two or more dimensions we must use the Laplacian Template:Math, which generalises the second derivative, obtaining the equation $\frac{\partial φ}{\partial t} = D Δ φ .$

Fick's second law has the same mathematical form as the Heat equation and its fundamental solution is the same as the Heat kernel, except switching thermal conductivity $k$ with diffusion coefficient $D$ : $φ (x, t) = \frac{1}{\sqrt{4 π D t}} \exp (- \frac{x^{2}}{4 D t}) .$

Derivation of Fick's second law

Fick's second law can be derived from Fick's first law and the mass conservation in absence of any chemical reactions: $\frac{\partial φ}{\partial t} + \frac{\partial}{\partial x} J = 0 \Rightarrow \frac{\partial φ}{\partial t} - \frac{\partial}{\partial x} (D \frac{\partial}{\partial x} φ) = 0 .$

Assuming the diffusion coefficient Template:Mvar to be a constant, one can exchange the orders of the differentiation and multiply by the constant: $\frac{\partial}{\partial x} (D \frac{\partial}{\partial x} φ) = D \frac{\partial}{\partial x} \frac{\partial}{\partial x} φ = D \frac{\partial^{2} φ}{\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 Fick's second law becomes $\frac{\partial φ}{\partial t} = D \nabla^{2} φ,$ 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 $\frac{\partial φ}{\partial t} = \nabla \cdot (D \nabla φ) .$

An important example is the case where Template:Math 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 Template:Mvar, the solution for the concentration will be a linear change of concentrations along Template:Mvar. In two or more dimensions we obtain $\nabla^{2} φ = 0,$ which is Laplace's equation, the solutions to which are referred to by mathematicians as harmonic functions.

Example solutions and generalization

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: $\frac{\partial φ}{\partial t} + \nabla \cdot 𝐣 = R,$ where Template:Math is the total flux and Template:Mvar is a net volumetric source for Template:Math. The only source of flux in this situation is assumed to be diffusive flux: $𝐣_{diffusion} = - D \nabla φ .$

Plugging the definition of diffusive flux to the continuity equation and assuming there is no source (Template:Math), we arrive at Fick's second law: $\frac{\partial φ}{\partial t} = D \frac{\partial^{2} φ}{\partial x^{2}} .$

If flux were the result of both diffusive flux and advective flux, the convection–diffusion equation is the result.

Example solution 1: constant concentration source and diffusion length

A simple case of diffusion with time Template:Mvar in one dimension (taken as the Template:Mvar-axis) from a boundary located at position Template:Math, where the concentration is maintained at a value Template:Math is $n (x, t) = n_{0} erfc (\frac{x}{2 \sqrt{D t}}),$ where Template:Math 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 – that is, the 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 Template:Math, Template:Math and that with Template:Math, Template:Math), then the solution is amended only with coefficient Template:Sfrac in front of Template:Math (as the diffusion now occurs in both directions). This case is valid when some solution with concentration Template:Math is put in contact with a layer of pure solvent. (Bokstein, 2005) The length $2 \sqrt{D t}$ is called the diffusion length and provides a measure of how far the concentration has propagated in the Template:Mvar-direction by diffusion in time Template:Mvar (Bird, 1976).

As a quick approximation of the error function, the first two terms of the Taylor series can be used: $n (x, t) = n_{0} [1 - 2 (\frac{x}{2 \sqrt{D t π}})] .$

If Template:Mvar is time-dependent, the diffusion length becomes $2 \sqrt{\int_{0}^{t} D (τ) d τ} .$ This idea is useful for estimating a diffusion length over a heating and cooling cycle, where Template:Mvar varies with temperature.

Example solution 2: Brownian particle and mean squared displacement

Another simple case of diffusion is the Brownian motion of one particle. The particle's Mean squared displacement from its original position is: $MSD \equiv ⟨ {(𝐱 - 𝐱_{0})}^{2} ⟩ = 2 n D t,$ where $n$ is the dimension of the particle's Brownian motion. For example, the diffusion of a molecule across a cell membrane 8 nm thick is 1-D diffusion because of the spherical symmetry; However, the diffusion of a molecule from the membrane to the center of a eukaryotic cell is a 3-D diffusion. For a cylindrical cactus, the diffusion from photosynthetic cells on its surface to its center (the axis of its cylindrical symmetry) is a 2-D diffusion.

The square root of MSD, $\sqrt{2 n D t}$ , is often used as a characterization of how far the particle has moved after time $t$ has elapsed. The MSD is symmetrically distributed over the 1D, 2D, and 3D space. Thus, the probability distribution of the magnitude of MSD in 1D is Gaussian and 3D is a Maxwell-Boltzmann distribution.

Generalizations

In non-homogeneous media, the diffusion coefficient varies in space, Template:Math. This dependence does not affect Fick's first law but the second law changes: $\begin{aligned} \frac{\partial φ (x, t)}{\partial t} & = \nabla \cdot (D (x) \nabla φ (x, t)) \\ = D (x) Δ φ (x, t) + \sum_{i = 1}^{3} \frac{\partial D (x)}{\partial x_{i}} \frac{\partial φ (x, t)}{\partial x_{i}} . \end{aligned}$
In anisotropic media, the diffusion coefficient depends on the direction. It is a symmetric tensor Template:Math. Fick's first law changes to $J = - D \nabla φ,$ it is the product of a tensor and a vector: $J_{i} = - \sum_{j = 1}^{3} D_{i j} \frac{\partial φ}{\partial x_{j}} .$ For the diffusion equation this formula gives $\frac{\partial φ (x, t)}{\partial t} = \nabla \cdot (D \nabla φ (x, t)) = \sum_{i = 1}^{3} \sum_{j = 1}^{3} D_{i j} \frac{\partial^{2} φ (x, t)}{\partial x_{i} \partial x_{j}} .$ The symmetric matrix of diffusion coefficients Template:Math should be positive definite. It is needed to make the right-hand side operator elliptic.
For inhomogeneous anisotropic media these two forms of the diffusion equation should be combined in $\begin{aligned} \frac{\partial φ (x, t)}{\partial t} & = \nabla \cdot (D (x) \nabla φ (x, t)) \\ = \sum_{i, j = 1}^{3} (D_{i j} (x) \frac{\partial^{2} φ (x, t)}{\partial x_{i} \partial x_{j}} + \frac{\partial D_{i j} (x)}{\partial x_{i}} \frac{\partial φ (x, t)}{\partial x_{i}}) . \end{aligned}$
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: $\frac{\partial φ_{i}}{\partial t} = \sum_{j} \nabla \cdot (D_{i j} \frac{φ_{i}}{φ_{j}} \nabla φ_{j}),$ where Template:Mvar are concentrations of the components and Template:Mvar is the matrix of coefficients. Here, indices Template:Mvar and Template:Mvar 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 Template:Math 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 Template:Math, where Template:Math refer to the components and Template:Math correspond to the space coordinates.

Applications

Equations based on Fick's law have been commonly used to model transport processes in foods, neurons, biopolymers, pharmaceuticals, porous soils, population dynamics, nuclear materials, plasma physics, and semiconductor doping processes. The theory of voltammetric methods is based on solutions of Fick's equation. On the other hand, in some cases a "Fickian (another common approximation of the transport equation is that of the diffusion theory)" description is inadequate. For example, in polymer science and food science a more general approach is required to describe transport of components in materials undergoing a glass transition. One more general framework is the Maxwell–Stefan diffusion equations^[11] of multi-component mass transfer, from which Fick's law can be obtained as a limiting case, 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).

Fick's flow in liquids

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 walks take 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.^[12]

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.

Sorption rate and collision frequency of diluted solute

Adsorption, absorption, and collision of molecules, particles, and surfaces are important problems in many fields. These fundamental processes regulate chemical, biological, and environmental reactions. Their rate can be calculated using the diffusion constant and Fick's laws of diffusion especially when these interactions happen in diluted solutions.

Typically, the diffusion constant of molecules and particles defined by Fick's equation can be calculated using the Stokes–Einstein equation. In the ultrashort time limit, in the order of the diffusion time a²/D, where a is the particle radius, the diffusion is described by the Langevin equation. At a longer time, the Langevin equation merges into the Stokes–Einstein equation. The latter is appropriate for the condition of the diluted solution, where long-range diffusion is considered. According to the fluctuation-dissipation theorem based on the Langevin equation in the long-time limit and when the particle is significantly denser than the surrounding fluid, the time-dependent diffusion constant is:^[13] $D (t) = μ k_{B} T (1 - e^{- t / (m μ)}),$ where (all in SI units)

k_B is the Boltzmann constant,
T is the absolute temperature,
μ is the mobility of the particle in the fluid or gas, which can be calculated using the Einstein relation (kinetic theory),
m is the mass of the particle,
t is time.

For a single molecule such as organic molecules or biomolecules (e.g. proteins) in water, the exponential term is negligible due to the small product of mμ in the ultrafast picosecond region, thus irrelevant to the relatively slower adsorption of diluted solute.

File:Diffusive sorption probability.png

Scheme of molecular diffusion in the solution. Orange dots are solute molecules, solvent molecules are not drawn, black arrow is an example random walk trajectory, and the red curve is the diffusive Gaussian broadening probability function from the Fick's law of diffusion.^[14]^{:Fig. 9}

The adsorption or absorption rate of a dilute solute to a surface or interface in a (gas or liquid) solution can be calculated using Fick's laws of diffusion. The accumulated number of molecules adsorbed on the surface is expressed by the Langmuir-Schaefer equation by integrating the diffusion flux equation over time as shown in the simulated molecular diffusion in the first section of this page:^[15] $Γ = 2 A C_{b} \sqrt{\frac{D t}{π}} .$

Template:Mvar is the surface area (m²).
$C_{b}$ is the number concentration of the adsorber molecules (solute) in the bulk solution (#/m³).
Template:Mvar is diffusion coefficient of the adsorber (m²/s).
Template:Mvar is elapsed time (s).
$Γ$ is the accumulated number of molecules in unit # molecules adsorbed during the time $t$ .

The equation is named after American chemists Irving Langmuir and Vincent Schaefer.

Briefly as explained in,^[16] the concentration gradient profile near a newly created (from $t = 0$ ) absorptive surface (placed at $x = 0$ ) in a once uniform bulk solution is solved in the above sections from Fick's equation, $\frac{\partial C}{\partial x} = \frac{C_{b}}{\sqrt{π D t}} exp (- \frac{x^{2}}{4 D t}),$ where Template:Mvar is the number concentration of adsorber molecules at $x, t$ (#/m³).

The concentration gradient at the subsurface at $x = 0$ is simplified to the pre-exponential factor of the distribution ${(\frac{\partial C}{\partial x})}_{x = 0} = \frac{C_{b}}{\sqrt{π D t}} .$ And the rate of diffusion (flux) across area $A .$ of the plane is ${(\frac{\partial Γ}{\partial t})}_{x = 0} = - \frac{D A C_{b}}{\sqrt{π D t}} .$ Integrating over time, $Γ = \int_{0}^{t} {(\frac{\partial Γ}{\partial t})}_{x = 0} = 2 A C_{b} \sqrt{\frac{D t}{π}} .$

The Langmuir–Schaefer equation can be extended to the Ward–Tordai Equation to account for the "back-diffusion" of rejected molecules from the surface:^[16] $Γ = 2 A C_{b} \sqrt{\frac{D t}{π}} - A \sqrt{\frac{D}{π}} \int_{0}^{\sqrt{t}} \frac{C (τ)}{\sqrt{t - τ}} d τ,$ where $C_{b}$ is the bulk concentration, $C$ is the sub-surface concentration (which is a function of time depending on the reaction model of the adsorption), and $τ$ is a dummy variable.

Monte Carlo simulations show that these two equations work to predict the adsorption rate of systems that form predictable concentration gradients near the surface but have troubles for systems without or with unpredictable concentration gradients, such as typical biosensing systems or when flow and convection are significant.^[17]

File:DiffusiveAdsorptionHistory.jpg

A brief history of the theories on diffusive adsorption.^[17]

A brief history of diffusive adsorption is shown in the right figure.^[17] A noticeable challenge of understanding the diffusive adsorption at the single-molecule level is the fractal nature of diffusion. Most computer simulations pick a time step for diffusion which ignores the fact that there are self-similar finer diffusion events (fractal) within each step. Simulating the fractal diffusion shows that a factor of two corrections should be introduced for the result of a fixed time-step adsorption simulation, bringing it to be consistent with the above two equations.^[17]

A more problematic result of the above equations is they predict the lower limit of adsorption under ideal situations but is very difficult to predict the actual adsorption rates. The equations are derived at the long-time-limit condition when a stable concentration gradient has been formed near the surface. But real adsorption is often done much faster than this infinite time limit i.e. the concentration gradient, decay of concentration at the sub-surface, is only partially formed before the surface has been saturated or flow is on to maintain a certain gradient, thus the adsorption rate measured is almost always faster than the equations have predicted for low or none energy barrier adsorption (unless there is a significant adsorption energy barrier that slows down the absorption significantly), for example, thousands to millions time faster in the self-assembly of monolayers at the water-air or water-substrate interfaces.^[15] As such, it is necessary to calculate the evolution of the concentration gradient near the surface and find out a proper time to stop the imagined infinite evolution for practical applications. While it is hard to predict when to stop but it is reasonably easy to calculate the shortest time that matters, the critical time when the first nearest neighbor from the substrate surface feels the building-up of the concentration gradient. This yields the upper limit of the adsorption rate under an ideal situation when there are no other factors than diffusion that affect the absorber dynamics:^[17] $⟨ r ⟩ = \frac{4}{π} A C_{b}^{4 / 3} D,$ where:

$⟨ r ⟩$ is the adsorption rate assuming under adsorption energy barrier-free situation, in unit #/s,
$A$ is the area of the surface of interest on an "infinite and flat" substrate (m²),
$C_{b}$ is the concentration of the absorber molecule in the bulk solution (#/m³),
$D$ is the diffusion constant of the absorber (solute) in the solution (m²/s) defined with Fick's law.

This equation can be used to predict the initial adsorption rate of any system; It can be used to predict the steady-state adsorption rate of a typical biosensing system when the binding site is just a very small fraction of the substrate surface and a near-surface concentration gradient is never formed; It can also be used to predict the adsorption rate of molecules on the surface when there is a significant flow to push the concentration gradient very shallowly in the sub-surface.

This critical time is significantly different from the first passenger arriving time or the mean free-path time. Using the average first-passenger time and Fick's law of diffusion to estimate the average binding rate will significantly over-estimate the concentration gradient because the first passenger usually comes from many layers of neighbors away from the target, thus its arriving time is significantly longer than the nearest neighbor diffusion time. Using the mean free path time plus the Langmuir equation will cause an artificial concentration gradient between the initial location of the first passenger and the target surface because the other neighbor layers have no change yet, thus significantly lower estimate the actual binding time, i.e., the actual first passenger arriving time itself, the inverse of the above rate, is difficult to calculate. If the system can be simplified to 1D diffusion, then the average first passenger time can be calculated using the same nearest neighbor critical diffusion time for the first neighbor distance to be the MSD,^[18] $L = \sqrt{2 D t},$ where:

$L = C_{b}^{- 1 / 3}$ (unit m) is the average nearest neighbor distance approximated as cubic packing, where $C_{b}$ is the solute concentration in the bulk solution (unit # molecule / m³),
$D$ is the diffusion coefficient defined by Fick's equation (unit m²/s),
$t$ is the critical time (unit s).

In this critical time, it is unlikely the first passenger has arrived and adsorbed. But it sets the speed of the layers of neighbors to arrive. At this speed with a concentration gradient that stops around the first neighbor layer, the gradient does not project virtually in the longer time when the actual first passenger arrives. Thus, the average first passenger coming rate (unit # molecule/s) for this 3D diffusion simplified in 1D problem, $⟨ r ⟩ = \frac{a}{t} = 2 a C_{b}^{2 / 3} D,$ where $a$ is a factor of converting the 3D diffusive adsorption problem into a 1D diffusion problem whose value depends on the system, e.g., a fraction of adsorption area $A$ over solute nearest neighbor sphere surface area $4 π L^{2} / 4$ assuming cubic packing each unit has 8 neighbors shared with other units. This example fraction converges the result to the 3D diffusive adsorption solution shown above with a slight difference in pre-factor due to different packing assumptions and ignoring other neighbors.

When the area of interest is the size of a molecule (specifically, a long cylindrical molecule such as DNA), the adsorption rate equation represents the collision frequency of two molecules in a diluted solution, with one molecule a specific side and the other no steric dependence, i.e., a molecule (random orientation) hit one side of the other. The diffusion constant need to be updated to the relative diffusion constant between two diffusing molecules. This estimation is especially useful in studying the interaction between a small molecule and a larger molecule such as a protein. The effective diffusion constant is dominated by the smaller one whose diffusion constant can be used instead.

The above hitting rate equation is also useful to predict the kinetics of molecular self-assembly on a surface. Molecules are randomly oriented in the bulk solution. Assuming 1/6 of the molecules has the right orientation to the surface binding sites, i.e. 1/2 of the z-direction in x, y, z three dimensions, thus the concentration of interest is just 1/6 of the bulk concentration. Put this value into the equation one should be able to calculate the theoretical adsorption kinetic curve using the Langmuir adsorption model. In a more rigid picture, 1/6 can be replaced by the steric factor of the binding geometry.

File:JChen2022JPCA.png

Comparing collision theory and diffusive collision theory.^[19]

The bimolecular collision frequency related to many reactions including protein coagulation/aggregation is initially described by Smoluchowski coagulation equation proposed by Marian Smoluchowski in a seminal 1916 publication,^[20] derived from Brownian motion and Fick's laws of diffusion. Under an idealized reaction condition for A + B → product in a diluted solution, Smoluchovski suggested that the molecular flux at the infinite time limit can be calculated from Fick's laws of diffusion yielding a fixed/stable concentration gradient from the target molecule, e.g. B is the target molecule holding fixed relatively, and A is the moving molecule that creates a concentration gradient near the target molecule B due to the coagulation reaction between A and B. Smoluchowski calculated the collision frequency between A and B in the solution with unit #/s/m³: $Z_{A B} = 4 π R D_{r} C_{A} C_{B},$ where:

$R$ is the radius of the collision,
$D_{r} = D_{A} + D_{B}$ is the relative diffusion constant between A and B (m²/s),
$C_{A}$ and $C_{B}$ are number concentrations of A and B respectively (#/m³).

The reaction order of this bimolecular reaction is 2 which is the analogy to the result from collision theory by replacing the moving speed of the molecule with diffusive flux. In the collision theory, the traveling time between A and B is proportional to the distance which is a similar relationship for the diffusion case if the flux is fixed.

However, under a practical condition, the concentration gradient near the target molecule is evolving over time with the molecular flux evolving as well,^[17] and on average the flux is much bigger than the infinite time limit flux Smoluchowski has proposed. Before the first passenger arrival time, Fick's equation predicts a concentration gradient over time which does not build up yet in reality. Thus, this Smoluchowski frequency represents the lower limit of the real collision frequency.

In 2022, Chen calculates the upper limit of the collision frequency between A and B in a solution assuming the bulk concentration of the moving molecule is fixed after the first nearest neighbor of the target molecule.^[19] Thus the concentration gradient evolution stops at the first nearest neighbor layer given a stop-time to calculate the actual flux. He named this the critical time and derived the diffusive collision frequency in unit #/s/m³:^[19] $Z_{A B} = \frac{8}{π} σ D_{r} C_{A} C_{B} \sqrt[3]{C_{A} + C_{B}},$ where:

$σ$ is the area of the cross-section of the collision (m²),
$D_{r} = D_{A} + D_{B}$ is the relative diffusion constant between A and B (m²/s),
$C_{A}$ and $C_{B}$ are number concentrations of A and B respectively (#/m³),
$\sqrt[3]{C_{A} + C_{B}}$ represents 1/⟨d⟩, where d is the average distance between two molecules.

This equation assumes the upper limit of a diffusive collision frequency between A and B is when the first neighbor layer starts to feel the evolution of the concentration gradient, whose reaction order is Template:Sfrac instead of 2. Both the Smoluchowski equation and the JChen equation satisfy dimensional checks with SI units. But the former is dependent on the radius and the latter is on the area of the collision sphere. From dimensional analysis, there should be an equation dependent on the volume of the collision sphere,^[21] e.g.,

$Z_{A B} = 4 V D_{r} C_{A} C_{B} (C_{A} + C_{B})^{2 / 3},$

V is the volume of the collision sphere

but eventually, all equations should converge to the same numerical rate of the collision that can be measured experimentally. The actual reaction order for a bimolecular unit reaction could be between 2 and Template:Sfrac, which makes sense because the diffusive collision time is squarely dependent on the distance between the two molecules.

These new equations also avoid the singularity on the adsorption rate at time zero for the Langmuir-Schaefer equation. The infinity rate is justifiable under ideal conditions because when you introduce target molecules magically in a solution of probe molecule or vice versa, there always be a probability of them overlapping at time zero, thus the rate of that two molecules association is infinity. It does not matter that other millions of molecules have to wait for their first mate to diffuse and arrive. The average rate is thus infinity. But statistically this argument is meaningless. The maximum rate of a molecule in a period of time larger than zero is 1, either meet or not, thus the infinite rate at time zero for that molecule pair really should just be one, making the average rate 1/millions or more and statistically negligible. This does not even count in reality no two molecules can magically meet at time zero.

Biological perspective

The first law gives rise to the following formula:^[22] $flux = - P (c_{2} - c_{1}),$ where

Template:Mvar is the permeability, an experimentally determined membrane "conductance" for a given gas at a given temperature,
Template:Math is the difference in concentration of the gas across the membrane for the direction of flow (from Template:Math to Template:Math).

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.

Under the condition of a diluted solution when diffusion takes control, the membrane permeability mentioned in the above section can be theoretically calculated for the solute using the equation mentioned in the last section (use with particular care because the equation is derived for dense solutes, while biological molecules are not denser than water. Also, this equation assumes ideal concentration gradient forms near the membrane and evolves):^[14] $P = 2 A_{p} η_{t m} \sqrt{\frac{D}{π t}},$ where:

$A_{P}$ is the total area of the pores on the membrane (unit m²),
$η_{t m}$ transmembrane efficiency (unitless), which can be calculated from the stochastic theory of chromatography,
D is the diffusion constant of the solute unit m²⋅s⁻¹,
t is time unit s,
c₂, c₁ concentration should use unit mol m⁻³, so flux unit becomes mol s⁻¹.

The flux is decay over the square root of time because a concentration gradient builds up near the membrane over time under ideal conditions. When there is flow and convection, the flux can be significantly different than the equation predicts and show an effective time t with a fixed value,^[17] which makes the flux stable instead of decay over time. A critical time has been estimated under idealized flow conditions when there is no gradient formed.^[17]^[19] This strategy is adopted in biology such as blood circulation.

Semiconductor fabrication applications

The semiconductor is a collective term for a series of devices. It mainly includes three categories：two-terminal devices, three-terminal devices, and four-terminal devices. The combination of the semiconductors is called an integrated circuit.

The relationship between Fick's law and semiconductors: the principle of the semiconductor is transferring chemicals or dopants from a layer to a layer. Fick's law can be used to control and predict the diffusion by knowing how much the concentration of the dopants or chemicals move per meter and second through mathematics.

Therefore, different types and levels of semiconductors can be fabricated.

Integrated circuit fabrication technologies, model processes like CVD, thermal oxidation, wet oxidation, doping, etc. use diffusion equations obtained from Fick's law.

CVD method of semiconductor fabrication

The wafer is a kind of semiconductor whose silicon substrate is coated with a layer of CVD-created polymer chain and films. This film contains n-type and p-type dopants and takes responsibility for dopant conductions. The principle of CVD relies on the gas phase and gas-solid chemical reaction to create thin films.

The viscous flow regime of CVD is driven by a pressure gradient. CVD also includes a diffusion component distinct from the surface diffusion of adatoms. In CVD, reactants and products must also diffuse through a boundary layer of stagnant gas that exists next to the substrate. The total number of steps required for CVD film growth are gas phase diffusion of reactants through the boundary layer, adsorption and surface diffusion of adatoms, reactions on the substrate, and gas phase diffusion of products away through the boundary layer.

The velocity profile for gas flow is: $δ (x) = (\frac{5 x}{{R e}^{1 / 2}}) R e = \frac{v ρ L}{η},$ where:

$δ$ is the thickness,
$R e$ is the Reynolds number,
Template:Mvar is the length of the substrate,
Template:Math at any surface,
$η$ is viscosity,
$ρ$ is density.

Integrated the Template:Mvar from Template:Math to Template:Mvar, it gives the average thickness: $δ = \frac{10 L}{3 {R e}^{1 / 2}} .$

To keep the reaction balanced, reactants must diffuse through the stagnant boundary layer to reach the substrate. So a thin boundary layer is desirable. According to the equations, increasing vo would result in more wasted reactants. The reactants will not reach the substrate uniformly if the flow becomes turbulent. Another option is to switch to a new carrier gas with lower viscosity or density.

The Fick's first law describes diffusion through the boundary layer. As a function of pressure (p) and temperature (T) in a gas, diffusion is determined.

$D = D_{0} (\frac{p_{0}}{p}) {(\frac{T}{T_{0}})}^{3 / 2},$ where:

$p_{0}$ is the standard pressure,
$T_{0}$ is the standard temperature,
$D_{0}$ is the standard diffusitivity.

The equation tells that increasing the temperature or decreasing the pressure can increase the diffusivity.

Fick's first law predicts the flux of the reactants to the substrate and product away from the substrate: $J = - D_{i} (\frac{d c_{i}}{d x}),$ where:

$x$ is the thickness $δ$ ,
$d c_{i}$ is the first reactant's concentration.

In ideal gas law $p V = n R T$ , the concentration of the gas is expressed by partial pressure.

$J = - D_{i} (\frac{p_{i} - p_{0}}{δ R T}),$ where

$R$ is the gas constant,
$\frac{p_{i} - p_{0}}{δ}$ is the partial pressure gradient.

As a result, Fick's first law tells us we can use a partial pressure gradient to control the diffusivity and control the growth of thin films of semiconductors.

In many realistic situations, the simple Fick's law is not an adequate formulation for the semiconductor problem. It only applies to certain conditions, for example, given the semiconductor boundary conditions: constant source concentration diffusion, limited source concentration, or moving boundary diffusion (where junction depth keeps moving into the substrate).

Invalidity of Fickian diffusion

Even though Fickian diffusion has been used to model diffusion processes in semiconductor manufacturing (including CVD reactors) in early days, it often fails to validate the diffusion in advanced semiconductor nodes (< 90 nm). This mostly stems from the inability of Fickian diffusion to model diffusion processes accurately at molecular level and smaller. In advanced semiconductor manufacturing, it is important to understand the movement at atomic scales, which is failed by continuum diffusion. Today, most semiconductor manufacturers use random walk to study and model diffusion processes. This allows us to study the effects of diffusion in a discrete manner to understand the movement of individual atoms, molecules, plasma etc.

In such a process, the movements of diffusing species (atoms, molecules, plasma etc.) are treated as a discrete entity, following a random walk through the CVD reactor, boundary layer, material structures etc. Sometimes, the movements might follow a biased-random walk depending on the processing conditions. Statistical analysis is done to understand variation/stochasticity arising from the random walk of the species, which in-turn affects the overall process and electrical variations.

Food production and cooking

The formulation of Fick's first law can explain a variety of complex phenomena in the context of food and cooking: Diffusion of molecules such as ethylene promotes plant growth and ripening, salt and sugar molecules promotes meat brining and marinating, and water molecules promote dehydration. Fick's first law can also be used to predict the changing moisture profiles across a spaghetti noodle as it hydrates during cooking. These phenomena are all about the spontaneous movement of particles of solutes driven by the concentration gradient. In different situations, there is different diffusivity which is a constant.^[23]

By controlling the concentration gradient, the cooking time, shape of the food, and salting can be controlled.

Citations

Template:Reflist

External links

Fick's equations, Boltzmann's transformation, etc. (with figures and animations)
Fick's Second Law on OpenStax

de:Diffusion#Erstes Fick'sches Gesetz

↑ Script error: No such module "citation/CS1".
↑
* Script error: No such module "Citation/CS1".
- Script error: No such module "Citation/CS1".
↑ Script error: No such module "Citation/CS1".
↑ Script error: No such module "citation/CS1".
↑ Script error: No such module "Citation/CS1".
↑ Script error: No such module "citation/CS1".
↑ Script error: No such module "citation/CS1".
↑ Script error: No such module "Citation/CS1".
↑ Script error: No such module "Citation/CS1".
↑ Script error: No such module "citation/CS1".
↑ Script error: No such module "citation/CS1".Template:Pn
↑ Script error: No such module "Citation/CS1".
↑ Script error: No such module "Citation/CS1".
↑ ^a ^b Script error: No such module "Citation/CS1".
↑ ^a ^b Script error: No such module "Citation/CS1".
↑ ^a ^b Script error: No such module "Citation/CS1".
↑ ^a ^b ^c ^d ^e ^f ^g ^h Cite error: Invalid <ref> tag; no text was provided for refs named JixinMCSimuAdsorption
↑ Script error: No such module "Citation/CS1".
↑ ^a ^b ^c ^d Script error: No such module "Citation/CS1".
↑ Script error: No such module "Citation/CS1".
↑ Script error: No such module "Citation/CS1".
↑ Script error: No such module "citation/CS1".
↑ Script error: No such module "Citation/CS1".

[1] Script error: No such module "citation/CS1".

[:0-2] * Script error: No such module "Citation/CS1".
Script error: No such module "Citation/CS1".

[3] Script error: No such module "Citation/CS1".

[3] Script error: No such module "Citation/CS1".

[4] Script error: No such module "citation/CS1".

[GorbanMMNP2011-5] Script error: No such module "Citation/CS1".

[6] Script error: No such module "citation/CS1".

[7] Script error: No such module "citation/CS1".

[8] Script error: No such module "Citation/CS1".

[9] Script error: No such module "Citation/CS1".

[10] Script error: No such module "citation/CS1".

[11] Script error: No such module "citation/CS1".Template:Pn

[12] Script error: No such module "Citation/CS1".

[13] Script error: No such module "Citation/CS1".

[Pyle-BJNano-14] Script error: No such module "Citation/CS1".

[LangmuirSchaefer1937JACS-15] Script error: No such module "Citation/CS1".

[WardTordai1946-16] Script error: No such module "Citation/CS1".

[JixinMCSimuAdsorption-17] ↑ ^a ^b ^c ^d ^e ^f ^g ^h Cite error: Invalid <ref> tag; no text was provided for refs named JixinMCSimuAdsorption

[Pandey-JPCB2024-18] Script error: No such module "Citation/CS1".

[JChen2022JPCA-19] Script error: No such module "Citation/CS1".

[Smoluchowski1916-20] Script error: No such module "Citation/CS1".

[Chen-AIP2024-21] Script error: No such module "Citation/CS1".

[22] Script error: No such module "citation/CS1".

[23] Script error: No such module "Citation/CS1".

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[19]

[20]

[21]

[22]

[23]

Fick's laws of diffusion: Difference between revisions

Latest revision as of 21:47, 13 November 2025

Contents

History

Fick's first law

Variations of the first law

Derivation of Fick's first law for gases

Fick's second law

Derivation of Fick's second law

Example solutions and generalization

Example solution 1: constant concentration source and diffusion length

Example solution 2: Brownian particle and mean squared displacement

Generalizations

Applications

Fick's flow in liquids

Sorption rate and collision frequency of diluted solute

Biological perspective

Semiconductor fabrication applications

CVD method of semiconductor fabrication

Invalidity of Fickian diffusion

Food production and cooking

See also

Citations

Further reading

External links

Navigation menu

@@ Line 2: / Line 2: @@
 {{for|the technique of measuring [[cardiac output]]|Fick principle}}
 {{Use dmy dates|date=July 2020}}
-[[File:DiffusionMicroMacro.gif|thumb|right|250px|[[Molecular diffusion]] from a microscopic and macroscopic point of view. Initially, there are [[solute]] molecules on the left side of a barrier (purple line) and none on the right. The barrier is removed, and the solute diffuses to fill the whole container. '''Top''': A single molecule moves around randomly. '''Middle''': With more molecules, there is a clear trend where the solute fills the container more and more uniformly. '''Bottom''': With an enormous number of solute molecules, randomness becomes undetectable: The solute appears to move smoothly and systematically from high-concentration areas to low-concentration areas. This smooth flow is described by Fick's laws.]]
+[[File:DiffusionMicroMacro.gif|thumb|[[Molecular diffusion]] from a microscopic and macroscopic point of view. Initially, there are [[solute]] molecules on the left side of a barrier (purple line) and none on the right. The barrier is removed, and the solute diffuses to fill the whole container. '''Top''': A single molecule moves around randomly. '''Middle''': With more molecules, there is a clear trend where the solute fills the container more and more uniformly. '''Bottom''': With an enormous number of solute molecules, randomness becomes undetectable: The solute appears to move smoothly and systematically from high-concentration areas to low-concentration areas. This smooth flow is described by Fick's laws.]]
 '''Fick's laws of diffusion''' describe [[diffusion]] and were first posited by [[Adolf Fick]] in 1855 on the basis of largely experimental results. They can be used to solve for the [[Mass diffusivity|diffusion coefficient]], {{mvar|D}}. Fick's first law can be used to derive his second law which in turn is identical to the [[diffusion equation]].
@@ Line 20: / Line 20: @@
 == Fick's first law ==
 '''Fick's first law''' relates the diffusive [[flux]] to the {{anchor|concentration gradient}}gradient of the concentration. 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 can be written in various forms, where the most common form (see<ref>{{cite book | vauthors = Atkins P, de Paula J |year=2006 |title=Physical Chemistry for the Life Science |url= https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Supplemental_Modules_(Physical_and_Theoretical_Chemistry)/Kinetics/Diffusion#Fick.E2.80.99s_First_Law_of_Diffusion }}</ref><ref>{{cite book |doi=10.1017/CBO9781139025614 |title=Essentials of Micro- and Nanofluidics |date=2013 |last1=Conlisk |first1=A. Terrence |isbn=978-0-521-88168-5 |page=43 }}</ref>) is in a molar basis:
-: <math>J = -D \frac{d \varphi}{d x}, </math>
+<math display="block">J = -D \frac{d \varphi}{d x}, </math>
 where
 * {{mvar|J}} is the '''diffusion flux''', of which the [[dimensional analysis|dimension]] is the [[amount of substance]] per unit area per unit time. {{mvar|J}} measures the amount of substance that will flow through a unit area during a unit time interval,
@@ Line 31: / Line 31: @@
 In two or more dimensions we must use {{math|∇}}, the [[del]] or [[gradient]] operator, which generalises the first derivative, obtaining
-: <math> \mathbf{J}=- D\nabla \varphi , </math>
+<math display="block"> \mathbf{J}=- D\nabla \varphi , </math>
 where {{math|'''J'''}} denotes the diffusion flux.
@@ Line 39: / Line 39: @@
 Another form for the first law is to write it with the primary variable as [[Mass fraction (chemistry)|mass fraction]] ({{mvar|y<sub>i</sub>}}, given for example in kg/kg), then the equation changes to
-: <math>\mathbf{J}_i = -\frac{\rho D}{M_i}\nabla y_i , </math>
+<math display="block">\mathbf{J}_i = -\frac{\rho D}{M_i}\nabla y_i , </math>
 where
 * the index {{mvar|i}} denotes the {{mvar|i}}<sup>th</sup> species,
@@ Line 47: / Line 47: @@
 The <math>\rho</math> is outside the [[gradient]] operator. This is because
-: <math>y_i = \frac{\rho_{si}}{\rho} , </math>
+<math display="block">y_i = \frac{\rho_{si}}{\rho} , </math>
 where {{mvar|ρ<sub>si</sub>}} is the partial density of the {{mvar|i}}th species.
 Beyond this, in chemical systems other than ideal solutions or mixtures, the driving force for the diffusion of each species is the gradient of [[chemical potential]] of this species. Then Fick's first law (one-dimensional case) can be written
-: <math>J_i = - \frac{D c_i}{RT} \frac{\partial \mu_i}{\partial x} , </math>
+<math display="block">J_i = - \frac{D c_i}{RT} \frac{\partial \mu_i}{\partial x} , </math>
 where
 * the index {{mvar|i}} denotes the {{mvar|i}}<sup>th</sup> species,
@@ Line 60: / Line 60: @@
 The driving force of Fick's law can be expressed as a [[fugacity]] difference:
-: <math>J_i = - \frac{D}{RT} \frac{\partial f_i}{\partial x} , </math>
+<math display="block">J_i = - \frac{D}{RT} \frac{\partial f_i}{\partial x} , </math>
 where <math> f_i </math> is the fugacity in Pa. <math> f_i </math> is a partial pressure of component {{math|''i''}} in a vapor <math> f_i^\text{G} </math> or liquid <math> f_i^\text{L} </math> phase. At vapor liquid equilibrium the evaporation flux is zero because <math> f_i^\text{G} = f_i^\text{L} </math>.
@@ Line 67: / Line 67: @@
 Four versions of Fick's law for binary gas mixtures are given below. These assume: thermal diffusion is negligible; the body force per unit mass is the same on both species; and either pressure is constant or both species have the same molar mass. Under these conditions, Ref.<ref>{{cite book| vauthors = Williams FA |year=1985 |chapter=Appendix E |title=Combustion Theory |publisher=Benjamin/Cummings}}</ref> shows in detail how the diffusion equation from the [[kinetic theory of gases]] reduces to this version of Fick's law:
-<math display="block"> \mathbf{V_i}=- D\nabla \ln{\left( y_i \right)} , </math>
+<math display="block"> \mathbf{V_i}=- D \, \nabla \ln y_i , </math>
 where {{math|'''V<sub>i</sub>'''}} is the diffusion velocity of species {{mvar|i}}. In terms of species flux this is
 <math display="block">\mathbf{J_i}=- \frac{\rho D}{M_i}\nabla y_i . </math>
@@ Line 80: / Line 80: @@
 == Fick's second law ==
 '''Fick's second law''' predicts how diffusion causes the concentration to change with respect to time. It is a [[partial differential equation]] which in one dimension reads
-: <math>\frac{\partial \varphi}{\partial t} = D\,\frac{\partial^2 \varphi}{\partial x^2},</math>
+<math display="block">\frac{\partial \varphi}{\partial t} = D\,\frac{\partial^2 \varphi}{\partial x^2},</math>
 where
 * {{mvar|φ}} is the concentration in dimensions of <math>[\mathsf{N}\mathsf{L}^{-3}]</math>, example mol/m<sup>3</sup>; {{math|1=''φ'' = ''φ''(''x'',''t'')}} is a function that depends on location {{mvar|x}} and time {{mvar|t}},
@@ Line 88: / Line 88: @@
 In two or more dimensions we must use the [[Laplacian]] {{math|1=Δ = ∇<sup>2</sup>}}, which generalises the second derivative, obtaining the equation
-: <math>\frac{\partial \varphi}{\partial t} = D\Delta \varphi . </math>
+<math display="block">\frac{\partial \varphi}{\partial t} = D\Delta \varphi . </math>
 Fick's second law has the same mathematical form as the [[Heat equation]] and its [[fundamental solution]] is the same as the [[Heat kernel]]<!-- This is true for the case of an initial Gaussian distribution. Other problem geometries will lead to different solutions. (e.g. diffusion with a fixed boundary concentration, inter-penetration of two solids, etc) -->, except switching thermal conductivity <math>k</math> with diffusion coefficient <math>D</math>:
@@ Line 101: / Line 101: @@
 At a given time step, half of the particles would move left and half would move right. Since half of the particles at point {{mvar|x}} move right and half of the particles at point {{math|''x'' + Δ''x''}} move left, the net movement to the right is:
-: <math>-\tfrac{1}{2}\bigl[N(x + \Delta x, t) - N(x, t)\bigr]</math>
+<math display="block">-\tfrac{1}{2}\bigl[N(x + \Delta x, t) - N(x, t)\bigr]</math>
 The flux, {{mvar|J}}, is this net movement of particles across some area element of area {{mvar|a}}, normal to the random walk during a time interval {{math|Δ''t''}}. Hence we may write:
-: <math>J = - \frac{1}{2} \left[\frac{ N(x + \Delta x, t)}{a \Delta t} - \frac{ N(x, t)}{a \Delta t}\right]</math>
+<math display="block">J = - \frac{1}{2} \left[\frac{ N(x + \Delta x, t)}{a \Delta t} - \frac{ N(x, t)}{a \Delta t}\right]</math>
 Multiplying the top and bottom of the right hand side by {{math|(Δ''x'')<sup>2</sup>}} and rewriting, one obtains:
-: <math> J = -\frac{\left(\Delta x\right)^2}{2 \Delta t}\left[\frac{N(x + \Delta x, t)}{a \left(\Delta x\right)^2} - \frac{N(x, t)}{a \left(\Delta x\right)^2}\right]</math>
+<math display="block"> J = -\frac{\left(\Delta x\right)^2}{2 \Delta t}\left[\frac{N(x + \Delta x, t)}{a \left(\Delta x\right)^2} - \frac{N(x, t)}{a \left(\Delta x\right)^2}\right]</math>
 Concentration is defined as particles per unit volume, and hence
-: <math>\varphi (x, t) = \frac{N(x, t)}{a \Delta x}.</math>
+<math display="block">\varphi (x, t) = \frac{N(x, t)}{a \Delta x}.</math>
 In addition, {{math|{{sfrac|(Δ''x'')<sup>2</sup>|2Δ''t''}}}} is the definition of the one-dimensional diffusion constant, {{mvar|D}}. Thus our expression simplifies to:
-: <math> J = -D \left[\frac{\varphi (x + \Delta x, t)}{\Delta x} - \frac{\varphi (x , t)}{\Delta x}\right]</math>
+<math display="block"> J = -D \left[\frac{\varphi (x + \Delta x, t)}{\Delta x} - \frac{\varphi (x , t)}{\Delta x}\right]</math>
 In the limit where {{math|Δ''x''}} is infinitesimal, the right-hand side becomes a space derivative:
-: <math> J = - D \frac{\partial \varphi}{\partial x}  </math> This is only the case for the initial condition of a very of a initial gaussian distribution. Other problem geometries will lead to different solutions
+<math display="block"> J = - D \frac{\partial \varphi}{\partial x}  </math> This is only the case for the initial condition of a very of a initial gaussian distribution. Other problem geometries will lead to different solutions
 -->
 Fick's second law can be derived from Fick's first law and the [[mass conservation]] in absence of any chemical reactions:
-: <math>\frac{\partial \varphi}{\partial t} + \frac{\partial}{\partial x}J = 0
+<math display="block">\frac{\partial \varphi}{\partial t} + \frac{\partial}{\partial x}J = 0
 \Rightarrow\frac{\partial \varphi}{\partial t} -\frac{\partial}{\partial x}\left(D\frac{\partial}{\partial x}\varphi\right)\,=0.</math>
 Assuming the diffusion coefficient {{mvar|D}} to be a constant, one can exchange the orders of the differentiation and multiply by the constant:
-:<math>\frac{\partial}{\partial x}\left(D\frac{\partial}{\partial x} \varphi\right) = D\frac{\partial}{\partial x} \frac{\partial}{\partial x} \varphi = D\frac{\partial^2\varphi}{\partial x^2},</math>
+<math display="block">\frac{\partial}{\partial x}\left(D\frac{\partial}{\partial x} \varphi\right) = D\frac{\partial}{\partial x} \frac{\partial}{\partial x} \varphi = D\frac{\partial^2\varphi}{\partial x^2},</math>
 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
-: <math>\frac{\partial \varphi}{\partial t} = D\,\nabla^2\varphi,</math>
+<math display="block">\frac{\partial \varphi}{\partial t} = D\,\nabla^2\varphi,</math>
 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
-: <math>\frac{\partial \varphi}{\partial t} =  \nabla \cdot (D\,\nabla\varphi).</math>
+<math display="block">\frac{\partial \varphi}{\partial t} =  \nabla \cdot (D\,\nabla\varphi).</math>
 An important example is the case where {{math|'''φ'''}} 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 {{mvar|D}}, the solution for the concentration will be a linear change of concentrations along {{mvar|x}}. In two or more dimensions we obtain
-: <math> \nabla^2\varphi = 0,</math>
+<math display="block"> \nabla^2\varphi = 0,</math>
 which is [[Laplace's equation]], the solutions to which are referred to by mathematicians as [[harmonic functions]].
 == Example solutions and generalization ==
 Fick's second law is a special case of the [[convection–diffusion equation]] in which there is no [[advection|advective flux]] and no net volumetric source. It can be derived from the [[Continuity equation#Differential form|continuity equation]]:
-: <math> \frac{\partial \varphi}{\partial t} + \nabla\cdot\mathbf{j} = R, </math>
+<math display="block"> \frac{\partial \varphi}{\partial t} + \nabla\cdot\mathbf{j} = R, </math>
 where {{math|'''j'''}} is the total [[flux]] and {{mvar|R}} is a net volumetric source for {{math|'''φ'''}}. The only source of flux in this situation is assumed to be ''diffusive flux'':
-: <math>\mathbf{j}_{\text{diffusion}} = -D \nabla \varphi . </math>
+<math display="block">\mathbf{j}_{\text{diffusion}} = -D \nabla \varphi . </math>
 Plugging the definition of diffusive flux to the continuity equation and assuming there is no source ({{math|1=''R'' = 0}}), we arrive at Fick's second law:
-: <math>\frac{\partial \varphi}{\partial t} = D\frac{\partial^2 \varphi}{\partial x^2} . </math>
+<math display="block">\frac{\partial \varphi}{\partial t} = D\frac{\partial^2 \varphi}{\partial x^2} . </math>
 If flux were the result of both diffusive flux and [[advection|advective flux]], the [[convection–diffusion equation]] is the result.
@@ Line 150: / Line 150: @@
 === Example solution 1: constant concentration source and diffusion length ===
 A simple case of diffusion with time {{mvar|t}} in one dimension (taken as the {{mvar|x}}-axis) from a boundary located at position {{math|1=''x'' = 0}}, where the concentration is maintained at a value {{math|''n''<sub>0</sub>}} is
-: <math>n \left(x,t \right)=n_0 \operatorname{erfc} \left( \frac{x}{2\sqrt{Dt}}\right) ,</math>
+<math display="block">n \left(x,t \right)=n_0 \operatorname{erfc} \left( \frac{x}{2\sqrt{Dt}}\right) ,</math>
-where {{math|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 – that is, the 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 {{math|1=''n''(''x'', 0) = 0}}, {{math|''x'' > 0}} and that with {{math|1=''n''(''x'', 0) = ''n''<sub>0</sub>}}, {{math|''x'' ≤ 0}}), then the solution is amended only with coefficient {{sfrac|2}} in front of {{math|''n''<sub>0</sub>}} (as the diffusion now occurs in both directions). This case is valid when some solution with concentration {{math|''n''<sub>0</sub>}} is put in contact with a layer of pure solvent. (Bokstein, 2005) The length {{math|2{{sqrt|''Dt''}}}} is called the ''diffusion length'' and provides a measure of how far the concentration has propagated in the {{mvar|x}}-direction by diffusion in time {{mvar|t}} (Bird, 1976).
+where {{math|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 – that is, the 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 {{math|1=''n''(''x'', 0) = 0}}, {{math|''x'' > 0}} and that with {{math|1=''n''(''x'', 0) = ''n''<sub>0</sub>}}, {{math|''x'' ≤ 0}}), then the solution is amended only with coefficient {{sfrac|2}} in front of {{math|''n''<sub>0</sub>}} (as the diffusion now occurs in both directions). This case is valid when some solution with concentration {{math|''n''<sub>0</sub>}} is put in contact with a layer of pure solvent. (Bokstein, 2005) The length <math>2\sqrt{Dt}</math> is called the ''diffusion length'' and provides a measure of how far the concentration has propagated in the {{mvar|x}}-direction by diffusion in time {{mvar|t}} (Bird, 1976).
 As a quick approximation of the error function, the first two terms of the [[Taylor series]] can be used:
-: <math>n(x,t)=n_0 \left[ 1 - 2 \left(\frac{x}{2\sqrt{Dt\pi}}\right) \right] . </math>
+<math display="block">n(x,t)=n_0 \left[ 1 - 2 \left(\frac{x}{2\sqrt{Dt\pi}}\right) \right] . </math>
 If {{mvar|D}} is time-dependent, the diffusion length becomes
-: <math> 2\sqrt{\int_0^t D( \tau ) \,d\tau}. </math>
+<math display="block"> 2\sqrt{\int_0^t D( \tau ) \,d\tau}. </math>
 This idea is useful for estimating a diffusion length over a heating and cooling cycle, where {{mvar|D}} varies with temperature.
 === Example solution 2: Brownian particle and mean squared displacement ===
 Another simple case of diffusion is the [[Brownian motion]] of one particle. The particle's [[Mean squared displacement]] from its original position is:
-<math display="block">\text{MSD} \equiv \left \langle (\mathbf{x}-\mathbf{x_0})^2 \right \rangle=2nDt , </math>
+<math display="block">\text{MSD} \equiv \left \langle \left(\mathbf{x} - \mathbf{x}_0\right)^2 \right \rangle=2nDt , </math>
 where <math>n</math> is the [[dimension]] of the particle's Brownian motion. For example, the diffusion of a molecule across a [[cell membrane]] 8&nbsp;nm thick is 1-D diffusion because of the spherical symmetry; However, the diffusion of a molecule from the membrane to the center of a [[Eukaryotic Cell|eukaryotic cell]] is a 3-D diffusion. For a cylindrical [[cactus]], the diffusion from photosynthetic cells on its surface to its center (the axis of its cylindrical symmetry) is a 2-D diffusion.
@@ Line 168: / Line 168: @@
 === Generalizations ===
-* In ''non-homogeneous media'', the diffusion coefficient varies in space, {{math|1=''D'' = ''D''(''x'')}}. This dependence does not affect Fick's first law but the second law changes: <math display="block">\frac{\partial \varphi(x,t)}{\partial t}=\nabla\cdot \bigl(D(x) \nabla \varphi(x,t)\bigr)=D(x) \Delta \varphi(x,t)+\sum_{i=1}^3 \frac{\partial D(x)}{\partial x_i} \frac{\partial \varphi(x,t)}{\partial x_i}. </math>
+* In ''non-homogeneous media'', the diffusion coefficient varies in space, {{math|1=''D'' = ''D''(''x'')}}. This dependence does not affect Fick's first law but the second law changes: <math display="block">\begin{align}
+\frac{\partial \varphi(x,t)}{\partial t}&= \nabla\cdot \bigl(D(x) \nabla \varphi(x,t)\bigr) \\
+&= D(x) \Delta \varphi(x,t) + \sum_{i=1}^3 \frac{\partial D(x)}{\partial x_i} \frac{\partial \varphi(x,t)}{\partial x_i}.
+\end{align} </math>
 * In ''[[anisotropic]] media'', the diffusion coefficient depends on the direction. It is a symmetric [[tensor]] {{math|1=''D<sub>ji</sub>'' = ''D<sub>ij</sub>''}}. Fick's first law changes to <math display="block">J=-D \nabla \varphi ,</math> it is the product of a tensor and a vector: <math display="block"> J_i=-\sum_{j=1}^3 D_{ij} \frac{\partial \varphi}{\partial x_j}.</math> For the diffusion equation this formula gives <math display="block">\frac{\partial \varphi(x,t)}{\partial t}=\nabla\cdot \bigl(D \nabla \varphi(x,t)\bigr)=\sum_{i=1}^3\sum_{j=1}^3D_{ij} \frac{\partial^2 \varphi(x,t)}{\partial x_i \partial x_j}. </math> The symmetric matrix of diffusion coefficients {{math|''D<sub>ij</sub>''}} should be [[Positive-definite matrix|positive definite]]. It is needed to make the right-hand side operator [[Elliptic operator|elliptic]].
-* For ''inhomogeneous anisotropic media'' these two forms of the diffusion equation should be combined in <math display="block">\frac{\partial \varphi(x,t)}{\partial t}=\nabla\cdot \bigl(D(x) \nabla \varphi(x,t)\bigr)=\sum_{i,j=1}^3\left(D_{ij}(x) \frac{\partial^2 \varphi(x,t)}{\partial x_i \partial x_j}+ \frac{\partial D_{ij}(x)}{\partial x_i }  \frac{\partial \varphi(x,t)}{\partial x_i}\right). </math>
+* For ''inhomogeneous anisotropic media'' these two forms of the diffusion equation should be combined in <math display="block">\begin{align}
+\frac{\partial \varphi(x,t)}{\partial t}&= \nabla\cdot \bigl(D(x) \nabla \varphi(x,t)\bigr) \\
+&= \sum_{i,j=1}^3\left(D_{ij}(x) \frac{\partial^2 \varphi(x,t)}{\partial x_i \partial x_j}+ \frac{\partial D_{ij}(x)}{\partial x_i }  \frac{\partial \varphi(x,t)}{\partial x_i}\right).
+\end{align} </math>
 * The approach based on [[Diffusion#Einstein's mobility and Teorell formula|Einstein's mobility and Teorell formula]] gives the following generalization of Fick's equation for the ''multicomponent diffusion'' of the perfect components: <math display="block">\frac{\partial \varphi_i}{\partial t} = \sum_j \nabla\cdot\left(D_{ij} \frac{\varphi_i}{\varphi_j} \nabla \, \varphi_j\right) ,</math> where {{mvar|φ<sub>i</sub>}} are concentrations of the components and {{mvar|D<sub>ij</sub>}} is the matrix of coefficients. Here, indices {{mvar|i}} and {{mvar|j}} are related to the various components and not to the space coordinates.
@@ Line 182: / Line 188: @@
 === Fick's flow in liquids ===
-When two [[miscibility|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 walk]]s take 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.<ref>{{cite journal | vauthors = Brogioli D, Vailati A | title = Diffusive mass transfer by nonequilibrium fluctuations: Fick's law revisited | journal = Physical Review E | volume = 63 | issue = 1 Pt 1 | pages = 012105 | date = January 2001 | pmid = 11304296 | doi = 10.1103/PhysRevE.63.012105 | bibcode = 2000PhRvE..63a2105B | arxiv = cond-mat/0006163 | s2cid = 1302913 }}</ref>
+When two [[miscibility|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 walk]]s take 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.<ref>{{cite journal | vauthors = Brogioli D, Vailati A | title = Diffusive mass transfer by nonequilibrium fluctuations: Fick's law revisited | journal = Physical Review E | volume = 63 | issue = 1 Pt 1 | article-number = 012105 | date = January 2001 | pmid = 11304296 | doi = 10.1103/PhysRevE.63.012105 | bibcode = 2000PhRvE..63a2105B | arxiv = cond-mat/0006163 | s2cid = 1302913 }}</ref>
 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 (logic)|tautology]], since the phenomena described by a lower order approximation is the result of a higher approximation: this problem is solved only by [[renormalization|renormalizing]] the fluctuating hydrodynamics equations.
@@ Line 190: / Line 196: @@
 Typically, the diffusion constant of molecules and particles defined by Fick's equation can be calculated using the [[Stokes–Einstein equation]]. In the ultrashort time limit, in the order of the diffusion time ''a''<sup>2</sup>/''D'', where ''a'' is the particle radius, the diffusion is described by the [[Langevin equation]]. At a longer time, the [[Langevin equation]] merges into the [[Stokes–Einstein equation]]. The latter is appropriate for the condition of the diluted solution, where long-range diffusion is considered. According to the [[fluctuation-dissipation theorem]] based on the [[Langevin equation]] in the long-time limit and when the particle is significantly denser than the surrounding fluid, the time-dependent diffusion constant is:<ref>{{cite journal | vauthors = Bian X, Kim C, Karniadakis GE | title = 111 years of Brownian motion | journal = Soft Matter | volume = 12 | issue = 30 | pages = 6331–6346 | date = August 2016 | pmid = 27396746 | pmc = 5476231 | doi = 10.1039/c6sm01153e | bibcode = 2016SMat...12.6331B }}</ref>
-: <math> D(t) = \mu \, k_{\rm B} T\left(1-e^{-t/(m\mu)}\right) , </math>
+<math display="block"> D(t) = \mu \, k_{\rm B} T\left(1-e^{-t/(m\mu)}\right) , </math>
 where (all in SI units)
 * ''k''<sub>B</sub> is the [[Boltzmann constant]],
@@ Line 199: / Line 205: @@
 For a single molecule such as organic molecules or [[biomolecule]]s (e.g. proteins) in water, the exponential term is negligible due to the small product of ''mμ'' in the ultrafast picosecond region, thus irrelevant to the relatively slower adsorption of diluted solute.
-[[File:Diffusive sorption probability.png|300px|thumb|right|Scheme of molecular diffusion in the solution. Orange dots are solute molecules, solvent molecules are not drawn, black arrow is an example random walk trajectory, and the red curve is the diffusive Gaussian broadening probability function from the Fick's law of diffusion.<ref name = "Pyle-BJNano">{{cite journal | vauthors = Pyle JR, Chen J | title = Photobleaching of YOYO-1 in super-resolution single DNA fluorescence imaging | journal = Beilstein Journal of Nanotechnology | volume = 8 | pages = 2296–2306 | date = 2017-11-02 | pmid = 29181286 | pmc = 5687005 | doi = 10.3762/bjnano.8.229 }}</ref><sup>:Fig. 9</sup>]]
+[[File:Diffusive sorption probability.png|thumb|Scheme of molecular diffusion in the solution. Orange dots are solute molecules, solvent molecules are not drawn, black arrow is an example random walk trajectory, and the red curve is the diffusive Gaussian broadening probability function from the Fick's law of diffusion.<ref name = "Pyle-BJNano">{{cite journal | vauthors = Pyle JR, Chen J | title = Photobleaching of YOYO-1 in super-resolution single DNA fluorescence imaging | journal = Beilstein Journal of Nanotechnology | volume = 8 | pages = 2296–2306 | date = 2017-11-02 | pmid = 29181286 | pmc = 5687005 | doi = 10.3762/bjnano.8.229 }}</ref><sup>:Fig. 9</sup>]]
 The [[adsorption]] or [[Absorption (chemistry)|absorption]] rate of a dilute solute to a surface or interface in a (gas or liquid) solution can be calculated using Fick's laws of diffusion. The accumulated number of molecules adsorbed on the surface is expressed by the Langmuir-Schaefer equation by integrating the diffusion flux equation over time as shown in the simulated molecular diffusion in the first section of this page:<ref name = "LangmuirSchaefer1937JACS">{{Cite journal| vauthors = Langmuir I, Schaefer VJ | date = 1937 | title = The Effect of Dissolved Salts on Insoluble Monolayers| journal = Journal of the American Chemical Society | volume = 29 | issue = 11 | pages = 2400–2414 | doi = 10.1021/ja01290a091| bibcode = 1937JAChS..59.2400L }}</ref>
-: <math> \Gamma= 2AC_b\sqrt{\frac{Dt}{\pi}}.</math>
+<math display="block"> \Gamma= 2AC_b\sqrt{\frac{Dt}{\pi}}.</math>
 * {{mvar|A}} is the surface area (m<sup>2</sup>).
 * <math>C_b</math> is the number concentration of the adsorber molecules (solute) in the bulk solution (#/m<sup>3</sup>).
@@ Line 212: / Line 218: @@
 Briefly as explained in,<ref name = "WardTordai1946">{{Cite journal| vauthors = Ward AF, Tordai L |date=1946| title = Time-dependence of Boundary Tensions of Solutions I. The Role of Diffusion in Time-effects| journal = Journal of Chemical Physics | volume = 14 | issue = 7| pages = 453–461 | doi = 10.1063/1.1724167| bibcode = 1946JChPh..14..453W}}</ref>
 the concentration gradient profile near a newly created (from <math>t=0</math>) absorptive surface (placed at <math>x=0</math>) in a once uniform bulk solution is solved in the above sections from Fick's equation,
-:<math> \frac{\partial C}{\partial x} = \frac{C_b}{\sqrt{\pi Dt}}\text{exp} \left (-\frac{x^2}{4Dt} \right ) , </math>
+<math display="block"> \frac{\partial C}{\partial x} = \frac{C_b}{\sqrt{\pi Dt}}\text{exp} \left (-\frac{x^2}{4Dt} \right ) , </math>
 where {{mvar|C}} is the number concentration of adsorber molecules at <math> x, t </math> (#/m<sup>3</sup>).
 The concentration gradient at the subsurface at <math>x = 0</math> is simplified to the pre-exponential factor of the distribution
-:<math> \left (\frac{\partial C}{\partial x} \right ) _{x = 0} = \frac{C_b}{\sqrt{\pi Dt}} . </math>
+<math display="block"> \left (\frac{\partial C}{\partial x} \right ) _{x = 0} = \frac{C_b}{\sqrt{\pi Dt}} . </math>
 And the rate of diffusion (flux) across area <math>A . </math> of the plane is
-:<math> \left (\frac{\partial \Gamma }{\partial t} \right ) _{x = 0} = -\frac{DAC_b}{\sqrt{\pi Dt}} . </math>
+<math display="block"> \left (\frac{\partial \Gamma }{\partial t} \right ) _{x = 0} = -\frac{DAC_b}{\sqrt{\pi Dt}} . </math>
 Integrating over time,
-: <math> \Gamma = \int_0^t \left( \frac{\partial \Gamma}{\partial t} \right) _{x = 0} = 2AC_b\sqrt{\frac{Dt}{\pi}} . </math>
+<math display="block"> \Gamma = \int_0^t \left( \frac{\partial \Gamma}{\partial t} \right) _{x = 0} = 2AC_b\sqrt{\frac{Dt}{\pi}} . </math>
 The Langmuir–Schaefer equation can be extended to the Ward–Tordai Equation to account for the "back-diffusion" of rejected molecules from the surface:<ref name = "WardTordai1946" />
-: <math> \Gamma= 2A{C_\text{b}}\sqrt{\frac{Dt}{\pi}} - A\sqrt{\frac{D}{\pi}}\int_0^\sqrt{t}\frac{C(\tau)}{\sqrt{t-\tau}} \, d\tau , </math>
+<math display="block"> \Gamma= 2A{C_\text{b}}\sqrt{\frac{Dt}{\pi}} - A\sqrt{\frac{D}{\pi}}\int_0^\sqrt{t}\frac{C(\tau)}{\sqrt{t-\tau}} \, d\tau , </math>
 where <math>C_b</math> is the bulk concentration, <math>C</math> is the sub-surface concentration (which is a function of time depending on the reaction model of the adsorption), and <math>\tau</math> is a dummy variable.
@@ Line 233: / Line 239: @@
 A more problematic result of the above equations is they predict the lower limit of adsorption under ideal situations but is very difficult to predict the actual adsorption rates. The equations are derived at the long-time-limit condition when a stable concentration gradient has been formed near the surface. But real adsorption is often done much faster than this infinite time limit i.e. the concentration gradient, decay of concentration at the sub-surface, is only partially formed before the surface has been saturated or flow is on to maintain a certain gradient, thus the adsorption rate measured is almost always faster than the equations have predicted for low or none energy barrier adsorption (unless there is a significant adsorption energy barrier that slows down the absorption significantly), for example, thousands to millions time faster in the self-assembly of monolayers at the water-air or water-substrate interfaces.<ref name = LangmuirSchaefer1937JACS/> As such, it is necessary to calculate the evolution of the concentration gradient near the surface and find out a proper time to stop the imagined infinite evolution for practical applications. While it is hard to predict when to stop but it is reasonably easy to calculate the shortest time that matters, the critical time when the first nearest neighbor from the substrate surface feels the building-up of the concentration gradient. This yields the upper limit of the adsorption rate under an ideal situation when there are no other factors than diffusion that affect the absorber dynamics:<ref name=JixinMCSimuAdsorption/>
-: <math> \langle r \rangle = \frac{4}{\pi}A C_b^{4/3}D , </math>
+<math display="block"> \langle r \rangle = \frac{4}{\pi}A C_b^{4/3}D , </math>
 where:
 * <math> \langle r \rangle </math> is the adsorption rate assuming under adsorption energy barrier-free situation, in unit #/s,
@@ Line 243: / Line 249: @@
 This critical time is significantly different from the first passenger arriving time or the mean free-path time. Using the average first-passenger time and Fick's law of diffusion to estimate the average binding rate will significantly over-estimate the concentration gradient because the first passenger usually comes from many layers of neighbors away from the target, thus its arriving time is significantly longer than the nearest neighbor diffusion time. Using the mean free path time plus the Langmuir equation will cause an artificial concentration gradient between the initial location of the first passenger and the target surface because the other neighbor layers have no change yet, thus significantly lower estimate the actual binding time, i.e., the actual first passenger arriving time itself, the inverse of the above rate, is difficult to calculate. If the system can be simplified to 1D diffusion, then the average first passenger time can be calculated using the same nearest neighbor critical diffusion time for the first neighbor distance to be the MSD,<ref name = "Pandey-JPCB2024">{{cite journal | vauthors = Pandey S, Gautam D, Chen J | title = Measuring the Adsorption Cross Section of YOYO-1 to Immobilized DNA Molecules | journal = Journal of Physical Chemistry B | volume =  128| pages = 7254–7262
-  | date = 2024-07-16 | issue = 29 | pmid =  39014882| pmc = 11286311| doi = 10.1021/acs.jpcb.4c03359 | pmc-embargo-date = July 25, 2025 }}</ref>
+  | date = 2024-07-16 | issue = 29 | pmid =  39014882| pmc = 11286311| doi = 10.1021/acs.jpcb.4c03359 }}</ref>
-:<math>L = \sqrt{2Dt} , </math>
+<math display="block">L = \sqrt{2Dt} , </math>
 where:
 *<math>L~=C_b^{-1/3} </math> (unit m) is the average nearest neighbor distance approximated as cubic packing, where <math>C_b</math> is the solute concentration in the bulk solution (unit # molecule / m<sup>3</sup>),
@@ Line 250: / Line 256: @@
 *<math>t</math> is the critical time (unit s).
 In this critical time, it is unlikely the first passenger has arrived and adsorbed. But it sets the speed of the layers of neighbors to arrive. At this speed with a concentration gradient that stops around the first neighbor layer, the gradient does not project virtually in the longer time when the actual first passenger arrives. Thus, the average first passenger coming rate (unit # molecule/s) for this 3D diffusion simplified in 1D problem,
-:<math> <r> = \frac{a}{t} = 2aC_b^{2/3}D , </math>
+<math display="block"> \langle r \rangle = \frac{a}{t} = 2aC_b^{2/3}D , </math>
 where <math> a</math> is a factor of converting the 3D diffusive adsorption problem into a 1D diffusion problem whose value depends on the system, e.g., a fraction of adsorption area <math>A</math> over solute nearest neighbor sphere surface area <math>4 \pi L^2 /4</math> assuming cubic packing each unit has 8 neighbors shared with other units. This example fraction converges the result to the 3D diffusive adsorption solution shown above with a slight difference in pre-factor due to different packing assumptions and ignoring other neighbors.
@@ Line 260: / Line 266: @@
 The bimolecular collision frequency related to many reactions including protein coagulation/aggregation is initially described by [[Smoluchowski coagulation equation]] proposed by [[Marian Smoluchowski]] in a seminal 1916 publication,<ref name=Smoluchowski1916>{{cite journal | vauthors = Smoluchowski M | title = Drei Vorträge über Diffusion, Brownsche Molekularbewegung und Koagulation von Kolloidteilchen | journal =  Zeitschrift für Physik| year = 1916 | volume = 17 | pages = 557–571, 585–599 | language = German | bibcode = 1916ZPhy...17..557S }}</ref> derived from [[Brownian motion]] and Fick's laws of diffusion. Under an idealized reaction condition for A + B → product in a diluted solution, Smoluchovski suggested that the molecular flux at the infinite time limit can be calculated from Fick's laws of diffusion yielding a fixed/stable concentration gradient from the target molecule, e.g. B is the target molecule holding fixed relatively, and A is the moving molecule that creates a concentration gradient near the target molecule B due to the coagulation reaction between A and B. Smoluchowski calculated the collision frequency between A and B in the solution with unit #/s/m<sup>3</sup>:
-: <math> Z_{AB} = 4{\pi}RD_rC_AC_B,</math>
+<math display="block"> Z_{AB} = 4{\pi}RD_rC_AC_B,</math>
 where:
 * <math>R</math> is the radius of the collision,
@@ Line 271: / Line 277: @@
 In 2022, Chen calculates the upper limit of the collision frequency between A and B in a solution assuming the bulk concentration of the moving molecule is fixed after the first nearest neighbor of the target molecule.<ref name=JChen2022JPCA>{{cite journal | vauthors = Chen J | title = Why Should the Reaction Order of a Bimolecular Reaction be 2.33 Instead of 2? | journal = The Journal of Physical Chemistry A | volume = 126 | issue = 51 | pages = 9719–9725 | date = December 2022 | pmid = 36520427 | pmc = 9805503 | doi = 10.1021/acs.jpca.2c07500 | bibcode = 2022JPCA..126.9719C }}</ref> Thus the concentration gradient evolution stops at the first nearest neighbor layer given a stop-time to calculate the actual flux. He named this the critical time and derived the diffusive collision frequency in unit #/s/m<sup>3</sup>:<ref name=JChen2022JPCA/>
-: <math> Z_{AB} = \frac{8}{\pi}{\sigma} D_rC_AC_B\sqrt[3]{C_A+C_B} , </math>
+<math display="block"> Z_{AB} = \frac{8}{\pi}{\sigma} D_rC_AC_B\sqrt[3]{C_A+C_B} , </math>
 where:
 * <math>{\sigma}</math> is the area of the cross-section of the collision (m<sup>2</sup>),
 * <math>D_r = D_A + D_B</math> is the relative diffusion constant between A and B (m<sup>2</sup>/s),
 * <math>C_A</math> and <math>C_B</math> are number concentrations of A and B respectively (#/m<sup>3</sup>),
-* <math>\sqrt[3]{C_A+C_B} </math> represents 1/<d>, where d is the average distance between two molecules.
+* <math>\sqrt[3]{C_A+C_B} </math> represents 1/⟨''d''⟩, where ''d'' is the average distance between two molecules.
-This equation assumes the upper limit of a diffusive collision frequency between A and B is when the first neighbor layer starts to feel the evolution of the concentration gradient, whose reaction order is {{sfrac|2|1|3}} instead of 2. Both the Smoluchowski equation and the JChen equation satisfy dimensional checks with SI units. But the former is dependent on the radius and the latter is on the area of the collision sphere. From dimensional analysis, there will be an equation dependent on the volume of the collision sphere but eventually, all equations should converge to the same numerical rate of the collision that can be measured experimentally. The actual reaction order for a bimolecular unit reaction could be between 2 and {{sfrac|2|1|3}}, which makes sense because the diffusive collision time is squarely dependent on the distance between the two molecules.
+This equation assumes the upper limit of a diffusive collision frequency between A and B is when the first neighbor layer starts to feel the evolution of the concentration gradient, whose reaction order is {{sfrac|2|1|3}} instead of 2. Both the Smoluchowski equation and the JChen equation satisfy dimensional checks with SI units. But the former is dependent on the radius and the latter is on the area of the collision sphere. From dimensional analysis, there should be an equation dependent on the volume of the collision sphere,<ref name = "Chen-AIP2024">{{cite journal | vauthors = Chen J | title = Dimensional analysis of diffusive association rate equations | journal = AIP Advances | volume =  14| article-number = 115218
+ | date = 2024-11-14 | issue = 11 | pmid =   39555209 | pmc = 11567696 | doi = 10.1063/5.0238119 | bibcode = 2024AIPA...14k5218C }}</ref> e.g.,
+<math display="block"> Z_{AB} = 4 VD_rC_AC_B(C_A+C_B)^{2/3}, </math>
+*V is the volume of the collision sphere
+but eventually, all equations should converge to the same numerical rate of the collision that can be measured experimentally. The actual reaction order for a bimolecular unit reaction could be between 2 and {{sfrac|2|2|3}}, which makes sense because the diffusive collision time is squarely dependent on the distance between the two molecules.
 These new equations also avoid the singularity on the adsorption rate at time zero for the Langmuir-Schaefer equation. The infinity rate is justifiable under ideal conditions because when you introduce target molecules magically in a solution of probe molecule or vice versa, there always be a probability of them overlapping at time zero, thus the rate of that two molecules association is infinity. It does not matter that other millions of molecules have to wait for their first mate to diffuse and arrive. The average rate is thus infinity. But statistically this argument is meaningless. The maximum rate of a molecule in a period of time larger than zero is 1, either meet or not, thus the infinite rate at time zero for that molecule pair really should just be one, making the average rate 1/millions or more and statistically negligible. This does not even count in reality no two molecules can magically meet at time zero.
@@ Line 284: / Line 296: @@
 === Biological perspective ===
 The first law gives rise to the following formula:<ref>{{cite book| title= Essentials of Human Physiology|  vauthors = Nosek TM | chapter=Section 3/3ch9/s3ch9_2 |chapter-url=http://humanphysiology.tuars.com/program/section3/3ch9/s3ch9_2.htm |archive-url=https://web.archive.org/web/20160324124828/http://humanphysiology.tuars.com/program/section3/3ch9/s3ch9_2.htm|archive-date=2016-03-24}}</ref>
-: <math>\text{flux} = {-P \left(c_2 - c_1\right)} , </math>
+<math display="block">\text{flux} = {-P \left(c_2 - c_1\right)} , </math>
 where
 * {{mvar|P}} is the permeability, an experimentally determined membrane "[[Electrical conductance|conductance]]" for a given gas at a given temperature,
@@ Line 294: / Line 306: @@
 Under the condition of a diluted solution when diffusion takes control, the membrane permeability mentioned in the above section can be theoretically calculated for the solute using the equation mentioned in the last section (use with particular care because the equation is derived for dense solutes, while biological molecules are not denser than water. Also, this equation assumes ideal concentration gradient forms near the membrane and evolves):<ref name = "Pyle-BJNano" />
-: <math> P= 2A_p\eta_{tm} \sqrt{ \frac{D}{\pi t}} , </math>
+<math display="block"> P= 2A_p\eta_{tm} \sqrt{ \frac{D}{\pi t}} , </math>
 where:
 * <math>A_P</math> is the total area of the pores on the membrane (unit m<sup>2</sup>),
@@ Line 313: / Line 325: @@
 [[Integrated circuit]] fabrication technologies, model processes like CVD, thermal oxidation, wet oxidation, doping, etc. use diffusion equations obtained from Fick's law.
-==== CVD method of fabricate semiconductor ====
+==== CVD method of semiconductor fabrication====
 The wafer is a kind of semiconductor whose silicon substrate is coated with a layer of CVD-created polymer chain and films. This film contains n-type and p-type dopants and takes responsibility for dopant conductions. The principle of CVD relies on the gas phase and gas-solid chemical reaction to create thin films.
@@ Line 371: / Line 383: @@
 == See also ==
+{{div col}}
 * [[Advection]]
 * [[Churchill–Bernstein equation]]
@@ Line 380: / Line 393: @@
 * [[Nernst–Planck equation]]
 * [[Osmosis]]
+{{div col end}}
 == Citations ==
@@ Line 385: / Line 399: @@
 <ref name=JixinMCSimuAdsorption>
-{{cite journal | vauthors = Chen J | title = Simulating stochastic adsorption of diluted solute molecules at interfaces | journal = AIP Advances | volume = 12 | issue = 1 | pages = 015318 | date = January 2022 | pmid = 35070490 | pmc = 8758205 | doi = 10.1063/5.0064140 | bibcode = 2022AIPA...12a5318C }}
+{{cite journal | vauthors = Chen J | title = Simulating stochastic adsorption of diluted solute molecules at interfaces | journal = AIP Advances | volume = 12 | issue = 1 | article-number = 015318 | date = January 2022 | pmid = 35070490 | pmc = 8758205 | doi = 10.1063/5.0064140 | bibcode = 2022AIPA...12a5318C }}
 </ref>
@@ Line 405: / Line 419: @@
 [[Category:Diffusion]]
-[[Category:Eponymous laws of physics]]
 [[Category:Mathematics in medicine]]
 [[Category:Physical chemistry]]

Fick's laws of diffusion: Difference between revisions

Latest revision as of 21:47, 13 November 2025

History

Fick's first law

Variations of the first law

Derivation of Fick's first law for gases

Fick's second law

Derivation of Fick's second law

Example solutions and generalization

Example solution 1: constant concentration source and diffusion length

Example solution 2: Brownian particle and mean squared displacement

Generalizations

Applications

Fick's flow in liquids

Sorption rate and collision frequency of diluted solute

Biological perspective

Semiconductor fabrication applications

CVD method of semiconductor fabrication

Invalidity of Fickian diffusion

Food production and cooking

See also

Citations

Further reading

External links

Navigation menu

Search