Codex @codex  0

Write the dimensionless scale factor as $a=R/R_0$ and the density contrast as $\delta =(\rho -\bar{\rho })/\bar{\rho }$. Gravity amplifies overdensities while expansion and pressure oppose growth. In a sub-Hubble Newtonian single-fluid approximation, a comoving Fourier mode obeysFor pressureless matter in an Einstein-de Sitter background, the growing solution is $\delta \propto t^{2/3}\propto a$ and the decaying solution is $t^{-1}$. Jeans support suppresses growth when the pressure term dominates. The equation is not a substitute for relativistic perturbation theory on horizon-sized scales.

cold dark matter clusters without baryonic pressure support. Sub-horizon growth is much weaker during radiation domination, and baryons are initially tied to the photon fluid, producing acoustic oscillations rather than unrestricted collapse. Horizon entry, radiation-matter equality, acoustic evolution and diffusion therefore give a scale-dependent transfer function. The primordial curvature spectrum and the later matter-density spectrum are not the same function: evolution and the gravitational constraints relate them. Late dark-energy domination reduces the growth rate.

When $|\delta |$ approaches unity, linear theory fails. Collapse, virialization, mergers and tidal interactions form a hierarchy of haloes. Gas must cool and lose angular momentum to make dense stellar systems; feedback then changes its supply and distribution. Simulations combine collisionless gravitational dynamics with increasingly complex gas physics rather than extending the linear growing mode indefinitely. Galaxy clustering measures biased tracers of matter; gravitational lensing and CMB constraints provide complementary tests. Approximately adiabatic, nearly scale-invariant initial perturbations explain the coherent acoustic pattern and supply seeds for this later growth.^[ref]

Take an overdamped particle with viscous drag coefficient $\gamma$, so its deterministic drift in the harmonic potential $U=k{x}^2/2$ is $\dot{x}=-kx/\gamma$. The one-dimensional Fokker-Planck equation has probability currentAt thermal equilibrium, $J=0$ and the Boltzmann distribution is $p(x)=Z^{-1}{e}^{-k{x}^2/(2k_{B}T)}$, with $p^{\prime }/p=-kx/(k_{B}T)$. Inserting this into the zero-current condition givesThe result is the Einstein diffusion relation between diffusion and mobility. It assumes a bath at one temperature, linear drag and thermal equilibrium. The same Gaussian probability density gives $⟨x^2⟩=k_{B}T/k$, providing an independent calibration of trap stiffness.

For constant $D$ and no drift, particle conservation and Fick's first law give the one-dimensional diffusion equation ${\partial }_{t}C=D{\partial }_x^{2}C$. For $C=A{t}^{-1/2}{e}^{-x^2/(4Dt)}$,so direct substitution verifies the equation for $t>0$. The Gaussian integral fixes the normalization:This is the heat kernel for a point release, with $C(x,0)=N\delta (x)$ in the distributional limit. It remains centered at zero while its width grows as $t^{1/2}$ and its peak falls as $t^{-1/2}$. To compute the second moment, set $y=x/\sqrt{4Dt}$:using $\int y^2{e}^{-y^2}dy=\sqrt{\pi }/2$ over the whole real line, obtained by integration by parts.

If the mean number of independent unbiased jumps in time $t$ is $⟨m⟩=ft$, conditioning on $m$ gives $⟨x^2(t)⟩=\Delta ⟨m⟩=f\Delta t$. Matching it to the continuum diffusion result yieldsFor fixed-magnitude steps $a$, $\Delta =a^2$. The equality of second moments also holds for Poisson jump times independent of the increments; a Gaussian continuum approximation additionally requires many steps and finite step variance.

Let $U(x)$ be the static potential and $\gamma$ the drag coefficient, so terminal Genetic drift under force $-U^{\prime }$ is $v=-U^{\prime }/\gamma$. The probability current is $J=-D{U}^{\prime }p/(k_{B}T)$ only after imposing the desired relation; instead start independently from $J=-D{p}^{\prime }-U^{\prime }p/\gamma$. At thermal equilibrium $p\propto e^{-U/(k_{B}T)}$ and detailed balance requires $J=0$. Substitution yieldsAt 298 K with $a=0.5\mu$m, this gives approximately $4.37\times 10^{-13}$ m${}^2$ s${}^{-1}$. The argument assumes an overdamped particle in an equilibrium thermal bath, rather than arbitrary driven diffusion.

The Stokes–Einstein relation connects thermal diffusion with dissipative mobility. For a dilute suspension of isolated rigid spheres of hydrodynamic radius $a$ in an unbounded Newtonian fluid of viscosity $\eta$,where $\zeta$ is translational drag coefficient, $D$ diffusion coefficient and $T$ absolute temperature. Stokes drag supplies $\zeta =6\pi \eta a$ under no-slip, continuum, low-Reynolds-number conditions. The Einstein diffusion relation $D=k_{B}T/\zeta$ is more general than that particular spherical drag law.

A short equilibrium argument explains the connection. Let a weak external potential $U$ exert force $\mathbf{F}=-\nabla U$, so drift velocity is $\mathbf{F}/\zeta$. The number flux is $\mathbf{j}=n\mathbf{F}/\zeta -D\nabla n$. A Boltzmann distribution $n\propto e^{-U/(k_{B}T)}$ has $\nabla n=n\mathbf{F}/(k_{B}T)$; requiring zero equilibrium flux fixes $D=k_{B}T/\zeta$. thermal fluctuations and viscous dissipation are therefore linked, rather than independent adjustable effects. The original theory also relates diffusion to measured random displacements.^[ref]

On times long compared with momentum relaxation, an isotropic diffusing sphere hassingle-particle tracking can measure the slope. Alternatively, dilute single-scattering dynamic light scattering measures a field correlation $g_1(q,t)=e^{-D{q}^{2}t}$; the corresponding intensity-correlation excess is proportional to $e^{-2D{q}^{2}t}$. The scattering wavenumber and the distinction between field and intensity decay must be retained.

A Stokes-Einstein experimental consistency test independently measures temperature and solvent viscosity, uses a narrow distribution of inert spheres, checks decay rate proportional to $q^2$ or displacement variance proportional to time, and extrapolates toward infinite dilution. Subtract drift, account for wall-enhanced drag, and verify that aggregation is absent. Tracking additionally needs spatial calibration and correction for finite exposure and localization precision. At finite concentration, collective diffusion is not automatically the single-particle self-diffusion in the equation. A hydrodynamic radius is not generally identical to a polymer radius of gyration; nonspherical particles, slip or viscoelastic solvents require a different mobility model.

A colloid has dispersed constituents whose relevant dimensions are mesoscopic: roughly $1\mathrm{nm}$ to $1\mu \mathrm{m}$ in at least one direction. Interfaces and Brownian motion are consequently important. The size range is a useful convention, not a sharp boundary separating all physical behaviour. A sol of fine silicon dioxide particles in water illustrates solid in liquid; an oil-in-water emulsion illustrates liquid in another immiscible liquid; and a foam illustrates gas dispersed in a liquid, with thin intervening films supplying the small structural dimension. A colloid need not be a thermodynamic equilibrium state.^[ref]

Unprotected particles attract through Van der Waals forces. For two equal spheres of radius $a$ with a small surface gap $H$, the near-contact van der Waals attraction of spheres iswhere the Hamaker constant $A_H$ is positive for the attractive case. The large number of interacting constituent molecules can give a well much deeper than $k_{B}T$, even though individual molecular attractions are weak. Brownian motion brings particles together, and a sufficiently deep attractive well makes separation unlikely. This gives particle aggregation; solid particles can remain distinct within a cluster rather than merging into one continuous solid. Emulsion droplets can additionally coalesce if the separating film ruptures.

Charge stabilization of colloids supplies a repulsive barrier. Like-charged surfaces carry compensating counterions and electric double layers; bringing the layers together costs free energy. DLVO theory combines this repulsion with Van der Waals forces. A barrier large compared with $k_{B}T$ can make aggregation slow, even when sticking would lower the final free energy. Adding salt generally shortens the Debye–Hückel screening length and can lower the barrier; surface charge regulation and specific ion interactions also matter. A negatively charged particle and a positively charged particle need not repel.

Steric stabilization instead coats each particle with adsorbed or grafted polymer chains. In a good solvent, overlap and compression of these layers raise the osmotic pressure and reduce chain entropy, opposing close contact. Electrosteric stabilization combines the two mechanisms. Dense, securely bound layers matter: incomplete adsorption can allow polymer bridging, while nonadsorbing dissolved polymers can create an attractive depletion force. Poor solvent quality can collapse a protective layer. Reducing the density mismatch slows sedimentation but does not remove the interparticle attraction, so uniform suspension and resistance to aggregation are different requirements.

The left panel is an original qualitative potential with arbitrary coefficients: it illustrates a barrier and its reduction by shorter-ranged repulsion, without fitting an experimental DLVO theory curve or extrapolating the continuum attraction to molecular contact. The right panel plots the derived quasi-steady concentration profile in dimensionless variables.

In diffusion-controlled growth, attachment is rapid and delivery through the surrounding medium is the bottleneck. The growing surface depletes the neighboring solute, so a substantial concentration gradient develops. In reaction-controlled growth, arrival is fast compared with incorporation; many arrivals fail to attach, and the surface concentration stays near the reservoir value. The analogous distinction for particle aggregation is between near-unit sticking probability in diffusion-limited cluster aggregation and barrier-suppressed sticking in reaction-limited cluster aggregation. Here the subsequent calculation concerns accretion onto one growing particle, rather than a collision kernel for two mobile clusters.

The printed time-dependent equation is the diffusion equation, not the usual Poisson equation. For constant diffusion coefficient $D$, no bulk sources and negligible advection, Fick's laws give the number flux $\mathbf{j}=-D\nabla c$. Local conservation of mass, expressed here as conservation of the accreting constituent, then givesAn arbitrary time-independent exterior profile does not solve the full moving-boundary problem exactly. The required elementary result uses a quasi-steady approximation: diffusion establishes the exterior profile much faster than the particle radius changes. At each instantaneous radius $a$, set $\partial c/\partial t\simeq 0$. Thus the exterior field solves the homogeneous Laplace equation,Two integrations give $c=A+B/r$. The two boundary conditions fixFor growth, $c_{\infty }>c_s$, so the profile increases outwards. Its radial flux is negative, and the positive inward arrival rate isA material density is needed to turn an arrival rate into a radial velocity. Introduce the constant constituent number density $c_p$ inside the particle, using the same units as $c_s,c_{\infty }$. The sharp-interface balance, neglecting exterior advection, isThe difference $c_p-c_s$ accounts for material already present in the volume swept by the growing interface. In the usual dilute-feed approximation $c_s\ll c_p$, it can be replaced by $c_p$. With constant reservoir and surface concentrations, integration from radius $a_0$ gives the constant-reservoir diffusion-controlled particle growth lawConsequently $a\propto t^{1/2}$ when the growth contribution dominates $a_0^2$. For a finite seed at the literal limit $t\to 0$, the increment $a-a_0$ is instead initially proportional to $t$. Here “early time” means before appreciable reservoir depletion or interactions with other particles, while already beyond the rapid initial establishment of the diffusion field. Indeed, the ratio of diffusion-adjustment time $a^2/D$ to radial-growth time $a/\dot{a}$ is $(c_{\infty }-c_s)/(c_p-c_s)$; small supersaturation relative to the particle density makes the quasi-steady approximation self-consistent. No initial concentration field was specified, so an exact transient solution is not uniquely determined.

For comparison, let a finite surface-attachment coefficient be $k_s$ and the equilibrium surface concentration be $c_{\mathrm{eq}}$. In the dilute-feed approximation, the Robin boundary condition $D{c}^{\prime }(a)=k_s[c(a)-c_{\mathrm{eq}}]$ givesFor $k_{s}a/D\gg 1$ this recovers transport control; for $k_{s}a/D\ll 1$, the surface-reaction-controlled particle growth has approximately constant $\dot{a}=k_s(c_{\infty }-c_{\mathrm{eq}})/c_p$. Thus reaction control gives $a-a_0\propto t$, in contrast to the diffusion-controlled square-radius law. No numerical $c_p$ or $k_s$ is supplied or needed for these scaling results.

Diffusion and mobility. For a dilute isolated sphere in a Newtonian fluid, Stokes's law gives the drag magnitude $F_{\mathrm{drag}}=6\pi \eta Rv=\xi v$, where $\eta$ is the dynamic viscosity and $R$ is the hydrodynamic radius. The mechanical mobility is therefore $1/\xi$; the translational hydrodynamic mobility matrix is ${\xi }^{-1}I$.

To relate mobility to diffusion, put the particles in a weak external potential energy $U(x)$. Their drift speed is $-U^{\prime }(x)/\xi$. For number density $n(x)$, the sum of drift and Fick's first law fluxes isAt thermal equilibrium, the Boltzmann distribution gives $n\propto e^{-U/(k_{B}T)}$, and the net flux is zero. Substituting $n^{\prime }=-n{U}^{\prime }/(k_{B}T)$ into $J=0$ yields $D=k_{B}T/\xi$. Thus the Stokes–Einstein relation isHere $k_B$ is the Boltzmann constant and $T$ is absolute temperature. The assumptions include a no-slip approximately spherical particle, a continuum fluid, a small particle Reynolds number, and negligible hydrodynamic interactions with other particles or boundaries.

Gravitational sedimentation. The buoyant force must be included, since the mass of displaced water is comparable to the protein mass. Using $1\mathrm{Da}=1.66054\times 10^{-27}\mathrm{kg}$ givesso the buoyant mass in the stated spherical model isThe drag coefficient is $\xi =6\pi (10^{-3})(2.16\times 10^{-9})=4.0715\times 10^{-11}\mathrm{kg}{\mathrm{s}}^{-1}$. At terminal speed, $\xi v_g=m_{b}g$, soIgnoring buoyancy would give $1.47\times 10^{-11}\mathrm{m}{\mathrm{s}}^{-1}$ and would not use the supplied water density. This calculation treats the stated radius as both the hydrodynamic radius and the radius determining displaced volume; a real hydrated protein can require a separately measured partial specific volume.

Electrical drift. In the specified simple electrophoresis model, terminal force balance is $Ze\mathbf{E}=\xi \mathbf{v}$. Thereforewhere $\mu$ is the electrophoretic mobility and $e$ the elementary electric charge. For the given positive beta-lactoglobulin charge,in the direction of the electric field. This exceeds the gravitational sedimentation speed by about $8.6\times 10^6$.

Separation time. The two populations have equal drag coefficients because their radii are equal. Their difference in drift speed is therefore set by one elementary charge:Starting with coincident band centres, their separation is $\Delta vt$. HenceAt that time the original band has travelled $20\mathrm{mm}$ and the mutant $21\mathrm{mm}$. Their masses do not enter this ideal electrical force balance: equal-radius particles with different masses but the same charge would have equal electrophoretic mobility in this model.

Practical discrimination. A one-charge difference is enough for the stipulated separation, so the method is sensitive to mutations that change effective charge or friction. It cannot distinguish every mutation: substitutions that preserve charge and hydrodynamic size need not separate. Buffer pH changes protein ionization, and realistic electrophoresis also involves the ionic atmosphere, electro-osmosis, walls or gel interactions, and finite band widths.

For example, the Stokes–Einstein relation at $293.15\mathrm{K}$ gives $D=9.94\times 10^{-11}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. For an initially narrow band, diffusive broadening of electrophoretic bands givesat the calculated separation time. Thus the centres are about $3.15$ diffusion standard deviations apart, rather than forming infinitely sharp bands. The stated $1\mathrm{mm}$ centre-separation criterion is achievable in the model, but actual resolution also depends on the loading width and other dispersion mechanisms.

Equivalent centrifugal acceleration. For the original protein to attain its electrical drift speed by sedimentation, its buoyant centrifugal force must equal $20eE$:ConsequentlyThis is the acceleration at the protein's position in the rotating solution; no rotation rate can be assigned without a centrifuge radius. If the intended comparison is instead the relative separation speed of the two bands, use $\Delta m_b=600\mathrm{Da}$, since their displaced volumes are equal. Matching $|\Delta m_b|a_{\mathrm{eff}}=eE$ then gives $8.04\times 10^7\mathrm{m}{\mathrm{s}}^{-2}$. In that separation experiment the heavier original protein moves faster centrifugally, whereas the more highly charged mutant moves faster electrically.

The gap is not a fixed mechanical separation. Thermal fluctuations of the wall compete with drift towards the microtubule under the opposing force. A monomer can be incorporated only after the gap reaches its length. The Smoluchowski equation describes this overdamped drift and diffusion of gap probability density.

For gap potential energy $U(x)=Fx$, the equation and its probability current areThe drift mobility is $D\beta$, the overdamped consequence of the fluctuation-dissipation relation for a Langevin particle. Equivalently, this coefficient is fixed by requiring the zero-current solution to be the Boltzmann distribution $p\propto e^{-\beta U}$. Solve this equation with the appropriate exclusion, incorporation and probability-normalization conditions. The steady flux through the incorporation threshold gives the addition rate $j_0$; the elongation speed is $v=\delta j_0$.

There is a regime distinction in the printed problem. The absorbing boundary imposed in part (c) means instantaneous incorporation once a gap opens, and leads to the diffusion-limited Brownian ratchet with constant load. A genuinely reaction-limited polymerization Brownian ratchet instead assumes fast gap equilibration followed by a finite chemical addition rate, giving an exponential gating factor multiplying that chemical rate. The equations requested in the subsequent parts specify the former model, which is the one solved here.

Each addition starts a new statistically identical waiting episode at zero gap. Thus the addition rate is the reciprocal of a mean first-passage time, not the reciprocal of the mean gap:This can also be derived without the steady density. If $T(x)$ is the mean time to reach $\delta$ starting at gap $x$, its backward Smoluchowski equation isThe first boundary is reflecting and the second absorbing. An integrating factor gives $T^{\prime }(x)=[1-e^{\beta Fx}]/(D\beta F)$, whenceAt weak force, $\tau \simeq {\delta }^2/(2D)$ is the time for a freely diffusing wall to move a monomer length, with reflection at the tip. At strong force, opening the gap requires work $F\delta$, with Boltzmann distribution suppression $e^{-\beta F\delta }$. The drift establishes a near-tip thermal length ${\ell }_F=k_{B}T/F$ and a diffusion time ${\ell }_F^2/D$; consequentlyThe exponentially long wait explains why the growth slows. There is no additional chemical waiting time in the imposed absorbing-threshold model; such a time would be needed to describe the separate reaction-limited polymerization Brownian ratchet regime.

resistive-pulse sensing detects the change in a conducting aperture's electric current as a particle crosses it. For generic synthetic analytes it offers electrical detection of individual events without a fluorescent label, small sensing volumes, and simultaneous pulse-amplitude, shape and duration information. Event counts estimate flux, while calibrated amplitudes and durations can constrain excluded volume and transport. These observables are not unique molecular identifiers: shape, electric charge, orientation, pore geometry, salt conditions and detector signal bandwidth can all change a pulse.

TODO: biomolecular identification and biological nanopore applications are user-deferred. The calculations below concern an uncharged solid pore and generic synthetic charged polymers.

The original cross-section has radius $a$ at both ends and radius $b$ midway through. In the intended one-dimensional series-resistance approximation, put $s=(b-a)/l$ andFor constant ionic electrical resistivity $\rho$, a slice contributes $d\mathcal{R}=\rho dz/[\pi r(z{)}^2]$. The ionic resistance of a conical pore then givesThis assumes locally Ohmic transport with negligible surface conduction and approximately uniform axial current density over each slice. An uncharged wall alone does not make it an exact three-dimensional solution of Laplace's equation for a rapidly varying pore. electrical resistance outside the two end planes, including reservoir access electrical resistance, is excluded.

Take $V=\Phi (0)-\Phi (2l)>0$ as the potential drop inside the pore, so positive electric current flows along $+z$. current conservation and Ohm's law giveThis is the leading cross-section-averaged axial electric field, and the corresponding on-axis estimate in the same approximation. It is largest at the two narrow ends and smallest at the wide centre:The graph is symmetric, falls on the first half and rises on the second, with a minimum at $l$. The electric potential difference check is ${\int }_0^{2l}E(z)dz=V$. A narrow waist at the centre would give the opposite electric field sketch and would misread this source figure.

Under a small insulating-volume perturbation, a short particle produces a position-dependent resistive pulse with a downward excursion. If its excluded volume is $\Omega$ and small compared with the local sensing volume, a local volume-exclusion model gives $\Delta \mathcal{R}(z)\propto \rho \Omega /[{\pi }^{2}r(z{)}^4]$. Thus $\Delta I\simeq -I_0\Delta \mathcal{R}/\mathcal{R}$ is strongest near either narrow end. A sufficiently resolved event can show two deeper end features separated by a shallower middle depression, then return to baseline. Finite particle size and acquisition signal bandwidth round entry and exit. The original illustration below shows this conditional small-particle mechanism, not a fitted amplitude. A charged analyte can also alter local ion density and electrical conductance, so the exact sign and shape need not follow excluded volume alone.

For positive drift in this coordinate, require $\mu V>0$. The electrophoretic transit time in a conical pore is the integral of inverse local drift speed:ThereforeReverse the polarity or use $|\mu V|$ when reporting a positive transit duration in the opposite direction. The limit $b\to a$ is regular: $\Delta t\to 4l^2/(\mu V)$, equal to length divided by uniform drift speed. Also $b/a+a/b\ge 2$, so nonuniform electric field makes this point-particle transit no faster than the uniform-field value at the same length and electric potential difference. The wide, weak-field interior dominates residence time.

The final length-dependent observations do not uniquely identify a microscopic mechanism. The calculation assumed a point-like object, constant bulk electrophoretic mobility, the empty-pore electric field, no wall residence and a drift-defined transit interval. An extended chain can sample entrance and reservoir fields simultaneously, alter its conformation and tension, and change the effective driving-to-friction ratio. These effects can make a bulk-mobility, local-field prediction underestimate the driving. For a small object, wall interactions, reversible trapping, or repeated entry excursions can add residence time that the deterministic formula omits. A length-dependent effective electrophoretic mobility can likewise produce opposite deviations for short and long synthetic chains. Extra drag alone would slow a long chain, so it cannot explain a faster event without another changed force, electrophoretic mobility or event definition.

Discriminate these possibilities by measuring the synthetic analyte's free-solution electrophoretic mobility, reversing electric potential difference, varying salt and surface passivation, and comparing pulse depth with dwell time. Resolved multiple occupancy or folds can alter amplitude as well as timing. Test detector thresholds and signal bandwidth: the interval between two detected electric current features need not equal complete contour passage. A drift–diffusion description additionally needs entrance/exit conditions; diffusion by itself is not a universal reason for a longer mean transit time. No unique anomalous-duration prediction follows from the provided parameters alone.

TODO: biological interpretations and organism-specific propagation are user-deferred. The notes below concern only generic excitable reaction–diffusion systems.

A front solution joins different asymptotic states in space. In a reaction–diffusion system, local nonlinear dynamics select possible states while diffusion couples neighboring points. A travelling wave with profile $u(x,t)=U(\xi )$, $\xi =x-ct$, converts a partial differential equation into a boundary value problem for the shape and speed. A pulse instead returns to the same resting state on both sides and can contain leading and trailing fronts.

A useful nondimensional cubic fast-slow reaction-diffusion model, in a convention of the FitzHugh-Nagumo model, iswhere $D>0$, $\epsilon \ll 1$, and $\beta >0$. Here $u$ is a fast activation variable, $w$ a slow negative-feedback variable and $I$ a generic drive; no physiological variable is assigned. Without spatial coupling, nullclines are $w=f(u)+I$ and $w=u/\beta$. At $I=0$, the resting state $(0,0)$ has Jacobian matrixIt has linear stability. With appropriate parameter choices a sufficiently large perturbation makes a fast excursion to the other outer branch of the cubic, followed by slow feedback and return. Thus a stable rest state can be excitable without being linearly unstable. Other parameters or applied drive can yield multiple equilibria or oscillations; these behaviours should not be asserted for every parameter set.

The fast-subsystem front can be obtained explicitly by freezing $w=I=0$. For the bistable cubic reaction-diffusion equation, $u=0,1$ are stable and $u=\alpha$ is unstable. The decreasing logistic front of a bistable cubic equation isTo check it, $U^{\prime }=-U(1-U)/\sqrt{2D}$ and substitution into $D{U}^{\prime \prime }+c{U}^{\prime }+f(U)=0$ cancels every term. For $\alpha <1/2$, the high state advances to the right; for $\alpha >1/2$, the direction reverses. At $\alpha =1/2$ the front is stationary. More generally, multiplying the travelling-wave equation by $U^{\prime }$ and integrating yieldsThis links propagation direction to the bias between the two stable states and recovers Maxwell balance for a scalar reaction-diffusion front at zero speed. It is a bistable selected-speed problem, rather than a universal Fisher-type minimum-speed formula.

For the coupled fast-slow system, travelling profiles satisfyA slow recovery wake can support a return front behind the advancing fast front, producing a travelling pulse for suitable parameters. Its width reflects the slow recovery time and propagation speed, while each sharp edge has the faster diffusion–reaction scale. The frozen-recovery speed above is a controlled subsystem result, not an exact speed for every coupled pulse. A definite pulse speed requires parameters and connecting boundary condition; none are specified by this brief-note question.

In gel sieving of synthetic polymers, gel electrophoresis combines drift under an electric field with size- and conformation-dependent obstruction by a porous matrix. In weak-field free solution, force balance gives $v=\mu E$ and, in a simple effective-charge/effective-friction model, $\mu =q_{\mathrm{eff}}/\zeta$. The mobility is motion relative to the medium, so fluid motion from electro-osmosis must be accounted for separately. This is the operational mobility convention, not a claim that the bare structural charge is always the effective driving charge.^[ref]

For a chain with free-draining polymer friction with approximately constant effective electric charge per segment, both charge and friction scale with $N$. The free-solution mobility can therefore become nearly independent of length, making size separation poor. A gel adds steric constraints. Coils smaller than the mesh can pass relatively easily; larger coils require deformation and pore-to-pore rearrangement. A schematic weak-sieving law is $\mu (c_g)\simeq \mu (0)e^{-K_r{c}_g}$, where matrix concentration $c_g$ and retardation coefficient $K_r$ are defined by a calibration; it is not a universal size law.

When a long flexible charged polymer is confined to a tube of available paths, reptation supplies another regime. An additive-friction model gives contour diffusion $D_c\propto N^{-1}$ and renewal time $\tau \sim N^2/D_c\propto N^3$. For ideal large-scale chain statistics, $R^2\propto N$, so centre-of-mass diffusion scales as $N^{-2}$. With total effective driving charge proportional to $N$, weak-field mobility can then scale as $N^{-1}$. This assumes screened/ideal large-scale conformations and linear response; strong-field alignment, trapping and nonideal matrices can alter the scaling.

Operationally, load narrow zones of independently characterized synthetic standards and sample into a controlled gel, apply a known field and detect migration positions after a measured time. Compare with a free-solution measurement, vary matrix concentration and field, and control temperature, ionic strength, pH and matrix charge. Larger gel pores or lower matrix concentration usually reduce sieving within a given calibrated material system. Excessive field causes Joule heating and can change conformation and mobility. adsorption or electro-osmotic drift can mimic a size effect; standards should match charge chemistry and topology as well as nominal mass.

Band centres move as $\mu Et$ while ordinary diffusion gives ${\sigma }_x^2(t)={\sigma }_x^2(0)+2Dt$. Resolution depends on separation of centres relative to band widths, not travel distance alone. A neutral tracer helps identify bulk flow, and multiple concentrations test adsorption or overloading. A mobility or band position does not uniquely determine size without a regime-specific calibration. TODO: biomolecular examples, sample preparation and biological applications are user deferred; the account here concerns synthetic charged polymers only.

velocity $U$, length $\ell$ and diffusivity $D$ have dimensions $L/T$, $L$, and $L^2/T$. Their dimensionless Péclet number isdiffusion dominates for $\mathrm{Pe}\ll 1$, advection for $\mathrm{Pe}\gg 1$. At fixed $U,D$, the crossover length is $D/U$; neglecting advection requires $\ell \ll D/U$.

TODO: user-deferred organism-specific swimming-speed and detection-scale estimates. No biological numerical substitution is made.

For the generic equilibrium transport derivation, let an external potential $V(x)$ give force $F=-V^{\prime }$ and drift $u=F/\zeta$. Boltzmann equilibrium has $c\propto e^{-V/(k_{B}T)}$, hence $c^{\prime }=Fc/(k_{B}T)$. Substituting in the zero-current condition givesThe Stokes–Einstein relation in the final equality is for an isolated no-slip sphere in a Newtonian liquid at low Reynolds number. Zero flux alone is insufficient: the equilibrium distribution and mechanical mobility law are also needed.

TODO: user-deferred organism-specific absorbing-surface calculation and its biological flux interpretation.

TODO: user-deferred receptor-number, receptor-density and receptor-competition calculations. These continuations have not been solved.

TODO: user-deferred biological kinetic mechanism and its biochemical pattern-forming applications. The remainder treats generic two-field reaction-diffusion physics only.

Let $\mathbf{u}=(u,v)$ obey ${\partial }_t\mathbf{u}=\mathbf{f}(\mathbf{u})+\mathrm{diag}(D_u,D_v){\nabla }^2\mathbf{u}$. At a uniform steady state, write its reaction Jacobian as $J=\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)$. Without diffusion stability requires $a+d<0$ and $\Delta =ad-bc>0$. A Fourier perturbation evolves under $J-q^2\mathrm{diag}(D_u,D_v)$, whose determinant isIts trace becomes more negative with increasing $q$. A stationary diffusion-driven instability therefore requireswith an allowed domain mode in the resulting negative-determinant band. Equal diffusivities cannot satisfy this while $a+d<0$. The physical mechanism is that transport can destabilize coupled fields even when the local kinetics relaxes a spatially uniform perturbation. Faster diffusion of an inhibitory response is one possible realization.

The selected wave numbers depend on growth rates and boundaries; nonlinear saturation determines the final stripes, spots or other states. Linear instability alone neither supplies a biochemical mechanism nor uniquely predicts a final finite-amplitude pattern.

Let the channel length be ${\ell }_c\gg R$ and assume uniform solution conductivity. Its Ohmic resistance scales as ${\ell }_c/(\sigma R^2)$, whereas access resistance of a reservoir mouth scales as $1/(\sigma R)$. The fraction of voltage in the entrance regions is therefore of order $R/{\ell }_c$. Most of the potential drop lies inside the long channel; the small entrance field region is neglected in the model, not asserted identically zero.

Write $x$ for inserted contour length and choose $qE>0$ to pull the chain inward. The charge already inside feels force $qEx$, givingConfinement removes polymer configurations. In the stiff-chain regime, the confinement deflection length obeys $R^2\sim {\lambda }^3/{\ell }_p$, hence $\lambda \sim (R^2{\ell }_p{)}^{1/3}$. A free-energy cost of order $k_{B}T$ per deflection gives a roughly constant resisting entropic force $f\sim k_{B}T/\lambda$. The insertion barrier of a charged polymer follows from $F(x)\simeq fx-qE{x}^2/2$, with an entry maximum at $x_{\ast }=f/(qE)$ when that point lies within the insertion interval. Its height is $f^2/(2qE)$. The field supplies work against this confinement entropy.

Put ${\zeta }_{\ell }=2\pi \eta /\ln (R/a)$, so $\gamma (x)={\zeta }_{\ell }x$. After capture, neglect the entry barrier, exterior-chain drag and thermal corrections to strong directed motion. ThenThe confined-polymer entry and recoil times have different length dependence; this entry graph is linear through the origin. Charge and drag both grow with inserted length, leaving a constant velocity. The same cancellation explains length-independent mobility in an ideal free-draining electrophoresis model. It does not include a capture waiting time or claim that every free-solution chain has identical hydrodynamic drag.

A measured superlinear entry time can arise when exterior-coil drag, tension propagation, wall adsorption or nonuniform friction adds length dependence absent from this model. For example $\gamma (x)\propto x^p$ with the same $qEx$ drive gives $\tau \propto L_0^p$ for $p>0$. Under the assumptions of uniform wall charge and polymer entry scaling, a homogeneous thin-double-layer electro-osmotic model adds a length-independent fluid velocity $u_{\mathrm{eo}}=-ϵ{\zeta }_{w}E/\eta$. It changes the drift coefficient, not its length exponent. This does not rule out charge heterogeneity, adsorption or near-cancellation followed by diffusion in a real channel; those are additional mechanisms.

With the field off, the confinement entropy drives outward motion. Assuming the constant force $f$ and the same local friction,The small final end region lies outside the long-confined-segment approximation but gives only a small correction. Thus ideal insertion scales as $L_0$, while recoil scales as $L_0^2$.

TODO: user-deferred DNA-specific numerical stopping-voltage estimate and its biological parameter-based inference. No biological quantities are evaluated.

For entropic-trap length separation, a reservoir/channel array can repeat length-dependent capture, confinement and relaxation delays, turning them into different mean residence times. The uniform inside-channel drift speed derived above is itself length independent, so that drift alone cannot separate lengths. Separation must retain length dependence of the barrier or of the capture/recoil dynamics and can be tuned with a forward-biased time-dependent field. For large bulk throughput, pressure-driven advection is preferable to diffusion alone; a pressure gradient can be combined with electric capture control. In the unobstructed continuum limit $Q=\pi R^4\Delta P/(8\eta {\ell }_c)$. Excessive driving erases trapping differences, so throughput and resolution compete; the question supplies no geometry or adsorption law selecting a unique optimal drive.

A generic chemical reaction–diffusion system couples local kinetics to spatial transport:A homogeneous stationary state satisfies $f=g=0$. Let its reaction Jacobian be $J=\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)$. Local well-mixed stability requires $a+d<0$ and $ad-bc>0$. A Fourier perturbation of wavevector magnitude $q$ evolves under $J-q^2\mathrm{diag}(D_u,D_v)$, whose determinant isSince its trace remains negative, a diffusion-driven stationary instability appears if the determinant becomes negative for some allowed wavevector. In a spatially large domain the conditions areTogether with reaction-only stability, these give a finite unstable wavelength range. Equal diffusion coefficients cannot generate this instability: they simply shift each reaction eigenvalue by $-D{q}^2$. On a finite domain, an admissible spatial eigenmode must actually lie in the unstable band. Nonlinearities then determine saturation and whether a predicted pattern is stable, so a linear unstable band alone does not calculate a final pattern.^[ref]

Other kinetics support travelling fronts or oscillatory waves rather than a stationary pattern. For instance, the scalar abstract equation ${\partial }_{t}u=D{\partial }_x^{2}u+ru(1-u/u_s)$ has a pulled-front speed $2\sqrt{Dr}$ for a sufficiently steep invading initial profile. Here it is used as an abstract kinetics-and-diffusion example, without a biological interpretation. The ratio of diffusion time $L^2/D$ to reaction time characterizes whether spatial mixing is rapid enough to justify a well-mixed model.

For a radial function $f(r)$ in three dimensions, the Laplacian is $f^{\prime \prime }+2f^{\prime }/r$, while the imposed radial drift contributes $-\alpha f^{\prime }/r$. Thus the radial generator of an SDE diffusion isThe hitting probability for the inner sphere is harmonic for this diffusion generator in the shell, with value one on the inner sphere and zero on the outer sphere. At $\alpha =1$, the boundary value problem is thereforeIntegrating $(r{p}^{\prime }{)}^{\prime }=0$ gives $p(r)=A\log r+B$. The two boundary conditions determineBoth arguments of the logarithms are dimensionless. The result decreases from one to zero across the shell and lies strictly between zero and one for $R_1<r_0<R_2$.

For a direct probability interpretation, stop at the first contact $T$ with either sphere. The stopped process $p(r(t\wedge T))$ is a bounded martingale by Itô's formula. The shell has positive inner radius, so its drift is bounded and its diffusion is nondegenerate; exit occurs almost surely. Optional stopping and dominated convergence therefore giveThe logarithm appears because the inward drift cancels one of the two geometric $f^{\prime }/r$ terms in the three-dimensional radial Laplacian, leaving the same radial diffusion generator as a planar Brownian motion with diffusion coefficient one.

biology is physically realized, but this does not make its explanations merely technological applications of physics. The strongest defensible position combines physicalism with qualified explanatory autonomy: biological mechanisms must respect physical constraints, while biological questions can require concepts, experiments and theories organized at other levels.

The contention associated with J. J. C. Smart should first be separated into different claims. Ontological reductionism concerns what living systems are made of; explanatory reductionism concerns whether explanations in terms of organisms, populations or regulatory organization can always be replaced without loss by descriptions of their physical constituents. A methodological preference for investigating smaller components is a further claim. Accepting the first does not logically establish the other two. The cited article supplies the historical setting for the contention; the distinctions here provide a way of assessing it.^[ref]

The reductionist case has substantial scientific force. DNA replication, enzyme catalysis and transport across a cell membrane involve material interactions, not an extra vital substance. thermodynamics constrains the energy available for metabolism; diffusion limits transport; chemical kinetics constrains the speeds of reactions. Describing a gene regulatory network without specifying how transcription factors bind and how their concentrations change can conceal an impossible mechanism. systems biology benefits from these constraints because a diagram of arrows alone does not guarantee that the proposed dynamics can occur. Nor should emergence be used as a name for an unexplained exemption from physical law.

Nevertheless, physical consistency leaves many explanatory questions open. Knowing the chemical composition of a population does not immediately explain why an allele frequency changes. In a simple haploid natural selection model with two alleles and reproductive outputs $w_A,w_a$, the next-generation frequency satisfiesThis relation identifies differential reproduction as the relevant comparison. It applies to many distinct molecular causes of the reproductive difference, provided the model's assumptions hold. An explanation of why an allele spreads can therefore be both physically implemented and usefully stated at the population level. genetic drift, migration or frequency dependence require changes to the model, rather than a conclusion that population-level explanation was never scientific.

Similarly, Mendelian segregation makes a conditional probabilistic claim about inheritance. For a heterozygote cross at one locus, equal segregation, random fertilization and complete dominance produce genotype probabilities $1/4,1/2,1/4$ and a phenotype ratio of $3:1$. A molecular description explains how chromosomes segregate and why exceptions occur; it does not remove the utility of the inheritance-level relation. Molecular detail and the higher-level relation answer connected questions at different resolutions. The experimental importance of tracking characters across controlled crosses is already clear in Gregor Mendel's original work.^[ref]

Organization also matters. The same components can form a stable negative feedback circuit or an amplifying positive feedback circuit, depending on how interactions are connected. Conversely, distinct biochemical implementations can produce the same adaptation behavior. This is multiple realizability: it supports explanations in terms of a relationship or control architecture that remains useful across several implementations. It does not establish that every implementation has identical dynamics, or that a sufficiently complete physical description could never recover the behavior. The claim is about what makes an explanation informative and generalizable, not a proof of nonphysical causation.

A further distinction concerns historical questions and biological function. To ask how a receptor transduces a signal is to ask for a present mechanistic explanation. To ask why a lineage retained that receptor may require evidence about ancestry, past environments and natural selection. Present physical laws alone do not supply those historical boundary conditions. A functional claim must still be tested: the existence of a useful feature does not prove that it evolved because of that use, and natural selection has no foresight. Explanations of evolutionary history complement rather than contradict physical accounts of current operation.

The objection that biology lacks exceptionless laws is also too quick. Many useful biological generalizations are conditional on explicitly stated assumptions, but physical models also idealize away interactions, specify boundary conditions and use statistical descriptions. stochasticity in a biological model need not mean absence of regularity: a probability distribution can be a precise prediction even when one individual outcome is uncertain. The scientific issue is whether assumptions, predictions and domains of applicability are clear enough to be tested.

There are limits to this defense of autonomy. A higher-level explanation should not be protected from evidence about its mechanism, and an attractive organizing principle should not substitute for quantitative testing. Conversely, merely listing every molecular event can fail to reveal which changes would alter the outcome. A productive investigation moves between levels: identify a population or cellular regularity, propose a mechanism, perturb its components, and test whether the higher-level prediction survives. biology's autonomy is therefore explanatory and methodological, not independence from the material world; calling it only technology understates its role in discovering and testing such regularities.

A simple model is especially valuable when it overturns a clearly stated expectation and explains why, but surprise is a means to understanding rather than a requirement that every useful model be counterintuitive. The expectation must be explicit enough to test, and the model should identify which assumptions make it fail. A surprising simulation without a mechanism or a plausible connection to observations is not sufficient.

One expectation is that a deterministic system with few variables must settle to a simple, predictable state. The logistic map,provides a compact counterexample. Its positive fixed point is $x_{\ast }=1-1/r$, and the derivative there is $f^{\prime }(x_{\ast })=2-r$. Thus the fixed point is locally stable for $1<r<3$; changing a single parameter can destroy that stability. At $r=3.2$, repeated iteration approaches a two-cycle near $0.513045$ and $0.799455$. For typical initial conditions at $r=4$, the Lyapunov exponent is $\log 2>0$, so nearby initial states separate exponentially on average until the bounded state space limits their separation. Robert M. May's analysis made the consequences of such simple difference equations particularly influential in population modeling.^[ref]

The broken expectation is not that determinism is false: the rule remains deterministic. Rather, an irregular trajectory need not be caused by a complicated mechanism or external stochasticity, and exact rules need not permit long-range prediction from imperfectly measured initial conditions. This does not show that every irregular population time series is chaotic. Real populations may experience demographic noise, environmental variation or overlapping generations, and the map's parameters and state must have a defensible biological interpretation. Evidence should distinguish those alternatives. Even within the mathematical example, special initial states can be periodic; the claim about $r=4$ concerns typical trajectories, not every point of the interval.

A second expectation is that diffusion can only erase spatial variation. In a coupled reaction-diffusion equation,diffusion smooths each species separately, yet can destabilize their coupled uniform equilibrium. Linearizing about that equilibrium gives the mode matrix $J-k^{2}D$. Local stability of $J$ therefore does not by itself imply stability of every spatial mode. This is the organizing insight of the Turing instability.^[ref]

For a concrete dimensionless linearization, takeThe well-mixed system has trace $-1$ and determinant $1$, so its perturbations decay. Writing $q=k^2$, the spatial mode has trace $-1-11q$ and determinantThis determinant is negative for $(4-\sqrt{6})/10<q<(4+\sqrt{6})/10$, creating one positive eigenvalue and a growing nonuniform mode. Equal diffusivities would instead shift both eigenvalues left and could not create this two-species instability. On a finite domain, at least one of the modes allowed by the boundary conditions must lie in the unstable band. The algebra exposes exactly why the intuition about smoothing fails: differential transport changes the balance of an interacting local system.

A Turing instability establishes a possible pattern-generating mechanism, not proof that a particular animal's markings use it. One must identify candidate interacting species, measure transport and response properties, and test perturbations that alter the predicted wavelength or stability. The linear example specifies onset only; nonlinear dynamics determine whether growth saturates into a stable pattern. Its value is to convert a vague statement about morphogenesis into discriminating questions, while leaving those experimental obligations visible.

A third expectation is that an infection with a basic reproduction number above one must produce a major outbreak. Consider the early phase as a continuous-time branching process, with independent per-infective transmission rate $\beta$ and recovery rate $\gamma$, while susceptible depletion is negligible. If $q$ is the probability of eventual extinction starting from one infective, conditioning on the first event givesTaking the smallest solution in $[0,1]$ yields $q=1$ when $\beta \le \gamma$, and $q=\gamma /\beta =1/R_0$ when $R_0>1$. Thus $R_0=4$ still gives an early-lineage extinction probability of $1/4$. A deterministic SIR model describes mean-field growth but omits this individual realization uncertainty. In a finite closed stochastic SIR model, every epidemic eventually ends; the branching calculation concerns early extinction before substantial susceptible depletion, not that eventual end or an arbitrary peak-size cutoff.

These examples show different benefits of simplicity: a stability calculation can separate predictable from sensitive dynamics, a mode calculation can uncover a new mechanism, and a probability recursion can expose what an average trajectory hides. The same discipline also explains why the word “should” needs qualification. A simple exponential-decay model that confirms a familiar expectation can still calibrate a rate, serve as a control, or expose a faulty measurement. A useful mathematical model must earn its simplifications through the question it answers. Breaking an expectation is a powerful route to discovery when the contradiction is precise, robust to reasonable assumptions and experimentally consequential; novelty alone is not a validation criterion.