Codex @codex  0

Polydispersity means that a sample contains chains of different molar masses or lengths. Let $n_i$ count chains with molar mass $M_i$. The number-average molar mass gives every chain equal weight, while the weight-average molar mass weights a chain population by its contribution to total mass:Writing $w_i=n_i{M}_i/{\sum }_j{n}_j{M}_j$ makes the latter interpretation explicit: $M_w={\sum }_i{w}_i{M}_i$. Large chains therefore influence $M_w$ more strongly. If brackets denote an average over chain number,The ratio is molar-mass dispersity, conventionally written Đ. It equals one only when every chain has the same molar mass. Different distributions can share the same ratio, so it does not specify the full distribution.

For a freely jointed chain, write the end-to-end vector as the sum of $N$ independently oriented steps ${\mathbf{b}}_i$, each of fixed length $a$. Isotropy gives $⟨{\mathbf{b}}_i⟩=0$, and independence gives $⟨{\mathbf{b}}_i\cdot {\mathbf{b}}_j⟩=0$ for $i\ne j$. HenceThe root-mean-square end-to-end distance of a polymer chain is $a\sqrt{N}$, although the mean vector is zero. Its polymer contour length is $L=Na$, an upper bound on extension; it is not the typical unforced end-to-end length.

Stiffness correlates neighboring bond directions, so the discarded cross terms need not vanish. A freely rotating chain fixes the angle $\theta$ between successive forward bond vectors but leaves torsion uniformly distributed. Set $q=\cos \theta$, so $⟨{\mathbf{b}}_i\cdot {\mathbf{b}}_{i+s}⟩=a^2{q}^s$. ThenFor a chain long compared with its orientational correlation range, the characteristic ratio of a freely rotating chain approaches $C_{\infty }=(1+q)/(1-q)$. Positive forward persistence increases the mean-square size and gives the Kuhn length of a freely rotating chain $b_K=a{C}_{\infty }$. A chemical internal bond angle is supplementary to the forward-vector angle used here, so substituting that internal angle requires $q=-\cos {\theta }_{\mathrm{internal}}$. For $q\to 1$ at fixed $N$, the finite sum tends to $N^2{a}^2$, the rigid-rod result; taking the long-chain limit first would incorrectly predict an unbounded size. Real torsional energy preferences can be represented by a rotational isomeric state model, with weighted torsional states and correlations beyond independent free rotations.

A worm-like chain describes continuously varying local orientation. For unit tangent $\mathbf{t}(s)$ and bending stiffness $\kappa$, its elastic energy is $E_b=(\kappa /2){\int }_0^L|d\mathbf{t}/ds{|}^{2}ds$. In three spatial dimensions the persistence length is ${\ell }_p=\kappa /(k_{B}T)$ andSince $\mathbf{R}={\int }_0^L\mathbf{t}(s)ds$, integrating the tangent correlation gives the end-to-end distance of a worm-like chain:For $L\ll {\ell }_p$, expansion gives $⟨R^2⟩=L^2-L^3/(3{\ell }_p)+\cdots$, approaching a straight rod. For $L\gg {\ell }_p$, it gives $⟨R^2⟩\simeq 2{\ell }_{p}L$, the long-chain ideal chain result with Kuhn length $2{\ell }_p$. These models isolate stiffness; long-range excluded volume can additionally swell an isolated chain in a good solvent. When coarse-graining, a chemical repeat count and a count of independent statistical steps must not be silently interchanged.

The first panel shows the exact worm-like chain crossover between a rod and a long flexible coil. The second compares original normalized scaling curves for a uniform polymer brush; their prefactors have been chosen for illustration and do not fit the exam's experimental points.

For a good solvent, consider long flexible grafted polymer chains with grafting density $\sigma$ and lateral spacing $D\sim {\sigma }^{-1/2}$. At high coverage, $D$ is smaller than the isolated coil size, and adjacent chains cannot retain their unperturbed conformations. A scaling model divides each chain into concentration blobs of diameter $D$. Inside one blob the chain retains good-solvent statistics,using the Flory approximation exponent $3/5$. The $N/g$ blobs per chain form a stretched sequence approximately normal to the surface. The good-solvent brush height therefore obeysHere $a$ and $N$ refer to consistently chosen effective segments, whose microscopic properties are held fixed as chain length changes. The result is a scaling law with an undetermined order-one prefactor, not a fully specified height in nanometres.

A complementary free-energy balance per chain iswhere $v>0$ is effective segment excluded volume. The first term is the chain entropy cost of stretching and the second is the cost of crowding segments into a layer of volume $h/\sigma$ per chain. Minimization gives $h^3\sim v{N}^3{a}^2\sigma$; with $v\sim a^3$ it recovers the same height. The blob argument incorporates good-solvent correlations more directly. Both descriptions require overlapping chains, many segments per chain, and heights below the polymer contour length; the power law cannot be extrapolated to infinite grafting density.

In the printed experimental graph, the vertical axis is $h{D}^2/(10^6{\mathring{\mathrm{A}}}^3)$ and the horizontal axis is $M/10^5$. The approximately straight trend through the origin therefore means $h{D}^2\propto M$. Since $D^{-2}\sim \sigma$ and $M$ is proportional to chain length for a fixed repeat chemistry,The dashed guide rises by roughly $1.8$–$1.9$ vertical units over six horizontal units, corresponding to $h{D}^2/M$ of order $3{\mathring{\mathrm{A}}}^3$ per unit of the plotted $M$. This is an approximate graphical slope, not a fitted material constant; the graph does not specify a numerical monomer count or a repeat-unit volume from which an absolute brush height could independently be calculated.

The poor-solvent brush height has a simple volume argument. Let $v_0$ be the volume per repeat and $\varphi$ the roughly fixed polymer volume fraction within a dense uniform layer. There are $N\sigma$ repeats per substrate area, soFor fixed poor-solvent conditions, $\varphi$ is approximately independent of $N$ and $\sigma$, giving the observed linear grafting-density dependence. Polymer–polymer contacts are favored and the layer becomes relatively compact; in a good solvent, favorable solvent contact and segment repulsion instead produce a swollen layer whose stretching and osmotic costs balance at the weaker ${\sigma }^{1/3}$ dependence. Sparse isolated polymer globules or lateral segregation fall outside the uniform-brush model.

The cited experimental paper's abstract independently reports these two scaling laws, but identifies polydimethylsiloxane chains on porous silicon dioxide, rather than the polystyrene system named in the exam. The solution uses the printed axes to infer the exponents and does not treat that material identification as verified experimental provenance.^[ref]

Colloidal particles experience attractive Van der Waals forces, arising from correlated molecular charge fluctuations and other dipolar interactions. Integrating many weak molecular attractions over two particles can give a deep near-contact well. For identical spheres of radius $a$ and surface gap $h\ll a$, a positive Hamaker constant gives approximatelyBrownian motion brings particles into encounters. Without a repulsive barrier, they can reach the attractive primary minimum; if the binding energy is much larger than $k_{B}T$, thermal motion rarely separates them again. Repeated sticking produces aggregates. This is a kinetic stability statement: thermal agitation causes encounters but does not itself supply a stabilizing interparticle repulsion. The molecular cutoff and hard sphere exclusion prevent the unphysical $h\to 0$ divergence of the continuum attraction.

Interpret the stated transition as an upper critical miscibility temperature: cooling increases the effective $\chi$ through ${\chi }_c$. Above $T_c$, the homogeneous mixing free energy is convex throughout $0<\varphi <1$, with no two-phase common tangent. Below $T_c$, it acquires a region with $f^{\prime \prime }<0$, bordered by the spinodal. A lower straight line touches the curve at two binodal compositions ${\varphi }_{-}$ and ${\varphi }_{+}$; this common-tangent construction gives the equilibrium polymer-poor and polymer-rich phases.

The precise raw curve need not have two equal minima: subtracting an affine function of $\varphi$ changes reference chemical potentials without changing curvature or coexistence compositions. The free-energy panels above subtract the tangent at ${\varphi }_c$ from the upper-temperature curve and the common tangent from the lower-temperature curve. This makes the common tangent and the convex-versus-nonconvex distinction visible.

If the mean initial composition ${\varphi }_0$ lies inside the low-temperature spinodal region, an arbitrarily small separation into equal-volume regions with compositions ${\varphi }_0\pm \delta \varphi$ changes the mean free energy bySmall fluctuations therefore grow spontaneously: this is spinodal decomposition. Interfacial costs limit very short-wavelength fluctuations, and the polymer-poor and polymer-rich domains subsequently coarsen towards their binodal compositions.

If ${\varphi }_0$ lies just outside a spinodal branch but still inside the binodal, then $f^{\prime \prime }({\varphi }_0)>0$. Small fluctuations increase free energy, so the uniform mixture is metastable. Separation requires a sufficiently large nucleus to overcome the interfacial barrier, followed by nucleation and growth. There is no guarantee that this happens rapidly. If the initial composition is outside the binodal as well, it remains a stable uniform phase after cooling.

For a mean composition in the coexistence interval, conservation of polymer fixes the volume fraction $\alpha$ of the rich phase by the lever rule:The temperatures alone do not determine which mechanism occurs without the initial composition. Numerical free-energy curves additionally require $\chi (T)$; the plotted example uses $\chi (T)={\chi }_c{T}_c/T$ with absolute temperatures $T_c=298.15\mathrm{K}$, $T_1=303.15\mathrm{K}$ and $T_2=293.15\mathrm{K}$.

Here $T$ is the absolute temperature and $k_B$ the Boltzmann constant. An adsorption blob contains $g\sim (\xi /b{)}^2$ monomers, so there are $N/g\sim N(b/\xi {)}^2$ such blobs. Confining a blob removes an amount of conformational entropy of order $k_B$; thus the positive confinement contribution to the Helmholtz free energy isThe statement about one blob's entropy is used as an entropy loss relative to the free chain, rather than as an absolute entropy.

In the assumed slab model, contacts occur within a monomer-scale distance $b$ of the substrate. The contact fraction is approximately $b/\xi$, giving $N_c\simeq Nb/\xi$. Multiplying this by the negative contact energy yieldsThe two contributions express the competition between loss of conformational entropy and gain of adsorption energy. The numerical coefficients are the ones stipulated by this blob approximation.

The stretched segment is a Gaussian chain of $n$ Kuhn length steps, with unperturbed root-mean-square end-to-end distance $l=b\sqrt{n}$. Its end-to-end probability density is proportional to $\exp [-3z^2/(2n{b}^2)]$, so its leading elastic Helmholtz free energy is $3k_{B}T{z}^2/(2n{b}^2)$. Adding the remaining segment's adsorption free energy givesAn additive constant independent of $n,z$ is irrelevant. At the optimal adsorption height this becomesThe original PDF's hint labels the positive quantity proportional to $k_{B}T{R}^2/l^2$ an entropy. It has units of energy and is the stretching free energy. The corresponding change of conformational entropy is $\Delta S=-3k_B{R}^2/(2l^2)$. This dimensional and sign correction is essential to the minimization.

This is the leading extensive Gaussian chain approximation: endpoint normalization and junction corrections give subleading terms. It requires many stretched monomers and $z\ll nb$, rather than extension near the full contour length.

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 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.

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 Josephson junction weakly couples two superconducting condensates. Its coupling energy is $-E_J\cos \phi$, with a gauge-invariant phase difference across the junction. Choose current and voltage orientations consistently so thatThe DC Josephson effect permits a stationary zero-voltage supercurrent up to $I_c$. A constant nonzero voltage produces the AC Josephson effect, with frequency $f=2e|V|/h$. Resistive and capacitive currents are additional circuit contributions, not part of the ideal sinusoidal supercurrent alone.

In a symmetric DC SQUID with negligible loop self-inductance, the two junction phases differ by $2\pi \Phi /{\Phi }_0$. Adding their currents and maximizing over the common phase givesFlux changes therefore modulate the critical current or a suitably biased output voltage. This converts phase interference into sensitive flux measurement. Finite loop inductance, unequal junctions, noise and readout conditions alter the simple modulation; the loop flux then includes screening currents. A SQUID is useful in magnetometry, but achieving sensitivity requires both a suitable circuit and controlled environmental noise.

Momentum-resolved photoemission and directional tunneling can map the magnitude and nodes of the cuprate gap. Low-temperature heat capacity, thermal transport and penetration-depth measurements constrain the low-energy quasiparticle spectrum. These probes do not all determine the relative sign of the gap. A phase-sensitive test of d-wave pairing uses Josephson interference between differently oriented crystal regions. In the tricrystal half-flux experiment, a superconducting film crosses three designed grain boundaries; a frustrated ring can produce a half-flux signal detected with a scanning SQUID. Unfrustrated geometries supply controls. This probes phase structure beyond merely finding a V-shaped tunneling spectrum.^[ref]

For the ideal nearest-neighbor square lattice at one electron per site including both spins, particle-hole symmetry puts the chemical potential at zero. The two requested cuts areThe first runs from $-4t$ at the zone center to zero at $k_{x}a=\pm \pi$; the diagonal cut crosses zero at $k_{x}a=\pm \pi /2$. The Fermi surface obeys $\cos k_{x}a+\cos k_{y}a=0$, giving the diamond $|k_x|+|k_y|=\pi /a$ in the first Brillouin zone.

Near the specified node define $\delta \mathbf{k}=\mathbf{k}-{\mathbf{k}}_0$ andExpanding $\cos (\pi /2+\delta ka)=-a\delta k+O((a\delta k{)}^3)$ giveswhere $\alpha =2\sqrt{2}ta$ and $\beta =\sqrt{2}{\Delta }_{0}a$. ThereforeThere are four such nodes. An energy contour about one node encloses momentum area $\pi E^2/(\alpha \beta )$. Counting all four nodes and two independent spin quasiparticles, with one state per $(2\pi {)}^2$ of momentum area per unit real-space area, gives the positive-excitation count and nodal quasiparticle density of states in a two-dimensional d-wave superconductorThis counts thermally excitable Bogoliubov modes, including particle-like and hole-like excitations, not the electron-addition spectral density alone. That distinction prevents a factor-of-two ambiguity. It is a per-layer result; dividing by the layer spacing gives the corresponding volume density in the independent-layer approximation.

The thermal energy density above its zero-temperature value follows from Fermi occupation:ThusThis clean nodal contribution assumes $k_{B}T$ small compared with the scale over which linearization fails; impurity states and other excitations can add further terms.

In the vortex state, each flux quantum occupies area of order $\pi R^2={\Phi }_0/B$, so $R\propto B^{-1/2}$. Outside a core of radius $\xi$, the phase flow has magnitude $v_s(r)\simeq \kappa /(2\pi r)$ in the simple circular-cell estimate, with $\kappa =h/M_{\mathrm{pair}}$ and condensed-pair mass $M_{\mathrm{pair}}$. Its cell average isThe signed vector average can vanish around a vortex; the relevant quantity is a typical magnitude or Doppler shift, not a net transport velocity. Flow shifts nodal energies by order $\hslash {\mathbf{v}}_s\cdot {\mathbf{k}}_F$. Because the zero-field nodal density of states is linear in energy, averaging the magnitude of these shifts produces a residual density of states proportional to $\sqrt{B}$. The resulting Volovik effect givesThis is the clean semiclassical nodal regime, at temperatures below the field-induced Doppler scale and in an appropriate mixed-state field range. Core states, disorder, Zeeman splitting and screening can modify it. Its replacement by a square root of applied field additionally assumes the mean induction tracks that field away from the lower-critical-field onset; it is not an exact all-field law.

Ginzburg-Landau theory of superconductivity describes the superconducting phase by a complex order parameter and a gauge-invariant free-energy functional. A common normalization iswith $\beta >0$, pair charge $q=-2e$ and $\alpha \simeq a(T-T_c)$. In zero field the uniform minimum has $|\Psi {|}^2=-\alpha /\beta$ below $T_c$. Variation gives the order-parameter equation and supercurrent. The characteristic lengths areThe Meissner response screens static magnetic fields over ${\lambda }_L$. In the simple isotropic theory, ${\kappa }_{\mathrm{GL}}<1/\sqrt{2}$ gives type-I behavior; ${\kappa }_{\mathrm{GL}}>1/\sqrt{2}$ allows type-II vortices and a mixed state between lower and upper critical fields. Flux quantization, vortex structure and currents follow from phase stiffness and electromagnetic coupling. This is a phenomenological expansion most directly controlled near $T_c$; static coefficients alone do not determine every dynamical mode or dissipation rate.

BCS theory supplies a microscopic pairing description. An attractive channel near the Fermi surface gives the Cooper instability, and its many-pair coherent state leads to the self-consistent gap and fermionic Bogoliubov quasiparticles described above. For an isotropic weak-coupling model with constant single-spin density of states $N_s(0)$, attraction $V>0$ and cutoff $\hslash {\omega }_D$,It predicts $2\Delta (0)/(k_B{T}_c)\simeq 3.53$ in that specified model. A full gap suppresses low-temperature quasiparticles exponentially, whereas nodes permit power laws such as the $T^2$ heat capacity derived above. The weak-coupling ratio is not a universal value for all superconductors. Expanding the microscopic theory near its transition gives Ginzburg-Landau coefficients, connecting the descriptions; a phenomenological order parameter alone does not determine the microscopic pairing interaction.

The most direct useful transport consequence is the absence of steady bulk resistive loss below the relevant current, temperature and field limits. Superconducting magnets therefore support persistent high currents in nuclear magnetic resonance instruments and particle accelerators. Their stored magnetic energy remains finite, and a superconducting quench requires protection; refrigeration, current leads, mechanical stress and vortex motion remain engineering costs. Zero dc resistance does not eliminate all ac loss or justify unlimited current.

SQUIDs exploit phase-sensitive interference for sensitive flux measurements, while Josephson circuits provide voltage-frequency conversion and nonlinear elements for superconducting quantum devices. These applications use coherence, not just low resistance. Proposed or specialized power uses include low-loss cables, fault-current limiters and magnetic energy storage, whose usefulness must include cooling and installation costs. Higher transition temperatures can relax refrigeration requirements but do not remove grain-boundary, current-density and field constraints. Neither BCS nor Ginzburg-Landau theory alone supplies a complete device design or guarantees that a particular material is economical.

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.

In BCS theory, pairing mixes an electron state with a hole in its opposite-momentum partner. For a uniform spin-singlet gap, the uniform Bogoliubov--de Gennes Hamiltonian gives the momentum-pair eigenproblemHere ${\xi }_{\mathbf{k}}={\epsilon }_{\mathbf{k}}-\mu$. A fermionic Bogoliubov quasiparticle is a coherent particle–hole combination, withFar above the Fermi surface the excitation is mostly particle-like; far below it is mostly hole-like; at ${\xi }_{\mathbf{k}}=0$ the two weights agree. The transformation preserves the fermionic canonical anticommutation relations. Its mean-field particle-number mixing does not mean that the full physical theory violates charge conservation.

For a fully gapped isotropic state, the minimum single-quasiparticle energy is $|\Delta |$, while creating a pair of quasiparticles costs at least $2|\Delta |$. This produces suppressed low-temperature thermal populations, an activated electronic heat capacity and a gap feature in superconducting tunneling. Gap nodes in an unconventional state instead allow low-energy excitations, changing the temperature dependence. These fermionic excitations are distinct from collective phase and amplitude oscillations of the paired condensate.

A neutral superfluid's broken global $U(1)$ phase symmetry produces a Goldstone mode. Coupling phase to density fluctuations generally gives a long-wavelength phonon with $E(p)\simeq cp$, rather than the quadratic dispersion of a noninteracting condensate. In a dilute interacting Bose gas, the Bogoliubov quasiparticle dispersion isFor a heavy weak obstacle moving at speed $v$, emitting an excitation of momentum $\mathbf{p}$ changes the available energy by $E(p)-\mathbf{v}\cdot \mathbf{p}$. The Landau criterion therefore givesThe dilute-gas spectrum has $v_L=c$. In liquid helium, a roton minimum can make the minimum ratio smaller than the sound speed; the slope of a line from the origin tangent to the full dispersion determines the kinematic bound. A neutral paired Fermi superfluid also has a fermionic pair-breaking bound of order $\Delta /p_F$ in weak coupling. The smallest relevant branch bound must be considered. Actual flow can fail earlier through vortex nucleation, surface defects or other instabilities: the excitation-energy criterion is not a complete dynamical guarantee of dissipation-free flow.

Use a complex pair order parameter $\Psi =|\Psi |e^{i\theta }$ with signed charge $q$, $|q|=2e$. A local isotropic Ginzburg-Landau theory of superconductivity has free energyFor $\alpha <0$, $\beta >0$, the uniform amplitude is $|{\Psi }_0{|}^2=-\alpha /\beta$. Variation gives the static field equation and supercurrentThe vector potential obeys Maxwell's equations, including $\nabla \times \mathbf{B}={\mu }_0({\mathbf{j}}_s+{\mathbf{j}}_n)+{\mu }_0{ϵ}_0{\partial }_t\mathbf{E}$ and $\nabla \times \mathbf{E}=-{\partial }_t\mathbf{B}$, with $\mathbf{B}=\nabla \times \mathbf{A}$ and $\mathbf{E}=-\nabla \varphi -{\partial }_t\mathbf{A}$. The static approximation drops the displacement current. A common dissipative near-transition dynamical extension is the time-dependent Ginzburg-Landau equationIt is a phenomenological relaxation law in its applicable regime, not a universal microscopic Schrödinger equation for every superconductor. Its covariant time derivative and the covariant spatial derivative ensure consistent phase and electromagnetic gauge transformations.

The amplitude and screening lengths areFor $\kappa >1/\sqrt{2}$ the ideal model is type II. An Abrikosov vortex has a $2\pi n$ phase winding, a core with suppressed $|\Psi |$ of size $\xi$, and circulating screening currents over $\lambda$. Single-valuedness implies fluxoid and flux quantization in a superconductor:The sign of $n$ depends on the chosen circulation and charge convention. Far from the core, where the contour current term is negligible, each singly wound vortex carries one flux quantum in magnitude.

The equilibrium lower critical magnetic field marks favorable entry of the first vortex; for $\kappa \gg 1$, ${\mu }_0{H}_{c1}\simeq [{\Phi }_0/(4\pi {\lambda }^2)](\ln \kappa +C)$ with an order-one core correction $C$. The upper critical magnetic field follows by linearizing the order-parameter equation in the field:Between these fields lies the superconducting mixed state, often with an Abrikosov vortex lattice. The ideal equilibrium entry and exit fields need not yield large irreversible loops. Flux pinning and surface barriers impede vortex motion; field sweeps then retain history-dependent flux and produce magnetic hysteresis and dissipative losses.

A superconducting weak link, such as a tunnel junction or narrow constriction, couples the phases of two superconductors. With consistent choices of phase and voltage orientation, the Josephson effect givesThe Gauge-invariant Josephson phase includes the line integral of $\mathbf{A}$ across the link. A zero voltage can sustain a phase-dependent steady current; a constant nonzero voltage makes that contribution oscillate. Microwave irradiation can lock this oscillation to an applied frequency, giving voltage steps at $V_n=nhf/(2e)$, an experimental test of coherent pair transport rather than an ordinary resistive current law.

A DC SQUID places two weak links in a loop. For identical junction critical currents $I_0$ and negligible self-inductance, flux quantization fixes their relative phase difference to $2\pi {\Phi }_{\mathrm{ext}}/{\Phi }_0$. Adding the two sinusoidal currents and maximizing over their common phase givesThe flux-periodic response provides sensitive magnetic measurements. Junction asymmetry or appreciable loop inductance changes the ideal contrast and can introduce additional hysteresis; the phase-coherent interference principle remains. The Jaklevic Josephson interference experiment is a motivating experimental realization of this link between a superconducting loop, weak links and magnetic flux.

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.

A bosonic coherent state is an annihilation-operator eigenstate,It has $⟨N⟩=\mathrm{Var}N=|\alpha {|}^2$, a selected complex amplitude, and Poisson-distributed particle number. A coherent-field condensate variational state extends this construction to a spatial field with $⟨\hat{\psi }(\mathbf{r})⟩=\Psi (\mathbf{r})$. Evaluating a normal-ordered interacting Bose Hamiltonian in that state gives a classical-field energy functional; its variational equation is the Gross–Pitaevskii equation. This connects a macroscopic phase to interference, quantized circulation and neutral Josephson effects. A fixed-number condensate instead has vanishing one-point field expectation; off-diagonal long-range order identifies the condensate without requiring exact number fluctuations.

The BCS ground state has the pair-coherent formIt is a product of Cooper-pair pseudospin coherent states: each paired momentum sector mixes its empty and doubly occupied configurations. The anomalous pair expectation is $⟨c_{-\mathbf{k}\downarrow }{c}_{\mathbf{k}\uparrow }⟩=u_{\mathbf{k}}^{\ast }{v}_{\mathbf{k}}$, while an ordinary single-electron expectation vanishes by fermion parity. Pair operators are not unrestricted canonical boson operators, so this state must not be identified with a naive oscillator coherent state of individual electrons. Number projection gives a fixed-even-number description with pair correlations. Phase differences between condensates enter the Josephson effect; in a charged condensate that phase difference must be gauge invariant.

The electron-phonon interaction lets one electron polarize the ionic lattice and another interact with that delayed deformation. Phonon-mediated electron attraction acts over energy transfers below a characteristic phonon scale. Pairing additionally requires a net attractive interaction in some allowed pairing channel after screened Coulomb repulsion is included. In the conventional weak-coupling situation, retardation separates the low phonon scale from the much larger electronic scale, reducing the effective repulsion. The Cooper instability then makes a Fermi surface unstable to an arbitrarily weak net attraction in the idealized channel; finite temperature cuts off its logarithmic enhancement.

In BCS theory, a momentum-independent pairing amplitude within the cutoff obeys a self-consistency equation containing $\tanh (E/2k_{B}T)/(2E)$, with quasiparticle energy $E=\sqrt{{\xi }^2+|\Delta {|}^2}$. At a continuous transition, take $\Delta \to 0$. A constant effective density of states and coupling combine into $\lambda$, while the band bottom truncates the negative-energy side at $-\sigma$, with $\sigma =\min (\mu ,\theta )$. This gives the supplied linearized gap equation. Its integrand is even, so withit becomes $1/\lambda =[I(\theta ,T_c)+I(\sigma ,T_c)]/2$. If both cutoffs are much larger than $k_B{T}_c$, the given large-cutoff integral gives $I(a,T_c)\simeq \ln (1.13a/(k_B{T}_c))$. Consequently the asymmetric-cutoff BCS transition temperature isFor $\mu >\theta$, this reduces to $1.13\theta e^{-1/\lambda }$; for $0<\mu <\theta$, it is $1.13\sqrt{\theta \mu }{e}^{-1/\lambda }$. The second limit requires $k_B{T}_c\ll \mu$ as well as $k_B{T}_c\ll \theta$. Weak $\lambda$ at fixed nonzero $\mu /\theta$ supplies that limit, but taking $\mu \to 0$ afterward need not preserve it.