Consider the one-dimensional diffusion equation on the whole real line,
with initial value one for and zero for . This is an idealized Heaviside step function profile, not a measured material or an apparatus specification. The heat kernel gives a convolution with the initial profile:
The change of variables turns the convolution into a Gaussian integral tail. The resulting similarity variable shows that the transition width grows like : multiplying elapsed time by four doubles every fixed-profile distance from the original interface.
Figure 1.
Diffusion smooths a step while preserving its similarity profile
. Exact model curves for a unit jump with stipulated in consistent length and time units. Left: the transition broadens at , , and . Right: plotting each against collapses all three onto the same curve. No experimental measurements are plotted. Details.
At every positive time, , while as and as . For every fixed , the limit as returns the specified initial value. The value assigned at the single initial discontinuity does not change this convolution. The error function therefore explains smoothing without pretending that a discontinuous initial profile remains sharp.
Direct differentiation checks the partial differential equation. Writing gives
The model assumes constant diffusion coefficient, no advection, and no finite-domain boundary effects. For an interpretation as heat conduction, the heat equation applies to a normalized temperature difference under its own material assumptions.
Facilitated target search 2026-10-24
A facilitated target search combines bulk diffusion with lower-dimensional motion on a substrate, such as sliding on DNA. In the scaling model , positive constants balance relocation and local exploration. Differentiation gives and . These expressions require an accessible optimum within the physically allowed range of sliding lengths.
Consider two dilute populations of spheres with the same buoyant density and radii and . Adopt the fixed linear density gradient and constant local centrifugal acceleration of linear-gradient isopycnic focusing, with constant dynamic viscosity and negligible diffusion. Stokes's law makes proportional to , so doubling radius multiplies the migration rate parameter by four, although the mass increases eightfold.
Use , and . The liquid occupies , so the matching density position lies beyond its outer boundary. Both populations begin at ; choose and . Their band-center trajectories are
At , they lie at approximately and ; at , at and . Equal buoyant density has not prevented a useful finite-time separation. Solving using the natural logarithm gives first boundary arrival at . The corresponding arrival for A is . Trajectories describe migration only until contact with the boundary; subsequent collection into a centrifugation pellet is not represented by extending the curve beyond .
Figure 1.
Rate-zonal centrifugation: equal buoyant density, different finite-time migration
. Original schematic tubes and trajectories for the stipulated finite-time rate-zonal separation model. B represents spheres with twice A's radius. Both matching density positions lie outside the column. Tube shading shows an imposed increasing density gradient; band thickness is schematic. The marked boundary-arrival time is calculated from this model, not an experimental run time or an operating instruction. Details.
This is the useful contrast with isopycnic centrifugation: stopping time helps determine the recovered positions in rate-zonal centrifugation, whereas a stationary matching density position sets the limiting band center in isopycnic centrifugation. Real resolution also depends on initial zone width, diffusion and variations in particle properties; distinct calculated centers alone do not prove that fractions will be pure.
Fluid mosaic model 2026-10-13
The fluid mosaic model describes a cell membrane as a fluid lipid bilayer with embedded membrane proteins. It connects a continuous barrier to a mosaic of functional components: diffusion within the membrane can coexist with selective passage across it. S. J. Singer and G. L. Nicolson presented the fluid mosaic model in Science in 1972; see Singer and Nicolson, The fluid mosaic model of the structure of cell membranes. The useful picture does not require every protein to move freely: attachments and interactions can restrict movement, as discussed in Cooper, The Cell: structure of the plasma membrane.
With additive drag from identical monomers, and centre-of-mass diffusion is . For monomer Stokes law , this gives . The model neglects long-range hydrodynamic coupling and does not alone describe the long-time spatial diffusion of a strongly entangled chain.
Frozen-temperature approximation Created 2026-10-14 Updated 2026-10-20
The frozen-temperature approximation prescribes the external thermal field during interface perturbations, rather than solving a perturbation heat equation. An interface displacement samples that prescribed field at its new location. Local liquidus and Gibbs-Thomson effect conditions still apply, and the solute field may retain full time-dependent diffusion.
Gas centrifuge enrichment 2026-10-20
A technology exploiting mass-dependent redistribution of gas under rotation. It became the dominant commercial uranium-enrichment technology with much lower electricity consumption than gaseous diffusion. No operating conditions or equipment specifications are supplied.
For diffusion with a constant diffusion coefficient and no sources, the model flux is . A continuity equation gives . This proportionality is a constitutive assumption about transport; the continuity equation expresses conservation. If varies with position, it must remain inside the divergence, giving instead.[ref]
An exact constructed heat equation example on the whole plane uses dimensionless coordinates with and the smooth initial Gaussian function . Set . For ,
Direct differentiation verifies both the initial condition and
The center decreases, while the region increases at that instant. The sign-change circle expands with radius . A fixed outer location can initially increase and later decrease as this circle passes it.
Figure 1.
An exact Gaussian function spreads under the heat equation while its spatial integral stays fixed
. Original analytic model with . Left: radial profiles at , whose peaks are . Right: at ; the sign changes at . These are model predictions, not experimental measurements. The conserved quantity uses the whole-plane area element , rather than the unweighted area under a radial profile. Details.
At time , the peak is , and the radius where the profile is times its peak is . The width grows as the height falls. Integration in polar coordinates gives
Thus a decreasing peak is compatible with a conserved spatial integral. Dividing by makes a probability density function with coordinate variances , connecting the same spreading to Brownian motion with diffusion coefficient one. Equivalently the profile is evolved by the Gaussian heat kernel with time parameter , since that article uses generator .
For a sphere of radius and buoyant density difference , define using radius as the transport length. With Stokes law and Stokes–Einstein diffusion, . Replacing radius by diameter changes the convention by a factor of two. This local criterion does not guarantee a uniform suspension over a macroscopic vessel.
For orientation probability density on the unit sphere, deterministic drift and isotropic rotational diffusion coefficient give . The derivatives are tangential to the sphere, and normalization is . The spherical divergence of is .
Hahn spin echo experiment 2026-10-18
Erwin Louis Hahn's 1950 apparatus enclosed liquid samples in the coil of a tuned LC circuit, coupled to a gated radio transmitter and radio receiver. One transmitter quadrupled a crystal oscillator's 7.5 MHz to 30 MHz. A separate proton resonance regulator stabilized the static magnetic field.[ref]
Waiting for thermal equilibrium prepared nuclear polarization. Oscilloscope photographs compared glycerol proton echoes at different field homogeneity: echo width changed; echo time followed twice the pulse separation. Receiver saturation obscured early free induction decay; transmitter leakage required care. The delayed signal tested refocusing; diffusion could still attenuate it.[ref]
Hahn analyzed equal pulses as well as noting improved recovery with a second pulse twice as long; the ideal sequence above isolates that latter case.[ref]
Heterogeneous catalysis Created 2026-10-14 Updated 2026-10-18
In heterogeneous catalysis, the catalyst and reactants occupy different thermodynamic phases; a solid surface can furnish the catalytic route.
For example, specify dissolved hydrogen peroxide as the reactant, water as the solvent, and solid manganese dioxide as the catalyst. The solid and liquid chemical solution are separate thermodynamic phases. A dispersion that looks thoroughly mixed still contains those separate phases; mixing distributes contact, rather than making the solid molecularly dissolved. The manganese dioxide peroxide contact demonstration makes the arrival of the liquid at the solid tangible.
A useful surface-route picture is approach from the liquid, adsorption at a catalytic surface site, conversion through surface-bound reaction intermediates, and release of reaction products. [ref] This sequence organizes possible roles; it is not a verified sequence of elementary reactions for every catalyst, or for the particular manganese dioxide demonstration. The net chemical equation supplies neither the binding geometry nor the identity of the working catalytic surface site.
This contact requirement makes surface area and diffusion practically important. More exposed solid can provide more catalytic surface sites, whereas diffusion through the surrounding liquid can limit their supply of reactants. Changing stirring or particle size can therefore change the observed reaction rate without identifying the microscopic rate constant. A recoverable solid is attractive for reuse, but separation alone does not establish retained activity: compare with catalyst deactivation.
Multiply the Hookean dumbbell model configuration equation by and integrate twice by parts, with sufficient decay at large separation. The linear velocity field contributes , the spring drift contributes , and diffusion contributes . Thus the dimensionless conformation tensor obeys , with polymer relaxation time . Its left side is an upper-convected derivative for this gradient convention.
Let be the unnormalized conformation tensor, the velocity gradient, and a symmetric positive hydrodynamic mobility matrix independent of the individual connector. The connector Fokker-Planck equation has drift and diffusion matrix . Integrating its second moment twice by parts gives the displayed upper-convected derivative, . If is a polynomial in , the two matrices commute and the relaxation term simplifies to . No Gaussian distribution assumption is required for this moment closure.
In a planar no-flow model with pure ice, a stable cold brine below freshwater can support an ice layer whose upper phase boundary freezes freshwater while its lower boundary melts into diluted brine. Heat conduction and salt diffusion obey different similarity scales. Local liquidus equilibrium, salt conservation and two Stefan conditions determine both boundary positions. The heat lost into the brine supports net ice growth even while some basal ice melts.
Instantaneous interface kinetics is the limit of negligible kinetic undercooling at a moving solidification interface. Its temperature is set by local liquidus equilibrium and curvature-induced melting-temperature depression, without a further velocity-dependent attachment term. It does not imply that diffusion or heat transport is instantaneous.
Isoelectric point Created 2026-10-13 Updated 2026-10-20
The isoelectric point is the pH at which a specified molecule's mean net electric charge is zero. Different protonation states can still carry nonzero electric charges; their contributions cancel in the average. For the two-step glycine speciation model, this gives . Additional ionizable amino acid side chains require their own chemical equilibrium calculations; averaging every pKa is not a general rule. [ref]
Below its isoelectric point, free glycine has positive mean electric charge and an electric field drives it toward the negative electrode; above its isoelectric point, the mean electric charge is negative and the direction reverses. At its isoelectric point the charge-driven contribution vanishes. Diffusion and bulk movement of the chemical solution can still move the sample. Gel electrophoresis applies this charge dependence to separation, although protein size, shape and interactions also affect migration. The pKa values of free glycine do not determine a protein's isoelectric point: its free N-terminus, C-terminus and ionizable amino acid side chains supply a different inventory. [ref]
Isopycnic centrifugation Created 2026-10-13 Updated 2026-10-20
In isopycnic centrifugation, a species concentrates near the position where the liquid's density matches its buoyant density. In a density gradient increasing outward, material displaced inward encounters less-dense liquid and tends to move outward; material displaced outward encounters denser liquid and tends to move inward. Diffusion broadens the resulting band. If the density gradient never reaches the species' buoyant density, a matching equilibrium position need not exist within the vessel.[ref]
The 1958 paper The replication of DNA in Escherichia coli by Matthew Meselson and Franklin W. Stahl (pp. 672–676, Figures 1–4) makes the distinction concrete. Its separation apparatus was a Spinco Model E ultracentrifuge with an Escherichia coli lysate in a cesium chloride density gradient. Ultraviolet wave absorption photographs of the rotating sample were scanned with a recording microdensitometer: peak position reported buoyant density, while signal above baseline reported DNA concentration.[ref]
Cultures received labeled nitrogen; chilled biological cells were collected as centrifugation pellets and disrupted before molecular banding.[ref]
Mixed lysates labeled fully with the nitrogen isotope and unlabeled lysates supplied reference bands; a mixture including intermediate material checked its position relative to those references. Figure 1 followed band formation over time, helping distinguish a migrating distribution from an approximately stationary one. These controls address resolution, calibration and approach to equilibrium. They do not make every ultraviolet-absorbing component identifiable as DNA, or turn a band into a direct image of individual strands. The DNA replication article develops the copying inference; here the key result is a separation by buoyant density, not a claim that heavier-labeled DNA is physically larger.[ref]
Kinetic theory Created 2026-10-06 Updated 2026-10-13
Kinetic theory describes many-particle systems through a distribution on phase space, retaining both position and velocity. In the number-density convention, is the expected number of particles in a small position–velocity region. For identical particles of mass , its integrals give the number density and bulk velocity :
This is not normalized to one: its integral over the whole phase space gives the total expected particle number. Transport moves particles along trajectories, while collision or diffusion terms change the distribution. The Boltzmann equation, Hamiltonian Liouville equation and Fokker-Planck equation describe different physical assumptions, rather than interchangeable collision laws.
An ideal gas is a motivating example: the pressure on a wall comes from microscopic momentum transfers, even when the bulk velocity is zero. The Maxwell-Boltzmann velocity distribution separates this random molecular motion from bulk motion; see ideal gas pressure from wall collisions. The attractive feature of kinetic theory is that a distribution can describe two opposing populations where one average velocity hides almost everything. It also tells us when averaging into a continuum approximation loses essential information: compare the intermolecular mean free path with the apparatus scale using the Knudsen number. [ref]
The Kuramoto-Sivashinsky equation combines destabilizing long-wave growth, stabilizing fourth-order diffusion, and nonlinear transport. Its potential form is , with . It arises as a long-wave phase amplitude equation near loss of the Benjamin-Feir stability condition. On a domain of period , the zero solution's Fourier modes grow at rate , .