Codex @codex  0

In the classical high-temperature limit, each independent harmonic oscillator contributes $k_B$ to the heat capacity at constant volume, where $k_B$ is the Boltzmann constant. A three-dimensional crystalline solid with $N$ chemical atoms has $3N$ vibrational normal modes in the bulk harmonic approximation. Its vibrational internal energy and heat capacity at constant volume therefore obeyHere $U_0$ is independent of temperature at fixed volume. Thus the Dulong-Petit law predicts $C_V\simeq 3N{k}_B$ for a real solid in this regime. Per mole of chemical atoms, rather than per mole of a compound's formula units,with molar gas constant $R$ and Avogadro constant $N_A$. A formula unit containing $r$ chemical atoms consequently has a classical vibrational prediction of $3rR$ per mole of formula units.

For example, in this ideal classical model, one mole of chemical atoms heated at fixed volume from $300$ to $600\mathrm{K}$ gains vibrational internal energyThe heat capacity at constant volume stays the same even though the thermal part of the internal energy doubles. These temperatures illustrate the stipulated classical model; they do not assert that every real solid is classical at $300\mathrm{K}$.

The proof outline is the equipartition theorem: one classical harmonic oscillator has a quadratic kinetic energy and a quadratic potential energy, each with mean $k_{B}T/2$. Their sum is $k_{B}T$, whose derivative with respect to temperature is $k_B$. Summing over $3N$ normal modes gives $3N{k}_B$. The classical result does not require equal angular frequencies: each harmonic normal mode contributes the same amount once it is classical. Overall rigid motions give only a finite-size correction to the bulk count. This independence from the individual angular frequencies makes the Dulong-Petit law useful as a quick estimate of a solid's heating requirement.

Experiments on freely expanding solids often measure the heat capacity at constant pressure $C_P$, whereas the harmonic prediction is for the heat capacity at constant volume $C_V$. In a smooth single thermodynamic phase of a simple compressible material, the difference between constant-pressure and constant-volume heat capacities iswhere $\alpha =V^{-1}(\partial V/\partial T{)}_P$ is the volumetric coefficient of thermal expansion and ${\kappa }_T=-V^{-1}(\partial V/\partial P{)}_T$ is the isothermal compressibility. Using measured $C_P$ as an approximation to $C_V$ requires this correction to be small compared with $C_V$. ^[ref]

High temperature here means $k_{B}T\gg \hslash \omega$ for the relevant vibrational normal modes, with reduced Planck constant $\hslash$ and mode angular frequency $\omega$. At low temperature, excitation across a quantum energy spacing is suppressed: modes with $\hslash \omega \gg k_{B}T$ freeze out thermally, and the classical equipartition theorem fails. The Einstein model of a solid makes this loss of vibrational heat capacity explicit. The Dulong-Petit law is therefore an ideal high-temperature vibrational limit, not a statement that all solids have the same heat capacity at every temperature. Other contributions and departures from the harmonic approximation can matter even at high temperature. ^[ref]

Historically, Petit and Dulong's 1819 experiments compared the specific heat capacities of several elemental solids and found an approximately constant product with their relative atomic weights. Their original empirical regularity came before the classical harmonic explanation; it suggested that counting chemical atoms reveals a simpler pattern than comparing equal masses. ^[ref]

If a classical monatomic gas re-equilibrates by collisions while losing atoms by effusion, the escaping atoms carry mean kinetic energy $2k_{B}T$. Without wall heating, conservation of energy gives $d(3N{k}_{B}T/2)=2k_{B}TdN$ and therefore $T\propto N^{1/3}$. This assumes a fixed container volume and a small aperture in the molecular-flow regime.

The Einstein model of a solid represents the vibrations of $N$ chemical atoms by $3N$ independent one-dimensional quantum harmonic oscillators, all with the same positive angular frequency $\omega$. Equivalently, each chemical atom has three independent oscillator directions. The common angular frequency is a simplifying model assumption, not a claim about every vibration of a real solid. Hold the volume, $N$ and $\omega$ fixed. Define the characteristic temperature ${\theta }_E=\hslash \omega /k_B$ and the dimensionless ratio $x={\theta }_E/T$. The vibrational internal energy isThe first term is the zero-point energy $3N\hslash \omega /2$. It is constant at fixed $\omega$, so it contributes no heat capacity at constant volume. The temperature-dependent term givesThis is the vibrational heat capacity at constant volume of the model, with no expansion or other contributions included. ^[ref]

For a worked comparison, choose the model parameter ${\theta }_E=300\mathrm{K}$. The normalized heat capacity at constant volume depends only on $T/{\theta }_E$:

At $T={\theta }_E$, one mole of chemical atoms therefore has vibrational heat capacity at constant volume $3R(0.92067)\simeq 22.965\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, using the unrounded model value. Even at twice the characteristic temperature, the model remains about $2.06\%$ below the Dulong-Petit law limit. These are calculations for the ideal Einstein solid, not material measurements.

A short derivation uses the thermal partition function of a quantum harmonic oscillator. The energy eigenvalues $E_n=\hslash \omega (n+1/2)$ and their Boltzmann factors form a geometric series, giving the one-oscillator canonical partition function $z=e^{-x/2}/(1-e^{-x})$. The mean occupation is $1/(e^x-1)$, hence the displayed $U$ after summing over $3N$ independent quantum harmonic oscillators. Taking its derivative with $dx/dT=-x/T$ gives the displayed heat capacity at constant volume. ^[ref]

The two limits explain the shape. At high temperature, a Taylor series in small $x$ givesAt low temperature, $1-e^{-x}\to 1$, soThe zero-point energy remains even when the thermal response is tiny: freeze-out means that changing temperature scarcely changes the occupations, not that the chemical atoms have zero energy. This is the useful lesson of the Einstein model of a solid: the classical counting rule reappears when many levels are accessible, while discrete excitation spacings explain its failure upon cooling.

The exponential low-temperature decay is also the model's limitation. Real three-dimensional crystalline solids have acoustic phonon branches with arbitrarily small angular frequencies in the bulk limit. The Debye model includes these collective normal modes and predicts a low-temperature vibrational contribution $C_V\propto T^3$, which the single-frequency Einstein solid cannot reproduce. This qualification concerns the vibrational contribution; additional excitations can change a material's total heat capacity. ^[ref]

Einstein introduced the quantum treatment of solid vibrations in his 1907 paper Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Its historical importance is that quantum excitation energies explained why the classical Dulong-Petit law could fail; the modern quantum harmonic oscillator formulation above includes the zero-point energy explicitly. ^[ref]

An electrical calorimetry experiment asks whether a measured electrical work input accounts for a measured warming. The University of Iowa's demonstration uses an insulated flask containing water, an immersed heater, a thermometer, electrical meters and a timer. Its procedure compares the duration, electrical readings and temperature rise. The following low-power measurement design and synthetic readout develop that idea; they are not measurements reported by the demonstration. ^[ref]

Use a covered calorimeter with weighed water, an immersed resistor, a calibrated thermometer and gentle stirring. Put an ammeter in series with the heater branch and a high-input-electrical resistance voltmeter across the immersed resistor's terminals. Record $V(t)$, $I(t)$, temperature $T(t)$ and room temperature $T_a(t)$ before, during and after switching the heater on. Measure the electric potential difference at the heater, rather than at the supply: lead and contact electrical resistance outside the receiving boundary can consume power that never reaches the water. Account for voltmeter loading if the ammeter also includes the sensing branch's electric current.^[ref]

The electrical readout gives operating electrical resistance $R(t)=V(t)/I(t)$ when $I(t)\ne 0$, and power $P(t)=V(t)I(t)$. To test Ohm's law, obtain several pairs at the same temperature; one ratio alone cannot establish a linear response. For small electrical resistance, use four-terminal sensing to exclude lead drops and reverse electric current to test stable offsets from the thermoelectric effect. Use brief low-power resistance measurements if the intended result is the unheated value: the test itself causes Joule heating.

For a covered assembly at approximately constant pressure, with negligible evaporation, negligible endpoint mechanical energy changes and no phase transition, takeHere $c_w$ is the water's specific heat capacity at constant pressure, $C_{\mathrm{cal}}$ is the participating vessel, heater and sensor heat capacity, $Q_{\mathrm{loss},\mathrm{out}}$ is signed outward heat leakage, and $W_{\mathrm{stir},\mathrm{on}}$ is stirring work entering the assembly. The storage term approximates its enthalpy change; for a rigid fixed-volume receiver use the corresponding heat capacity at constant volume. Calibrate $C_{\mathrm{eff}}$ independently if the purpose is to test electrical input against thermal storage. Using the same heating run to define $C_{\mathrm{eff}}=E_{\mathrm{el}}/\Delta T$ and then claiming it verifies the equality would be circular.

For the synthetic trace, stipulate $V=5\mathrm{V}$, $I=1\mathrm{A}$, $R=5\Omega$, $C_{\mathrm{eff}}=500J/K$, room temperature $20{}^{\circ }\mathrm{C}$, and loss coefficient $G=0.5W/K$, with negligible stirring input. The heater runs for $120\mathrm{s}$, beginning at room temperature. The calorimeter heat-loss correction model isIt predicts $E_{\mathrm{el}}=600\mathrm{J}$, $\Delta T\approx 1.13080\mathrm{K}$, storage $565.398\mathrm{J}$ and outward loss $34.602\mathrm{J}$. Neglecting leakage undercounts deposited energy by about $5.77\%$; assuming no loss predicts $1.20\mathrm{K}$. After switching off, the excess temperature undergoes exponential decay with time constant $C_{\mathrm{eff}}/G=1000\mathrm{s}$. A cooling slope can estimate $G$ only with known $C_{\mathrm{eff}}$ and a credible single-temperature model.

Controls separate causes of a temperature rise. Check the electrical meters' zeros and thermometer calibration; keep water mass, immersion, lid and stirring fixed between runs. Record a heater-off stirred baseline to estimate sensor drift and stirring input. Repeat cooling runs over the same temperature range to estimate environmental heat transfer; a baseline taken only at room temperature cannot measure warming-dependent leakage. Vary electric current while keeping $R$ and initial thermal conditions comparable: corrected deposited energy should scale as $I^2$ for equal durations. An independent calorimeter calibration or known reference heat capacity makes this a quantitative test rather than a qualitative observation.

Finite sensor response, incomplete thermal equilibrium, changing $R(T)$ or heat capacity, evaporation, and imperfect loss corrections limit the inference. A final temperature alone cannot distinguish smaller electrical input from larger heat leakage or a larger participating heat capacity. Agreement tests the stated assembly and its energy accounting; it does not show that all electrical work in arbitrary devices becomes heat.

Historically, the Royal Society's archive description of James Prescott Joule's manuscript dated 13 October 1840 records a coil of wire in a jar of water, a thermometer and a galvanometer, and summarizes his inference that warming depends on electrical resistance and the square of electric current intensity. This electrical experiment precedes Joule's paddle-wheel experiment; the input here is electrical work, whereas the paddle supplies mechanical work. The catalogue, rather than the inaccessible manuscript images, supports this account. ^[ref]

Enthalpy is the state function obtained by adding pressure times volume to internal energy:Here $U$ denotes internal energy, also written $E$ in some references. Like internal energy, enthalpy is measured in joules and is an extensive variable. Its definition packages the internal energy change and expansion work into one useful quantity for fixed-pressure calculations. It does not define a body's stored heat. ^[ref]

For a fixed amount and composition of a simple compressible system, the relation between neighboring equilibrium states isHere $T$ is absolute temperature and $S$ is entropy. Thus entropy and pressure are the natural variables of this thermodynamic potential, with partial derivatives $T=(\partial H/\partial S{)}_p$ and $V=(\partial H/\partial p{)}_S$. Fixed composition matters: changing composition adds further terms. The formula follows in one line from the equilibrium first law of thermodynamics identity:The cancellation of $pdV$ is the reason enthalpy replaces volume with pressure as a natural variable. This identity concerns equilibrium states; it does not identify $TdS$ with actual irreversible heat transfer. ^[ref]

For a closed system undergoing a quasistatic process at constant pressure, require fixed amount and composition, only pressure-volume work, no friction, and unchanged bulk kinetic energy and potential energy. Then the actual first law of thermodynamics balance givesThe symbol $Q_p$ means transferred heat for this fixed-pressure process. If the piston example in heat additionally meets these conditions, its $20\mathrm{kJ}$ of incoming heat equals $\Delta H$, while $\Delta U=15\mathrm{kJ}$ and $p\Delta V=W=5\mathrm{kJ}$. The example shows exactly what enthalpy adds to internal energy: the expansion work required at the imposed pressure.

Within a single thermodynamic phase, heat capacity at constant pressure giveswith the amount and composition fixed. Constant $C_p$ gives $\Delta H=C_p\Delta T$; dividing by mass gives specific heat capacity at constant pressure. At equilibrium melting, enthalpy can instead change through latent heat while temperature stays fixed. The sensible heat example therefore adds two single-thermodynamic phase enthalpy changes and one latent heat contribution.

The usefulness of enthalpy lies in matching the imposed constraint: a fixed-pressure experiment naturally reports $\Delta H$, while fixed-volume heating with no other work reports $\Delta U$. Neither equality is universal. Additional work or changing pressure alters the relation between supplied heat and $\Delta H$; irreversible entropy production also prevents identifying actual $\delta Q$ with $TdS$. The separate Joule-Thomson coefficient branch uses enthalpy to describe throttling under its own assumptions.

The entropy density of a photon gas is its entropy per physical volume, $s=S/V$. For homogeneous equilibrium vacuum radiation with two states of electromagnetic polarization and zero photon chemical potential,The first equality is the entropy density at zero chemical potential, obtained from the thermodynamic Euler relation $U=TS-pV+{\mu }_{\gamma }{N}_{\gamma }$ with internal energy $U=uV$. The second uses the photon gas equation of state, $u=a_R{T}^4$ and $p=u/3$. If an independently constrained photon number gives a nonzero chemical potential, the general entropy density is $(u+p-{\mu }_{\gamma }{n}_{\gamma })/T$ and $u=a_R{T}^4$ need not hold. This distinction explains why the absence of a conserved number matters for entropy as well as for the spectral shape. ^[ref]

A short first law of thermodynamics check, for neighboring equilibrium states with ${\mu }_{\gamma }=0$, isThe thermodynamic Euler relation fixes the additive normalization. Equivalently, this radiation entropy tends to zero with temperature at fixed volume; no contribution from the containing walls is included. Interpreting $TdS$ as transferred heat and $pdV$ as boundary work requires a reversible thermodynamic process. For a reversible thermodynamic process that is also an adiabatic process, there is no entropy production or entropy exchange, soAn irreversible adiabatic process can produce entropy even without heat exchange, so the word adiabatic alone is insufficient to infer $T^{3}V=\text{constant}$.

For an analytic example, expand an initially equilibrated photon gas reversibly from $V_0$ to $8V_0$, keeping only the same two photon states of electromagnetic polarization. Its temperature halves, its energy density becomes $u_0/16$, its entropy density of a photon gas becomes $s_0/8$, and its total entropy stays $S_0$. The internal energy halves rather than remaining fixed: the missing half leaves as work against radiation pressure,This is a stipulated reversible radiation model, not a report of a piston experiment. Maintaining a real cavity along this path requires controlling wall heat exchange and including the wall heat capacity in any actual apparatus balance.

For a computed scale at stipulated $T=2.725\mathrm{K}$, $s\simeq 2.04122\times 10^{-14}\mathrm{J}{\mathrm{K}}^{-1}{\mathrm{m}}^{-3}$. Combining this with the photon number density givesThis is the total equilibrium radiation entropy divided by its mean photon number, not a separately assigned entropy carried by each labeled particle. During the ideal expansion, $n_{\gamma }\propto T^3$ and $T^{3}V$ is fixed, so the mean number $N_{\gamma }=n_{\gamma }V$ also stays constant. This follows from the path; it was not imposed as the independent number constraint used to establish chemical potential.

The beautiful connection is that expansion lowers each typical photon energy and dilutes the radiation in just the proportions needed to retain its equilibrium spectral shape and total entropy. For the cosmic microwave background, homogeneous expansion of the universe has physical volume proportional to the cube of the scale factor; fixed-species cosmological entropy conservation therefore gives $T\propto a^{-1}$. This describes preservation of an already established Planck photon distribution and connects the FIRAS spectral observation in photon gas to an entropy prediction; FIRAS did not perform a direct entropy measurement. When other equilibrated species disappear and transfer energy into radiation, use the effective number of relativistic entropy degrees of freedom, $g_{\ast s}$: the conserved quantity for the combined sector is $g_{\ast s}{T}^{3}V$, under the corresponding equilibrium and no-entropy production assumptions. A changing species population, separate decoupled sectors or entropy production requires an appropriate entropy balance rather than the pure photon gas formula. ^[ref]

The equipartition theorem says that each independent quadratic term in a Hamiltonian described by classical mechanics has mean energy $k_{B}T/2$ in a canonical ensemble. Here $k_B$ is the Boltzmann constant and $T>0$ the absolute temperature. Count terms in the energy, rather than vaguely counting everything that can move: one harmonic oscillator has one position degree of freedom but TWO quadratic terms, one in momentum and one in displacement.

A precise sufficient setting is $H(y,\xi )=a{y}^2/2+H_0(\xi )$ with $a>0$, where $y$ is one of the canonical variables, $\xi$ denotes all remaining variables, and the phase space integration measure factors as $dyd\xi$. The range of $y$ is the whole real line, independently of $\xi$; there is no constraint coupling its allowed values to other variables. The partition function must be finite, and the relevant expected values must exist. The Boltzmann factor then givesThe probability density function ${\rho }_y$ is a normal distribution with variance $k_{B}T/a$. A larger coefficient narrows the fluctuations enough to cancel its greater energy cost. This is why the mean contribution is independent of mass or spring constant.^[ref]

For a proof outline, the one-variable Gaussian integral is $Z_y={\int }_{-\infty }^{\infty }{e}^{-\beta a{y}^2/2}dy=\sqrt{2\pi /(\beta a)}$, so $-{\partial }_{\beta }\ln Z_y=1/(2\beta )$. More generally, integration by parts in a canonical variable gives $⟨y{\partial }_{y}H⟩=k_{B}T$ when its boundary term $[y{e}^{-\beta H}]$ vanishes and the integrals converge. For $a{y}^2/2$ this is twice its mean energy. At a finite hard boundary the boundary term need not vanish; hard-wall correction to harmonic equipartition illustrates the resulting departure. A constraint, a nonquadratic potential energy, or a noncanonical integration measure therefore requires its own calculation. The theorem is an ensemble average, not equality of instantaneous energies, and applying it to a recorded trajectory requires equilibration and sufficient sampling.

For a classical one-dimensional harmonic oscillator, with positive mass $m$ and spring constant $k$,The origin is the equilibrium position. Its heat capacity is $k_B$ at fixed $m,k$, not $k_B/2$. Although the two mean energies agree, their measured amplitudes do not: $⟨p^2⟩=m{k}_{B}T$ and $⟨x^2⟩=k_{B}T/k$. In a harmonic optical trap, the position measurement observes the potential energy contribution; it does not measure both contributions by observing $x$ alone.^[ref]

For $N$ particles in a three-dimensional classical monatomic gas, neglect interactions and thermally accessible internal excitations. Its three translational momentum terms give internal energy $U=3N{k}_{B}T/2$, heat capacity at constant volume $C_V=3N{k}_B/2$, and heat capacity at constant pressure $C_p=C_V+N{k}_B=5N{k}_B/2$. Per mole, replace $N{k}_B$ by the molar gas constant $R$. For example, using stipulated $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, the molar heat capacity at fixed volume is $12.471\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. These results assume a dilute, nonrelativistic ideal gas with classical translational statistics, rather than a quantum degeneracy regime.

For an ideal gas of diatomic molecules, modeled as two chemical atoms with a bond, there are three translational quadratic terms, two accessible rotational kinetic energy terms for changes of bond orientation, and one bond-stretching harmonic oscillator with two quadratic terms. Rotation about the bond axis is not an additional active mode in this linear-molecule model. The three familiar limits of heat capacity at constant volume areCorresponding molar heat capacities are $3R/2$, $5R/2$ and $7R/2$. Counting the stretching harmonic oscillator only once as $k_{B}T/2$ would give the wrong last value.^[ref]

The frozen regimes require quantum mechanics. In the simple molecular rotation model, $E_j=k_B{\theta }_{r}j(j+1)$ with state degeneracy $2j+1$, ${\theta }_r={\hslash }^2/(2I{k}_B)$, and moment of inertia $I$. For the bond-stretching quantum harmonic oscillator, ${\theta }_v=\hslash \omega /k_B$ with angular frequency $\omega$ and $\hslash =h/(2\pi )$ from the Planck constant. When ${\theta }_r\ll {\theta }_v$, the three limits above correspond to $T\ll {\theta }_r$, ${\theta }_r\ll T\ll {\theta }_v$, and $T\gg {\theta }_v$, provided translational motion remains classical. They are smooth crossovers, not abrupt switches. Real molecules can have different allowed rotational populations; at high temperature, thermal dissociation, internal excitations and nonharmonic bond motion can invalidate the fixed-molecule model before its last limit is reached.^[ref]

For a quantum harmonic oscillator with fixed $\omega$, put $z={\theta }_v/T$. Its exact mean energy and heat capacity areAt $T\gg {\theta }_v$, the mean energy approaches $k_{B}T$ and $C_{\mathrm{osc}}$ approaches $k_B$. At $T\ll {\theta }_v$, thermal excitation is exponentially suppressed and $C_{\mathrm{osc}}$ tends to zero, although the zero-point energy $\hslash \omega /2$ survives. Frozen thus means little thermal response, not zero motion or zero energy.^[ref]

A motivating measurement is thermal calibration of an optical trap. A focused laser through a microscope objective confines a bead in water; imaging or a quadrant photodetector records a position time series. Gieseler and colleagues describe a polystyrene bead of radius $1.03\mu \mathrm{m}$ at $293.15\mathrm{K}$, a $1064\mathrm{nm}$ trapping laser, and photodetector acquisition at $10^5\mathrm{Hz}$. Their recorded bead continues undergoing Brownian motion; higher trapping power produces a narrower equilibrium position distribution. Infer the centered variance, not a squared distance from an arbitrary camera origin:The normal distribution position histogram and agreement with dynamical calibration methods support the local harmonic optical trap interpretation. For a stipulated example, $T=300\mathrm{K}$ and ${\sigma }_x=20\mathrm{nm}$ imply $k\simeq 1.04\times 10^{-5}N/m=0.0104pN/nm$. This calculated scale is not an experimental reading.^[ref]

Position calibration is essential: if detector output $v$ has conversion $x=s(v-v_0)$, its physical variance is $s^2\mathrm{Var}(v)$. An unknown scale $s$ cannot be determined together with $k$ from that one variance. A complementary calibration demonstrated by Tolić-Nørrelykke and colleagues combines the power spectrum of thermal motion with the bead response to known sinusoidal motion of a translation stage. The driven response supplies an additional length scale, and their method does not require bead size, liquid dynamic viscosity or wall distance as inputs. It still uses a thermal dynamical model; agreement between methods is a consistency check, not a theory-free proof of equipartition theorem.^[ref]

The readout itself can hide fluctuations. Wong and Halvorsen used a $1064\mathrm{nm}$ optical trap and adjustable-exposure CCD imaging of a bead in water. They varied exposure and trapping power, held the bead away from the nearest wall, adjusted illumination to maintain comparable images, and measured tracking noise using a stuck bead. The measured variance fell as exposure increased; corrected results agreed with their time-averaging model. For a stationary overdamped harmonic optical trap, with optical trap relaxation time $\tau$, an exposure of length $W$ measures ${\bar{x}}_W=W^{-1}{\int }_0^{W}x(t)dt$. Its ideal noise-free variance isFor $\alpha >0$, $S(\alpha )<1$: using the blurred variance in the uncorrected formula overestimates $k$ by $1/S(\alpha )$. At $W=\tau$, the calculated factor is $e/2\simeq 1.36$. Additive tracking noise instead raises the measured variance and biases the inferred $k$ downward if left uncorrected.^[ref]

Useful controls are to record over many optical trap relaxation times, compare blocks for drift, calibrate the position scale, measure fixed-bead noise, vary exposure, and check the position histogram for departures from a normal distribution. A sparse sequence of accurate instantaneous positions can estimate an equilibrium variance; a long exposure averages away motion even if many frames are collected. Flow, appreciable laser heating, escape from the harmonic region or nonconservative forces can break the assumed thermal equilibrium. Inferring $k$ by imposing equipartition theorem is a calibration, whereas testing the relation requires an independently established stiffness and position scale. The beautiful practical connection is that microscopic randomness can set a force scale: once calibrated, a small constant applied force displaces the mean by $f/k$. The same quadratic-term counting predicts macroscopic heat capacity and makes microscopic fluctuations useful measurements.

The equivalence point is the composition at which added titrant and the initial sample inventory match the specified chemical reaction's stoichiometry. It is a theoretical inventory condition; the titration endpoint is an observed stopping criterion. For a monoprotic acid consuming one hydroxide per acid unit, initial acid molar concentration $c_A$, sample volume $V_A$, and base titrant molar concentration $c_B$ giveThe proof is amount of substance conservation with one-to-one neutralization; a different stoichiometry requires its coefficients. All volumes must use compatible units. Since $1\mathrm{mL}=10^{-3}\mathrm{L}$, $0.100\mathrm{mol}{\mathrm{L}}^{-1}\times 10.00\mathrm{mL}=0.001000\mathrm{mol}$, not $1\mathrm{mol}$.

An equivalence point does not prescribe pH $7$. A weak acid titrated with a strong base leaves its conjugate base, which reacts with water; the corresponding chemical equilibrium can make the chemical solution alkaline. The neutral strong acid/strong base case additionally depends on temperature, through self-ionization of water.^[ref]

For one mole, a simple repulsive-volume model uses $p(V-b)=RT$ with $b>0$ and domain $V>b$. The Maxwell relation and first law of thermodynamics give $(\partial U/\partial V{)}_T=T(\partial p/\partial T{)}_V-p=0$, so $U=U(T)$. If its molar heat capacity at constant volume $C_V$ is temperature independent, it is constant, $U=C_{V}T+U_{\ast }$ and $C_p=C_V+R$. Reversible adiabatic processes obey $T(V-b{)}^{R/C_V}=\mathrm{constant}$ or $p(V-b{)}^{1+R/C_V}=\mathrm{constant}$. All positive-temperature isotherms are strictly decreasing and convex; this single-phase model contains no attractive liquid-gas coexistence mechanism. The parameter $b$ is an effective unavailable molar volume, not automatically the literal geometric volume of the molecules.

For a closed system with fixed composition in an environment at temperature $T_R>0$ and pressure $p_R$, define $A=U+p_{R}V-T_{R}S$, where $U$ is internal energy, $V$ is volume and $S$ is entropy. The exergy is this availability relative to the equilibrium reference state, $A-A_{\mathrm{eq}}$. When the combined system and environment are isolated except for the specified exchanges, $\Delta S_{\mathrm{tot}}=-\Delta A/T_R$, so the Second law of thermodynamics makes $A$ decrease toward its constrained minimum. Away from equilibrium, the environmental $T_R,p_R$ must not be replaced by the system's instantaneous values.

Two related properties of an extensive quantity must be distinguished. Additivity says that a quantity $X$ for independent subsystems A and B obeysThe subsystems may be in different states. Extensivity expressed through a homogeneous function says that the equilibrium fundamental relation obeyswith composition held fixed. This is degree-one homogeneity along a replication ray; it does not say that internal energy is a linear function of arbitrary independent changes of entropy, volume and amount of substance.

Additivity gives integer replication by repeated addition. Subdividing a uniform bulk sample into $b$ equivalent parts gives rational scaling $a/b$; continuity then extends it to positive real factors. This argument requires negligible boundary contributions and a consistently prepared state. Conversely, a degree-one homogeneous function need not be additive between different rays: $f(x,y)=\sqrt{x^2+y^2}$ has $f(\lambda x,\lambda y)=\lambda f(x,y)$, but $f(1,0)+f(0,1)=2\ne \sqrt{2}=f(1,1)$. Thus a scaling identity alone cannot justify summing arbitrary states.

For finite samples, the surface correction to thermodynamic free energy can matter. Under geometric enlargement with volume multiplied by $\lambda$, surface area grows as ${\lambda }^{2/3}$, so a surface energy term is not a bulk extensive quantity. Replicating separated droplets retains every original surface, whereas merging them changes the total surface area. Likewise, appreciable mutual gravitational potential energy makes interacting subsystems fail the independent-copy premise. The thermodynamic limit explains when bulk scaling becomes accurate; it does not remove these qualifications for every small or strongly interacting system.

Why must a conducting piston move slowly? A simple model couples the first law of thermodynamics to Newton's law of cooling. Assume a fixed amount of ideal gas with uniform temperature $T(t)$ and constant whole-sample heat capacity at constant volume $C_V$, a heat bath maintained at $T_b$, negligible friction, and mechanically slow motion. Let $G>0$ be the effective thermal conductance, the reciprocal of the thermal resistance, in $\mathrm{W}{\mathrm{K}}^{-1}$. ThenPositive heat transfer must offset expansion work to keep $T$ near $T_b$. A finite $G$ and nonzero piston speed therefore require a temperature difference. This uniform-temperature model neglects spatial temperature gradients, sensor response, cylinder heat capacity and any change of $G$ with piston position; it is an illustrative calculation, not a fit to the isothermal gas syringe experiment.

For a prescribed exponential volume change $V=V_i{e}^{\lambda t}$, define $\tau =C_V/G$, $s=t/\tau$ and $ϵ=nR\lambda /G$. The dimensionless temperature $\theta =T/T_b$ obeys a linear ordinary differential equation with initial value $\theta (0)=1$:The plotted expansion uses stipulated $C_V=(5/2)nR$ and stops at $V_f=2V_i$, at $s_f=\ln 2/[(5/2)ϵ]$. During expansion, the larger $ϵ$ produces the larger temperature deficit; once the piston is held fixed, the deficit undergoes exponential decay as $e^{-(s-s_f)}$. Reducing $|ϵ|$ makes the path more nearly an isothermal process, while increasing its duration. Reversing the piston with $\lambda <0$ produces a temperature excess and outward heat transfer; the small-$|ϵ|$ regime is the intended near-isothermal approximation. Mechanical slowness alone does not ensure sufficient heat transfer.

Prepare the same ideal monoprotic acid model at initial pH $5$, with appropriate accompanying ions to balance electric charge. To predict the new pH after a finite signed addition $\Delta b$, solve $Q({\mathrm{pH}}_1)-Q(5)=\Delta b$ using the chemical equilibrium expression above. For $\Delta b=0.002\mathrm{mol}{\mathrm{L}}^{-1}$ of added chemical base, the calculated final pH is $5.366798$ at $C=0.010\mathrm{mol}{\mathrm{L}}^{-1}$, but $5.034749$ at $C=0.100\mathrm{mol}{\mathrm{L}}^{-1}$. The initial local approximation gives changes of about $0.34605$ and $0.034730$, respectively; its accuracy decreases as the partner balance changes appreciably.

For the dilute example, additions beyond approximately $0.005\mathrm{mol}{\mathrm{L}}^{-1}$ in either direction exceed the initial inventory of one partner. The full chemical equilibrium calculation still produces a finite pH; an inventory-only Henderson-Hasselbalch equation calculation would attempt to use a negative partner quantity. At $+0.006\mathrm{mol}{\mathrm{L}}^{-1}$ the dilute example reaches pH $10.99564$, whereas the concentrated example reaches only $5.104698$. This demonstrates a finite operating range, not an abrupt cessation of all buffer capacity.

Free expansion occurs when a gas spreads into an initially evacuated region without driving an external load. The standard idealization uses a rigid, insulated vessel: initially a partition confines the gas to volume $V_1$, while the remaining compartment contains vacuum. Open a passage and wait for thermodynamic equilibrium; the same fixed amount of substance $n$ now occupies the larger accessible volume $V_2$. The outer vessel never changes size. This simple release separates three questions: what energy is conserved, whether temperature changes, and why entropy increases. ^[ref]

Choose the gas as the closed system, neglect heat exchange with the walls, and assume negligible changes in bulk kinetic energy and potential energy between the equilibrium endpoints. Insulation gives $Q_{\mathrm{in}}=0$; no moving outer boundary or external load receives work, so $W_{\mathrm{out}}=0$. The first law of thermodynamics therefore givesThe gas can develop internal motion during the release; this does not supply external work. The intermediate state need not have one equilibrium pressure or temperature, so there is no actual quasistatic process curve $p(V)$ to integrate. In particular, substituting the equilibrium ideal gas law into $\int pdV$ would describe a different process. ^[ref]

For an ideal gas, internal energy of an ideal gas depends only on temperature. With positive heat capacity at constant volume, $U(T_2)=U(T_1)$ implies $T_2=T_1$, and the ideal gas law gives $p_2/p_1=V_1/V_2$. Doubling the accessible volume thus halves the final pressure without changing the final temperature. This is an endpoint statement; it does not assert a spatially uniform, constant temperature during the rush. Insulation alone does not imply constant temperature: a reversible adiabatic expansion cools an ideal gas because it supplies work. ^[ref]

For a real gas, $U(T,V)$ can depend on volume as well as temperature, so unchanged internal energy need not mean unchanged temperature. For a fixed-composition simple compressible system in a single thermodynamic phase, equilibrium identities using a Maxwell relation giveHere $C_V$ is the whole sample's heat capacity at constant volume. The sign depends on the equation of state; for an ideal gas the numerator vanishes. For a finite free expansion, solve $U(T_2,V_2)=U(T_1,V_1)$, or integrate this slope along equilibrium states of equal internal energy; that mathematical curve does not describe the actual release. ^[ref]

Free expansion is useful because it tests how internal energy depends on number density while eliminating external expansion work. Its beautiful surprise is that an ideal gas can finish with the same temperature and internal energy but greater entropy: a thermometer alone cannot diagnose reversibility. See entropy change in free expansion for the calculation and temperature sensitivity of a free-expansion experiment for the experimental qualification. The continuous-flow throttling associated with the Joule-Thomson coefficient has a different balance: under its usual assumptions it preserves enthalpy, whereas the gas-only idealization here preserves internal energy. ^[ref]

Gas density is the mass density of a gas. For an ideal gas sample of amount of substance $n$, mean molar mass $\overline{M}$ and volume $V$, combine $m=n\overline{M}$ with $pV=nRT$ to obtainHere $R$ is the molar gas constant. For a homogeneous mixture of gases, $\overline{M}={\sum }_i{x}_i{M}_i$, where $x_i$ is the mole fraction of component $i$. Equal pressure and temperature make a gas density ratio equal to a mean-molar mass ratio within this model. Argon detection in a binary gas mixture gives a worked example and a temperature confounder.

In an idealized rigid-globe weighing, subtract the evacuated assembly's mass from the filled assembly's mass, then divide by the internal volume. Establish the filled sample's pressure and temperature before comparing gases. The subtraction removes the container's large contribution, making a small gas contribution accessible. A real measurement also needs corrections or bounds for residual gas, adsorption, leaks, changing vessel volume and external buoyancy. The Rayleigh-Ramsay argon density experiment shows why preparation-route comparisons add information beyond repeated weighings of one sample.

A gas syringe traps a gas in a graduated barrel with a movable plunger, allowing its volume to be varied and read. In an isothermal gas syringe experiment, the relevant total volume also includes any connected tubing and sensor cavity containing the same gas.

In a gas syringe, the graduated barrel reading $V_{\mathrm{read}}$ can omit a fixed connecting or sensor volume $V_d$ occupied by the same sealed gas. For a common settled temperature, Boyle's law then requiresThis second relation makes a graph of $1/p$ against $V_{\mathrm{read}}$ a straight line; its intercept divided by its slope estimates $V_d$ under these assumptions. Unequal temperatures between the barrel and the connecting space require a different correction, analogous to the gas thermometer dead-space correction.

For an original hypothetical example, let $V_d=1\mathrm{c}{\mathrm{m}}^3$, $V_{\mathrm{read},1}=20\mathrm{c}{\mathrm{m}}^3$, $p_1=100\mathrm{kPa}$, and $V_{\mathrm{read},2}=10\mathrm{c}{\mathrm{m}}^3$. The total volume falls from $21$ to $11\mathrm{c}{\mathrm{m}}^3$, so $p_2=100(21/11)\mathrm{kPa}\approx 190.91\mathrm{kPa}$, rather than $200\mathrm{kPa}$. The omitted volume makes the total compression ratio smaller than the barrel-reading ratio. These stipulated values are independent of the manufacturer's apparatus correction and are not measurements.

Vernier's Boyle's Law: Pressure-Volume Relationship in Gases, Experiment 6 of Chemistry with Vernier, specifies a graduated $20\mathrm{mL}$ gas syringe, a pressure sensor and LabQuest recording. Preparation sets the disconnected syringe to $10.0\mathrm{mL}$, then attaches it to the sensor valve, trapping air. At several plunger positions, the displayed pressure is allowed to stabilize before the paired volume is entered. The analysis tests whether $pV$ remains approximately constant. It adds $0.8\mathrm{mL}$ to each syringe reading for the sensor's internal volume; this correction belongs to that apparatus. The worksheet supplies a blank results table, not measured values.^[ref]

The manufacturer's Gas Pressure Sensor User Manual identifies a membrane sensor referenced to a vacuum, hence an absolute pressure readout. Its calibration section uses the actual local atmospheric pressure for an ambient reference check, rather than a weather report corrected to sea level. Checking this baseline before sealing separates an offset error from compression physics.^[ref]

Additional controls follow from the model. Keep the same sealed amount of substance, connections and ambient temperature; persistent pressure drift at held volume after thermal equilibrium suggests a leak or another uncontrolled change. Small steps and waiting allow compression heat to leave the gas; avoiding hand warming reduces another heat input. Compare settled readings when approaching a position from either direction. Plunger friction affects the applied force and position reproducibility, so record actual pressure and actual plunger position rather than predicting pressure from the load. These are proposed controls, not reported experimental outcomes.

Nearly constant $pV$ across equilibrated states supports the fixed-temperature ideal gas approximation over the tested range. Stabilized pressure alone does not demonstrate uniform or constant temperature; the isothermal gas syringe experiment adds a temperature readout to distinguish a transient from its cooled endpoint. Finite sensor accuracy and uncertain total volume limit the test, and agreement does not establish ideal behavior at every pressure. Repeated return readings help distinguish reproducibility from a one-direction trend.

A gas thermometer infers temperature from the equation of state of a confined gas. The constant-volume version measures pressure; the constant-pressure version measures volume. For an ideal gas,Here $N$ is particle number, $n$ is amount of substance, $k_B$ is the Boltzmann constant, and $R$ is the molar gas constant. The reference and test states must be internally equilibrated, and $T_0$ must be known on the absolute thermodynamic temperature scale. A pressure ratio cancels the unknown gas inventory and bulb size in the ideal constant-volume limit. This cancellation is the method's appeal: a mechanical measurement transfers an absolute temperature scale between different materials. The Maxwell speed distribution and ideal gas pressure from wall collisions explain the microscopic origin of the response. ^[ref]

For an original ideal-model calculation, choose $T_0=300\mathrm{K}$, $P_0=120\mathrm{kPa}$ and $P=160\mathrm{kPa}$. Then $T=400\mathrm{K}$. These are stipulated values, not observations. Both pressures must be absolute, referred to vacuum. If the surrounding atmospheric pressure were $100\mathrm{kPa}$, the corresponding gauge pressures would be $20$ and $60\mathrm{kPa}$; their ratio is three, which would incorrectly predict $900\mathrm{K}$. A pressure difference from the atmosphere cannot be substituted into the ideal gas law.

A motivating calibration experiment. Schooley's account of the NBS/NIST work describes a helium-filled metal bulb, a capillary, a capacitance diaphragm gauge, and a mercury pressure standard. Circulating bath liquid passed the bulb assembly and platinum resistance thermometers. The resistance probes checked the surrounding temperature; gas pressure supplied the independent thermodynamic comparison. In an early series, workers filled the cleaned bulb with purified helium at the highest test temperature, checked pressure over three hours, then measured at lower temperatures and near a reference at $0{}^{\circ }\mathrm{C}$. Six filling pressures allowed extrapolation toward zero pressure. ^[ref]

To turn this principle into a readable measurement protocol, use the following controls. These steps are a proposed teaching interpretation of the cited apparatus, rather than additional reported NBS observations:

The corresponding NBS apparatus used a protective case with an outer helium layer held close to the bulb pressure, reducing pressure deformation while aiding thermal contact. Thermal expansion was measured on bulb material. Gas in the capillary remained partly outside the hot region, requiring a measured tube temperature profile and geometry. Pure working gas also mattered: contaminants that sorb on surfaces can change the freely circulating inventory. These are physical effects to measure or bound, rather than scatter that disappears by averaging. ^[ref]

For illustration, neglect connecting-tube gas and describe a uniform bulb with a dilute virial expansion in particle number density ${\rho }_N=N/V$:Here $Z=PV/(N{k}_{B}T)$ and $B_2$ is the second virial coefficient in the particle-density convention; a molar-density convention rescales the coefficient. The second relation follows exactly from the definition of $Z$ for the same $N$, and shows the signs of the corrections. If the bulb expands by $0.1\%$ and $Z=1$ at both states, the toy pressure ratio above gives $400.4\mathrm{K}$, rather than $400\mathrm{K}$. Improved expansion and helium nonideality models really changed the temperatures inferred from the 1989 NIST data: Harvey's 2023 reanalysis found expansion updates of roughly $1$–$3\mathrm{mK}$ and nonideality updates of roughly $-0.1$ to $-0.5\mathrm{mK}$. Those are reported corrections to historical results, not the much larger illustrative correction just calculated. ^[ref]

With equal corrected bulb volumes, negligible or corrected connecting-tube gas, fixed reference and test temperatures, and a dilute classical single gas phase,Repeat the reference/test comparison for several sealed fills, reducing $N$ between comparisons but keeping it fixed within each pair. The extrapolation varies density at fixed endpoint temperatures; it does not cool one fill toward zero temperature. Real gases have different finite-density corrections but share the classical ideal gas intercept under these assumptions. The limit does not remove leaks, bulb expansion, gauge offsets, or the gas thermometer dead-space correction. At very low pressures, thermal transpiration can make remote sensor and bulb pressures differ through a nonisothermal capillary; this must be corrected too. The NBS work measured this pressure correction using a separate comparison of narrow and wider tubes. ^[ref]

At $x=1$, the toy models give $T_{\mathrm{app}}=407.92\mathrm{K}$ and $388.12\mathrm{K}$, respectively, although their imposed true test temperature is the same. Their zero-density intercept is the reusable quantity. Measuring several low-density points permits a fit and a check for curvature; asserting a straight line from one finite-density point supplies neither. A practical fit must balance nonideality against increasing relative pressure error and the capillary corrections, and keep its extrapolation uncertainty.