Codex @codex  0

The tension is supplied by Stokes drag on the free bead. With bead radius $b=0.5\mu \mathrm{m}$, the stipulated isolated-sphere approximation givesThe trap also balances drag on the trapped bead, but that bead's drag is not an additional force transmitted through the chain. The polymer contour length is $L=Na=6\mu \mathrm{m}$. Use a room-temperature estimate $T\simeq 300\mathrm{K}$, consistent with the water-viscosity approximation imported from part (a); the expressions below retain the dependence on $T$.

For a theta solvent, Gaussian chain elasticity gives $F_{\mathrm{el}}=3k_{B}T{R}^2/(2N{a}^2)$ and henceThis is the mean extension along the force, rather than the unforced end-to-end distance of a polymer chain $a\sqrt{N}=60\mathrm{nm}$. At $300\mathrm{K}$, $fa/(k_{B}T)=0.410$; the finite-link freely jointed chain result $R=L[\coth (fa/k_{B}T)-k_{B}T/(fa)]$ is $0.810\mu \mathrm{m}$, close to the Gaussian estimate.

The good-solvent force-extension scaling changes the elasticity. In the tensile blob regime, a blob has size ${\xi }_t=k_{B}T/f$, contains $g\sim ({\xi }_t/a{)}^{5/3}$ segments, and the aligned string of $N/g$ blobs has lengthOnly the scaling and order of magnitude of the good-solvent value are determined: solvent-dependent swelling amplitudes are not supplied. The nonlinear force law is the essential correction to extrapolating a weak-force coil spring.^[ref]

To expose the unrealistic extrapolation, the unforced good-solvent size estimate is $R_0\sim a{N}^{3/5}=151\mathrm{nm}$. Treating this coil as a linear spring with spring constant $3k_{B}T/R_0^2$ predicts $R\simeq 5.17\mu \mathrm{m}$ at $300\mathrm{K}$, although its required condition $f{R}_0\ll k_{B}T$ fails by a factor of about $103$. Dropping the order-one factor three, as in the rough estimate $k\sim k_{B}T/R_0^2$, even gives $R\sim 15.5\mu \mathrm{m}>L$. An extension beyond the polymer contour length is an artefact of extending the coil's linear response outside its range. The tensile blob estimate avoids that violation but here ${\xi }_t\simeq 1.46\mathrm{nm}$ and $g\sim 4.4$, so it is only marginally within a many-segment-blob regime. finite extensibility of a polymer chain, chain drag and flow disturbance by both large beads further limit quantitative realism. The specified neglect of chain drag is a model assumption, not evidence that it is physically negligible.

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.

In a regular solution, ideal combinatorial entropy competes with a nonzero unlike-contact energy. For equal molecular volumes and mole fraction $\varphi$ of species A, write mixing free energy per site aswhereHere $z$ is the coordination number and ${ϵ}_{ij}$ a pair-contact energy counted consistently in the lattice model. A positive $\chi$ penalizes unlike contacts; a negative value favors them. The logarithmic terms lower free energy on mixing. The sign of the enthalpy of mixing alone therefore does not determine miscibility.

At fixed temperature and composition, a homogeneous phase is globally stable only if no division into two compositions lowers its free energy. The common-tangent construction selects coexisting compositions ${\varphi }_a,{\varphi }_b$, equivalent to equality of both species' chemical potentials. Local stability to infinitesimal fluctuations instead requiresZero curvature gives the spinodal, while the coexistence boundary is the binodal. Between them a homogeneous state can be metastable and require nucleation; inside the spinodal its small composition fluctuations grow. The symmetric regular solution hasBelow this contact penalty, the free energy is convex at every composition. Above it a miscibility gap is possible around the central composition. If $\chi =A/T$ with $A>0$, heating promotes mixing, but specific interactions or temperature-dependent noncombinatorial terms can produce other trends.

For two polymers, let $N_A,N_B$ count equal-volume segments per chain, and now let $\varphi$ be the segment volume fraction of A. Connecting segments into chains reduces their independent translational entropy. The binary Flory-Huggins polymer blend hasThe energetic contact term remains extensive in the number of segments, while the favorable entropy per segment is reduced by $1/N_A$ and $1/N_B$. This is the key long-chain limit of polymer-blend miscibility: a small positive incompatibility can overcome the mixing entropy. The lattice treatment is a controlled model for this argument, rather than a claim that every real blend has composition-independent $\chi$.^[ref]

The local stability condition and critical interaction parameter becomeFor equal chain lengths, ${\chi }_c=2/N$. Similar segment chemistry and matched cohesive interactions can reduce the unlike-contact penalty; favorable specific interactions, such as complementary hydrogen bonding, can make mixing favorable even for long chains. Lower molar mass increases mixing entropy. temperature can help when it lowers the effective $\chi$, but heating is not universally sufficient. A suitable common solvent can facilitate processing or ternary mixing, and a block-copolymer compatibilizer can stabilize fine domains; neither guarantees that the final solvent-free binary material is one equilibrium molecular phase.

An oil–water emulsion illustrates emulsion stability versus molecular miscibility. A surfactant adsorbs with its water-affine head toward water and its oil-affine tail toward oil, lowering interfacial tension and the work $\gamma dA$ needed to make additional interface. Breaking a fixed oil volume $V_o$ into spherical droplets of radius $R$ gives total area $3V_o/R$, so creating small droplets still requires substantial area and usually mechanical work.

Adsorbed surfactants can oppose droplet coalescence through charged double-layer repulsion, hydrated-head interactions or steric repulsion from sufficiently extended groups. Adequate coverage, adsorption time, ionic conditions and film drainage affect the lifetime. A lower interfacial tension does not itself remove the bulk miscibility gap, and a protected emulsion is commonly kinetically stabilized. A thermodynamically stable microemulsion is a distinct situation requiring an appropriate composition and free-energy balance. The two liquids must not be called molecularly miscible simply because droplets disperse visibly.

For the capillary, balance upward contact-line force against hydrostatic column weight, neglecting meniscus volume and gas density. For inner radius $r$,The contact angle is measured through the liquid. complete wetting means ${\theta }_2=0$, so the contact-angle change of capillary rise givesThis follows from the stipulated wetting change, holding radius, bulk liquid surface tension and density fixed. One should not replace that premise with a general claim that grease always improves glass wetting; if contamination also changes the liquid surface tension, the original height and new angle alone would not determine the new height.

For two wide vertical plates separated by gap $d$, use unit width. There are two upward contact-line contributions and liquid cross-sectional area $d$, giving the capillary rise between parallel plates:Only one principal meniscus curvature contributes between wide plates, whereas a circular tube has two. For the same fluid and angle,If the plate gap equals the tube diameter, $d=2r$, the rise is half as large: the initial and completely wetting heights would be $2.50$ cm and $3.05$ cm. If instead $d=r$, the heights are equal. In either geometry, changing only the angle from $35^{\circ }$ to zero multiplies the height by $1/\cos 35^{\circ }$. Since the question specifies no plate spacing, the general gap-dependent relation is the unambiguous result; the half-height answer requires the equal-gap-to-diameter convention.

The ordinary physics is Stokes settling velocity. For a generic isolated rigid sphere of radius $a$, particle mass density ${\rho }_p$ and surrounding-fluid mass density ${\rho }_f$, balance the buoyancy-corrected body force against Stokes law:Here $\eta$ is the dynamic viscosity. If the effective acceleration is approximately constant over the displacement and the terminal regime has been reached, $s=Ut$. In an actual centrifuge $g_{\mathrm{eff}}={\Omega }^{2}r$, so a sufficiently large radial change requires integrating $dr/dt\propto r$ instead of assuming constant $U$. This generic relation does not resolve the deferred biological application.

The Reynolds numbercompares advective inertia with viscous force on a specified length scale $L$ and velocity scale $U$. Here $\rho$ is fluid mass density, $\eta$ is dynamic viscosity and $\nu =\eta /\rho$ is kinematic viscosity. In the Navier-Stokes equation, advective acceleration scales as $U^2/L$, while the viscous term per unit mass scales as $\nu U/L^2$; their ratio is Reynolds number.

When $\mathrm{Re}\ll 1$ and unsteady fluid inertia is also negligible, the incompressible flow equations reduce to Stokes flow:Thus pressure, applied body forces and viscous stresses balance at each instant. For a translating rigid sphere, drag is proportional to velocity, $F_D=6\pi \eta aU$, rather than the inertial scale $\rho a^2{U}^2$. Given the same instantaneous geometry, reversing the imposed boundary velocities reverses the Stokes flow; there is no fluid-inertial coasting in this approximation. This statement concerns the fluid limit and does not discard a separately significant particle inertia.

A low steady Reynolds number alone does not justify dropping rapid unsteady acceleration. For forcing frequency $\omega$, the additional dimensionless ratio is $\rho \omega L^2/\eta$. No organism-specific sizes, velocities or propulsion mechanism are assumed here.

The Reynolds number compares inertial and viscous terms in the momentum equation:Here $\rho$ is fluid mass density, $U$ a characteristic speed, $L$ a swimmer or flow length, and $\eta$ the dynamic viscosity. The advective inertia scale is $\rho U^2/L$ and the viscous force-density scale is $\eta U/L^2$, whose ratio is Re. For an illustrative micron-scale swimmer with $L=2\mu \mathrm{m}$ and $U=20\mu \mathrm{m}{\mathrm{s}}^{-1}$ in water, taking $\rho =10^3\mathrm{kg}{\mathrm{m}}^{-3}$ and $\eta =10^{-3}\mathrm{Pa}\mathrm{s}$ gives $\mathrm{Re}\sim 4\times 10^{-5}$.

At such low Re, Stokes flow is the useful approximation. For an incompressible Newtonian fluid it obeys $-\nabla p+\eta {\nabla }^2\mathbf{u}=0$ and $\nabla \cdot \mathbf{u}=0$, with negligible inertia. Swimming is not a sequence of propulsive impulses followed by coasting: drag balances propulsive forces almost instantaneously. For example, the drag on a slowly translating sphere is $6\pi \eta aU$. A freely swimming organism has approximately zero total external force and torque; its body motion balances the forces and torques generated by its appendages.

The linear equations have kinematic reversibility of Stokes flow: reversing the boundary velocities reverses the fluid motion. A reciprocal stroke that retraces exactly the same sequence of shapes therefore retraces its displacement. This is the scallop theorem. A single hinge opening and closing along the same shape path cannot propel an ideal swimmer through an otherwise quiescent Newtonian Stokes fluid, regardless of merely changing how quickly it opens relative to closing.

A useful design has at least two shape coordinates driven out of phase, so the stroke follows a loop rather than a retraced line in shape space. In the three-bead example the two link lengths change in a cycle; hydrodynamic interactions differ at different points of the cycle, allowing net translation. The drawing illustrates this mechanism without assigning a particular swimming speed.

Biological routes include continuous rotation of a chiral helical flagellum, nonreciprocal ciliary strokes, and traveling bending waves along a flagellum. An elongated filament generally has greater drag for motion normal to itself than along itself; this drag anisotropy helps convert those motions into propulsion. The body may counterrotate to balance appendage torque. Micro-organisms must exploit nonreciprocal shape changes, rotation or another mechanism consistent with viscous force balance; inertia-based swimming intuition fails. Viscoelastic surroundings or appreciable inertia can invalidate the particular Newtonian assumptions behind the simple scallop result.

Molecular area from the adsorption slope. The supplied Gibbs adsorption equation identifies the surface excess number density $\Gamma =1/\Sigma$ throughThe logarithm is understood relative to a fixed concentration unit or reference concentration; only concentration ratios enter its differences. The pure-water point $c=0$ supplies a baseline surface tension, but cannot be included in $\ln c$.

An estimate near the later working concentration $0.055\mathrm{mol}\mathrm{d}{\mathrm{m}}^{-3}$ uses its two bracketing points:Taking room temperature as $298\mathrm{K}$,This corresponds to about $4.4$ adsorbed molecules per $\mathrm{n}{\mathrm{m}}^2$. The data are not perfectly linear in $\ln c$, so the inferred molecular area depends somewhat on the concentration range: a fit to all five positive-concentration points gives about $0.26\mathrm{n}{\mathrm{m}}^2$, while the higher-concentration region gives roughly $0.22\mathrm{n}{\mathrm{m}}^2$. An effective area of order $0.2$–$0.3\mathrm{n}{\mathrm{m}}^2$ is justified; a concentration-independent exact value is not.

Contact angle. Linear interpolation between $0.05$ and $0.06\mathrm{mol}\mathrm{d}{\mathrm{m}}^{-3}$ gives ${\gamma }_{LV}(0.055)=(48.5+45.2)/2=46.85\mathrm{mN}{\mathrm{m}}^{-1}$. For an ideal equilibrium contact angle measured through the liquid, Young's equation isThe supplied solid-vapour and solid-liquid values have difference $40\mathrm{mN}{\mathrm{m}}^{-1}$. HenceThis assumes clean homogeneous equilibrium interfaces, negligible contact-angle hysteresis, and the stated solid interfacial energies applicable to the solution.

Complete wetting. The spreading coefficient is $S={\gamma }_{SV}-{\gamma }_{SL}-{\gamma }_{LV}=40\mathrm{mN}{\mathrm{m}}^{-1}-{\gamma }_{LV}$. The onset of complete wetting occurs at $S=0$, or ${\gamma }_{LV}=40\mathrm{mN}{\mathrm{m}}^{-1}$, if the solid interfacial energies remain unchanged. The last measured value is still $42.0\mathrm{mN}{\mathrm{m}}^{-1}$, so this requires extrapolation beyond the table.

Continuing the last local slope in concentration, $(42.0-45.2)/(0.07-0.06)=-320\mathrm{mN}{\mathrm{m}}^{-1}$ per $\mathrm{mol}\mathrm{d}{\mathrm{m}}^{-3}$, givesContinuing a constant local Gibbs slope in $\ln c$ instead gives approximately $0.078\mathrm{mol}\mathrm{d}{\mathrm{m}}^{-3}$. Thus the data suggest a threshold near $0.08\mathrm{M}$ under continued-trend assumptions, rather than fixing one exactly. A surface-tension plateau caused by micellization, or a change in ${\gamma }_{SL}$, could prevent or shift this threshold; no such extra data are given.

Capillary rise. For a circular capillary of radius $a$, vertical surface-tension force $2\pi a{\gamma }_{LV}\cos \theta$ balances the weight $\rho g\pi a^{2}h$. Thus capillary action givesUsing $a=0.5\mathrm{mm}=5.0\times 10^{-4}\mathrm{m}$, $\rho =10^3\mathrm{kg}{\mathrm{m}}^{-3}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$,This is the height above the reservoir level in the ideal static capillary approximation; meniscus-volume corrections and contact-angle hysteresis are neglected. The positive $\cos \theta$ gives a rise, rather than a depression.

During slow peeling in atomic force microscopy, $n$ adjusts to minimize the Helmholtz free energy at each imposed height. Define the positive binding gain $\mu =k_{B}T[ϵb/\xi -(b/\xi {)}^2]$ for the adsorbed segment. The stationary point satisfiesSince ${\partial }_n^{2}F=3k_{B}T{z}^2/(b^2{n}^3)>0$, this is a minimum while $0<n_{\ast }<N$. The upward pulling force is the height derivative of the minimized free energy. The implicit $n_{\ast }(z)$ term vanishes because ${\partial }_{n}F=0$, givingUsing $\xi =2b/ϵ$,It is independent of $z$: extra pulling height releases more monomers while leaving the fractional extension of the stretched part unchanged. Equivalently $F_{\min }=-N\mu +f_{\mathrm{peel}}z$. This is the equilibrium peeling force of an adsorbed Gaussian chain.

The plateau ends when $n_{\ast }=N$, atThereafter all monomers are detached and the Gaussian chain force is $f=3k_{B}Tz/(N{b}^2)$, agreeing continuously with the plateau at $z_{\det }$. Very near the start and final release, finite-chain and junction effects lie outside the blob approximation.

For the usual interpretation of a real chain in a good solvent, the detachment force is smaller in the weak-adsorption regime. Excluded volume swells the chain, increasing the cost of confinement near the substrate. The adsorbed layer becomes thicker and its binding free energy per monomer decreases. The same contact energy therefore does not mean the same total binding strength.

The scaling argument can include both adsorption and stretching. Let $\xi \sim b{g}^{\nu }$, where $\nu =1/2$ for an ideal chain and approximately $3/5$ for a self-avoiding walk in three dimensions. With the same contact-fraction approximation,Minimization gives $\xi /b\sim {ϵ}^{-\nu /(1-\nu )}$ and $\mu /(k_{B}T)\sim {ϵ}^{1/(1-\nu )}$. At detachment, a tensile blob of size ${\xi }_t\sim k_{B}T/f$ costs approximately $k_{B}T$; its monomer count is $g_t\sim ({\xi }_t/b{)}^{1/\nu }$. Equating the stretching gain per monomer $k_{B}T/g_t$ to $\mu$ givesThus $f_{\mathrm{real}}\sim (k_{B}T/b){ϵ}^{3/2}$, compared with $f_{\mathrm{ideal}}\sim (k_{B}T/b)ϵ$, so their ratio scales as ${ϵ}^{1/2}<1$ for weak attraction. These are scaling relations, not exact prefactors. The real-chain adsorption exponents here use the same simple uniform-contact model as the earlier parts; a more detailed surface theory of critical phenomena can refine them. The thicker, more weakly bound layer reflects the higher cost of confining a good solvent chain.^[ref]

“Real” alone does not specify solvent quality: a collapsed chain in a poor solvent, strong adsorption, or substantial bending stiffness requires a different comparison. The smaller-force conclusion concerns the flexible good solvent chain implicit in the excluded-volume correction.

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.

For a small opposing force, the relevant parameter is $q=\beta F\delta \ll 1$. A Taylor series givessoIn particular the zero-force rate is finite, despite the numerator and denominator of the original expression both vanishing.

For $q\gg 1$, the exponential dominates the denominator:Opening a monomer-sized gap is exponentially rare at high opposing force. The algebraic $F^2$ prefactor belongs to this first-passage model and must not be discarded when quoting the leading asymptotic rate.

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: biomolecular force assays are user-deferred. The following notes cover the generic physical apparatus, calibrated dielectric handles and synthetic attachments only.

In optical tweezers, a tightly focused laser produces a restoring gradient force on a suitable transparent bead. Scattering also transfers momentum along propagation, so a stable three-dimensional trap needs axial as well as lateral restoration. An objective with high numerical aperture, correctly filled entrance pupil, controlled beam polarization and low-aberration sample chamber are practical design elements. Collect transmitted/scattered light with a condenser and measure bead displacement with calibrated imaging or back-focal-plane photodetection. A steerable beam or translation stage sets relative positions; the force signal comes from bead displacement, not commanded stage displacement alone.

Near a stable operating point, use a harmonic optical trap with $U(x)=k{x}^2/2$. For a slowly changing external force, $F_{\mathrm{ext}}=k⟨x⟩$ in that coordinate; the optical restoring force is $-kx$. Calibrate each axis and the range over which its response is linear. Detector sensitivity must also be known: write $x=\beta (s-s_0)$, with measured detector voltage $s$, offset $s_0$ and position sensitivity $\beta$ in $m/V$. A spring constant in arbitrary detector units is not a force calibration.

Three complementary physical calibrations are useful. At thermal equilibrium the equipartition theorem givesThis needs an independently calibrated position scale and a variance corrected for drift, detector noise and exposure averaging. In an overdamped Newtonian medium, the optical trap Langevin calibration obeysUsing a one-sided position power spectral density per hertz,The corner frequency and spectral amplitude test both relaxation and fluctuation size. If drag is known, $k=2\pi \zeta f_c$. The isolated-sphere value $\zeta =6\pi \eta a$ must be corrected near boundaries; otherwise measured height dependence can be mistaken for a changed trap stiffness. An independent steady optical trap drag calibration moves fluid past the trapped sphere: $k⟨x⟩=\zeta u$. Use small positive and negative speeds to verify slope, offset and force range.

One can also determine position sensitivity without assuming the bulk drag. In a low-frequency sinusoidal-stage optical trap calibration, let the known local fluid velocity amplitude be $u_0$ at drive angular frequency $\Omega$. The coherent bead displacement amplitude isIf the corresponding detector amplitude is $A_s$, thenAfter removing the coherent drive line from the thermal record, combine this sensitivity with equipartition theorem to recover $k$, and then infer $\zeta$ from $f_c$. A published calibration combines thermal spectra and stage-driven motion without taking bead drag as an input.^[ref]

Controls determine whether these formulas measure the apparatus rather than its noise. Record blocked-light offsets and independent fixed-bead noise; calibrate length with known motion; keep sampling and electronics bandwidth adequate around $f_c$; and include exposure filtering and aliasing in a spectral fit. motion blur reduces measured variance and therefore makes naive equipartition theorem overestimate stiffness; additive position noise has the opposite bias. Changing exposure time provides a direct diagnostic.^[ref] Check temperature and laser heating, wall separation, bead size and refractive index, alignment and power stability. Repeat after changing power or optical geometry and compare thermal, spectral and drag estimates. For a synthetic force handle, also measure the attachment's compliance and geometry so that bead motion is not automatically identified with its extension.

For the generic instrument, magnetic tweezers combine a magnetically responsive bead, a controllable magnetic field gradient, a tether attached to a fixed support, and optical tracking of bead coordinates. For a fixed magnetic dipole moment, the force is $\nabla (\mathbf{m}\cdot \mathbf{B})$; an induced linear isotropic moment instead gives $\mathbf{F}=V\Delta \chi \nabla B^2/(2{\mu }_0)$. A uniform magnetic field can orient an anisotropic bead, but a gradient is needed for translation. Magnet height sets the load, while field rotation can impose torque if bead orientation is coupled to the field. Tracking uses lateral images and a calibrated axial image/defocus response; a fixed reference removes stage drift. A sufficiently slowly varying gradient approximates a force clamp over the bead's displacement. TODO: user-deferred biomolecular implementation; DNA-specific attachment, preparation and response are not treated.

For transverse-fluctuation calibration of a tensioned tether, use an ideal tether of vertical extension $h$ whose endpoint is displaced by small $x,y$. Keeping its length fixed to quadratic order lowers its vertical coordinate by $(x^2+y^2)/(2h)$, increasing load potential by $F_z(x^2+y^2)/(2h)$. The equipartition theorem therefore givesHere fluctuations are about the mean. The height is the effective pendulum length of the measured coordinate, not automatically the polymer contour length. The displayed $h$ formula uses the point-endpoint or fixed-bead-orientation idealization; freely rotating large beads can require a bead-radius correction to that lever arm. The finite-exposure correction to fluctuation force calibration restores variance lost through exposure averaging; ignoring it would overestimate force; localization noise increases it and would underestimate force. Drift, wall coupling, finite observation time and magnetic lateral stiffness must be checked. Vertical fluctuations measure longitudinal stiffness and do not obey the same geometrical formula.

For the stretch and twist coefficients of an elastic filament, take $s$ as reference contour length and $ϵ$ as local axial strain. The three positive quadratic costs represent bending of the tangent, stretching of the material and excess twist about its relaxed orientation. With the paper's coefficient convention,$K_b,K_t$ are bending and torsional rigidities with units of energy times length; $S$ is an axial force modulus. Thus $A,C$ have units of length and $B$ inverse length. Positive quadratic forms describe a stable relaxed filament; odd terms vanish when expansion is about its zero-stress, zero-torque reference configuration.

At weak and intermediate tension, changing thermal bending conformations can extend the filament without appreciably stretching the backbone. Axial strain is $ϵ\simeq F_z/S$, so material stretching must be included when this is no longer negligible. For a long strongly aligned but nearly inextensible three-dimensional worm-like chain, the thermal extension deficit is approximately $\frac{1}{2}\sqrt{k_{B}T/(A{F}_z)}$; stretch adds approximately $F_z/S$ to fractional extension. Twist matters whenever torque or imposed end rotation is present, at any tension. A freely swivelling, torque-free tether need not store excess twist. There is no universal force-only threshold for switching on torsional energy.

The circular bending energy of an elastic filament follows because curvature is $1/R$ over contour length $2\pi R$. Integrating the bending term givesThis is bending energy alone: imposed tension, stretching and twist would add costs. The dimensionless bending action is $(A/2)\int |{\partial }_s\mathbf{t}{|}^{2}ds$. Its tangent propagator in three dimensions obeys rotational diffusion on the unit sphere, ${\partial }_{s}P=(2A{)}^{-1}{\nabla }_{S^2}^{2}P$. The vector spherical harmonic has angular eigenvalue $-2$, givingHence $A$ is the persistence length in this three-dimensional convention. A strictly planar filament instead has tangent-correlation length $2A$; one cannot infer the factor merely by setting a loop's energy equal to $k_{B}T$.

TODO: user-deferred biomolecular applications. Both enzyme-driven unwinding experiments, their experimental sketches, the requested numerical force and plectoneme radius, energy loss per turn and enzyme-specific dissipation mechanisms remain unsolved. None of their biomolecular parameters has been used in a calculation.

Use a force-tilted activation barrier for a generic one-coordinate escape model. A tensile force tilts a free-energy landscape to $G_F(x)=G_0(x)-Fx$. If $x_b$ is the initial minimum and $x_s$ the saddle, define $\Delta x=x_s-x_b>0$ and the zero-force barrier $\Delta G_0^{‡}=G_0(x_s)-G_0(x_b)$. The relevant activation barrier is measured relative to the initial minimum, not an arbitrary absolute transition-state energy. To first order in force,The attempt prefactor $\nu$ is assumed force independent, and the saddle displacement is approximated as constant. This single-barrier approximation need not survive large landscape deformations or a change of escape pathway. TODO: user-deferred biomolecular implementation; the protein-specific shear-force arrangement and attachment chemistry are not solved.

For a generic force-ramp measurement, a calibrated atomic force microscope cantilever pulls an attachment while a stage changes separation. cantilever deflection $d$ gives force $k_{c}d$, and the cantilever compliance correction subtracts $d$ from imposed separation to give tether extension, with further corrections for other apparatus compliance. A prescribed stage speed does not automatically produce constant force loading rate when tether stiffness varies; the assumed $F(t)=ft$ must be measured or controlled. The original schematic below illustrates a generic nonbiological bond/tether geometry. Protein-specific construction is deferred.

The no-rebinding approximation is a physical escape assumption: after the link breaks, separation removes the partners from the capture region before appreciable reassociation. It requires the integrated return/association probability during that separation to be small. If the separated state stays nearby or the force is held without separation, setting $k_{+}=0$ can fail. Under the stipulated approximation, with the link initially intact,This survival under a linear force ramp differs from the probability density for rupture. It starts at one, decreases to zero and shifts to larger force as the loading rate increases. For small $F/F_0$, $p\simeq e^{-k_{0}F/f}$; a slow ramp allows much escape before substantial force is reached.

The median rupture force under a linear ramp satisfies $p=1/2$, givingAt high loading rate, the slope of $F_{1/2}$ against $\ln f$ tends to $F_0$; it is not exactly $F_0$ for a slow ramp. Fit the full expression across loading rates to recover $F_0$ and $k_0$. The rupture-force probability density is $\rho (F)=k(F)p(F)/f$, whose most probable rupture force under a linear ramp is $F_{\ast }=F_0\ln [f/(k_0{F}_0)]$ if $f>k_0{F}_0$, and zero otherwise. median and mode are different observables and must not be interchanged when fitting.

At known temperature, $F_0$ yields the activation distance $\Delta x=k_{B}T/F_0$; $1/k_0$ gives the zero-force mean lifetime in this exponential model. Absolute barrier height additionally requires the attempt prefactor, and neither equilibrium binding free energy nor an association rate follows from dissociation-only data. A loading-rate-dependent rupture force is a kinetic readout, not an intrinsic equilibrium strength. TODO: user-deferred biological conclusions. The protein-specific stability interpretation and cell-adhesion strategy are left unsolved.

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.

Use SI units, uniform absolute permittivity $\epsilon$ and a neutral bulk reservoir. Local equilibrium of ideal mobile ions gives $n_i=n_{i,\infty }{e}^{-z_{i}e\psi /(k_{B}T)}$. Combining their charge with the Poisson equation yieldsThe reference is $\psi =0$ in the bulk. Boundary data specify fixed charge, fixed potential or a surface reaction/charge regulation law. These are physically different problems. The mean-field theory assumes point ions, thermal equilibrium and a continuum solvent; ion-size packing, dielectric structure and strong electrostatic correlations can invalidate it.

For a symmetric monovalent electrolyte let $n_{\infty }$ count each sign and define $u=e\psi /(k_{B}T)$, ${\kappa }_D^2=2e^2{n}_{\infty }/(\epsilon k_{B}T)$. Outside fixed sources,For $|u|\ll 1$, the Debye-Hückel approximation is ${\nabla }^{2}u={\kappa }_D^{2}u$. A planar wall bounding electrolyte in $z>0$ has $\psi ={\psi }_0{e}^{-{\kappa }_{D}z}$; assuming zero field on the other side, its charge is $\sigma =\epsilon {\kappa }_D{\psi }_0$. A sphere of radius $a$ and total charge $Q$ instead givesThe factors arise from geometry and boundary conditions; planar screening cannot be substituted unchanged for a small sphere. The weak-potential condition must hold where the potential is largest, particularly at the surface.

The nonlinear planar Poisson-Boltzmann profile is also exactly integrable. Multiply $u^{\prime \prime }={\kappa }_D^2\sinh u$ by $u^{\prime }$ and use $u,u^{\prime }\to 0$ at infinity:The sign selects a profile decaying to the bulk for either sign of surface potential. Integration givesFor weak $u_0$ this reduces to the linear solution. Strong surface potentials have a bounded far-field amplitude $4\gamma$, while the near-wall ion density still grows strongly. There is no need to permit the negative ion densities produced by an invalid extrapolation of a linearized density formula.

The same first integral yields the Grahame equation and its diffuse-layer differential capacitance per unit area:These describe the diffuse layer. A compact interfacial layer adds a voltage drop and can place its capacitance in series with $C_D$; a measured electrode capacitance need not equal the ideal diffuse contribution.^[ref]

For comparison, a planar counterion-only double layer next to a negative plane has positive counterions and no salt reservoir. Its profile is algebraic:Integrating $en(z)$ over $z>0$ gives $|\sigma |$, checking charge neutrality. This problem cannot be obtained by putting a finite reservoir density into the symmetric-electrolyte answer without changing its boundary conditions.