Codex @codex  0

The thermodynamic limit describes a sequence of larger systems, not a single enormous specimen: particle number $N$ and volume $V$ tend to infinity with number density $n=N/V$ fixed. Hold the chosen intensive and extensive thermodynamic quantities such as temperature $T$ and an applied magnetic field fixed along that sequence. In the canonical ensemble the controls are $T,N,V$; in the grand canonical ensemble one instead fixes $T$ and chemical potential, and the limiting number density is a result, possibly with a jump at phase coexistence. These are different specifications of a limiting experiment.

The useful limit is usually an intensive and extensive thermodynamic quantities normalization of the Helmholtz free energy, not a finite limit of the total extensive quantity itself:Here $Z_{N,V}$ is the canonical partition function, $k_B$ is the Boltzmann constant, $f$ is free-energy density, and $a$ is Helmholtz free energy per particle. Thus $F_{N,V}=Vf+o(V)$ when the limit exists; finite systems need not be exactly extensive quantities because of the surface correction to thermodynamic free energy.

Specify the interactions and the growing containers too. For a conventional fluid, use a fixed stable pair potential with a tempered pair potential tail, and containers whose boundary layer of any fixed thickness occupies a vanishing fraction of the volume; uniformly enlarged cubes work, a tube of fixed cross-section does not give the same bulk limit. For a finite-range Ising model, keep the lattice spacing and bond coupling fixed while enlarging all intended spatial directions. Ordinary bounded wall interactions or periodic boundary conditions can then give the same bulk free-energy density, although boundary conditions can select different spontaneous symmetry breaking states. This is a set of hypotheses, not a theorem that every interaction has a boundary-independent limit.^[ref]

The Curie–Weiss model below deliberately uses a different limiting prescription: its all-to-all coupling is divided by $N$. Without that factor its interaction energy would grow as $N^2$, rather than as an extensive quantity. The model has no spatial wall or correlation length to compare with a container; it is an exactly summable illustration of finite-size rounding, not a quantitative model of the fluid experiment below.

A useful original model treats the thermopile absorber as one approximately uniform temperature $T_h(t)$ with heat capacity $C_{\mathrm{th}}$, connected to a fixed reference temperature $T_c$ by thermal conductance $G_{\mathrm{th}}$. Both parameters are positive constants here. Let $P_{\mathrm{abs}}(t)$ be an added absorbed power above the initial background. Linearize the outgoing heat transfer about that background. With $\Delta T=T_h-T_c$ measured from an initial state having $\Delta T=0$, conservation of energy gives the first-order linear differential equationThe left side is the rate of stored energy; the two terms on the right are added power and outgoing heat transfer. The effective thermal conductance can include several paths, rather than just heat conduction along the wires. The model assumes a stable reference, negligible electrical loading, a common absorber temperature, and changes small enough for constant parameters.

For a step from zero to constant $P_0$ at $t=0$, subtract the steady solution $\Delta T_{\infty }=P_0/G_{\mathrm{th}}$. The difference obeys exponential decay, so the initial condition fixesThe time constant has units $(\mathrm{J}{\mathrm{K}}^{-1})/(\mathrm{W}{\mathrm{K}}^{-1})=\mathrm{s}$. With the same constructed $N=20$ and $\alpha =40\mu \mathrm{V}{\mathrm{K}}^{-1}$, choose $C_{\mathrm{th}}=0.02\mathrm{J}{\mathrm{K}}^{-1}$, $G_{\mathrm{th}}=0.01\mathrm{W}{\mathrm{K}}^{-1}$ and $P_0=0.03\mathrm{W}$. ThenAt one time constant, the response is $1-e^{-1}\approx 0.632$ of its final value: $V_{\mathrm{oc}}(2\mathrm{s})\approx 1.52\mathrm{mV}$. Reaching $95\%$ requires $t=-\tau \ln (0.05)\approx 5.99\mathrm{s}$. A voltmeter read too soon therefore underestimates the steady signal even when its electrical calibration is perfect.

The steady ratio $V_{\infty }/P_0=N\alpha /G_{\mathrm{th}}$ is $0.080\mathrm{V}{\mathrm{W}}^{-1}$ in this constructed model. Halving thermal conductance doubles this ratio but also doubles the time constant if heat capacity stays fixed. Reducing heat capacity alone speeds the response without changing its steady value. Increasing $N$ gives a larger signal only at fixed individual temperature differences: adding legs can also change thermal conductance, heat capacity and electrical resistance. The algebra identifies coupled design choices, rather than predicting a particular manufactured detector.

A finite electrical resistance produces Thermal electrical noise. For nearly uniform absolute temperature $T$ and effective measurement signal bandwidth $B$, its root mean square open-circuit voltage isHere $k_{\mathrm{B}}$ is the Boltzmann constant. Thermal electrical noise in the sensing element is an important contribution, and the signal bandwidth enters its voltage variance.^[ref] If identical pairs have fixed electrical resistance $r$ and fixed $\Delta T$, then $R_{\mathrm{pile}}=Nr$: the signal grows as $N$ while this noise grows as $\sqrt{N}$. The resulting signal-to-noise ratio grows only as $\sqrt{N}$ under those assumptions. Reference temperature drift, environmental electromagnetic radiation, and meter noise remain separate limits; reducing signal bandwidth trades fast information for less integrated electrical noise.

The torsion constant, or torsional spring constant, is the positive restoring torque per unit twist angle of a complete suspension in its linear elasticity range:Its SI unit is $\mathrm{N}m/rad$. The minus sign describes restoration toward the untwisted orientation. The derivative specifies a local stiffness if the response is not globally linear; using one constant over a finite range additionally requires a straight calibration curve. Return to the same zero after unloading tests recoverable elastic deformation; repeat loading in both directions to check whether the response depends on its previous loading. Staying below permanent deformation is necessary, but does not alone prove proportionality between torque and twist.

For an ideal uniform circular wire of length $\ell$, radius $r$, and shear modulus $G_s$, the filament torsional stiffness is $C=G_s{J}_p$, with polar second moment of area $J_p=\pi r^4/2$. Uniform twist gives $\theta /\ell$ as the twist per unit length, soHere $U$ is stored elastic energy. Doubling the wire's radius multiplies the ideal torsion constant by $16$; doubling its length halves it. This explains why a long thin fiber is useful for weak torque measurement. The formula assumes linear elasticity with uniform material and circular section; it does not calibrate arbitrary ribbons, clamps, or a compound suspension. Within this ideal model, the equipartition theorem predicts an equilibrium thermal mean-square angular displacement $⟨{\theta }^2⟩=k_B\mathcal{T}/\kappa$, where $k_B$ is the Boltzmann constant and $\mathcal{T}$ the uniform absolute temperature. A softer suspension magnifies both the signal and its thermal fluctuations.

Static calibration supplies known torques without assuming the unknown interaction being measured. Apply a known perpendicular force $F$ at lever arm $a$, record the change from the unloaded twist, and useSeveral positive and negative settings test proportionality; fit the slope of applied torque against twist and retain a possible intercept. A nonzero intercept must not be silently absorbed into the slope. Convert degrees to radians before reporting $\mathrm{N}m/rad$. PASCO's procedure turns the Coulomb balance on its side and adds small known masses to generate a gravitational weight at the sphere's lever arm, returning the rod to its index mark by turning the dial. Its tabulated coefficient is force per degree: multiply by the lever arm and $180/\pi$ to obtain our torsion constant. ^[ref]

The figure's ideal static calibration uses $\kappa =1.0\times 10^{-8}\mathrm{N}m/rad$: an applied torque of $2.0\times 10^{-10}\mathrm{N}\mathrm{m}$ gives $0.020\mathrm{rad}$. Both signs lie on the same straight line. The plotted points are calculated examples, not experimental observations; they cannot demonstrate a real fiber's linearity.

Dynamic calibration instead measures the oscillation period $T_0$ with known total suspended moment of inertia $I$. For negligible damping, small free oscillations about a fixed reference obey simple harmonic motion:Include the rod, mirrors, and any attachments in $I$ unless their contributions are demonstrably negligible. Keep the applied torque absent or independent of angle. If an external torque has derivative $g_{\tau }=d{\tau }_{\mathrm{ext}}/d\theta$ at equilibrium, linearization gives effective stiffness $\kappa -g_{\tau }$: timing alone then measures that combination, and stability requires $\kappa -g_{\tau }>0$. A constant applied torque merely shifts equilibrium. Non-negligible damping changes the measured period; damped torsional oscillation gives a correction under a specified damping model. Agreement between static and dynamic calibrations is a useful check on the assumed moment of inertia, reference stability, and restoring law, rather than a reason to ignore their assumptions.

Consider an original ideal calculation with no heating, decoherence, detuning or spectator modes. In the ordered basis $(|g,n⟩,|e,n-1⟩)$ for $n\ge 1$, the red motional sideband has $H_r=\hslash \eta \Omega \sqrt{n}{\sigma }_x/2$. Here ${\sigma }_x$ is the Pauli matrix exchanging the two basis states. Since ${\sigma }_x^2=I$, its unitary time evolution isThus $U_n^{†}{U}_n=I$, and excitation from $|g,n⟩$ has probability $P_r(n,t)={\sin }^2(\theta \sqrt{n}/2)$. For $n=0$ it vanishes. The blue motional sideband similarly gives $P_b(n,t)={\sin }^2(\theta \sqrt{n+1}/2)$. A red Rabi pi pulse calibrated for $n=1$ has $\theta =\pi$; it is not a Rabi pi pulse for every occupied $n$.

Choose $\eta =0.08$, $\Omega /(2\pi )=100\mathrm{kHz}$ and $\nu /(2\pi )=2\mathrm{MHz}$. The $n=1$ red Rabi pi pulse lasts $t_{\pi }=\pi /(\eta \Omega )=62.5\mu \mathrm{s}$. Let motion have the Boltzmann distribution of a quantum harmonic oscillator, with chosen mean occupation $\bar{n}=0.25$:Here $\sum p_n=1$, $\sum n{p}_n=\bar{n}$ and $p_0=0.8$. Reindex the red sum using $n=m+1$ and $p_{m+1}=r{p}_m$: $P_r(t)=r{P}_b(t)$ for equal resonant pulse parameters, even outside the short-pulse limit. At times with $P_b>0$, the ratio $R=P_r/P_b$ therefore gives $\bar{n}=R/(1-R)$. A nonthermal distribution, unequal drive strengths, detuning or state-detection bias can invalidate that thermometer. Ground-state population cannot generally be deduced from the mean alone without specifying the distribution.

The chosen pulse gives $P_r\simeq 0.181393$ and $P_b\simeq 0.906967$. These are internal excitation probabilities, not quantum gate fidelities or cooling-cycle efficiencies. An independent matrix exponentiation of each two-state block agrees with the sine formula and preserves normalization. Summing $n=0,\dots ,79$ omits probability $r^{80}<1.21\times 10^{-56}$; because every sine squared lies between zero and one, this also bounds each omitted response. The thermal relation supplies a separate dimensional check:Here $k_B$ is the Boltzmann constant. Substitution recovers $\bar{n}=0.25$; the thermal root-mean-square expansion parameter is $\eta \sqrt{2\bar{n}+1}\simeq 0.098$. This checks the populated part of the Lamb-Dicke regime, while exponentially rare large $n$ still exceed any fixed small-parameter criterion. The pulse is a leading-order model; it does not bound corrections from the full coupling or experimental noise.

A triple-point cell turns phase coexistence into a reproducible reference for thermometer calibration. The NIST water-cell demonstration uses a sealed glass vessel containing water of specified isotopic composition and an open central thermometer well. Cooling this well with solid carbon dioxide and alcohol freezes an ice mantle while leaving liquid water and its vapor. The coolant is outside the sealed sample; the sample's vapor space does not communicate with room air.^[ref]

In the actual NIST cell-comparison experiment, newly frozen mantles were compared with an aged reference mantle in a controlled maintenance bath. A warm aluminium rod melted a thin inner layer around the well, allowing the mantle to rotate freely and establishing an ice–water interface close to the sensor. A standard platinum resistance thermometer and a resistance bridge measured the electrical resistance ratio to a calibrated, temperature-controlled reference resistor. Readings at two currents permitted extrapolation to zero electrical power, correcting self-heating; measured immersion depths supplied hydrostatic pressure corrections. Comparing preparations over time tested stabilization rather than assuming freshly frozen ice was ready.^[ref]

The inference is a reproducible fixed-point response for the thermometer, not a direct reading of three chemical potentials. The reference cell anchors the comparison; a water point alone does not determine a thermometer's response at every other temperature. The ideal $F=0$ argument excludes impurity and isotope variation, finite interface effects and unintended heat transfer; real cells also require adequate immersion and stable electrical references. Repeatability is evidence that these effects are controlled, not a reason to omit their uncertainty.

Since 20 May 2019 the kelvin has been defined through a fixed numerical value of the Boltzmann constant. The thermodynamic temperature of the water triple point is consequently a measured quantity with uncertainty, rather than an exact value imposed by the kelvin definition. Its continued use as an assigned reference on the International Temperature Scale of 1990 is a separate convention.^[ref][ref]

The Van der Waals equation is a simple real gas equation of state that can also describe a liquid-gas critical point. With particle count $N$ fixed, write the volume per particle as $v=V/N$:or equivalently $(p+a{N}^2/V^2)(V-Nb)=N{k}_{B}T$. Here $p$ is pressure, $T$ is absolute temperature, and $k_B$ is the Boltzmann constant. In SI, $v$ and $b$ have units ${\mathrm{m}}^3$ per particle, $a$ has units $\mathrm{Pa}{\mathrm{m}}^6=\mathrm{J}{\mathrm{m}}^3$, and $k_B$ has units $\mathrm{J}{\mathrm{K}}^{-1}$; the particle count is dimensionless. The pole at $v=b$ excludes smaller volume values from this model.

The $b$ correction represents excluded volume; the attraction term lowers pressure at fixed temperature and particle number density $1/v$. This is a mean-field theory, replacing detailed intermolecular structure by two parameters. At $v\gg b$ with $a/v^2\ll k_{B}T/v$, it approaches the ideal gas law $pv=k_{B}T$. It supplies a useful qualitative route from a dilute gas to phase coexistence, without a claim of quantitative accuracy for every fluid.^[ref]

Do not mix this particle convention with the molar volume convention elsewhere in this article. With amount of substance $n_{\mathrm{mol}}=N/N_A$, Avogadro constant $N_A$, and molar gas constant $R=N_A{k}_B$, one has $v_m=N_{A}v$, $b_m=N_{A}b$, and $a_m=N_A^{2}a$, so $p=RT/(v_m-b_m)-a_m/v_m^2$. The molar parameters have units ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$ and $\mathrm{Pa}{\mathrm{m}}^6\mathrm{mo}{\mathrm{l}}^{-2}$, respectively.

Johannes Diderik van der Waals introduced this equation of state in 1873, connecting molecular attraction with departures from the dilute ideal gas law. Its historical significance is the unified treatment of gases and liquids, rather than a universally exact fit.^[ref]

The virial expansion describes the first departures from the ideal gas law as a power series in number density. For a homogeneous, single-species gas in three dimensions at thermal equilibrium, write $n=N/V$, with $V$ the volume, and define the dimensionless compressibility factor $Z$ byHere $p$ is absolute pressure, $T$ is absolute thermodynamic temperature, and $k_B$ is the Boltzmann constant. The second virial coefficient determines the initial slope away from $Z=1$; higher coefficients describe further corrections. The symbol $Z$ here is an equation of state ratio, not a canonical partition function. ^[ref]

The left panel sketches the measurement principle; the right panel uses deliberately chosen polynomials, $Z=1+bx+6x^2$ with $x=n{\ell }^3$ and $b=\pm 2$, where $\ell$ is a reference length. These curves are not measured data. Their initial slopes have opposite signs, but both bend because of the quadratic term. Using a finite secant from $x=0$ to $x=h$ gives $(Z(h)-1)/h=b+6h$, rather than $b$: reducing the fitting range reduces this particular truncation bias. This algebra explains why measuring one nonideal pressure is insufficient to extract a reliable second virial coefficient.

Waxman and Hastings's 1971 Burnett experiment used sample and expansion chambers in a thermostatic bath, with a valve between them and a pressure transducer linked to a piston gauge. They measured successive isothermal expansions: evacuate the expansion chamber, open the valve, record the equilibrated pressure, isolate and evacuate that chamber again, then repeat. This avoids separately measuring the absolute gas amount and chamber volumes. For carbon dioxide at $50{}^{\circ }\mathrm{C}$, their measurements below $35$ bar gave a molar second virial coefficient of approximately $-102.2\mathrm{c}{\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$. ^[ref]

An idealized reduction makes the signal explicit. Let $q=(V_A+V_B)/V_A$ be the calibrated chamber volume ratio, $p_i$ the pre-expansion pressure, and $p_f$ the post-expansion pressure. Particle conservation during one valve opening gives $n_f=n_i/q$. At the same temperature, division of the two equations of state yieldsThus the intercept determines $q$, while the dilute-limit slope determines the second virial coefficient. In the actual experiment, helium runs calibrated the cell constant and chamber deformation corrections were included. Agreement of the carbon dioxide low-pressure intercept with the helium calibration, and the absence of an anomalous low-pressure upturn, checked for gas adsorption on chamber walls. ^[ref]

For a limited dilute-gas fit, hold temperature fixed, allow thermal equilibrium after every expansion, calibrate absolute pressure and the chamber ratio, and check leaks and sample composition. Repeat the analysis with a smaller upper pressure cutoff and an added curvature term; compare the resulting second virial coefficients with their uncertainties. Too wide a range mixes in higher coefficients; too narrow a range makes the nonideal signal hard to distinguish from instrument error. The displayed chamber equation neglects deformation, dead spaces and adsorption, so it is a teaching reduction rather than the complete experimental analysis.

Units matter. With number density in ${\mathrm{m}}^{-3}$, $B_2$ has units ${\mathrm{m}}^3$ and $B_3$ has units ${\mathrm{m}}^6$. For molar density $c=n/N_{\mathrm{A}}$ in $\mathrm{mol}{\mathrm{m}}^{-3}$, the same equation of state isHere $N_{\mathrm{A}}$ is the Avogadro constant and $R=N_{\mathrm{A}}{k}_B$ is the molar gas constant; the molar coefficients have units ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$ and ${\mathrm{m}}^6\mathrm{mo}{\mathrm{l}}^{-2}$. A coefficient tabulated per mole must therefore multiply molar density, not number density.

For the classical formulas below, use a nondegenerate gas with thermal de Broglie wavelength ${\lambda }_{\mathrm{th}}$ satisfying $n{\lambda }_{\mathrm{th}}^3\ll 1$, negligible quantum corrections to relative motion, and a stable, sufficiently short-ranged pair potential. Diluteness must also control the omitted interaction terms: $|B_{2}n|\ll 1$ alone is not enough when positive and negative contributions cancel. A finite truncation is a local low-number density approximation, not an equation of state for a dense liquid or across a phase transition. ^{[ref] [ref]}

The Wien displacement constant is the constant of proportionality $b$ in ${\lambda }_{\max }T=b$. From the exact International System of Units values of the Planck constant, speed of light and Boltzmann constant,The units follow from $hc/k_B$: $(\mathrm{J}\mathrm{s})(\mathrm{m}{\mathrm{s}}^{-1})/(\mathrm{J}{\mathrm{K}}^{-1})=\mathrm{m}\mathrm{K}$. The defining constants are exact; the displayed decimal is rounded. ^[ref]