Codex @codex  0

For $\alpha =1/2$, the preceding density of states hasWhen $T\gg T_D$, the Bose-Einstein distribution has the classical approximationEach of the $2n$ modes therefore contributes $k_{B}T$ to the energy, so the Dulong-Petit law gives

When $T\ll T_D$, the upper integration limit tends to infinity. The Bose integralthen yieldsHence the heat capacity at constant volume isso $d=4$. The high-temperature constant is the classical equipartition theorem result, while the low-temperature power law makes $C_V\to 0$ as required by the third law of thermodynamics.

Interpret the given noise correlation as the three-dimensional dot-product correlation. Each Cartesian component then has $⟨{\xi }_i(t){\xi }_j(t^{\prime })⟩=2\gamma k_{B}T{\delta }_{ij}\delta (t-t^{\prime })$. For the usual Gaussian thermal-noise idealization, each coordinate is an Ornstein-Uhlenbeck process. With $\tau =\gamma /k$,For a particle initially at ${\mathbf{x}}_0=0$,Hence its long-time displacement about the trap/elastic equilibrium hasThis also follows from equipartition theorem: each of three quadratic coordinates contributes $k_{B}T/2$.

A different common experimental quantity is the stationary time-lag displacement. It obeyswith long-lag RMS $\sqrt{6k_{B}T/k}$. The factor of two distinguishes fluctuations about a mean position from differences between two decorrelated positions. The observation convention must be stated rather than silently interchanged.

The Debye model treats a solid's low-frequency lattice vibrations as acoustic phonons in an elastic continuum. There are three polarizations per wavevector, one longitudinal and two transverse, with approximately linear dispersion relations. A cutoff is chosen to preserve the total $3N$ vibrational modes of $N$ atoms. The phonon number itself is not fixed; mode occupations follow the Bose-Einstein distribution with zero chemical potential.

For one polarization of speed $v_s$, counting wavevector states in a spherical shell gives $g_s(\omega )=V{\omega }^2/(2{\pi }^2{v}_s^3)$. Define an average Debye speed by $3/v_D^3=1/v_L^3+2/v_T^3$. The total density of states can then be writtenThe Debye wavevector is $k_D=(6{\pi }^{2}N/V{)}^{1/3}$, the Debye frequency is ${\omega }_D=v_D{k}_D$, and the Debye temperature is ${\Theta }_D=\hslash {\omega }_D/k_B$.

Omitting the temperature-independent zero-point energy, the vibrational energy and resulting Debye heat capacity areAt low temperature the upper limit can be extended to infinity. The integral is $4{\pi }^4/15$, givingOnly the long-wavelength acoustic modes are appreciably populated, explaining the Debye T-cubed test. At high temperature every mode contributes $k_B$ by equipartition theorem, so $\boxed{C_V\to 3N{k}_B}$, the Dulong-Petit law.

The model succeeds especially at low temperature, where long-wavelength elasticity is reliable. It replaces real anisotropic branches by a simple continuum and cannot reproduce detailed optical phonon branches or every intermediate-temperature feature. Metals also have an electronic heat capacity approximately proportional to $T$; defects, magnetic excitations or a phase transition can add other terms. Experiments often measure heat capacity at constant pressure, which need not equal the harmonic model's $C_V$ exactly.

The photoelectric effect, black-body radiation and the Davisson–Germer experiment challenge different parts of a fully classical picture. Classical electromagnetic waves successfully describe propagation, interference and diffraction, while classical particles follow trajectories determined by forces. The difficulty is combining that picture with the observed exchange of energy between radiation and matter, the equilibrium distribution of radiation, and electron scattering. The three experiments expose the roles of discrete energy exchange and coherent matter waves.

In a photoelectric effect experiment, monochromatic light illuminates a metal cathode in vacuum and a second electrode collects emitted electrons. Varying the collecting voltage measures a photocurrent; reversing it until collection ceases determines the stopping potential. The observations distinguish two controls: at fixed frequency above threshold, increasing light intensity increases the emission rate, whereas the largest emitted kinetic energy is determined mainly by frequency. There is a material-dependent photoelectric threshold frequency, and emission above threshold does not require the long charging-up delay predicted by a simple picture of gradual classical energy accumulation at low intensity.

A classical wave carries continuously distributed energy, with intensity proportional to squared field amplitude. On a fully classical absorption picture, increasing intensity would be expected to increase the energy delivered to an electron, and sufficiently long exposure could accumulate the escape energy even below the observed threshold. It does not naturally produce the measured universal straight-line relation between stopping potential and light frequency. In 1905, Albert Einstein instead treated light-energy transfer as localized quanta. His argument distinguished the fastest escaping electrons from electrons losing energy before escape.^[ref]

In the single-photon regime, conservation of energy gives the Einstein photoelectric equationwhere $\Phi$ is the surface work function. The photon supplies energy $h\nu$; the escape cost is paid first, and the remainder sets the upper electron-energy edge. Increasing intensity at fixed frequency increases the number of available photons, so it raises current without raising this edge. Below ${\nu }_0$, one photon cannot pay the escape cost. In an actual measurement, stable surface preparation, suppression of unwanted higher-frequency light, and the photoelectric contact-potential correction are needed to test the slope $h/e$ reliably. Robert Andrews Millikan's experiment used selected mercury spectral lines, freshly exposed metal surfaces and an electrometer to make such a quantitative test.^[ref] The interpretation is for weak single-photon illumination; multiphoton photoemission at sufficiently intense fields is a different regime.

Black-body radiation tests equilibrium physics rather than the escape of individual electrons. A small aperture in an isothermal absorbing cavity approximates an ideal blackbody: radiation entering the hole undergoes repeated interactions with the walls, and the equilibrium spectrum depends on temperature rather than wall material. A spectrally resolving detector outside the aperture measures how emitted power varies with frequency or wavelength. Experimentally the spectrum has a finite peak, is strongly suppressed at short wavelengths, shifts toward shorter wavelengths when the cavity is heated, and has integrated emission proportional to $T^4$.

For a large cavity, the electromagnetic mode density in a box, including two independent transverse polarizations, is $g(\nu )=8\pi {\nu }^2/c^3$ per volume per unit frequency. The classical equipartition theorem assigns mean energy $k_{B}T$ to each mode, leading to the Rayleigh-Jeans lawThis agrees at low frequencies but diverges as frequency increases; integrating it over all frequencies gives infinite energy density at any nonzero temperature. This is the ultraviolet catastrophe, a failure of classical equipartition applied to the electromagnetic modes.

Max Planck's essential modification was that energy at a frequency is exchanged in units $h\nu$. In a modern oscillator description, the thermal excitation levels are $nh\nu$, $n=0,1,\dots$. Using the Boltzmann distribution and summing a geometric series givesMultiplying by the mode density yields Planck's law for the spectral energy density,When $h\nu \ll k_{B}T$, expansion of the denominator recovers the classical law. When $h\nu \gg k_{B}T$, the large excitation cost suppresses occupancy exponentially, removing the divergence. Its integral gives the Stefan–Boltzmann law, while its peak gives Wien's displacement law. Thus quantization repairs the spectral shape, total emission and temperature dependence together. The displayed thermal energy excludes the oscillator's temperature-independent zero-point energy.

The Davisson–Germer experiment supplies the complementary evidence that matter cannot always be treated as classical particles with definite scattering trajectories. In 1927, Clinton J. Davisson and Lester Halbert Germer sent a narrow, accelerated electron beam normally onto a single nickel crystal in vacuum. A movable collector measured the current of scattered electrons as direction and accelerating voltage were varied; a retarding potential selected electrons retaining approximately the incident energy. Sharp angular beams depended on crystal orientation and beam energy, with a prominent feature near $54\mathrm{V}$ and detector colatitude $50^{\circ }$. These are scattering directions and currents, not photographs of electron paths.^[ref]

The de Broglie wavelength predictsAt $U=54\mathrm{V}$, this gives $\lambda \simeq 0.167\mathrm{nm}$, comparable to atomic spacings. A crystalline arrangement supplies regularly spaced scattering sites, so scattered probability amplitudes can add constructively in selected directions. For a suitable family of crystal planes, Bragg's law expresses the condition as $2d\sin {\theta }_B=n\lambda$, with ${\theta }_B$ measured relative to those planes. The detector colatitude is not automatically this angle. The measured angular structure and its variation with electron momentum are described by electron diffraction with the de Broglie wavelength. A classical point-particle deflection model has no phase wavelength $h/p$ and does not supply that interference condition.

The common quantum account preserves the successes of wave propagation and particle detection while changing what is quantized. Photons exchange energy in frequency-dependent units; a thermal distribution of those excitations gives the blackbody spectrum. Electrons have a wavefunction whose amplitudes interfere on scattering from a crystal, while detection transfers localized charge. The Born rule turns squared amplitudes into event probabilities. These results invalidate a fully classical description of matter and radiation; the photoelectric energy law alone does not establish every property of a quantized electromagnetic field, but the three phenomena together provide complementary evidence for the quantum description.

The Debye model approximates a three-dimensional crystal as an elastic continuum supporting quantized normal modes, or phonons. At long wavelengths, the acoustic phonon branches have approximately linear dispersion relations ${\omega }_s(k)=v_{s}k$: one longitudinal branch and two transverse branches. Replace their complicated actual spectra by these linear forms, then introduce a cutoff chosen to preserve the total $3N$ vibrational degrees of freedom for $N$ atoms. This is an approximation to the spectrum, not a claim that a real crystal has no non-acoustic branches.

For periodic boundary conditions in volume $V$, each allowed wavevector occupies $(2\pi {)}^3/V$ in the space of wavevectors. Counting a spherical shell for each branch gives the density of states per unit angular frequencyAn effective speed of sound $v_D$ can be defined by $3/v_D^3=1/v_l^3+2/v_t^3$. In the single-speed Debye model, choose ${\omega }_D=v_D{k}_D$ so thatThus the Debye wavevector, Debye frequency and Debye temperature satisfyHere $g$ has units of seconds when angular frequency is measured in ${\mathrm{s}}^{-1}$; its integral is the dimensionless number of normal modes. The spherical cutoff replaces the real Brillouin zone by an equal-count continuum construction.

Each phonon mode is a quantum harmonic oscillator. Since phonon number is not conserved at thermal equilibrium, the chemical potential is zero and its Bose-Einstein distribution gives mean thermal energy $\hslash \omega /[e^{\hslash \omega /(k_{B}T)}-1]$. The zero-point energy is independent of temperature in the harmonic, fixed-volume approximation and does not contribute to the heat capacity at constant volume. Excluding that constant term,Differentiating the first integral at fixed volume avoids losing the temperature-dependent upper limit in the second one. The resulting Debye heat capacity isFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{r\ge 1}{e}^{-rx}$ and integrating term by term gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\zeta (4)={\pi }^4/15$. ThereforeThis $T^3$ law follows from the three-dimensional ${\omega }^2$ density of states and the availability only of low-frequency modes at low temperature. It is a central success of the Debye model.

For $T\gg {\Theta }_D$, the Bose-Einstein distribution reduces to the classical equipartition theorem: every mode contributes thermal energy $k_{B}T$. Consequentlythe Dulong-Petit law, or $3R$ per mole of atoms. The simple model omits detailed phonon dispersion relation, distinct optical phonon branches, and anharmonic effects such as thermal expansion. In a metal, the low-temperature total heat capacity additionally contains a term approximately proportional to $T$ from electrons, so the lattice $T^3$ term should not be mistaken for the whole measurement.

The classical equipartition theorem assigns mean energy $k_{B}T/2$ to each independent quadratic term in the Hamiltonian, provided the canonical distribution and its relevant boundary integrals exist. For a quadratic term $a{q}^2$, the Boltzmann distribution is proportional to $e^{-\beta a{q}^2}$, with $\beta =1/(k_{B}T)$. Integrating a total derivative givesAfter dividing by the normalization,The same argument applies to quadratic momentum terms. Other independent canonical coordinates factor out of the ratio. It is a classical statement: widely spaced quantum energy levels can freeze out a contribution at low temperature.

The confining quartic potential energy $U(x)=Q{x}^4$, with $Q>0$, is not quadratic, so the $k_{B}T/2$ rule cannot be assigned to it. This does not invalidate the canonical Boltzmann distribution or its generalized integration-by-parts identity. Indeed,Hence equipartition for a homogeneous potential givesEquivalently, $Z_x=\Gamma (1/4)(\beta Q{)}^{-1/4}/2$ and $-\partial \ln Z_x/\partial \beta =1/(4\beta )$. This is the potential contribution only; an additional standard kinetic term $p^2/(2m)$ would contribute a further $k_{B}T/2$ to energy and $k_B/2$ to heat capacity.

For an intact diatomic molecule, the classical degrees of freedom comprise three translational quadratic terms, two rotational quadratic terms, and two terms for its one bond-stretching vibration: one kinetic and one potential. Thus complete classical equipartition predictsThe heat capacity of a diatomic gas has lower-temperature regimes because rotational and vibrational levels are quantized. In a dilute molecular gas, the expected progression isThe rotational energy scale is normally much smaller than the vibrational scale, so these steps are separated. The high-temperature value treats the vibration as harmonic and excludes dissociation and electronic excitation; the lowest-temperature gas regime is meaningful only while the sample remains dilute and gaseous.^[ref]

The plotted example uses a quantum rigid rotor and harmonic vibration, with an illustrative ratio of characteristic temperatures. It is a model calculation, not measured data for a particular molecular species. For the vibrational part,which tends to zero on freezing out and to one when many vibrational levels are accessible. The classical equipartition theorem is recovered when the relevant level spacings are small compared with $k_{B}T$.

Using $Z_{\mathrm{rot}}\simeq CT=C/(k_B\beta )$ givesThe internal-energy identity from (b), followed by differentiation with respect to temperature at fixed moment of inertia, yieldsFor one mole of independent rotors, the rotational energy is $RT$ and the rotational molar heat capacity is $R\simeq 8.31\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$.

This agrees with the equipartition theorem. The given spectrum describes a linear rotor with two independent rotational degrees of freedom, each contributing $k_{B}T/2$ to the high-temperature mean energy. Rotation about the molecular axis is not an additional degree of freedom of this ideal linear-rotor model, so the result is $k_B$ rather than $3k_B/2$ per rotor. These are rotational contributions alone; translational and vibrational heat capacities have not been included. At low temperature the discrete excitation gap freezes rotational excitation out, so the rotational heat capacity tends to zero instead of remaining at its classical value.

For a diatomic ideal gas in the usual regime with three translational and two rotational quadratic degrees of freedom, but frozen vibrational excitation, the equipartition theorem givesThe Brayton cycle efficiency is thereforeThe constant-pressure expansion has $T_3/T_2=V_3/V_2=2$, while the compression has $T_2/T_1=10^{2/7}=1.93070$. The same relations give $T_4/T_1=2$. Thus, in units of the minimum temperature,The maximum and minimum temperatures are $T_3$ and $T_1$. A Carnot cycle between those extremes hasThe difference is about $25.9$ percentage points. In the Brayton cycle, heat is supplied and rejected over ranges of temperature; the whole input is not supplied at $T_3$, nor is the whole rejection made at $T_1$. This explains the lower efficiency even for mechanically ideal stages. Absolute state temperatures are not supplied, so $\gamma =7/5$ is an explicit standard diatomic-gas approximation rather than a conclusion about all possible rotational and vibrational excitation regimes.

For a unit mass of ideal gas, the first law of thermodynamics in enthalpy form is $\delta Q=c_{p}dT-vdp$. A reversible adiabatic displacement has $\delta Q=0$, and the ideal gas law gives $v=RT/p$. ThusThe atmospheric profile approaches this dry adiabat when convection efficiently mixes potential temperature and there is no condensation or other strong diabatic source. The first law by itself determines a parcel path, not every possible environmental temperature profile; the convective mixing assumption is needed. Hydrostatic balance gives the corresponding height lapse rate $dT/dz=-g/c_p$.

For rigid linear molecules with translation and rotation active but vibration frozen, the molecular heat capacities are $C_V=(5/2)R_u$ and $C_P=(7/2)R_u$, hence $R/c_p=2/7$. Earth is dominated by the diatomic gases nitrogen and oxygen. Mars is dominated by carbon dioxide, which is triatomic but linear and has the same two rotational degrees of freedom. In the rigid-molecule approximation it therefore has the same ratio. The value is not exact for carbon dioxide when its bending vibrations contribute appreciably to heat capacity, and variable heat capacities give a nonconstant dry-adiabatic exponent. Here $\gamma$ denotes $R/c_p$, not the distinct adiabatic index $c_p/c_v$.

The fluctuation-dissipation theorem relates spontaneous fluctuations in thermal equilibrium to the dissipative part of the linear response to a weak conjugate force. For an overdamped trapped coordinate, the Langevin equation isThe fluctuation-dissipation relation for a Langevin particle fixes the noise strength from the same drag coefficient $\gamma$ that dissipates motion. With a drive proportional to $e^{-i\omega t}$, the linear response function is $\chi (\omega )=1/(K-i\gamma \omega )$. The two-sided angular-frequency power spectral density then obeysIts integral $\int S_{xx}d\omega /(2\pi )=k_{B}T/K$ agrees with the equipartition theorem. Without a restoring force, the same noise-drag relation gives the diffusion-mobility relation $D=k_{B}T/\gamma$.

In thermal calibration of an optical trap, record a bead's spontaneous displacement fluctuations and compare their spectrum with its response to a weak calibrated oscillatory force. The overdamped harmonic-trap power spectrum and the dissipative response must be consistent with the same temperature and drag. The equilibrium and weak-drive conditions are essential: an actively driven bath need not obey this equilibrium relation.

In the three-dimensional Debye model, the lattice has $3N$ normal modes for $N$ atoms. Replacing their dispersion by a linear acoustic dispersion up to a cutoff gives the Debye heat capacityAt low temperature the upper limit can be replaced by infinity. Since the integral then equals $4{\pi }^4/15$,At high temperature each harmonic mode contributes $k_B$ by the equipartition theorem, giving the Dulong-Petit lawFor one mole replace $N{k}_B$ by $R$. The curve starts with zero slope, rises as $T^3$, and levels off toward the classical limit. This is the lattice contribution; a metal's electronic heat capacity is a separate low-temperature term.

The Debye wavevector follows by counting three polarizations inside the cutoff sphere:With the supplied speed of sound treated as the effective Debye speed, ${\omega }_D=v{k}_D$ andHere $k_D\simeq 1.52\times 10^{10}{\mathrm{m}}^{-1}$. Real crystals have different longitudinal and transverse sound speeds; the one-speed estimate uses the stated value, rather than inferring those unprovided branches.

Put $x=h\nu /(k_{B}T)=hc/(\lambda k_{B}T)$. In the long-wavelength limit $x\ll 1$, the Taylor expansion $e^x-1=x+O(x^2)$ in Planck's law gives the Rayleigh-Jeans law,This is the classical equipartition theorem result: each electromagnetic mode carries mean energy $k_{B}T$. Extending it to arbitrarily high frequency would make the total energy density diverge, the ultraviolet catastrophe.

In the short-wavelength limit $x\gg 1$, $e^x-1\simeq e^x$, giving the Wien approximation,Large photon energies are thermally suppressed, removing the classical high-frequency divergence.

If the requested “form” is expressed per unit wavelength, the change of variables includes its Jacobian determinant: $U_{\lambda }=U_{\nu }|d\nu /d\lambda |=U_{\nu }c/{\lambda }^2$. ThusSubstituting $\nu =c/\lambda$ without the Jacobian determinant would give an incorrect spectral density.

A synthetic fluid membrane is a thin, laterally fluid sheet whose components can rearrange in its plane while the sheet resists changes of area and curvature. Synthetic amphiphiles can assemble into two opposed layers with solvent-compatible groups facing outward. Their fluid in-plane response gives no equilibrium static shear modulus, but does not imply zero membrane bending rigidity or zero area-compression modulus. A thin membrane can therefore bend easily while remaining difficult to stretch. These notes cover generic synthetic interfaces.

TODO: biological, biofluid and biomolecular membrane applications are user-deferred.

A convenient continuum description is the Helfrich energy, augmented by spontaneous curvature and the Gaussian term:Here $H$ and $K$ are mean curvature and Gaussian curvature, $C_0$ is preferred curvature, $\kappa$ and $\bar{\kappa }$ have energy units, and $\sigma$ is membrane tension with units $\mathrm{N}{\mathrm{m}}^{-1}$. Leaflet asymmetry can favor nonzero $C_0$. At fixed topology the integrated Gaussian term is constant for a closed surface by the Gauss-Bonnet theorem. A symmetric tensionless spherical sheet has mean-curvature bending energy $8\pi \kappa$, independent of its radius in this ideal model. Edge energy can favor closure, while fixed area and enclosed volume constrain the resulting shapes.

For a nearly flat sheet, its small-slope elastic membrane energy becomesWith $h_{\mathbf{q}}=A^{-1}\int h(\mathbf{x})e^{-i\mathbf{q}\cdot \mathbf{x}}{d}^{2}x$, the equipartition theorem gives the membrane fluctuation spectrum under tension:Bending suppresses short wavelengths, whereas positive tension suppresses the longest undulations. Their crossover is $q\sim \sqrt{\sigma /\kappa }$. Measuring the spectrum can separate these two contributions, provided finite image resolution, projected versus true area, and external confinement are accounted for. The drawn modes below are synthetic calculations, not an image of a sample.

Write the membrane height above its mean plane as a sum of Fourier modes, $h(\mathbf{r})={\sum }_{\mathbf{q}}{h}_{\mathbf{q}}{e}^{i\mathbf{q}\cdot \mathbf{r}}$. The supplied average is a thermal mean-squared modal height; for complex Fourier coefficients it is understood as $⟨|h_{\mathbf{q}}{|}^2⟩$, with $h_{-\mathbf{q}}=h_{\mathbf{q}}^{\ast }$ for a real height field. The terms are:

In the small-slope approximation, Helfrich membrane energy containsThe equipartition theorem, with the independent real modes counted consistently, gives the stated membrane fluctuation spectrum under tension. Its units are length squared: $\sigma q^2$ and $\kappa q^4$ both have units $\mathrm{J}{\mathrm{m}}^{-4}$, so $k_{B}T/[A(\sigma q^2+\kappa q^4)]$ has units ${\mathrm{m}}^2$.

At small length scales, $q$ is large and the bending term dominates:The fluctuation amplitude is controlled mainly by the membrane bending modulus and thermal energy, rather than tension. The crossover occurs when $\sigma q_c^2=\kappa q_c^4$, henceThe bending-dominated regime extends up to the membrane elastic length when scale means $1/q$, or up to ${\lambda }_c$ when it means wavelength. At still larger scales tension gives a $q^{-2}$ spectrum. At the opposite end, the continuum formula requires wavelengths larger than the microscopic membrane thickness or molecular cutoff. Without numerical $\kappa$ and $\sigma$, no numerical crossover scale is determined.

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.

Structure through pair correlations. For a homogeneous isotropic monatomic liquid of number density $n=N/V$, the radial distribution function $g(r)$ is defined byIt compares the conditional density around a tagged particle with the bulk density. A repulsive core makes close neighbours unlikely; a first peak identifies a preferred coordination shell, and damped peaks at larger radii describe short-range order. Far away, $g(r)\to 1$ in the bulk limit. Unlike a crystal, an ordinary liquid has no indefinitely repeated positional order. An ideal uncorrelated gas has $g(r)=1$ away from the tagged particle. The shell factor $4\pi r^2$ must not be confused with $g(r)$ itself.

The same pair structure enters the static structure factor measured through scattering. In the bulk isotropic convention,Thus real-space coordination and reciprocal-space scattering peaks describe the same correlations.

Energy from the structure. If interactions are sums of a central pair potential $v(r)$, the equipartition theorem gives the classical kinetic energy $3N{k}_{B}T/2$. Count the interaction of a tagged particle with its neighbours and divide by two to avoid counting each pair twice. The pair-correlation expression for internal energy isThis makes structure thermodynamically useful: short-distance repulsion and the attractive region contribute with different weights. Additional internal molecular modes or many-body interactions require extra terms. The formula is a classical bulk pair-potential result, not a universal energy formula for every liquid.

Pressure and the dilute expansion. The mechanical virial theorem relates the pair forces to pressure. The virial pressure from pair correlations isAt low density, $g(r)\simeq e^{-\beta v(r)}$, giving the virial expansionThe second virial coefficient has units of volume per particle; higher coefficients encode additional cluster correlations. A hard-sphere exclusion diameter $d$ gives $B_2=2\pi d^3/3>0$. Attractions contribute negative terms and can make $B_2$ negative. Although the pair-correlation formulas can describe a dense liquid when its actual $g(r)$ is known, truncating the virial series after a few terms is normally a dilute-fluid approximation. It should not be silently extrapolated through condensation or into a dense liquid.

Canonical probabilities. A canonical ensemble describes a system of fixed particle number and fixed external parameters such as volume, weakly coupled to a much larger thermal reservoir. They exchange energy and reach thermal equilibrium, while the reservoir's temperature stays effectively constant. The interaction energy is negligible in the state counting. The combined isolated system is described by equal probabilities for its accessible microstates.

If its total energy is $E_{\mathrm{tot}}$, the probability of a particular system microstate $i$ is proportional to the reservoir multiplicity at $E_{\mathrm{tot}}-E_i$. Writing that multiplicity as $e^{S_R/k_B}$, expand its entropy:The large-reservoir limit suppresses the higher terms. Normalization therefore givesThe sum counts individual microstates; when grouping them by energy, their state degeneracies must be included. This $Z$ is the canonical partition function.

Classical orientational heat capacity. Hold the magnetic field fixed. Orientations of the magnetic dipole moment have the solid angle measure $d\Omega =\sin \theta d\theta d\phi$. With $x=\beta mB$, the angular partition function isChoosing a normalized angular measure would remove $4\pi$, with no effect on the energy derivatives. Differentiating at fixed $B$,The angular average $⟨\cos \theta ⟩=\coth x-1/x$ is the Langevin function. Since $dx/dT=-x/T$,This is the heat capacity of a classical magnetic dipole, in $\mathrm{J}{\mathrm{K}}^{-1}$ for one magnetic moment. The word specific here refers to that single moment; $N$ independent magnetic moments have $N$ times this heat capacity. At $B=0$ the energy is independent of orientation and $C_B=0$, obtained by the continuous $x\to 0$ limit. At fixed nonzero $B$, the classical low-temperature limit instead gives $C_B\to k_B$.