Diffusion ($\mathbf{j}=-D\nabla C$)

Diffusion is net transport associated with random microscopic motion. A localized patch of a dilute solute can spread through a stationary solvent: individual molecules move in both directions, but more cross from the crowded side of an imaginary surface than from the less crowded side. This imbalance gives a net diffusive flux down the concentration gradient; the molecules do not need to sense which direction has lower concentration. At uniform concentration, microscopic exchanges continue even though the net diffusive flux vanishes. ^[ref]

The distinction from advection is easiest to see in a dye patch. Uniform bulk motion carries its centre downstream; diffusion broadens it relative to that centre. Both can happen together. A moving or spreading patch by itself does not identify the mechanism: a nonuniform velocity field can also stretch a patch. ^[ref]

For a dilute nonreacting solute in a uniform medium at fixed temperature, Fick's first law models the relative diffusive flux as $\mathbf{j}=-D\nabla C$, with scalar diffusion coefficient $D>0$ of units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Here $C$ is the amount of substance per unit volume, so $\mathbf{j}$ has units $\mathrm{mol}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$. With carrier velocity field $\mathbf{u}$, the total flux is $C\mathbf{u}+\mathbf{j}$. The continuity equation, with no sources or chemical reactions, givesFor constant $D$ and incompressible flow, this reduces to $C_t+\mathbf{u}\cdot \nabla C=D{\nabla }^{2}C$, where ${\nabla }^2$ is the Laplacian. In a stationary carrier it is the diffusion equation. These assumptions matter: variable $D$ must remain inside the divergence, and changing sources require extra terms. ^[ref]

The useful surprise is the common mathematics: thermal conduction smooths temperature through the analogous heat equation, while molecular diffusion smooths composition. This analogy already motivated the original formulation of Fick's laws. ^[ref] However, the thermal diffusivity and the solutal diffusivity need not be equal. A diffusive time scale $t_D\sim L^2/D$ means that crossing ten times the distance requires roughly one hundred times the time. This is why molecular diffusion is especially effective across thin layers, and why chaotic advection can assist mixing by making those layers thinner. Physical diffusion is a transport mechanism, distinct from the broader probabilistic class of Markov diffusion.

Relaxation of a density or composition fluctuation involving many constituents. Its coefficient can depend on osmotic compressibility and collective mobility, and need not equal the tagged-particle self-diffusion coefficient at nonzero concentration.

Random displacement of a tagged constituent, even in a material of uniform composition. The displayed long-time relation defines ordinary isotropic diffusion in $d$ dimensions.

The scale $\sqrt{Dt}$ characterizes diffusion over time $t$. For unrestricted isotropic diffusion, the rms displacement in one coordinate is $\sqrt{2Dt}$ and the total three-dimensional rms displacement is $\sqrt{6Dt}$. Stating the convention avoids spurious disagreement over a numerical prefactor.

For translational drag coefficient $\zeta$ and mobility $1/\zeta$, thermal equilibrium requires $D=k_{B}T/\zeta$. This connects random displacement to dissipative response and is a form of the fluctuation-dissipation theorem. Combining it with Stokes drag law gives the Stokes–Einstein relation.

A transport closure in which a current is proportional to the negative gradient of a scalar field. It discards detailed angular dependence and becomes unreliable where a nearly isotropic local distribution cannot be assumed. Effective coefficients and boundary modelling must match the physical application.

A diffusion approximation for a scalar neutron flux, with one prescribed energy group. Its terms balance current loss, absorption, nuclear fission production and a source. A solution also requires appropriate boundary conditions.

The positive separation parameter in $-{\nabla }^2\varphi =B_g^2\varphi$ for a specified fundamental spatial mode and boundary conditions. The corresponding material parameter must match it for a source-free steady solution.

The coefficient in $\mathbf{J}=-D\nabla \varphi$ for the scalar neutron flux convention. It has dimensions of length, unlike a particle-density diffusion coefficient with dimensions length squared per time. For a single fixed speed $v$, the corresponding density diffusion coefficient is $vD$.

Movement of material by diffusion, convection or migration.

diffusion inferred by tracking labelled particles in an equilibrium material. Local jump correlations and carrier concentration must be accounted for when relating it to electrical conductivity.

The factor $f_D$ correcting an independent-jump mean-square-displacement estimate. If ${\tau }_j$ is the total waiting time between jumps of length $\ell$, then $D=f_D{\ell }^2/(6{\tau }_j)$ in three dimensions.

diffusion by hopping of an interstitial particle between interstitial sites without displacing a host lattice particle. Direction sums over neighbour vectors determine the conductivity tensor in a dilute independent-jump model.

Migration through successive exchanges of atoms or ions with vacant lattice sites. The diffusivity depends on both defect abundance and jump kinetics; mobile-defect diffusivity and atomic tracer diffusivity are different quantities.

Molecular diffusion is the component of material transport caused by the random thermal motion of molecules, measured relative to bulk motion of their surroundings. A diffusion coefficient is a property of a particular substance in a particular medium under specified conditions, not a universal speed of spreading. Its dependence on the solute, solvent and temperature is precisely what makes a spreading experiment informative. A random walk idealizes this motion; the Brownian motion description provides a continuous approximation over sufficiently long times. ^[ref]

A motivating experiment: watch dye enter a gel. Samprovalaki, Robbins and Fryer studied rhodamine 6G and methylene blue entering agar gel from two side reservoirs. A covered slab, $10\times 5\times 1\mathrm{mm}$, sat on a controlled thermal stage; the dye was brought to the same temperature. They photographed the gel every ten minutes, averaged image intensity across a strip, and fitted the changing concentration to a diffusion equation model. The cover limited evaporation and stopped dye travelling over the gel surface. Runs at $30$, $50$ and $70{}^{\circ }\mathrm{C}$ were repeated in triplicate. The measured spatial colour change therefore supplied a quantitative transport test, rather than merely showing that dye eventually spreads. Changing lighting also changes image intensity, so the optical response had to be checked. Importantly for interpretation, their constant-diffusion coefficient fit was imperfect and the gel coefficient is an effective transport coefficient, not automatically the value in free water. ^[ref]

A finite-pulse version to test width growth. The following is a proposed teaching variant, not that paper's apparatus or a report of measurements. Make a long, narrow, transparent gel channel of uniform cross-section. Join a thin dyed gel band between two undyed lengths of the same gel, with good contact and no air gaps; the dye should initially fill the cross-section so the dominant concentration variation is along the channel coordinate $x$. After assembly, keep the covered channel horizontal at uniform temperature, away from vibration. The dye is a finite initial pulse: no reservoir supplies more dye. Choose a dye whose image response can be calibrated and whose interaction with the gel is small enough for an approximately constant diffusion coefficient to be plausible.

Photograph the whole pulse at recorded times with fixed lighting, focus and exposure, including a length scale. Use uniformly dyed standards of the same gel thickness and a blank gel to convert image intensity to concentration $C(x,t)$; do not assume raw brightness is proportional to concentration. Average across the channel width and examine several parallel strips to check the one-dimensional approximation. Measure the initial profile instead of treating the band as infinitely thin. For a finite pulse with negligible end effects, define its spatial integral $Q$, centre $\bar{x}$ and width $\sigma$ byHere ${\sigma }^2$ is the variance of the normalized profile $C/Q$, and $\sigma$ is its standard deviation. If $C$ is measured in $\mathrm{mol}{\mathrm{m}}^{-3}$, then $Q$ has units $\mathrm{mol}{\mathrm{m}}^{-2}$; multiplying by the channel cross-section gives the total amount of substance. For the unbounded one-dimensional diffusion equation with constant $D$, conserved $Q$, finite variance and sufficiently decaying tails,Thus plot width squared against time: its slope determines the diffusion coefficient, while its intercept includes the finite starting width. This moment law holds even for a non-Gaussian distribution initially; it follows by applying integration by parts twice to the diffusion equation. ^[ref]

For the particularly readable case of an initially Gaussian distribution centred at zero, uniform advection at speed $u$ and constant diffusion coefficient $D$ give the one-dimensional concentration profileIts centre moves as $\bar{x}=ut$ and its variance grows as ${\sigma }_0^2+2Dt$. Our illustrative choice ${\sigma }_0=0.50\mathrm{mm}$ and $D=0.20\mathrm{m}{\mathrm{m}}^2{\mathrm{h}}^{-1}$ gives widths $0.50$, $0.81$ and $1.36\mathrm{mm}$ at zero, one and four hours. These are calculated values, not a measured dye coefficient. Once $2Dt$ dominates the initial variance, width grows as $\sqrt{2Dt}$: four times the elapsed time gives approximately twice the width. A rectangular starting band will not immediately have a Gaussian distribution, although the same width-squared law still applies under the stated assumptions.

Controls and limits of the interpretation. In the proposed experiment, check that $Q$ stays constant and that the centre remains stationary while the width increases. Include a blank and a uniformly dyed gel photographed over the same period: an apparent signal change there exposes lighting drift, dye fading or evaporation. Repeat with a lower starting dye concentration and a second initial width; agreement of the slopes is a check on constant $D$ and on preparation effects. Keep the pulse far from the channel ends and record its faint tails: truncation and background errors can strongly bias a variance. Poor contact, dye binding, chemical reactions or transport in a liquid film can invalidate the simple model. A gel suppresses bulk liquid motion, but its inferred diffusivity need not equal that in the pure solvent. A stationary centre alone does not exclude symmetric flows, and linear width-squared growth alone does not prove a molecular mechanism: turbulence can also produce an effective diffusive description. These are reasons to combine the profile test with the controlled setup.

Rotational diffusion is the random change of a particle's rigid body orientation through Brownian motion. A suspended ellipsoid can turn while its centre hardly moves; translational diffusion instead spreads its position. Random molecular collisions supply fluctuating torques as well as forces. The rotational diffusion coefficient $D_r$ measures how quickly angular uncertainty grows, rather than a steady spin speed. Its units are ${\mathrm{s}}^{-1}$, often written $\mathrm{ra}{\mathrm{d}}^2{\mathrm{s}}^{-1}$ to distinguish angular spreading from translational spreading, whose diffusion coefficient $D_t$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. The radian is dimensionless, but writing angular variance in $\mathrm{ra}{\mathrm{d}}^2$ remains useful.

First identify what the image can reveal. A marked end of an elongated particle defines an arrow: turning it through $\pi$ reverses that arrow. An unmarked ellipsoid with identical ends reveals only its long axis, a director: angles $\varphi$ and $\varphi +\pi$ describe the same visible state. An unmarked perfect sphere has no visible rigid body orientation at all. Tracking its centre cannot measure its rotational diffusion; a surface mark or another identifiable internal direction is needed.

This distinction makes rotational diffusion useful: it measures how long a particle remembers its pointing direction. For an elongated particle, turning also changes which laboratory direction is easiest for translational diffusion. Randomization of rigid body orientation can therefore erase directional preferences in the observed displacement, even though the particle's shape has not changed.

For free planar turning with constant $D_r$, no imposed torque and negligible inertia, model an unwrapped angle by the stochastic differential equationHere $W$ is standard one-dimensional Brownian motion, with increment variance $dt$, and brackets denote expected values over repeated realizations. The angular variance is conventionally reported in $\mathrm{ra}{\mathrm{d}}^2$. “Unwrapped” means recording accumulated turns continuously, rather than resetting the angle at each revolution. Directly subtracting wrapped camera angles creates artificial jumps; for an unmarked ellipsoid the wrapping period is $\pi$, and resolving successive turns additionally requires sufficiently frequent images.

For a planar arrow $\mathbf{u}=(\cos \varphi ,\sin \varphi )$, the orientational autocorrelation isFor an unoriented planar axis, use an autocorrelation unchanged by replacing either endpoint angle by that angle plus $\pi$:The same $D_r$ thus gives arrow and axis correlation times $1/D_r$ and $1/(4D_r)$, respectively, where a correlation time here is the lag at which the relevant curve reaches $1/e$. A short explanation is enough: $\Delta \varphi$ has a Gaussian distribution of variance $2D_{r}t$, so its characteristic function gives $⟨\cos (n\Delta \varphi )⟩=e^{-n^2{D}_{r}t}$. The integer $n$ selects what the observable can distinguish; it does not change the underlying angular spreading.

For an original numerical example, choose $D_r=0.25{\mathrm{s}}^{-1}$. After $2\mathrm{s}$, the angular variance is $1\mathrm{ra}{\mathrm{d}}^2$ and its root mean square displacement is $1\mathrm{rad}$, although the mean displacement is zero. The planar arrow autocorrelation is $e^{-0.5}\simeq 0.607$, while the planar axis autocorrelation is $e^{-2}\simeq 0.135$. Calling both “the orientation correlation” without specifying the observable would hide a factor of four in decay rate.

For isotropic turning in three dimensions, the tip of a marked unit vector $\mathbf{u}$ wanders over the unit sphere. Its probability density function $f$ obeys the orientational Fokker-Planck equationHere ${\nabla }_{\mathbf{u}}^2$ is the Laplace-Beltrami operator on the unit sphere. The factor two follows from the sphere surface derivative identities: ${\nabla }_{\mathbf{u}}^2{u}_i=-2u_i$, so each mean vector component decays at rate $2D_r$. This is a different orientation space from the circle of planar rotational diffusion, whose corresponding rate is $D_r$. For the same illustrative $D_r=0.25{\mathrm{s}}^{-1}$ and $t=2\mathrm{s}$, this three-dimensional autocorrelation is $e^{-1}\simeq 0.368$.

One unit vector records a pointing direction, not every aspect of a rigid body's rigid body orientation: rotation about that vector remains unspecified. Also, $C_1$ requires an identifiable arrow. An unmarked ellipsoid's long axis alone does not supply one, and projecting a three-dimensional trajectory onto a camera plane does not automatically produce the planar model.

For a sphere of radius $a$ in an otherwise stationary, unbounded Newtonian fluid with dynamic viscosity $\eta$, impose the no-slip boundary condition and the slow-flow assumptions of Stokes flow. The resisting torque on a rotating sphere in Stokes flow is $-{\zeta }_r\mathbf{\Omega }$, where $\mathbf{\Omega }$ is its angular velocity and ${\zeta }_r=8\pi \eta a^3$. Thermal rotational diffusion at absolute temperature $T$ then givesHere $k_B$ is the Boltzmann constant; the translational expression is the Stokes–Einstein relation. The ratio has units of inverse length squared, so comparing the numerical values of $D_r$ and $D_t$ without their units is meaningless. At fixed temperature and dynamic viscosity, decreasing radius tenfold increases $D_r$ a thousandfold but $D_t$ only tenfold. This strong size dependence explains why small particles can lose directional memory extremely quickly.

Choose $T=298\mathrm{K}$ and $\eta =0.001\mathrm{Pa}\mathrm{s}$ as illustrative model inputs. Using $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$ gives $k_{B}T=4.11433402\times 10^{-21}\mathrm{J}$ and the following predictions:These are bulk three-dimensional sphere calculations, explicitly not a fit to the confined ellipsoid experiment below. Walls, nonspherical shape and departure from the no-slip boundary condition can change rotational resistance. A visible surface mark is still necessary to read the sphere's turning in an image, even when its rotational diffusion coefficient is well defined physically.

In 2006, Yilong Han and colleagues dispersed plastic ellipsoids in water inside a thin glass cell, using a very dilute suspension to limit particle interactions. A camera on an optical microscope recorded 30 frames per second. Image analysis extracted centre positions and long-axis angles, with reported angular resolution about $1^{\circ }$. Their shape supplied the observable axis. Reflected interference colours calibrated cell thickness; nearly constant particle interference colour checked vertical wandering. One tracked particle had dimensions $2.4\times 0.3\times 0.3\mu \mathrm{m}$ in an $846\mathrm{nm}$ gap. Angular displacement variance grew linearly with lag; displacement resolved along the particle's axes and in the laboratory showed how turning washes out translational anisotropy. These are observations of confined particles: walls alter drag, and slight out-of-plane tumbling affected positional resolution. They do not establish the bulk sphere formula. The paper's Figure 1A–B shows simulations; Figure 1F shows a measured trajectory.^[ref]

For interpreting such images, the director distinction is essential: identical ends give an angle modulo $\pi$. A continuously unwrapped trace can follow small successive turns, but it does not make the ends physically distinguishable. The invariant planar axis autocorrelation is $⟨\cos (2\Delta \varphi )⟩$, and its decay determines $D_r$ through $e^{-4D_{r}t}$. Separately measuring position and rigid body orientation is what makes the connection between translational anisotropy and loss of angular memory experimentally accessible.

Diffusion of particle position, as distinguished from rotational diffusion of its orientation. Its diffusion coefficient has units of length squared per time.