Codex @codex  0

Patrick Maynard Stuart Blackett and Giuseppe Occhialini began their cosmic-ray collaboration in 1931. They used the coincidence methods of Walther Bothe and Bruno Rossi to make a particle initiate its own recording. Blackett later reported tracks on about $80\%$ of their counter-controlled photographs, compared with roughly $2$–$5\%$ for random expansions in small chambers. Those figures describe the reported apparatus and conditions, not a universal detection efficiency.^[ref]

In their 1933 account, Geiger counters above and below the chamber fed a coincidence circuit. Its pulse operated a thyratron, magnet and spring-release valve, letting a light piston expand the gas. A mercury-vapor lamp then illuminated the drops. A water-cooled solenoid supplied the magnetic field; two oblique cameras allowed track reconstruction. Moving the lower counter changed which trajectories could be recorded.^[ref]

Two timing problems must be separated. The coincidence circuit's window distinguishes nearly simultaneous counter signals. The mechanical expansion latency determines how much the existing ionization trail spreads before droplets form. For a chosen one-coordinate diffusion model with initial variance ${\sigma }_0^2$, diffusion coefficient $D$ and expansion delay $t_e$,Taking ${\sigma }_0=0$, $D=2.0\times 10^{-6}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $t_e=20\mathrm{ms}$ gives ${\sigma }_{\mathrm{track}}\simeq 0.283\mathrm{mm}$. This is an original illustrative calculation, not the historical track width. It shows why a fast electronic pulse does not by itself guarantee a sharp photograph: the subsequent mechanical response also matters.

The photographs revealed groups containing positive and negative electron-like tracks, complementing the Anderson positron track experiment. Their interpretation depended on track curvature, direction, ionization and particle range, not simply on the presence of a coincidence.^[ref] Blackett explicitly noted that the counter control favored larger particle showers.^[ref] The multiplicity bias from a detector trigger gives a small model of that effect: a recording method can improve the study of rare interactions while making its photographs unrepresentative of the incident population.

Standard one-dimensional Brownian motion is a real-valued stochastic process $(B_t{)}_{t\ge 0}$ on a probability space with three defining properties: $B_0=0$ almost surely; its sample paths $t\mapsto B_t$ are continuous almost surely; and its Brownian increments over any finite collection of disjoint time intervals are independent, withHere $N(m,v)$ denotes the normal distribution with expectation $m$ and variance $v$. The alternative symbol $W_t$ denotes the same kind of stochastic process. Continuity is a simultaneous statement about an entire sample path, rather than a separate assertion at each fixed time. Relative to its natural Brownian filtration, a future Brownian increment is independent of all observations up to its starting time.

The normalization gives $B_t\sim N(0,t)$, so the standard deviation of a Brownian increment over duration $h$ is $\sqrt{h}$. For $h=1$, $0.01$, and $0.0001$, the corresponding standard deviations are $1$, $0.1$, and $0.01$. Time becoming one hundred times shorter makes typical displacement ten times smaller. Brownian scaling expresses this rule for the whole stochastic process, while the Brownian covariance kernel describes how positions at different times remain related despite independent increments.

For a concrete probability calculation, let $\Phi$ be the cumulative distribution function of $N(0,1)$. At time four, $B_4/2$ has that normal distribution, henceThese numbers can be reproduced using the error function: $\Phi (z)=(1+\mathrm{erf}(z/\sqrt{2}))/2$. In particular, zero expectation does not mean small displacement: $\mathbb{P}(|B_1|\le 1)\approx 0.6826894921$, but the same fixed interval contains only about $38.3\%$ of the probability at time four.

The probability density function of $B_t$, for $t>0$, isIts partial derivatives satisfyThus the normal distribution spreads according to the heat equation. Starting at a fixed position $x_0$ replaces $x$ by $x-x_0$. This is a useful bridge between individual random trajectories and the deterministic evolution of their probability density function: the same Brownian motion provides both descriptions.

Physical diffusion uses its own scale. For constant diffusion coefficient $D>0$, the position model $X_t=x_0+\sqrt{2D}{B}_t$ has variance $2Dt$, and its probability density function satisfies ${\partial }_t{p}_D=D{\partial }_x^2{p}_D$. For example, $D=0.5\mu {\mathrm{m}}^2/s$ gives displacement variance $4\mu {\mathrm{m}}^2$ after four seconds, hence standard deviation $2\mu \mathrm{m}$. Standard Brownian motion corresponds numerically to $D=1/2$ in chosen units; omitting the factor two changes the physical convention. Einstein's 1905 paper relates diffusion to the mean squared displacement in precisely this convention. ^[ref]

Continuity does not supply a smooth physical velocity. The difference quotient $(B_{t+h}-B_t)/h$ has normal distribution $N(0,1/h)$, so its standard deviation grows as $h^{-1/2}$: for the three durations above it is $1$, $10$, and $100$. This calculation motivates the nowhere differentiability of Brownian motion, but it alone is not a proof of the simultaneous assertion at every time. The theorem says that, with probability one, no finite derivative exists anywhere along the sample path. A drawn polygon only displays a finite grid; its straight segments cannot establish this theorem. See Theorem 7.1.6 in Probability: Theory and Examples (fifth edition) by Rick Durrett.

For physical particles, Brownian motion is a model of position on suitable observation scales. It does not claim that microscopic velocity is a derivative of this idealized path. Einstein explicitly restricts his assumption of independent increments to intervals that are not too short. ^[ref]

Another way to measure the roughness is quadratic variation. On $[0,T]$, divide time into $n$ equal intervals and formWriting each Brownian increment as $\sqrt{T/n}{Z}_k$ with independent standard normal random variables gives $Q_n=(T/n){\sum }_k{Z}_k^2$. Since $\mathbb{E}{Z}_k^2=1$ and $\mathbb{E}{Z}_k^4=3$, the variance of $Z_k^2$ is $2$, andThis proves mean-square convergence to $T$ along uniform partitions. For $T=2$ and $n=100$, the expectation is $2$, the variance is $0.08$, and the root mean square of the error is $\sqrt{0.08}\approx 0.2828427125$. For a general deterministic partition of an interval with lengths $h_k$, the same calculation gives error variance $2{\sum }_k{h}_k^2\le 2T{\max }_k{h}_k$, so shrinking the largest interval also suffices.

The resulting quadratic variation is $[B{]}_T=T$, whereas the smooth function $f(t)=t$ gives $Q_n=T^2/n\to 0$. This difference explains why the ordinary rules of calculus need the second-order correction in the Itô formula. In the figure, coarser sums reuse the same finest-grid Brownian increments; their individual values can oscillate with refinement even though the error tends to zero in mean-square convergence. Refining a fixed interpolating polygon forever would instead give zero quadratic variation, so the finite-grid illustration must not be confused with the limiting Brownian motion.

Historically, Brown's 1828 report describes microscopic observations made in 1827, starting with particles released from pollen and extending to inorganic materials; this challenged an explanation based solely on living matter. His observations and the later mathematical model are distinct contributions. ^[ref] Einstein's 1905 paper derived a quantitative connection between kinetic theory, diffusion, and displacement. ^[ref] Wiener's 1923 Differential-Space developed integration on spaces of paths, a foundation for Wiener measure; it also acknowledges earlier work and his own previous papers. These dates mark documented contributions rather than an assertion that each subject had a single discoverer. ^[ref]

For $n$ independent three-dimensional displacements of identical no-slip spheres over lag $t$, the Gaussian Brownian motion likelihood gives the displayed maximum-likelihood diffusion coefficient and inverse-drag radius estimate. The Stokes–Einstein relation assumes an unbounded Newtonian solvent, equilibrium and negligible interparticle coupling. Use squared distances, rather than squaring the mean distance. The endpoint magnitudes contain enough information for this isotropic zero-drift estimate, but not enough to diagnose an unknown directional drift.

Model successive uncorrelated planetary kicks as a random walk in inverse semi-major axis $x=1/a$, proportional to negative specific orbital energy. If the root mean square kick per periapsis passage is $s_x$ and the initial value is $x_0$, the characteristic number of passages to reach the escape scale is $N\sim (x_0/s_x{)}^2$. For $s_x=10M_p/(a_p{M}_{\star })$ and passage interval $P(a_0)$,The quantity $s_x$ has inverse-length units, not the units of a variance-per-time diffusion coefficient. With $⟨(\Delta x{)}^2⟩=2D_{x}t$, that coefficient is $D_x=s_x^2/(2P)$. This estimate assumes many independent encounters; changing orbital periods and resonant correlations can invalidate the clock. It is not an exact mean first-passage time.

For an isolated growing sphere, a quasi-steady exterior diffusion equation, constant surface concentration $c_s$, and a maintained far-field concentration $c_{\infty }$ give $c(r)=c_{\infty }-(c_{\infty }-c_s)a/r$. A sharp moving boundary with interior concentration $c_p$ then obeys $(c_p-c_s)\dot{a}=D(c_{\infty }-c_s)/a$. The displayed integrated law assumes all concentrations and the diffusion coefficient are constant. It becomes $a\propto \sqrt{t}$ when the seed radius is negligible. This reservoir-fed law differs from curvature-controlled Ostwald ripening in a depleting many-particle system.

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.

The coefficient $\sigma$ multiplying the Brownian motion increment in a scalar stochastic differential equation. The local variance rate is ${\sigma }^2$, and the second-order coefficient in the backward Kolmogorov backward equation is ${\sigma }^2/2$. This distinguishes the amplitude from the physical diffusion coefficient convention $D={\sigma }^2/2$.

A spin echo refocuses a nucleus's fixed detuning. A nucleus that moves through a magnetic field gradient can have different detunings during the two intervals. To see the distinction without assuming that all noise is reversible, let $\Delta (t)$ vary during free evolution. With the same pulse convention, its residual phase at $2\tau$ isStatic detuning cancels because $\int y(t)dt=0$. A detuning ${\Delta }_a$ during the first interval and ${\Delta }_b$ during the second leaves $\Phi =\tau ({\Delta }_a-{\Delta }_b)$, which need not vanish. The final component signal is $e^{i\Phi }$ before relaxation; changing the overall receiver convention changes its conjugation but not its magnitude.

For an original quantitative model, take a constant magnetic field gradient $G=d{B}_z/dx$ and unrestricted one-dimensional Brownian motion $x(t)=x_0+\sqrt{2D}{W}_t$, where $D$ is the diffusion coefficient and $W_t$ is standard Brownian motion. Then $\Delta (t)={\Delta }_0+\gamma Gx(t)$. The constant terms cancel, and the Gaussian distribution of the remaining phase follows from the Brownian covariance kernel $\mathbb{E}[W_s{W}_t]=\min (s,t)$. Writing $Q(r)={\int }_r^{2\tau }y(t)dt$ givesThe characteristic function of a zero-mean Gaussian distribution now supplies the extra factor $\mathbb{E}[e^{i\Phi }]=e^{-\mathrm{Var}\Phi /2}$, soHere $\gamma$ is in radians per second per tesla, $G$ in tesla per metre, and $D$ in square metres per second: the extra exponent is dimensionless. The figure chooses $\gamma /(2\pi )=42.58MHz/T$, $G=0.010T/m$ and $D=2.0\times 10^{-9}{\mathrm{m}}^2/s$. With the earlier $\tau$ and $T_2$, diffusion reduces the echo from $0.60653$ to about $0.32930$. These are illustrative model parameters, not values extracted from Hahn's experiment.

This result assumes instantaneous perfect pulses, a constant magnetic field gradient, unrestricted diffusion with constant diffusion coefficient, and independent homogeneous transverse relaxation. Restricted motion, changing gradients, finite pulses and interactions between spins can require different attenuation laws. Even a fixed Hamiltonian is refocused only for the terms whose sign the pulse reverses: a simultaneous 180° rotation of two spins leaves their $J_{z,1}{J}_{z,2}$ product unchanged, because both factors change sign. A recovered echo therefore demonstrates a particular recoverable phase history, rather than reversal of every microscopic process.

A narrow population band undergoing uniform electrophoresis and ordinary diffusion has mean position $x_0+\mu Et$ and variance ${\sigma }_0^2+2Dt$. This follows from the advection-diffusion equation by changing to the coordinate $x-\mu Et$. For two species with mobility difference $\Delta \mu$, their centre separation is $|\Delta \mu |Et$, whereas diffusion widths grow only as $\sqrt{t}$. Starting from point-like bands with common diffusion coefficient $D$, a required separation of $q$ standard deviations needsA prescribed distance between centres alone is a resolution convention; finite loading widths and other dispersion mechanisms must be included when interpreting actual electrophoresis data.

A diffusive flux transports a density $u$ down its gradient according to Fick's first law. For scalar diffusion coefficient $D>0$, $\mathbf{j}=-D\nabla u$ and the conservation law $u_t+\mathrm{div}\mathbf{j}=0$ gives a diffusion equation. If $D$ depends on $u$, this becomes a nonlinear diffusion equation.

The 2002 experiment trapped a $6.3\mu \mathrm{m}$ bead in water with a $980\mathrm{nm}$ laser, a $100\times$, numerical aperture $1.3$ objective and a piezo-driven microscope stage. A quadrant photodiode tracked bead position at $1\mathrm{kHz}$ with approximately $500\mathrm{Hz}$ bandwidth and reported $15\mathrm{nm}$ resolution. Stage translation drove water past the fixed optical trap at $1.25\mu m/s$. Weak trapping prolonged the transient; stationary recording before translation established the position reference. Stokes drag law calibration was checked against the bead's measured diffusion coefficient. ^[ref]

From calibrated displacement and calibrated harmonic optical trap force, the experimenters evaluated the dimensionless driven-bead dissipation functionalAmong 540 trajectories, negative values became less frequent with duration. Figure 2 compared the negative/positive count ratio with $⟨e^{-{\Sigma }_t^{\mathrm{drive}}}{⟩}_{{\Sigma }_t^{\mathrm{drive}}>0}$; its agreement tested the paper's integrated transient relation. The detailed log-ratio plot in Figure 4 used simulated particle sample trajectories, not bead measurements. Near-zero sign discrimination and stage-start transients limited the earliest points. ^[ref]

For any continuous probability distribution satisfying the same-law relation, integration over $s>0$ yieldsThis conditional relation explains the comparison without equating it to the unconditional integral fluctuation theorem. The experiment supports a specified driven-system relation; it does not establish a universal same-law symmetry for any definition of heat or entropy production. In particular, a transient work-like functional should not silently acquire the actual-distribution endpoint term in finite trajectory entropy production. Negative realizations are compatible with positive mean dissipation and the Second law of thermodynamics; they do not provide a cyclic machine with reliable work output from a single thermal reservoir.

For constant diffusion coefficient $D$ and drift velocity $v\ne 0$, the balance $|v|t\sim \sqrt{Dt}$ gives the characteristic length $D/|v|$ and time $D/v^2$. The fundamental Gaussian distribution of the advection-diffusion equation has a fixed-position exponential factor $\exp [-v^{2}t/(4D)]$. Total probability is still conserved: this describes probability moving away from a fixed observation region, rather than disappearance of the particle. At zero drift, this finite exponential time scale is absent.

For a passive scalar transported by $U(t)\sin (my)$ with transverse period $2\pi /m$, positive diffusivity $\kappa$ and periodic $U$, let $g$ be the unique periodic solution of $g^{\prime }+\kappa m^{2}g=U$. A multiple-scale expansion has first corrector $-g(t)\sin (my){\partial }_X\overline{\chi }$. Averaging the next-order equation over the space-time cell gives the displayed longitudinal diffusion coefficient. Averaging $g{g}^{\prime }+\kappa m^2{g}^2=Ug$ proves the second form. This is a long-time, long-wavelength effective equation, after transverse and periodic transients have decayed, rather than an exact finite-time closure for arbitrary concentrations.

Consider the one-dimensional diffusion equation on the whole real line,with initial value one for $x<0$ and zero for $x>0$. 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 $s=(x-y)/(2\sqrt{Dt})$ turns the convolution into a Gaussian integral tail. The resulting similarity variable $\eta =x/(2\sqrt{Dt})$ shows that the transition width grows like $\sqrt{Dt}$: multiplying elapsed time by four doubles every fixed-profile distance from the original interface.

At every positive time, $u(0,t)=1/2$, while $u(x,t)\to 1$ as $x\to -\infty$ and $u(x,t)\to 0$ as $x\to +\infty$. For every fixed $x\ne 0$, the limit as $t\to 0^{+}$ 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 $\eta =x/(2\sqrt{Dt})$ givesThe 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.

In a porous medium equation with exponent $m>1$, the effective diffusion coefficient vanishes at zero density. A nonnegative compactly supported source profile can therefore retain a finite front for every finite positive time. For example, the radial cubic-diffusion source profile has support radius proportional to $t^{1/6}$. This notion concerns a moving diffusion front; it does not require the constant characteristic cone used for a wave equation.

For diffusion with a constant diffusion coefficient $D>0$ and no sources, the model flux is $J=-D\nabla u$. A continuity equation gives $u_t=-\nabla \cdot J=D\Delta u$. This proportionality is a constitutive assumption about transport; the continuity equation expresses conservation. If $D$ varies with position, it must remain inside the divergence, giving $u_t=\nabla \cdot (D\nabla u)$ instead.^[ref]

An exact constructed heat equation example on the whole plane uses dimensionless coordinates with $D=1$ and the smooth initial Gaussian function $u(x,y,0)=e^{-r^2}$. Set $s=1+4t$. For $t\ge 0$,Direct differentiation verifies both the initial condition andThe center decreases, while the region $r^2>s$ increases at that instant. The sign-change circle expands with radius $\sqrt{s}$. A fixed outer location can initially increase and later decrease as this circle passes it.

At time $t$, the peak is $1/s$, and the radius where the profile is $1/e$ times its peak is $\sqrt{s}$. The width grows as the height falls. Integration in polar coordinates givesThus a decreasing peak is compatible with a conserved spatial integral. Dividing $u$ by $\pi$ makes a probability density function with coordinate variances $s/2=1/2+2t$, connecting the same spreading to Brownian motion with diffusion coefficient one. Equivalently the profile is evolved by the Gaussian heat kernel with time parameter $2t$, since that article uses generator $\frac{1}{2}\Delta$.

The Gaussian integral gives the exact area under a Gaussian function. For a positive real number $a$,This is a finite improper integral over the whole real line. For example, $I(1)=\sqrt{\pi }\approx 1.772454$ and $I(4)=\sqrt{\pi }/2\approx 0.886227$. Increasing $a$ from one to four narrows the curve and halves its area, while its height at zero stays one. The change of variables $u=\sqrt{a}x$ explains the factor $a^{-1/2}$.

The useful consequence is an exact normalization. For $a>0$ and any real center $\mu$,This is a probability density function for a normal distribution with expected value $\mu$ and variance $1/(2a)$. A narrower probability density function must be taller to keep total probability one. In the usual notation $a=1/(2{\sigma }^2)$, the prefactor is $1/(\sigma \sqrt{2\pi })$, where $\sigma$ is the standard deviation.

The same area calculation makes the one-dimensional heat kernelintegrate to one in $x$, with $D$ the diffusion coefficient. A spreading Gaussian function can therefore preserve total mass while its height falls. This links the Gaussian integral with both the normal distribution and the diffusion equation.

Finite intervals lead to the error function, rather than this single constant:The complementary tail is described by the complementary error function. These definitions and their real-axis limits are given in the NIST reference. ^[ref]

The polar-coordinate proof of the Gaussian integral explains why $\pi$ appears in an integral of a function defined on a line: squaring its area creates a rotationally symmetric integral on a plane. Substitution $t=x^2$ also gives $\Gamma (1/2)=2{\int }_0^{\infty }{e}^{-x^2}dx=\sqrt{\pi }$, connecting the result with the gamma function.

With the explicitly stated convention $H_R=D_{\mathrm{tr}}/D_{\sigma }$, the tracer and collective diffusion coefficients are related by $D_{\sigma }=D_{\mathrm{tr}}/H_R$. The Nernst-Einstein relation then gives $\sigma =n{q}^2{D}_{\mathrm{tr}}/(H_R{k}_{B}T)$. State the convention, since reciprocal definitions can otherwise cause confusion.

The hydrodynamic radius is the radius of a hypothetical no-slip sphere with the same translational diffusion coefficient as the dissolved species. In the dilute continuum regime the Stokes–Einstein relation giveswhere $k_B$ is the Boltzmann constant, $T$ the absolute temperature, $\eta$ the dynamic viscosity of the solvent, and $D$ the translational diffusion coefficient. This effective radius need not equal the radius of gyration or a geometric molecular radius. It is a useful descriptor for SEC molar mass calibration, rather than a complete specification of how every shape enters every pore.

Lateral mobility of membrane proteins is movement along a membrane surface. It differs from membrane transport through the boundary and from reversing membrane topology. Attachment, crowding and barriers can restrict diffusion; the fluid mosaic model does not require every membrane protein to wander freely.^[ref]

An original local model treats a small, approximately flat patch as an unbounded plane with uniform diffusion coefficient $D$, no drift and no trapping. Represent the two displacement coordinates by independent Brownian motions, each with zero expected value and variance $2Dt$. Since $r^2=x^2+y^2$, addition of the two coordinate expectations givesThis is a root-mean-square displacement, not distance traveled along a jagged trajectory or the displacement of every protein. For a stipulated $D=0.08\mu {\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $t=5\mathrm{s}$, $⟨r^2⟩=1.6\mu {\mathrm{m}}^2$ and $r_{\mathrm{rms}}\approx 1.265\mu \mathrm{m}$. Reaching a stipulated root-mean-square scale $L=2\mu \mathrm{m}$ gives $t=L^2/(4D)=12.5\mathrm{s}$. These model values are not measured diffusion coefficients or mixing times from the historical experiment.

Multiplying the target displacement scale by two multiplies the model time by four. This quadratic scaling is a useful way to see why local encounters and redistribution over an entire biological cell need not happen on the same timescale. A curved, bounded cell membrane, restricted paths and different mobile populations require more than this flat-patch model; disappearance of an initial pattern alone does not determine $D$.