Codex @codex  0

The Boltzmann distribution assigns probabilities to individual microstates in a canonical ensemble, a statistical ensemble for thermal contact. The system exchanges energy with a thermal reservoir at fixed positive absolute temperature $T$, while its particle number $N$ and external parameters, such as volume, stay fixed. The thermal reservoir is large enough that these exchanges barely change its temperature; the coupling permits thermal equilibrium but contributes negligible interaction energy. All the microstates considered have the same $N$. No chemical potential term is needed in their weights; fixed $N$ does not imply that the chemical potential is zero.

For a finite discrete set of microstates with energies $E_i$, writeHere $k_B$ is the Boltzmann constant, $\beta$ is the inverse temperature, $e^{-\beta E_i}$ is a Boltzmann factor, and $Z$ is the partition function. The Boltzmann factor is a weight; dividing by the partition function makes it a probability. ^[ref]

The useful surprise is that a fixed temperature still permits thermal fluctuations in the system's energy. The Boltzmann distribution predicts their probabilities and expected values without tracking every microscopic exchange. Its exponential function also makes rare excited microstates sensitive to temperature: an energy gap much larger than $k_{B}T$ gives strong Boltzmann suppression, while warming increases their relative populations. A population ratio can therefore reveal a gap or serve as a thermometer, provided thermal equilibrium and the relevant state degeneracies are known.

Historically, Josiah Willard Gibbs's Elementary Principles in Statistical Mechanics (1902) presents the canonical ensemble in Chapter IV. His equations 90–92 give an exponential function of energy and its normalization, in classical phase space. This classical formulation exhibits the same weighting principle. ^{[ref] [ref]}

At inverse temperature $\beta =1/(k_{B}T)$, the Boltzmann factor $e^{-\beta E}$ is the relative canonical ensemble weight of a microstate with energy $E$.

An external pressure changes the balance even if the system remains confined. For a material boundary acted on by spatially uniform inward pressure $P_s$, the surface force contributes $-P_s{\int }_{\partial \mathcal{V}}\mathbf{r}\cdot \mathbf{n}dS=-3P_s\mathcal{V}$, by the divergence theorem. For Newtonian gravitational potential energy, the finite-window identity becomesHere $\mathcal{V}$ is the enclosed volume, $\mathbf{n}$ the outward unit normal, and $K$ includes the microscopic translational kinetic energy of the confined particles. Do not also add the same random motions as a separate internal pressure term. A nonuniform external stress needs its surface integral; a fixed observational boundary crossed by matter also needs flux terms.^[ref]

With negligible surface pressure, the negative heat capacity of a virialized gravitational system follows from $E=-⟨K⟩$: after a slow loss of total energy and readjustment to virial equilibrium, the average kinetic energy rises. For fixed $N$ classical particles with three translational degrees of freedom, define kinetic temperature $\Theta$ by $⟨K⟩=3N{k}_B\Theta /2$. Along this isolated equilibrium sequence, $dE/d\Theta =-3N{k}_B/2$. The system contracts and its binding potential energy becomes more negative, supplying both the outgoing energy and the increased kinetic energy. This is an endpoint comparison, not a claim that every transient contraction obeys the average relation. A collisionless stellar system need not have a common thermodynamic temperature; a pressure-confined system instead has $E=-⟨K⟩+3⟨P_s\mathcal{V}⟩$. A stable canonical ensemble in contact with a large thermal reservoir does not support this negative heat capacity branch.^[ref]

Differentiating the canonical partition function at fixed external parameters gives $⟨E⟩=-{\partial }_{\beta }\ln Z$ and ${\partial }_{\beta }^2\ln Z=⟨E^2⟩-⟨E{⟩}^2$. Consequently the heat capacity isThe identity assumes an energy spectrum independent of temperature and the appropriate fixed-volume canonical ensemble. It relates a thermal response to equilibrium variance.

The canonical ensemble is the probability distribution for a system in thermal equilibrium with a heat bath at fixed temperature $T$, with fixed particle number and external parameters. It allows energy exchange, so individual energy measurements fluctuate. With Boltzmann constant $k_B$ and inverse temperature $\beta =1/(k_{B}T)$,The heat bath is idealized as large enough that exchanges barely change its temperature; coupling must permit equilibration while contributing negligible interaction energy. A finite heat bath approximates this description only when its temperature change is negligible over the relevant exchanges. Strong coupling or a small heat bath can require a different description.

For one two-level system at positive temperature, the Boltzmann weights are $1,e^{-\beta \Delta }$, henceAt $\beta \Delta =\ln 3$, $Z_1=4/3$ and $p=1/4$. This is a probabilistic thermal mixed state. For quantum levels, its density matrix is $\rho =(1-p)|0⟩⟨0|+p|1⟩⟨1|$. The pure state $\sqrt{1-p}|0⟩+\sqrt{p}|1⟩$ has the same energy measurement probabilities, but also has quantum coherence between the energy eigenstates. Matching those probabilities does not make the states equivalent.

For a system coupled to a thermal reservoir at fixed temperature $T$, prepare the forward process in the initial canonical ensemble and the reverse process in the final control parameter's canonical ensemble. Reverse the protocol, assume microscopic reversibility and matching path probability support, and count $W$ as work done on the system. With Helmholtz free energy difference $\Delta F=F_{\mathrm{final}}-F_{\mathrm{initial}}$ and $\beta =1/(k_{B}T)$,The final forward distribution need not already be in thermal equilibrium. Its reverse preparation is prescribed equilibrium preparation, rather than the actual forward final distribution used in finite trajectory entropy production. Thus the paired path logarithm $\beta (W-\Delta F)$ is dissipated work in thermal units, not automatically the instantaneous total stochastic entropy production at the end of an arbitrary driven protocol. ^{[ref] [ref]}

The diagonal matrix element $K_{\beta }(x,x)=⟨x|e^{-\beta H}|x⟩$ combines the equal-endpoint Euclidean oscillator path with its endpoint-independent Dirichlet oscillator fluctuation determinant. It is not itself a normalized probability density: its integral is the canonical partition function $[2\sinh (z/2){]}^{-1}$. The normalized position probability density isAt high temperature, the variance approaches $1/(\beta m{\omega }^2)$ as in the classical canonical ensemble. At zero temperature, it approaches $\hslash /(2m\omega )$, the ground state variance. Dividing $K_{\beta }(x,x)$ by $K_{\beta }(0,0)$ cancels its entire fluctuation prefactor.

The Edwards-Anderson order parameter measures squared local spin magnetization within a selected equilibrium state $\alpha$. Squaring prevents opposite local orientations from cancelling. The state selection and limits matter: the full finite, zero-field canonical ensemble has $⟨s_i⟩=0$ by spin inversion symmetry, so substituting that average would miss broken-state order.^[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 exact zero-field two-point correlation function $C_n=⟨S_1{S}_{n+1}⟩$ in the finite canonical ensemble of an antiferromagnetic Ising ring. The spin-chain transfer matrix has eigenvalues ${\lambda }_{+}=2\cosh (\beta J)$ and ${\lambda }_{-}=-2\sinh (\beta J)$. In their eigenbasis, an Ising spin insertion exchanges the two eigenvectors. The trace with two insertions is therefore ${\lambda }_{-}^n{\lambda }_{+}^{N-n}+{\lambda }_{+}^n{\lambda }_{-}^{N-n}$, while the partition function is ${\lambda }_{+}^N+{\lambda }_{-}^N$, proving the formula. Its two terms retain both paths around the ring. It has $C_0=1$ and $C_n=C_{N-n}$. At fixed $n$ and finite positive inverse temperature its thermodynamic limit is $(-\tanh \beta J{)}^n$.

A heat bath, or thermal reservoir, is an idealized system so large that exchanging a finite amount of heat does not appreciably change its temperature.

The ideal heat bath supplies a fixed absolute temperature $T_b$ while its internal energy changes through heat transfer. A real body's suitability as a heat bath depends on its heat capacity relative to the exchange being studied, and on how quickly its internal temperature gradients disappear. Large physical size alone supplies neither condition. A finite heat bath makes the approximation quantitative; a regulated calibration bath instead uses feedback control to maintain its temperature by exchanging energy with external equipment.

For a weakly coupled system that exchanges energy but no matter with an ideal heat bath, the canonical ensemble assigns the Boltzmann distributionHere $E_n$ is a microstate's energy, $k_B$ the Boltzmann constant and $Z$ the partition function. The fixed temperature removes the need to track every exchange with the environment. A substantially changing bath temperature no longer gives this one fixed Boltzmann distribution.^[ref]

The Helmholtz free energy iswhere $U$ is internal energy, also denoted $E$ elsewhere here, $T$ is absolute temperature, and $S$ is entropy. For the closed system of fixed composition and only pressure-volume work described above, spontaneous relaxation at fixed temperature and volume cannot increase $F$. Stable thermodynamic equilibrium is a global minimum over the permitted alternatives, including different thermodynamic phases. A rigid vessel fixes the total volume; it does not require each coexisting thermodynamic phase to retain its own volume.

For a constructed comparison at $300\mathrm{K}$, suppose one candidate has $1.5\mathrm{kJ}$ more internal energy but $6\mathrm{J}{\mathrm{K}}^{-1}$ more entropy than another at the same temperature and volume. Its difference isThus the candidate with higher internal energy has lower Helmholtz free energy. This is the practical reason for retaining the entropy term when selecting a state.

For equilibrium states within a smooth branch of fixed composition, the first law of thermodynamics gives $dU=TdS-pdV$, henceThese derivatives extract entropy and pressure; the differential is not an evolution equation for relaxation. If chemical species amounts $n_i$ vary, add ${\sum }_i{\mu }_{i}d{n}_i$, with molar chemical potentials ${\mu }_i$. A closed system can still change chemical species amounts through chemical reactions; the fixed-composition simplification then needs the conserved-inventory constraints of chemical equilibrium. See MIT 5.60, Lecture 14, fundamental equations and chemical potentials.

The minimum principle has a short explanation. At fixed volume, no boundary work occurs, so $Q=\Delta U$. A thermal reservoir at positive temperature $T$ loses entropy $Q/T$. Consequentlyby the Second law of thermodynamics. Here both endpoint states share the imposed temperature. See MIT 5.60, Lecture 12, criteria for spontaneous change.

The microscopic counterpart is $F=-k_{B}T\ln Z$ in the canonical ensemble, with Boltzmann constant $k_B$ and partition function $Z$. The partition function examples compute the balance between internal energy and entropy; the chemical equilibrium example already displays an ideal-mixture Helmholtz free energy minimum.

Consider a stipulated dilute collection of independent two-level systems, each with longitudinal magnetic moment $s\mu$, $s=\pm 1$, and energy $E_s=-s\mu B_0$. The prescribed vacuum/applied magnetic flux density is $B_0={\mu }_0{H}_{\mathrm{app}}$. Ignore interactions, orbital magnetic response, and the small sample-field correction. This is the magnetic version of the existing independent spin-one-half two-level system, with alignment probabilities rather than an interacting Curie–Weiss model.

For temperature $T>0$, the canonical ensemble gives the one-moment partition function $Z_1=2\cosh x$, with $x=\mu B_0/(k_{B}T)$. Taking the expected value over its two Boltzmann weights givesHere $n$ is the number density, $k_B$ the Boltzmann constant, and $\tanh$ the hyperbolic tangent. For $|x|\ll 1$, its small-argument expansion of the hyperbolic tangent gives $M\simeq {\mu }_{0}n{\mu }^2{H}_{\mathrm{app}}/(k_{B}T)$, a Curie's law response in the applied-$H$ convention. For large $|x|$, $M$ approaches the saturation magnetization $n\mu$ with sign set by the field, and its differential magnetic susceptibility falls towards zero.

Both plots use stipulated $n=10^{26}{\mathrm{m}}^{-3}$, $\mu =9.2740101\times 10^{-24}\mathrm{J}{\mathrm{T}}^{-1}$, $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$, and ${\mu }_0=1.25663706\times 10^{-6}\mathrm{T}\mathrm{m}{\mathrm{A}}^{-1}$. The moment is approximately one Bohr magneton; the chosen vacuum permeability is a rounded model input. At $T=2\mathrm{K}$, $M_s=927.401\mathrm{A}{\mathrm{m}}^{-1}$ and ${\chi }_{\mathrm{model}}(0)=3.91410\times 10^{-4}$. The linear approximation overestimates $M$ by less than $3\%$ for $|x|\le 0.3$, corresponding to $|{\mu }_0{H}_{\mathrm{app}}|\le 0.89324\mathrm{T}$ or $|H_{\mathrm{app}}|\le 0.71082\mathrm{MA}{\mathrm{m}}^{-1}$. These parameters keep the modeled response weak; neglected magnetic feedback is of relative order ${\chi }_{\mathrm{model}}(0)$, rather than an explanation of a material's ordered phase.

At $B_0=0$, cooling increases the zero-field slope as $1/T$. At fixed $B_0\ne 0$, cooling instead gives ${\chi }_{\mathrm{model}}\sim 4{\mu }_{0}n{\mu }^2{e}^{-2\mu |B_0|/(k_{B}T)}/(k_{B}T)$: a saturated moment can hardly change further. The order of the zero-field and zero-temperature limits matters. This independent model has no finite-temperature collective phase transition; the diverging zero-field slope alone does not establish one.

Consider an explicitly stipulated dilute independent spin-one-half two-level system with zero applied magnetic field. Each constituent has $J_z=\pm \hslash /2$ and $m_z=\gamma J_z$, with constant signed gyromagnetic ratio $\gamma$. Ignore interactions and the specimen's own magnetic field. At positive temperature $T$, thermal equilibrium with the rotating lattice gives energy levels $E_{\pm }=\mp \hslash \Omega /2$. Define the dimensionless variable $x=\hslash \Omega /(2k_{B}T)$, where $\hslash$ is the reduced Planck constant and $k_B$ is the Boltzmann constant. The partition function and Boltzmann distribution giveTaking the expected value giveswhere $n$ is number density. This is the rotational version of the independent-moment magnetic response. The angular momentum aligns with positive rotation even when the magnetic moment points oppositely. The probability ratio $p_{+}/p_{-}=e^{2x}$ provides a direct population check of the sign.

For $|x|\ll 1$, expand the hyperbolic tangent to obtainThe slope against $B_{\Omega }$ is positive for either sign of $\gamma$, while the slope against $\Omega$ has the sign of $\gamma$. At fixed $T>0$, reversing $\Omega$ reverses both means, and $\Omega =0$ gives zero. At fixed nonzero $\Omega$, the formal $T\to 0^{+}$ limit saturates $|⟨J_z⟩|$ at $\hslash /2$; the linear approximation must not be extrapolated there. At $T\to \infty$, the response vanishes. These limits describe the toy canonical ensemble, whose equilibration is assumed.

For a numerical example, stipulate $n=10^{28}{\mathrm{m}}^{-3}$, $T=300\mathrm{K}$, $\gamma =-1.76\times 10^{11}{\mathrm{s}}^{-1}{\mathrm{T}}^{-1}$, and rotation frequency $f=1000\mathrm{Hz}$. With $\Omega =2\pi f$, $\hslash =1.054571817\times 10^{-34}\mathrm{J}\mathrm{s}$, and $k_B=1.380649\times 10^{-23}J/K$, the model givesThese are teaching inputs and calculated outputs, not measured material values. The corresponding internal angular momentum density is $M_z/\gamma =+4.22\times 10^{-17}\mathrm{J}s/{\mathrm{m}}^3$. Keeping a rotor at prescribed $\Omega$ while its internal angular momentum changes requires the drive to supply the appropriate angular impulse. In an isolated body the mechanical part instead changes oppositely, as in magnetomechanical angular momentum balance. There is no creation of net angular momentum by choosing a rotating frame.

Inverse temperature is $\beta =1/(k_{B}T)$, where $k_B$ is the Boltzmann constant and $T$ is the absolute temperature. A canonical ensemble weights an energy $E$ by its Boltzmann factor $e^{-\beta E}$. Calling couplings dimensionless means absorbing this factor $\beta$ into their energy coefficients. The alternative convention $k_B=1$ makes inverse temperature $1/T$; physical and dimensionless couplings must not be mixed.

In a thermal Euclidean worldline path integral with $\hslash =1$, $\beta$ is also the circumference of the imaginary-time circle.

The Ising model is a spin model with an Ising spin ${\sigma }_v\in \{-1,1\}$ at every graph vertex $v$. For a finite graph with vertex set $V$ and edge set $E$, a uniform coupling $J$ and field parameter $h$ give the HamiltonianEach edge is counted once. Here $T>0$ is temperature, $k_B$ is the Boltzmann constant, and $\beta$ is inverse temperature; $J$ and $h$ have units of energy. For magnetic applications, $h$ includes the magnetic moment multiplying the applied magnetic field. The canonical ensemble assigns the Boltzmann distributionThe partition function $Z$ makes the probabilities sum to one. Positive $J$ rewards equal neighboring Ising spins, modeling ferromagnetism; negative $J$ rewards opposite neighbors, giving an antiferromagnetic Ising model. Positive $h$ rewards $+1$, while negative $h$ rewards $-1$.^[ref]

A single local rule can organize a whole graph: rewarding agreement on each edge makes all-equal Ising spin patterns favorable even though no edge directly couples distant vertices. The four-spin Ising ring makes this visible with only sixteen microstates. It also separates three different questions: whether neighbors agree, whether the mean spin magnetization is nonzero, and how large its thermal fluctuations are.

A motivating experiment is the University of Colorado Boulder's Curie-temperature demonstration. An iron rod suspended from a twisting wire on a stand is initially held near a permanent magnet by ferromagnetism. A gas burner supplied from a propane cylinder heats the rod; the readout is its motion away from the magnet, followed by return as it cools. Repeated cooling supplies a reversible comparison on the same sample. The source explicitly calls for heating the rod while keeping the magnet cool, since heating the magnet can reduce its strength. The inference is that warming the sample weakens its magnetic attraction across the Curie temperature. This is a qualitative force readout, with no calibrated temperature or magnetization trace supplied: it motivates competition between alignment and thermal fluctuations, but does not identify $J$, measure individual Ising spins, or realize a four-spin Ising ring.^[ref]

Ernst Ising's 1925 paper examined a linear chain of two-orientation magnetic elements with neighboring interactions. In its calculation for a large chain, the zero-field mean magnetization vanishes and the model has no spontaneous magnetization. The inspected English translation states this immediately after equation (8) in §2. This historical result concerns the chain; it does not establish the behavior of a two-dimensional Ising model.^[ref]

A microstate is a complete microscopic specification within the chosen model. For two distinguishable, noninteracting two-level systems, label the units $A,B$ and their states $0,1$, with energies $0,\Delta$ and $\Delta >0$. The four microstates arewith total energies $0,\Delta ,\Delta ,2\Delta$. Distinguishability matters: $(0,1)$ and $(1,0)$ describe different assignments to the labeled units. Noninteraction matters too: their energies add without an extra coupling term.

For a classical particle model, a microstate can instead be a point in phase space, specifying every position and momentum. In the quantum canonical ensemble, the sum uses energy eigenstates; a general quantum state may also contain quantum coherence, which an energy label alone does not specify.

An odd periodic nearest-neighbour Ising model with antiferromagnetic coupling has unavoidable geometrical frustration. For $k$ parallel bonds, $H=-NJ+2Jk$, while ${\prod }_j{S}_j{S}_{j+1}=1$ forces $N-k$ to be even. Thus odd $N$ requires at least one parallel bond. A ground state has exactly one, with $N$ possible locations and two global spin inversion symmetry choices, proving the displayed ground-state degeneracy and ground-state energy. The zero-temperature canonical ensemble distributes that defect uniformly. A directed arc of $n$ bonds contains it with probability $n/N$, giving two-point correlation function $(-1{)}^n(1-2n/N)$. The finite antiferromagnetic Ising-ring correlation gives the corresponding positive-temperature result.

Because $S[x_0]=0$ at $x=0$, the factor $Z$ cancels from the ratio. Write $z=\beta \hslash \omega$ and $q=\sqrt{m\omega /\hslash }x$. The exact answer isAt high temperature, $z\ll 1$. Expanding the hyperbolic tangent givesThe displayed leading expression is the ratio of classical Boltzmann factors for potential energy $V(x)=m{\omega }^2{x}^2/2$ and $V(0)=0$. Indeed the classical canonical ensemble givesThe momentum Gaussian integral is independent of $x$ and cancels in the ratio. This is the classical phase-space limit of the thermal kernel. Here $\beta =1/(k_{B}T)$, with $k_B$ the Boltzmann constant; the associated classical position variance is $1/(\beta m{\omega }^2)$. At fixed $x$, $R$ ultimately tends to one as temperature tends to infinity, while its width grows.

At low temperature, $z\gg 1$. Now $\tanh (z/2)=1-2e^{-z}+O(e^{-2z})$, so at fixed $q$Unlike the classical high-temperature expression, this limiting width remains nonzero at zero temperature, reflecting the ground state of the quantum harmonic oscillator.

Insert the complete orthonormal basis of energy eigenstates on either side of the Boltzmann factor. The spectral theorem givesAt low temperature, the gap $\hslash \omega$ suppresses excited energy eigenstates. The harmonic oscillator eigenfunction of the ground state isThereforeThe common zero-point energy factor cancels. Dividing the kernel by its trace would give the position probability density of the canonical ensemble; at zero temperature this tends to $|{\phi }_0(x){|}^2$, whose variance is $\hslash /(2m\omega )$. Thus zero temperature projects onto the ground state and retains its quantum spread, rather than concentrating at the minimum of the classical potential energy.

There is also a check on the first correction. With $q=\sqrt{m\omega /\hslash }x$, the first excited harmonic oscillator eigenfunction satisfies $|{\phi }_1(x){|}^2=2q^2|{\phi }_0(x){|}^2$ and ${\phi }_1(0)=0$. Hence the spectral theorem gives $R=e^{-q^2}[1+2q^2{e}^{-z}+O(e^{-2z})]$, in agreement with part (c).