Codex @codex  0

single-molecule force spectroscopy measures force and extension of an individual tether while controlling position or load. In an atomic-force arrangement, a synthetic chain connects a cantilever probe to a translated substrate. For the cantilever compliance correction, a deflection sensor gives $d$, a calibrated spring constant gives $F=k_{c}d$, and substrate displacement $Z$ gives chain extension $x=Z-d$ after geometric offsets are removed. Additional linkers and instrument compliance must be included if appreciable. Thermal deflection calibration uses $k_c⟨d^2⟩=k_{B}T$ only after isolating the cantilever's calibrated thermal contribution.

An optical-tweezers implementation detects a tethered bead's displacement from a calibrated trap centre, with $F=k_{\mathrm{trap}}\Delta x$ in the harmonic regime. Detector calibration, drag and trap linearity determine the usable force and frequency range. These instruments measure individual force–extension curves and distributions; averaging many curves can hide multiple tether populations.

For a synthetic freely jointed chain with Kuhn length $b$,A strongly stretched worm-like chain instead has $1-x/L\simeq \sqrt{k_{B}T/(4F{\ell }_p)}$ before backbone stretching becomes important. Fits can estimate contour and persistence lengths only in their appropriate regimes. Hysteresis and pulling-rate dependence test whether an apparent mechanical transition is quasistatic or activated.

Control single versus multiple attachments, subtract nonspecific adhesion, vary pulling speed, calibrate temperature and stiffness, and correct drift. A rupture peak is not a universal bond strength independent of loading history. A primary nonbiological example measured individual polydimethylsiloxane chains by atomic-force methods; its flexibility interpretation combines those measurements with molecular calculations.^[ref]

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.

Write $\epsilon ={\epsilon }_0{\epsilon }_r$ and take $n_0$ to be the reservoir density of each monovalent ion species. Equal electrochemical potentials in equilibrium give the Boltzmann ion densities in a symmetric electrolyte,Combining this with the electrostatic Poisson equation gives the planar Poisson-Boltzmann equation,The Debye length is ${\kappa }^{-1}$. If the symbol $n_0$ instead denotes the combined density of both species, the factor two is absent; the per-species convention above makes the normalization explicit. At $|e\psi |\ll k_{B}T$, the linearized planar screening solution obeys ${\psi }^{\prime \prime }={\kappa }^2\psi$. With $\psi (0)=\zeta <0$ and $\psi (\infty )=0$,

The positive ions are counterions and the negative ions are co-ions. At weak potential,They approach $n_0$ from above and below respectively, with equal first-order excess and deficit. For stronger negative potential, interpret the source's strong-potential comparison as a statement about $|e\psi |$, since $\psi$ itself is negative. Counterions become strongly concentrated near the wall, while co-ions are exponentially depleted and remain nonnegative. Linearization misses this asymmetry and would eventually predict unphysical negative co-ion density. Point-ion mean-field theory itself eventually needs finite-size and correlation corrections; it provides no steric cap at arbitrarily large counterion concentration.

For steady tangential electro-osmosis, with no imposed pressure gradient and negligible inertia, the Stokes equation is $\eta v^{\prime \prime }+{\rho }_{e}E=0$. Since ${\rho }_e=-\epsilon {\psi }^{\prime \prime }$,The conditions $v(0)=0$ and $v^{\prime }(\infty )=0$ giveThe printed constant speed is the bulk velocity outside the double layer, not the velocity at the no-slip wall. With negative $\zeta$ the positive counterion excess drives fluid along the applied field. The far-field value is the Helmholtz-Smoluchowski electro-osmotic velocity.^[ref]

For the finite channel, first model the added wall as an ideally uncharged insulating boundary with no normal field on its outer side: ${\psi }^{\prime }(d)=0$, while $\psi (0)=\zeta$. The one-charged-wall electro-osmotic channel then hasIntegrating the same force balance and enforcing $v(0)=v(d)=0$ givesFor $U>0$, the maximum occurs atIn that limit $v/U\simeq 1-x/d-e^{-\kappa x}$: the velocity rises rapidly from zero at the charged wall, reaches a maximum a few screening lengths away, then falls approximately linearly to zero at the neutral wall. One charged wall does not give a symmetric plug across a channel with no-slip at the other wall.

An alternative electrode condition $\psi (d)=0$ means a wall held at reservoir potential; it is not the same as zero surface charge. Under that convention,with $\cosh [\kappa (d-x_{\max })]=\sinh (\kappa d)/(\kappa d)$. Both boundary conventions have the same leading thin-double-layer profile above. An externally imposed pressure gradient would add a pressure-driven profile; none is specified here.

For dilute nonadsorbing uncharged polymers in equilibrium or sufficiently weak flow, the entropic polymer depletion near nonadsorbing walls suppresses chain configurations near each wall. A depletion layer of order the coil radius of gyration is expected at both walls, with a bulk-like central region when $d$ is much larger than the coil size. Strong confinement rounds this into a central enrichment. Neutrality removes a direct electrophoretic force, not excluded-volume entropy. Adsorption, concentration and strong-flow migration were not specified, so this is a conditional equilibrium sketch rather than a universal flowing-polymer profile.

The generic instrument uses a magnetically responsive bead, a controllable external field and optical tracking of bead position or orientation. A field gradient supplies translational force; rotating an appropriate anisotropic field or bead moment can supply torque. For a permanent dipole,For an isotropic unsaturated induced moment with effective susceptibility contrast $\Delta \chi$ and bead volume $V$, the induced magnetic bead force instead follows fromThe factor one-half accounts for inducing the moment. A spatially uniform field alone supplies no translational force; an isotropic induced moment parallel to the field supplies no orienting torque either. Saturation or magnetic anisotropy requires the corresponding magnetization model.

Position tracking needs a length calibration, and force calibration needs measured field gradients, a known reference load, or calibrated thermal fluctuations. For a generic tensioned tether of extension $\ell$, small transverse displacement $x$ costs energy $F{x}^2/(2\ell )$, so $F=k_{B}T\ell /⟨x^2⟩$ under equilibrium, small-angle and single-mode assumptions. Attachment geometry, bead rotation, wall drag and finite acquisition bandwidth can affect this calibration. These equations describe generic magnetic tweezers mechanics, without assigning a biological response.

TODO: biological single-molecule applications are user-deferred; no biological systems or protocols have been researched or solved.

Fluctuation-induced forces arise because separation changes the fluctuations available to a system and therefore changes its free energy. The force is $F=-\partial \mathcal{F}/\partial d$ at fixed relevant control variables, not necessarily a force from a permanent charge or dipole.

For two neutral polarizable particles, correlated electromagnetic fluctuations give a short-distance dispersion attraction of the form $-C_6/r^6$. Retardation changes the asymptotic distance dependence. Boundaries produce a related Casimir force; for parallel ideal conducting plates at zero temperature,The negative pressure denotes attraction. Material response, nonplanar geometry and finite temperature alter the ideal result. The effect can be described through the separation-dependent electromagnetic mode spectrum rather than by treating a divergent absolute vacuum energy as observable.^[ref]

In a classical depletion force, smaller particles or polymer coils are excluded from a region near larger surfaces. Overlap of excluded regions increases their available volume. For ideal dilute depletants of osmotic pressure $\Pi$, $U(d)=-\Pi V_{\mathrm{overlap}}(d)$; the resulting attraction is entropic and scales with temperature. Confined interface fluctuations can instead generate repulsion when close approach removes entropy. Near a continuous phase transition, a large correlation length makes boundary-conditioned order-parameter fluctuations produce a critical Casimir force. Its sign depends on the boundary preferences, so fluctuation forces are not universally attractive.

The useful common idea is a free-energy difference with explicit boundary conditions; quantum vacuum, thermal field and excluded-volume mechanisms have different physical fluctuations and different ranges.

Only the separable instrument and mechanical tasks are treated here; biological operation and stepping energetics remain deferred.

For generic optical tweezers, a tightly focused laser and objective with high numerical aperture create an intensity gradient that attracts a suitable dielectric bead toward the focus. Momentum transfer also creates a scattering force; stable three-dimensional trapping requires sufficient restoring gradient force. A condenser and photodetector or camera track displacement. Near a stable focus, $F_x=-kx$ defines the stiffness and requires both a position calibration and a force calibration.^[ref]

At thermal equilibrium in a harmonic well, the equipartition trap calibration is $k⟨(x-⟨x⟩{)}^2⟩=k_{B}T$. Alternatively, fitting the overdamped harmonic-trap power spectrum gives $f_c=k/(2\pi \gamma )$, hence $k=2\pi \gamma f_c$. Detector noise, finite sampling and exposure time, wall corrections to drag, and drift must be controlled rather than attributed to stiffness.

For an isolated sphere in a Newtonian fluid, define $\mathrm{Re}={\rho }_{f}v(2a)/\eta$ using its diameter. The supplied viscosity converts to $10^{-3}\mathrm{Pa}\mathrm{s}$. The generic Stokes drag calculation givesHere $\mathrm{Re}=2\times 10^{-6}[{\rho }_f/(10^3\mathrm{kg}{\mathrm{m}}^{-3})]$; this retains the unspecified fluid density as a parameter rather than treating a density as supplied data. Stokes drag assumes an isolated no-slip sphere; nearby surfaces change the coefficient. TODO: the requested assessment of biological motor operation is deferred.

TODO: the biological stepping-energy and step-detection estimate is deferred in its entirety; no biological energy source or motor mechanism has been researched.

Ignoring biological coupling, the mechanical thermal model is the Langevin equation for a harmonic trap,For independent repeats with the specified initial conditions, the ensemble mean $X=⟨x⟩$ obeys the homogeneous damped-oscillator equation. An average over $N$ independent zero-mean force samples reduces its variance by $N$, with a finite detector bandwidth needed to define samples of white noise. High bandwidth alone does not remove random forces; averaging independent records or sufficiently separated samples is the additional operation. A single trajectory retains thermal fluctuations. For stationary overdamped position fluctuations, a long time average over duration $\mathcal{T}\gg {\tau }_k=\gamma /k$ has approximately $\mathrm{Var}(\bar{x})\simeq 2(k_{B}T/k){\tau }_k/\mathcal{T}$.

The initial-value solution of a damped harmonic trap has characteristic roots $r_{\pm }=[-\gamma \pm \sqrt{{\gamma }^2-4mk}]/(2m)$. For distinct roots,For critical damping, ${\gamma }^2=4mk$,For the underdamped case the same real solution can be written $e^{-\gamma t/(2m)}[x_0\cos {\omega }_{d}t+(v_0+\gamma x_0/(2m))\sin {\omega }_{d}t/{\omega }_d]$, with ${\omega }_d^2=k/m-{\gamma }^2/(4m^2)$. These describe the deterministic response or conditional mean, not the entire noisy trajectory.

For overdamped trap relaxation after an inertial boundary layer, let $ϵ=km/{\gamma }^2\ll 1$. The slow and fast roots areAfter $t\gg {\tau }_m$, the fast exponential is negligible and the surviving amplitude is $x_0+(m/\gamma )v_0+O(ϵx_0,ϵm{v}_0/\gamma )$. If also $|(m/\gamma )v_0|\ll |x_0|$, then to leading overdamped orderSmall Reynolds number justifies the linear viscous force; it does not alone prove the separate small inertial-trap parameter $ϵ$. Neglecting $v_0$ requires both rapid momentum relaxation and a negligible slow-mode amplitude shift. If $x_0=0$, a nonzero initial velocity can leave a small slow displacement even after its own fast transient decays. The result concerns generic bead mechanics only.

We first derive the stress and energy facts used in both parts. The Cauchy traction $s=Tn$ is force per unit area in the deformed configuration. The first Piola-Kirchhoff traction $s_R=T_{R}N$ is force per unit reference area. Their force vectors are measured in the same spatial coordinates, but the area measures and the normals used to specify the surface are different.

Parametrize a reference surface by $x=x(u,v)$. Its oriented area vector is $(x_u\times x_v)dudv$. Under deformation, the two tangents become $F{x}_u,F{x}_v$. The cross-product identitytherefore proves Nanson's formula $nda=J{F}^{-T}NdA$. Equality of the physical force on corresponding surface patches now givesSince this holds for every reference normal,The first Piola-Kirchhoff stress tensor is generally not symmetric, whereas the Cauchy stress tensor is symmetric in the absence of couple stresses.

For the stated energy on orientation-preserving gradients, differentiation givesA stress-free state satisfies $F{F}^T=J^{-1}I$. Taking determinants yields $J^2=J^{-3}$, so $J=1$ and $F{F}^T=I$: precisely $F\in SO(3)$. These states are also the global energy minimizers. If $v_i>0$ are the singular values, the arithmetic-geometric mean inequality gives $|F{|}^2\ge 3J^{2/3}$, with equality only when the singular values are equal. The function $3J^{2/3}/2+J^{-1}$ decreases up to $J=1$ and increases thereafter, with minimum $5/2$. Thus $W\ge 5\mu /2$, with equality precisely at rotations. A constant shift of the energy can set this minimum to zero without changing the stress.

For the given simple shear,For a unit normal $N$, $|T_{R}N{|}^2={\mu }^2(N_1^2+N_2^2)={\mu }^2(1-N_3^2)$. Hence every reference plane whose normal lies in the $e_1,e_2$ plane is maximizing:There is a family of maximizing planes, rather than one distinguished normal.

A dead-load traction is prescribed force per unit reference area whose vector is independent of the deformation. It is imposed through $T_{R}N=t_R$ on the reference boundary. It does not rotate or rescale with the current surface as a follower load would.

Elastic incompressibility imposes $\det F=1$ and adds a scalar Lagrange multiplier to the constitutive law. With the pressure sign convention used here,The pressure field is determined by equilibrium and boundary conditions, rather than by the strain energy density alone.

For the diagonal homogeneous deformation, $F=\mathrm{diag}(v_1,v_2,v_3)$ and $v_1{v}_2{v}_3=1$. Isotropy makes $DW(F)$ diagonal here, with entries ${\Phi }_i$, soThe six face conditions are exactly $T_R{e}_i=\tau e_i$. Thus ${\Phi }_i-p/v_i=\tau$ for each $i$, orNecessity follows from the boundary conditions. Conversely these equations give a constant pressure and constant stress satisfying those conditions, while $\mathrm{Div}{T}_R=0$ because the stress is constant and there is no body force. This proves sufficiency as well.

For the particular energy, ${\Phi }_i=\mu (v_i-v_i^{-2})$. Put $S=v_1+v_2+v_3$ and $Q=v_1{v}_2+v_1{v}_3+v_2{v}_3$. Elastic incompressibility gives $v_i^{-1}=v_j{v}_k$, henceEquality for every pair is therefore equivalent to $(\mu S-\tau )(v_i-v_j)=0$. If all stretches are equal, their positive product one gives $v_1=v_2=v_3=1$, with $p=-\tau$. Otherwise $S=\tau /\mu$. Conversely any positive triple with this sum and product one solves the equations, with $p=-\mu Q$. Thus the complete solutions within the stated homogeneous diagonal ansatz are the unit triple and the stated sum-product family.

By the arithmetic-geometric mean inequality, $S\ge 3$, with equality only at the unit triple. Consequently $\tau \le 3\mu$ admits only that triple, including equality, when the two descriptions coincide.

To prove infinitely many for $s=\tau /\mu >3$, choose $v_1=t$ near one and setAt $t=1$ the discriminant is $(s-1{)}^2-4>0$; it stays positive on an interval about one, and both roots remain positive. Every such $t$ gives a different admissible triple. HenceThe degeneracy has an energy explanation: under the product constraint the strain energy density is $\mu S^2/2$, and the total potential under dead loading is $\mu S^2/2-\tau S$. It depends on the three stretches only through their sum. Its minimum over attainable $S\ge 3$ moves from $S=3$ to the whole level set $S=\tau /\mu$ when the load exceeds $3\mu$.

For a hyperelastic material $T_R(F)=DW(F)$. Material frame indifference of this first Piola-Kirchhoff law meansThe spatial force components rotate under a superposed rigid rotation, whereas the reference normal does not.

If $W(QF)=W(F)$, differentiate with respect to $F$ in direction $H$ to obtain $DW(QF):(QH)=DW(F):H$, hence the displayed stress transformation. Conversely, assume that stress transformation. For fixed $Q$, the derivative of $W(QF)-W(F)$ with respect to $F$ vanishes. The positive-determinant matrix domain is connected, as follows from polar decomposition and the connectedness of $SO(3)$ and of the positive-definite cone. Thus $W(QF)-W(F)=c(Q)$, independent of $F$. Iterating givesThe assumed lower bound on $W$ forces $c(Q)=0$: either nonzero sign would tend to minus infinity for one direction of integer $n$. This proves the required equivalence and explains the use of the lower-bound hypothesis.

In the frame-indifferent setting, isotropy means invariance under reference rotations as well: $W(FQ)=W(F)$ for $Q\in SO(3)$. The distinction matters: left rotations encode changes of observer, and right rotations encode material symmetry.

An orientation-preserving singular value decomposition is $F=R\mathrm{diag}(v_1,v_2,v_3)S$, with $R,S\in SO(3)$ and $v_i>0$. Material frame indifference and isotropy remove $R,S$, leaving a function $\Phi$ of the stretches. Permuting the stretches using signed permutation rotations does not change the value, so $\Phi$ is symmetric. Conversely a symmetric stretch function is unchanged by both left and right rotations, in particular is isotropic.

The eigenvalues of $B=F{F}^T$ are $v_i^2$. Its three principal invariants areThey determine the unordered positive eigenvalues as the roots of ${\lambda }^3-I{\lambda }^2+II\lambda -III=0$. Therefore a symmetric stretch function defines a well-defined $h(I,II,III)$ on the admissible invariant triples. Conversely a function of these invariants is a symmetric stretch function, proving the equivalence of all three descriptions.

The printed domain of $\Phi$ must be $(0,\infty {)}^3$, since it takes three arguments. Also isotropy alone as right invariance would not imply the spectral representation without material frame indifference: $\mathrm{tr}(AF{F}^T)$ for a fixed nonscalar symmetric positive-definite $A$ is right invariant but not generally objective. The equivalence here uses the frame-indifferent constitutive setting.

For the specific energy, Newton identities for the three eigenvalues gives ${\sum }_i{v}_i^6=I^3-3III+3III$, while $v_1{v}_2{v}_3=J=\sqrt{III}$. ThusTo differentiate it, useThese follow from $I=|F{|}^2$, the displayed trace formula for $II$, and $III=J^2$. The three partial derivatives of $h$ are $h_I=\alpha +3\beta I^2-3\beta II$, $h_{II}=-3\beta I$ and $h_{III}=3\beta -\gamma II{I}^{-3/2}/2$. The chain rule therefore givesFinally $-6\beta II=-3\beta (\mathrm{tr}B{)}^2+3\beta \mathrm{tr}(B^2)$, and $BF=F{F}^{T}F$. This is precisely the requested expression, with the coefficient of $F$ separated from the $BF$ and cofactor terms. The positive sign of the $6\beta J$ term comes from the $3\beta III$ contribution in the energy.

In the classical smooth setting, a deformation in continuum mechanics is a sufficiently regular, injective map $y:\Omega \to {\mathbb{R}}^3$ with $\det Dy>0$. Injectivity excludes material interpenetration, and a positive determinant preserves orientation and local volume. We use $C^1$ regularity for these kinematics and additional differentiability where stress-divergence identities are applied.

Cauchy's stress theorem states that the contact force on an oriented surface element at a current point depends on the point and its unit normal, rather than on the shape of the chosen subbody. Thus the contact force on the boundary of a current subbody $v$ has the form ${\int }_{\partial v}t(z,n)da$. Static force balance, for a body force $b$ per unit current volume, isfor every subbody. A shrinking pillbox gives $t(z,-n)=-t(z,n)$; a shrinking tetrahedron then gives linear dependence on the normal, $t(z,n)=S(z)n$. Applying the divergence theorem and using arbitrary subbodies givesIf body force is specified per unit mass instead, replace $b$ by $\rho b$. Symmetry of $S$ follows additionally from conservation of angular momentum in the absence of body couples; it is not a consequence of force balance alone.

The surface element relation is Nanson's formula. The force on a corresponding pair of current and reference surface elements must agree:Therefore the first Piola-Kirchhoff stress tensor isIt measures force per reference area, with a current force direction and a reference normal direction.

In components, differentiate $T_{iA}=S_{ij}(y(x))(\mathrm{cof}F{)}_{jA}$. The Piola identity says ${\partial }_A(\mathrm{cof}Dy{)}_{jA}=0$. One proof writes the cofactor matrix as $\frac{1}{2}{\epsilon }_{jkl}{\epsilon }_{ABC}{y}_{k,B}{y}_{l,C}$; differentiating in $A$ produces symmetric second derivatives contracted against antisymmetric epsilon symbols, hence zero. The chain rule consequently givessince $F_{kA}(\mathrm{cof}F{)}_{jA}=J{\delta }_{kj}$. Thus, with the coordinate composition made explicit,The distinction between reference divergence and current divergence is essential to this identity.

The Green-Lagrange strain tensor has eigenvalues $({\lambda }_j^2-1)/2$. Its Frobenius norm and trace are invariant under the orthogonal principal-axis change, soFor the homogeneous dilation $F=\rho I$, put $K=\mu +3\lambda /2$. Then $E=({\rho }^2-1)I/2$. Differentiating the energy through $E$ gives the second Piola-Kirchhoff stress tensor $2\mu E+\lambda \mathrm{tr}(E)I=K({\rho }^2-1)I$, and multiplication by $F$ gives the first Piola-Kirchhoff stress tensor. ThusThe normals of the reference and deformed spheres point in the same radial direction. As $\rho \downarrow 0$, the reference traction magnitude $|TN|$ tends to zero, while the current traction magnitude $|Sn|$ tends to infinity, provided $K>0$, as for physically stable Lamé parameters. The negative signs of the stress coefficients for $0<\rho <1$ indicate compression.

There is no contradiction: the current area is ${\rho }^2$ times the reference area, so the same force divided by the vanishing current area has the larger stress. Indeed $TN={\rho }^{2}Sn$. The limiting map at $\rho =0$ is not itself an admissible deformation in continuum mechanics; the statement concerns the approach through positive $\rho$. A degenerate choice $K=0$ would not have the asserted growth without bound.

The displacement field is $u(x)=y(x)-x$, and its gradient is $Du=Dy-I$. Linear elasticity expands about an undeformed stress-free reference state, taking $y=x+ϵv$ with $|Du|\ll 1$ and external loads of the same small order. A rigid translation need not itself be small; it is the deformation gradient that controls the infinitesimal strain tensor.

For $W\in C^2$ with $DW(I)=0$, a Taylor series givesThe constant $W(I)$ can be set to zero. The tensor $\mathbb{A}$ is the elasticity tensor. Commutation of continuous mixed partial derivatives proves its major symmetry:Material frame indifference at a stress-free reference also makes the quadratic energy depend only on the infinitesimal strain tensor $e=(H+H^T)/2$, leading to the usual minor symmetries after spatial and reference indices are identified.

In this scaling, an external force of order $ϵ$ may be evaluated on the reference geometry: its change due to the order-$ϵ$ displacement contributes only at order ${ϵ}^2$. Thus the leading forces are fixed dead loads. Likewise $F=I+O(ϵ)$, $J=1+O(ϵ)$ and $S=O(ϵ)$, so the Piola transformation gives $T=JS{F}^{-T}=S+O({ϵ}^2)$. Hence Cauchy and first Piola stresses agree to linear order. This reasoning uses the usual small-load, stress-free linearization; an order-one prestress or leading follower-load effect requires a different expansion.

Taking traces of the constitutive equation givesThe assumed finite limits of stress and strain make the right side bounded. Therefore $\mathrm{tr}{e}_{\lambda }=O({\lambda }^{-1})$, and the limiting strain satisfies $\mathrm{tr}e=0$. The permitted interchange of derivatives and limits gives $\mathrm{tr}e=\mathrm{div}u$. DefineThen the limit stress is $S=-pI+2\mu e$. Passing to the limit in force balance and using $2\mathrm{Div}e=\Delta u+\nabla \mathrm{div}u$ givesThe pressure acts as the Lagrange multiplier enforcing the volume constraint. The differential equations determine it up to an additive constant, while a specified absolute limit stress can fix that constant.

Nonlinear elastic incompressibility is $\det Dy=1$. Expanding the determinant for $Dy=I+Du$ gives $\det (I+Du)=1+\mathrm{div}u+O(|Du{|}^2)$, so its first-order form is precisely $\mathrm{div}u=0$. The volumetric energy term $\lambda (\mathrm{div}u{)}^2/2$ becomes increasingly expensive as $\lambda$ grows, while the fixed positive $\mu$ leaves shear strain possible. Thus the infinite-bulk-stiffness limit is the incompressible limit of isotropic linear elasticity, with the same equilibrium operator as the steady Stokes equations.

The first equation is the local balance of linear momentum. Use $\mathsf{M},\mathsf{L},\mathsf{T}$ for the dimensions of mass, length, and time, to distinguish dimensional time from the Cauchy stress tensor $\mathbf{T}$.

Here $t$ has dimension $\mathsf{T}$. The hypothesis ${\partial }_t\mathbf{b}=0$ means that the prescribed spatial body force field is the same for every member of the static family and the resulting motion.

The electron-phonon interaction can produce a retarded attraction. A carrier distorts the ionic lattice; a second carrier can lower its energy in that distortion. Integrating out a phonon gives a frequency-dependent interaction with an attractive low-frequency region. A schematic Matsubara-frequency contribution is $V_{\mathrm{ph}}(\mathbf{q},i{\omega }_n)\propto -|g_{\mathbf{q}}{|}^2{\omega }_{\mathbf{q}}/({\omega }_n^2+{\omega }_{\mathbf{q}}^2)$. It must be combined with the screened Coulomb repulsion, not treated as a statement that the bare electron-electron force is attractive.

Retardation separates the electronic and phonon energy scales and reduces the repulsive contribution in the pairing channel to a Coulomb pseudopotential. In the simplest weak-coupling approximation, a positive net dimensionless pairing strength givesThe formula assumes a phonon cutoff and simplified coupling; strong-coupling calculations need the spectral interaction rather than this one-number replacement. The phonon frequency's dependence on ionic mass motivates an isotope effect, whose actual exponent can be modified by other physics.

An interaction need not be attractive in every momentum channel to support pairing. Spin fluctuations or other electronic correlations can favor a sign-changing gap: in the gap equation, repulsive scattering between two regions can contribute constructively when their gap signs differ. Such mechanisms are candidates for unconventional superconductivity. Observing a d-wave gap identifies symmetry information but does not, by itself, establish a unique microscopic pairing mechanism in a cuprate.

A neutral superfluid has a complex order parameter whose phase describes coherent particle transport. For a Bose condensate it is represented by a field amplitude; for a paired Fermi superfluid it is a pair amplitude. In a symmetry-breaking description, selecting its phase breaks the global particle-number $U(1)$ symmetry. A common phase shift leaves the energy unchanged, while a phase gradient costs energy and produces superfluid velocity. Long-wavelength phase oscillations give a gapless Goldstone mode, generally a sound wave when density fluctuations provide the conjugate dynamics. A finite isolated state with exact particle number need not have a nonzero field expectation: the phase-selected language describes the ordered limit and relative-phase coherence.

In a superconductor, the paired order parameter carries charge of magnitude $2e$ and couples to the electromagnetic field. The meaningful spatial phase combination is $\hslash \nabla \theta -q\mathbf{A}$, with signed pair charge $q$. Gauge symmetry is a redundancy, so “breaking local gauge symmetry” is not a literal breaking of a measurable physical symmetry. The gauge-invariant electromagnetic response includes the Meissner effect and a finite London penetration depth. Long-range Coulomb forces move the would-be neutral phase mode to a plasma scale through the Anderson-Higgs mechanism in a superconductor.

Additional symmetries can be broken by particular pairing states. Conventional isotropic spin-singlet pairing need not break spatial rotations or time reversal. A d-wave superconductor has a gap transforming nontrivially under crystal symmetries, and some triplet or multicomponent states can also break spin-rotation or time-reversal symmetries. These additional choices are not required by superconductivity itself. A nonzero pairing amplitude, phase coherence and electromagnetic screening describe related but distinct aspects of the ordered state.

A neutral superfluid can select a phase of its complex order parameter, breaking a global $U(1)$ symmetry. A spatially uniform change of that phase costs no energy. Slowly varying changes therefore produce a gapless Goldstone mode. For a dilute weakly repulsive Bose gas, the Bogoliubov quasiparticle dispersion isso $E_k\simeq \hslash c_{s}k$ at long wavelengths, with $c_s=\sqrt{g_B{n}_0/m_B}$. Density and phase fluctuations are canonically coupled components of this sound branch, rather than two independent positive-frequency modes. Superfluid helium and dilute atomic condensates supply neutral examples, although the dilute-gas formula is not quantitatively valid for strongly interacting liquid helium.

In a superconductor, the pair field is charged and its phase couples to the electromagnetic gauge field. Only combinations such as $\nabla \theta -q\mathbf{A}/\hslash$ are physically meaningful. The Anderson-Higgs mechanism in a superconductor produces electromagnetic screening and removes the neutral sound mode as an independent gapless charged-density oscillation. In a homogeneous three-dimensional continuum, long-range Coulomb forces give a longitudinal plasma gap, approximately ${\omega }^2={\omega }_p^2+c_s^2{k}^2$ at small $k$. The static transverse response gives the Meissner effect andwhen $n_s$ counts superconducting electrons of effective mass $m$. Counting pairs instead uses density $n_s/2$, charge $-2e$ and mass $2m$, giving the same stiffness. The condensate amplitude may also have a gapped Higgs mode, but that amplitude oscillation is distinct from the electromagnetic Higgs mechanism. Local gauge transformations remain a redundancy; this is not spontaneous breaking of an observable local symmetry in the same sense as the neutral global phase symmetry.

The superconducting critical current depends on temperature, magnetic field, sample geometry and the chosen dissipation criterion. For a magnet, a large superconducting critical current density matters because it permits a useful winding current with manageable conductor volume. The intrinsic depairing limit is an upper bound, but many practical type-II conductors become dissipative first through vortex motion.

A current exerts a transverse force per unit vortex length ${\mathbf{f}}_L=\mathbf{J}\times {\mathbf{\Phi }}_0$. Moving flux produces an electric field and dissipation. Flux pinning holds vortices at energetically favorable defects, allowing a larger transport current before motion. Useful conductors therefore need pinning as well as a sufficiently high transition temperature and upper critical field. Simply increasing disorder indefinitely is not a prescription: it can damage connectivity or superconductivity. High-angle grain boundaries in some copper-oxide materials also act as weak links, motivating controlled texture and connected current paths.

Conventional magnet conductors include niobium-titanium and niobium-tin; high-temperature coated conductors offer other operating regimes. Their performance must be compared on a superconducting critical surface in temperature, field and current density, with strain and thermal constraints also relevant. A high self-field current measured in a short sample does not by itself establish the current available throughout a high-field coil.

A superconducting magnet produces a strong, stable field without the normal winding's steady resistive power loss. Applications include magnetic resonance imaging, nuclear magnetic resonance and charged-particle beam steering or focusing in accelerators. A persistent-current loop is valuable when field stability matters. Refrigeration still requires power, and the magnet stores magnetic energy, approximately $L{I}^2/2$ in a linear inductance model. A local transition to a resistive state can release that stored energy as heat, so superconducting quench detection and protection are essential to operation. These requirements connect the microscopic vortex physics with useful engineering performance.

Let $m$ denote the constant rest mass, set $x(0)=0$, and put $a=f_x/m$. The relativistic three-force is $f_x=d{p}_x/dt$, soIt follows thatIntegrating the velocity with the initial position fixed givesThese formulas allow either sign of $a$; their $a\to 0$ limit is rest at the origin. For $a\ne 0$, the trajectory obeysThis is hyperbolic motion. Its coordinate acceleration decreases while its proper acceleration has constant magnitude $|a|$. In particular, a constant three-force does not give constant coordinate acceleration at relativistic speed.