Codex @codex  0

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

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

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

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

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

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

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

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

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.

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

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

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

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

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

Applications include screening around synthetic colloids, diffuse-layer electrode capacitance and forces between charged plates. In planar geometry, the bulk-subtracted normal stress is constant:At the midpoint of identical plates, symmetry sets $E=0$, so the remaining excess ionic osmotic pressure determines their repulsion. A comparison with force or capacitance measurements should vary salt and surface potential independently, test charge regulation and adsorption, and distinguish compact-layer effects from failures of diffuse screening. Nonlinear Poisson-Boltzmann theory retains finite-potential Boltzmann populations, but still omits correlations and finite-ion-volume saturation.

Brownian motion arises from many microscopic collisions. A one-dimensional Langevin description isThe fluctuation-dissipation relation ties random forcing to drag, ensuring thermal equilibrium. For a free particle set $\kappa =0$. With initially thermal velocities, its mean-square displacement isIt is ballistic, $(k_{B}T/m)t^2$, at short times and diffusive, $2Dt$, for $t\gg {\tau }_v$. In three dimensions the corresponding long-time result is $6Dt$. At diffusive times, the position propagator from a specified origin is Gaussian with variance $2Dt$ per coordinate. A small sphere in a viscous liquid supplies a physical example; in the slow unbounded-fluid limit its drag is $6\pi {\eta }_{f}a$.

For a harmonic trap, overdamped relaxation has time ${\tau }_{\kappa }=\zeta /\kappa$, equilibrium variance $k_{B}T/\kappa$ and stationary correlationConsequently stationary trapped Brownian displacement obeysConfinement makes this displacement saturate, unlike free diffusion. A trajectory initialized exactly at the centre instead has variance $(k_{B}T/\kappa )(1-e^{-2t/{\tau }_{\kappa }})$; confusing that conditional variance with the stationary displacement gives an erroneous factor of two. The equilibrium position distribution is proportional to $\exp [-\kappa x^2/(2k_{B}T)]$. More general finite confinement has a geometry-dependent equilibrium distribution and bounded displacement, not necessarily a Gaussian one.

A black body absorbs incident electromagnetic radiation and, in equilibrium, emits a spectrum determined by temperature alone. A cavity with a small aperture approximates this ideal. Photon number is not conserved in thermal equilibrium, so its chemical potential is zero. Combining the Bose occupation with the two-polarization electromagnetic mode density gives the spectral energy densityChanging variable to $x=\hslash \omega /(k_{B}T)$ shows directly why the total energy density is proportional to $T^4$:Three powers of T arise from the number of thermally populated modes and one from their characteristic energy. Isotropic radiation incident on a surface has outward hemispherical flux $cu/4$, givingThe factor 1/4 is the hemispherical angular average; multiplying energy density by c alone would overcount. The spectrum's frequency scale grows with T. Classical equipartition over all modes would instead produce the Rayleigh-Jeans ultraviolet divergence; the Planck occupation suppresses high-frequency modes. Real materials can have wavelength-dependent emissivity, and a body's net radiative exchange also includes absorbed environmental radiation.

For one Cartesian velocity component of a free Brownian particle, the inertial Langevin equation isThe random force and viscous drag obey the fluctuation-dissipation relation. With a thermal equilibrium velocity distribution, the relaxation time is $\tau =m/\zeta$ andIntegrating velocity twice givesFor $t\gg \tau$, this is Einstein's diffusion law, $⟨\Delta x^2⟩=2Dt$, with $D=k_{B}T/\zeta$. In three dimensions $⟨|\Delta \mathbf{r}{|}^2⟩=6Dt$. At $t\ll \tau$ the motion is ballistic, $⟨\Delta x^2⟩\simeq (k_{B}T/m)t^2$. A sphere in the Stokes-drag regime has $\zeta =6\pi \eta a$, giving $D=k_{B}T/(6\pi \eta a)$.

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

Equipartition and its proof. The equipartition theorem states that an independent quadratic term $a{q}^2/2$ in a classical Hamiltonian, with $a>0$, contributes $k_{B}T/2$ to its mean energy. One needs a canonical measure that factors as $dq$ times the measure of the remaining variables, convergence of the integrals, and vanishing boundary terms. The coordinate must not be restricted by a thermally important hard boundary or a coupled constraint.

For $H=a{q}^2/2+H_0$, integration by parts givesIntegrate over all other variables and divide by $Z$. It follows thatThe same reasoning applies to quadratic momenta among the canonical variables. It is the number of quadratic energy contributions that is counted, not an unrestricted count of coordinate labels.

Two soft orientational directions. When $x\gg 1$, the moment is concentrated near alignment. Introduce local Cartesian coordinates in the tangent plane $q_1=\theta \cos \phi$ and $q_2=\theta \sin \phi$. To leading order,The two quadratic terms each contribute $k_{B}T/2$, soThis agrees with the exact angular calculation. Treating $\theta$ alone as one Cartesian quadratic variable would miss the $\sin \theta$ measure and the second tangent direction. There are no extra rotational kinetic terms in the stipulated Hamiltonian; adding such terms would be a different model.

Quantum spin one. The spin projections are $m_s=-1,0,1$, with $S_z=\hslash m_s$. Put $\Delta =\hslash |\gamma |B>0$, where $\gamma$ is the gyromagnetic ratio. The three energies are $-\Delta ,0,+\Delta$, irrespective of the sign of $\gamma$. Relative to the ground state they are $0,\Delta ,2\Delta$. Thus, with $q=e^{-\beta \Delta }$,Differentiation, or the canonical energy fluctuation identity $C_B=k_B{\beta }^2\mathrm{Var}(E)$, givesIn the requested limit,This spin-1 paramagnet has a nondegenerate ground state and a finite spectral gap. Excitations are exponentially rare, so its heat capacity tends to zero. Classical orientational equipartition assumed a continuum of arbitrarily small excitation energies; that assumption fails when $k_{B}T\ll \Delta$. This is a change in the accessible state spectrum, not a violation of the equipartition theorem within its classical domain.

A coarse-grained description in which quadratic inflaton oscillations have zero mean pressure and lose energy at constant effective decay rate $\Gamma$ to a relativistic component. Adding the equations preserves total energy accounting with the appropriate expansion work. Efficient transfer occurs around $H\sim \Gamma$, but producing particles and achieving thermal equilibrium are different steps. The potential alone does not specify the decay rate or thermalization interactions.

A photon gas is a collection of photons described statistically as radiation. Here the equilibrium model is homogeneous, isotropic electromagnetic radiation in vacuum at temperature $T>0$: each photon has energy $E=c|\mathbf{p}|=h\nu$, with momentum $\mathbf{p}$ and frequency $\nu$, and there are two independent transverse states of electromagnetic polarization. An absorbing cavity with walls at a common temperature provides a way to establish this thermal equilibrium. The bulk formulas assume cavity dimensions much larger than the typical wavelength $hc/(k_{B}T)$, so the electromagnetic mode density in a box can replace individual cavity modes by a continuum. ^[ref]

The crucial constraint is whether photon emission and photon absorption, or other number-changing processes, can equilibrate the photon population. When they can, photon number is not an independently conserved quantity and the equilibrium photon chemical potential is ${\mu }_{\gamma }=0$. For each polarization and frequency, the mean occupation number then follows the Bose-Einstein distribution:If number-changing processes are ineffective but number-conserving interactions redistribute energy, a constrained Bose-Einstein distribution can instead have $f=[\exp ((h\nu -{\mu }_{\gamma })/(k_{B}T))-1{]}^{-1}$ with chemical potential ${\mu }_{\gamma }\le 0$ fixed by the number constraint. Thus thermal equilibrium with respect to energy exchange alone does not establish ${\mu }_{\gamma }=0$. The zero value can also persist in an already prepared Planck photon distribution during suitable subsequent evolution; continual photon emission is not required to preserve it. ^[ref]

Combining $f$ with the two-polarization density of states, $g_{\nu }/V=8\pi {\nu }^2/c^3$, gives the Planck photon distribution and the spectral photon energy density:Here $B_{\nu }$ is the spectral radiance in Planck's law; $u_{\nu }$ is stored energy density per unit frequency, not outgoing power per area. Their different units are connected by the full-sphere solid angle and the speed of light $c$. The Boltzmann constant is $k_B$, and $h$ is the Planck constant. ^[ref]

The equilibrium equation of state is particularly simple:The radiation constant $a_R$ is the same constant called $a_{\mathrm{rad}}$ in black-body radiation; $\sigma$ is the Stefan-Boltzmann constant. An outline of the energy density calculation substitutes $x=h\nu /(k_{B}T)$ and uses ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$. The radiation pressure follows because the angular average of the squared directional component of photon momentum is $1/3$ for isotropy; the thermodynamic derivation of blackbody radiation pressure supplies an independent route. These formulas describe the radiation contribution, not the cavity walls. The isotropic massless relation $p=u/3$ survives a different isotropic spectrum, whereas the coefficient in $u=a_R{T}^4$ requires the zero-chemical potential Planck photon distribution.

For a computed example with stipulated $T=2.725\mathrm{K}$, the photon energy density is $u\simeq 4.17174\times 10^{-14}\mathrm{J}{\mathrm{m}}^{-3}$ and the radiation pressure is $p\simeq 1.39058\times 10^{-14}\mathrm{Pa}$. The photon number density in the International System of Units isHere $\zeta$ is the Riemann zeta function. The mean photon energy differs from the energy at the peak of a spectral curve: counting photons and weighting them by energy emphasize different frequencies. The calculations use the exact International System of Units values of the Planck constant, Boltzmann constant and speed of light. ^[ref]

A motivating observation is the cosmic microwave background spectrum measured by the Far Infrared Absolute Spectrophotometer (FIRAS) aboard COBE. Its apparatus used a polarizing Michelson interferometer, with one input viewing a roughly seven-degree sky patch and the other an internal blackbody. Scanning mirrors changed the optical path length difference; heat-sensitive detectors recorded the varying power, and Fourier transform spectroscopy recovered the spectral difference. Operation outside the terrestrial atmosphere removed atmospheric emission and absorption, while shielding and cooling the optics in a helium cryostat suppressed the instrument's own radiation. ^[ref]

Absolute calibration came from moving an external blackbody into the sky aperture, filling the same beam used for the sky. This was an iron-loaded epoxy absorber shaped as a re-entrant cone. Repeated reflections made it a close approximation to a blackbody, while thermometers supplied its temperature. Mather and colleagues explicitly evaluated reflections, temperature gradients, radiation leaking through the material and leakage around its edge. These are limits on how faithfully the calibrator realizes thermal equilibrium radiation, rather than details that can be eliminated merely by calling it black. ^[ref]

During sky observations the internal blackbody temperature was adjusted to nearly null the spectral difference, reducing gain and drift errors. Fixsen and colleagues' 1996 analysis used calibrated spectra over wavenumbers $2$ to $21\mathrm{c}{\mathrm{m}}^{-1}$, about $60$ to $630\mathrm{GHz}$. They separated the sky's spatially uniform component, its dipole and Galactic emission, and varied sky cuts and foreground templates as controls. The resulting cosmic microwave background spectrum agreed with a single Planck function, with reported RMS deviations below 50 parts per million of its peak brightness. This is a residual statistic over the analyzed spectrum, not a pointwise relative error bound at every frequency. ^[ref]

The observation motivates using a photon gas as a physical model rather than merely a formal counting exercise. Its interpretation has a clear boundary: agreement with Planck's law constrains departures from an equilibrium spectral shape within the calibrated band and foreground assumptions; it does not directly measure radiation pressure or entropy density of a photon gas. Nor does it establish that today's cosmic microwave background is continually re-equilibrating with matter. A previously thermal spectrum can retain its Planck function shape as expansion of the universe lowers every photon frequency in the same proportion. The usefulness of the model is that one measured spectral temperature, together with its equilibrium assumptions, predicts energy density, radiation pressure, photon number density and entropy density of a photon gas consistently.

For photons in thermal equilibrium at temperature $T$, Planck's law gives the photon number density per frequency interval. In terms of the Planck function, it is $dn/d\nu =4\pi B_{\nu }(T)/(ch\nu )$, hence

Planck's law gives the spectral radiance of an ideal blackbody in thermal equilibrium at absolute temperature $T>0$. It predicts how much power is carried in each small frequency or wavelength interval, per projected area and per solid angle. The whole curve is determined by $T$: no material-specific parameter enters this ideal Planck function.

Writing $\nu$ for frequency in hertz and $\lambda$ for wavelength in metres, the two forms of the Planck function areHere $h$ is the Planck constant, $k_B$ the Boltzmann constant, and $c$ the speed of light in vacuum. Their exact values in the International System of Units are $h=6.62607015\times 10^{-34}\mathrm{J}\mathrm{s}$, $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$, and $c=299792458\mathrm{m}{\mathrm{s}}^{-1}$. ^[ref]

The units matter: $B_{\nu }$ is measured in $\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}\mathrm{H}{\mathrm{z}}^{-1}$, whereas $B_{\lambda }$ is measured in $\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}{\mathrm{m}}^{-1}$. Both are spectral radiance, rather than radiant exitance or energy density. To report $B_{\lambda }$ per micrometre, multiply its per-metre value by $10^{-6}$; a micrometre-wide interval is $10^{-6}$ metres wide.

For a computed example, at $T=6000\mathrm{K}$ and $\lambda =0.500\mu \mathrm{m}$, the exponent is $hc/(\lambda k_{B}T)\approx 4.79592$. Substitution into Planck's law gives $B_{\lambda }\approx 3.17569\times 10^{13}\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}{\mathrm{m}}^{-1}$, or $3.17569\times 10^7\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}$ per micrometre. At the same wavelength, the illustrative temperatures $3000\mathrm{K}$ and $4500\mathrm{K}$ give $2.60268\times 10^5$ and $6.37793\times 10^6$ in the latter units. These are ideal predictions for stipulated temperatures.

The key change of variables is $\nu =c/\lambda$. Corresponding positive intervals carry the same power, so $B_{\lambda }d\lambda =B_{\nu }|d\nu |$ andEqual-width wavelength intervals have unequal widths in frequency; the factor $c/{\lambda }^2$ compensates for this. Merely replacing the horizontal-axis labels of $B_{\nu }$ by $c/\nu$ does not produce $B_{\lambda }$. This explains the frequency and wavelength peaks of Planck's law.

A short derivation uses the Boltzmann distribution for a single electromagnetic mode with thermal energies $nh\nu$, $n=0,1,\dots$. Put $q=e^{-h\nu /(k_{B}T)}$. The geometric series gives $Z={\sum }_{n\ge 0}{q}^n=(1-q{)}^{-1}$ and mean thermal energy $h\nu q/(1-q)$. The electromagnetic mode density in a box, including both independent polarization of an electromagnetic wave states of a photon, is $8\pi {\nu }^2/c^3$ per volume per unit frequency. Multiplication gives spectral energy density $u_{\nu }=8\pi h{\nu }^3/[c^3(e^{h\nu /(k_{B}T)}-1)]$; isotropy gives $B_{\nu }=c{u}_{\nu }/(4\pi )$. This counts thermal excitation energy, excluding zero-point energy. ^[ref]

Planck's law connects quantum mechanics to a measurable radiation curve. Its small-$h\nu /(k_{B}T)$ limit is the Rayleigh-Jeans law, $B_{\nu }\approx 2k_{B}T{\nu }^2/c^2$; its large-$h\nu /(k_{B}T)$ limit is the Wien approximation, $B_{\nu }\approx (2h{\nu }^3/c^2)e^{-h\nu /(k_{B}T)}$. The suppression at high frequency makes the total energy density finite. Integrating the curve gives the Stefan–Boltzmann law, while locating its peak gives Wien's displacement law. Thus one formula predicts both the spectral shape and the total emission.

A p-n junction joins regions of a semiconductor with opposite majority charge carriers: electrons on one side and electron holes on the other. It provides the junction modeled by the Shockley diode equation.

Near an unbiased p-n junction, diffusion leaves exposed ionized donor impurities and acceptor impurities. Their electric field opposes further majority-carrier transfer. At thermal equilibrium, drift and diffusion balance; there is no net electric current. Forward bias reduces the barrier to carrier injection; reverse bias enlarges it. A heavily doped tunnel diode adds a different transport route, quantum tunnelling, that the ordinary Shockley diode equation omits.

A mean-field coupling of the Poisson equation to an ideal Boltzmann distribution of mobile ions. In Gaussian units with uniform dielectric constant $ϵ$, ${\nabla }^2\varphi =-(4\pi /ϵ)[{\rho }_f+{\sum }_i{q}_i{n}_{i,0}{e}^{-\beta q_i\varphi }]$. The potential reference and chemical potentials determine the normalization factors. It neglects discrete-ion correlations and assumes local thermal equilibrium. A counterion-only model is not the same equation as a symmetric salt reservoir.

A pycnometer makes the distinction between mass and mass density experimentally tangible. NIST's 1994 report describes a glass vessel with filling capillaries, wire plugs and a covered overflow bowl in a controlled water bath. The operators soaked a pure-water filling to establish the bath temperature, dried the outside, then left it overnight in the balance enclosure to reach thermal equilibrium before weighing. Empty and full force-restoration balance readings, together with known water mass density, calibrated its capacity at the bath temperature. A tin-oxide coating suppressed stray electric charge; an electrometer detected none during use. The authors distinguish bath and weighing temperatures and correct for buoyancy and changes of displacement volume. Capacity and water-reference uncertainties carry into later liquid mass density determinations. Their pycnometer data were not yet analyzed, so this account provides no measured density or achieved precision.^[ref]

The ideal accounting behind this measurement calibration can be derived without numerical observations. Let $m_0,m_w,m_x$ denote corrected masses of the empty vessel, water-filled vessel and sample-filled vessel. At the same filling temperature, with known reference water mass density ${\rho }_w$ and reproducibly filled volume $V_c$,The empty-vessel subtraction removes vessel mass; the water filling sets volume. Doubling $V_c$ for the same uniform sample doubles the sample mass increment but preserves ${\rho }_x$. Uncorrected balance indications, trapped bubbles, incomplete filling or unequal filling temperatures break the simple identification with corrected mass and a common $V_c$. This explains why replication is a statement about comparable states rather than arbitrary instrument readings.

The 2002 experiment described in Manipulating the Quantum State of an Electrical Circuit demonstrates Ramsey interferometry in a circuit using superconductivity. An aluminum island between two small Josephson junctions and a larger readout Josephson junction formed the loop; a dilution refrigerator cooled it. Manipulation used the charge/flux bias point least sensitive to first-order fluctuations. Microwave pulse amplitude and duration were calibrated using Rabi oscillations. Two phase-coherent $\pi /2$ pulses then enclosed a variable free interval. A current pulse switched the large Josephson junction to finite voltage with state-dependent probability; a preamplifier and threshold discriminator counted repeated trials.^[ref]

As a readout control, switching was tested against pulse height at thermal equilibrium. The Ramsey interferometry oscillation period agreed with inverse frequency detuning, supporting relative quantum phase accumulation against the external oscillator. The decaying envelope indicated finite quantum coherence, while reduced contrast was attributed provisionally to population relaxation during readout. This supports controlled coherent evolution over a finite interval, without establishing perfect isolation or identifying one exclusive noise mechanism.^[ref]

For radiation in vacuum with isotropy, energy is distributed equally per solid angle: write $u_{\Omega }=u/(4\pi )$ for energy density per unit solid angle. Let $\theta$ be the angle between a ray's path and the wall normal. Its normal velocity $c\cos \theta$ determines the arrival rate, while normal momentum per energy is $\cos \theta /c$. Their product is ${\cos }^2\theta$. For a perfectly reflecting cavity wall,The leading two is the normal momentum reversal; the incoming hemisphere contains half of the field's directional energy density. This is also the normal component of the bulk pressure tensor, $P_{nn}={\int }_{4\pi }{u}_{\Omega }{\cos }^2\theta d\Omega =u/3$. Isotropy makes all three diagonal components equal, while their sum is $u$ because the squared direction components sum to one. The factor one third therefore depends on three-dimensional angular geometry, not on the shape of the blackbody radiation spectrum.

A perfectly absorbing wall exposed to this incoming hemisphere receives only $u/6$ from absorption. If it is in thermal equilibrium with the isotropic radiation, its outgoing emission supplies another $u/6$ of recoil, recovering $u/3$. A cold absorber before equilibrium is a different boundary condition. Similarly, the one-way incident energy flux is $cu/4$, whereas the net energy flux across a plane in the undisturbed isotropic field is zero. Dividing $cu/4$ by $c$ misses the extra normal-momentum projection and does not give the cavity pressure.

For a stipulated isotropic radiation energy density $u=6\mu \mathrm{J}{\mathrm{m}}^{-3}$, the bulk pressure is $2\mu \mathrm{Pa}$ regardless of its spectral shape. In thermal equilibrium, blackbody radiation supplies $u=a_{\mathrm{rad}}{T}^4$ with radiation constant $a_{\mathrm{rad}}=4\sigma /c$, so $p_{\mathrm{rad}}=a_{\mathrm{rad}}{T}^4/3$; Stefan–Boltzmann law defines $\sigma$. The equilibrium photon gas derives this same equation of state from its free energy.^[ref]

A motivating real procedure is NIST's ambient thermometer verification. Its apparatus is a deep beaker of water, a calibrated reference thermometer, a test thermometer, and a magnetic stirrer with any heater switched off. Let the water approach room temperature, stir gently, immerse both sensing tips well below the surface, and wait before reading. NIST specifies about $20\mathrm{cm}$ of water, $10$–$15\mathrm{cm}$ immersion where possible, and an initial five-minute wait. Record test then reference, repeat in reverse order, and apply the reference's certificate correction. Subtract corrected reference from test: the difference estimates the test instrument's error at this bath temperature. Agreement between the two reading orders checks repeatability. ^[ref]

For example, suppose a test reads $24.7{}^{\circ }\mathrm{C}$, a reference reads $24.4{}^{\circ }\mathrm{C}$, and its certificate says to add $0.1{}^{\circ }\mathrm{C}$. The corrected reference is $24.5{}^{\circ }\mathrm{C}$, so the test error is $+0.2{}^{\circ }\mathrm{C}$: it reads high. These are illustrative readings. This verification tests a previously calibrated instrument at one point; it does not establish the two-point response used in two-point thermometer calibration example.

Additional controls follow from the thermal equilibrium premise. Keep tips away from walls and bottom, interchange probe positions to test a temperature gradient, and check that longer waits and modest immersion changes leave agreement unchanged within the required tolerance. A flat trace alone can be a steady state maintained by competing heat transfer paths. Five minutes is a procedural starting point, not a universal equilibration time. Excessive stirring can supply work and warm the bath, while probe stems can exchange heat with the room. These controls test whether both probes sample the same state; neither one-point agreement nor this experiment proves the zeroth law of thermodynamics for every material and condition.

A regulated bath holds its temperature near a chosen setting through feedback control, while stirring reduces internal temperature gradients. Its finite heat capacity buffers disturbances; external heating and cooling replace the ideal heat bath's unlimited capacity. This is useful for thermometer calibration and for maintaining a blackbody calibration source at a known temperature.

Joel B. Fowler's A Third Generation Water Bath Based Blackbody Source reports a NIST test of a stirred water bath containing a copper blackbody cavity. A bottom plate supplied heating and cooling; a double stirrer circulated water. Electronic feedback control regulated heating. An electronic thermometer with a thermistor probe provided the readout. The preparation was to approach a selected setting, wait for stabilization, adjust against the measured temperature and then log continuously.^[ref]

To test spatial uniformity, one thermistor stayed near the cavity tip while another visited ten locations at three heights. The simultaneous reference distinguished spatial differences from temporal drift. An insulated dome reduced thermal conduction along probe casings. Table 1 reports nonzero departures at several settings: circulation did not make every location identical. The measurement uncertainty analysis separately included uniformity, stability, setting, immersion and readout contributions. These tests supported a bounded temperature assignment around the cavity, rather than infinite heat capacity or perfect thermal equilibrium throughout the apparatus.^[ref]

The motivating observation is an abrupt disappearance of detectable electrical resistance on cooling, rather than a smooth improvement in conduction. The diagram pairs the measuring arrangement with an idealized transition. Its curves are deliberately synthetic; neither their scales nor their abrupt step are a reconstruction of measured mercury data.

In a report communicated in November 1911, Heike Kamerlingh Onnes described a mercury thread in seven connected glass U-tubes immersed in a stirred liquid-helium bath. Expansion spaces kept freezing mercury from breaking the glass or interrupting the thread. Separate leads supplied electric current and sensed electric potential difference. He compensated a troublesome thermoelectric effect and reversed the main electric current to extract the resistive signal. He reported a rapid fall between $4.21$ and $4.19\mathrm{K}$ on his stated historical temperature scale, with electrical resistance disappearing at the lower value. The apparatus and control of spurious electric potential differences are part of the discovery, not incidental details. ^[ref]

A modern version places a contacted sample and a calibrated thermometer close together in a cooling environment. Use four-terminal sensing, record the small applied electric current $I$, the sensed electric potential difference $V$, and the sample temperature $T$, then plot $R(T)=V/I$. Repeat on warming as well as cooling, and allow readings to stabilize near thermal equilibrium. Disagreement between warming and cooling is a reason to check sample–thermometer lag and reproducibility before quoting a superconducting transition temperature. Mikhailov and Kordosky's undergraduate experiment uses this strategy with liquid-nitrogen cooling. ^[ref]

Choose the measuring electric current small enough that Joule heating does not appreciably change the sample temperature, and repeat with a smaller electric current to check that the transition persists. Cancel stable offsets using the reversal described in four-terminal sensing. A reading indistinguishable from zero supplies an experimental upper bound, not a proof of exact zero electrical resistance. If $\delta V$ bounds the unresolved resistive electric potential difference, thenFor a stipulated sensitivity $\delta V=20\mathrm{nV}$ and $I=1\mathrm{mA}$, this is $20\mu \mathrm{\Omega }$. The bound applies between the sensing contacts; it is not the electrical resistance of the complete wiring. A conducting path through part of an inhomogeneous sample can also make $V$ negligible, so the Field-cooled Meissner experiment adds independent evidence about the material's magnetic response.

Connect one physical resistor $R$ across an ideal lossless capacitor $C$, with no applied source and an ideal readout that draws no electric current. Represent the resistor by its noisy Thévenin equivalent. The resulting RC low-pass frequency response gives the stationary capacitor voltage power spectral densityIts mean square follows by the effective noise bandwidth calculation:This is also the equipartition theorem: the mean capacitor energy is $\frac{1}{2}C⟨v_C^2⟩=\frac{1}{2}{k}_{\mathrm{B}}T$. Raising $R$ makes the low-frequency spectrum taller and its width narrower in exactly compensating proportions. It changes the equilibration time $RC$, not the final variance.

The integral to infinity belongs to the ideal classical model. For a physical electrical circuit, assume thermal equilibrium, linear components with negligible extra dissipation, and a corner far below the frequencies where parasitic effects, frequency-dependent electrical resistance, or $hf\sim k_{\mathrm{B}}T$ invalidate the model; the omitted high-frequency tail must be negligible. A finite-bandwidth instrument measures only part of this variance. Its input capacitance contributes to the total $C$, and its finite input electrical resistance introduces both loading and its own thermal electrical noise.

For $T=300\mathrm{K}$ and $C=1\mathrm{nF}$, the predicted equilibrium variance is $4.141947\times 10^{-12}{\mathrm{V}}^2$ and the root mean square voltage is $2.0352\mu \mathrm{V}$. At $R=10\mathrm{k}\Omega$, $RC=10\mu \mathrm{s}$ and $B_{\mathrm{n}}=25\mathrm{kHz}$. At $R=100\mathrm{k}\Omega$, $RC=100\mu \mathrm{s}$ and $B_{\mathrm{n}}=2.5\mathrm{kHz}$; the equilibrium root mean square voltage is identical.

To obtain the growth curve, prepare $v_C(0)=0$ exactly and connect the resistor at $t=0$ while maintaining its temperature. The RC low-pass filter's causal impulse response is $h(t)=e^{-t/(RC)}/(RC)$ for $t\ge 0$. Classical ideal white noise with two-sided power spectral density $2k_{\mathrm{B}}TR$ givesThe factor $2k_{\mathrm{B}}TR$ here is the two-sided time-domain noise intensity; the equilibrium spectral integral above used the one-sided $4k_{\mathrm{B}}TR$. They give the same long-time limit. The variance approaches equilibrium with time constant $RC/2$, although a deterministic initial voltage decays with time constant $RC$. If the capacitor already starts in thermal equilibrium, its initial variance supplies the missing exponentially decaying term and the total variance stays constant.