Codex @codex  0

A small increase in temperature can produce a large increase in a reaction rate. In the simplest Arrhenius equation model, warming from $300$ to $310\mathrm{K}$ increases a rate constant with $E_a=80\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$ by about $2.81$ times, even though the absolute thermodynamic temperature rises by only about $3.33\%$. This sensitivity helps explain why cooling can slow unwanted chemical reactions and why a hot process can become difficult to control. The comparison assumes the same rate law and reactant concentrations; the activation energy characterizes $k(T)$, rather than every change in the measured reaction rate.

The physical motivation is the distribution of molecular energies. Warming changes the Boltzmann distribution, increasing the relative population of energetically costly configurations. A transition state can therefore become much more accessible than the modest change in average kinetic energy might suggest. The exponential function in an Arrhenius equation captures such sensitivity in a simple model; it is not a universal formula for the fraction above every barrier in a Maxwell-Boltzmann distribution.^[ref]

Activation energy quantifies how strongly a rate constant depends on thermodynamic temperature. For a positive rate constant $k(T)$, the empirical definition iswhere $R$ is the molar gas constant and $k_{\mathrm{ref}}>0$ is a fixed reference with the same units as $k$. The ratio is a dimensionless quantity, keeping the natural logarithm meaningful without changing its derivative. Thus $E_a$ is a molar energy, commonly in $\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}$ or $\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$. The Arrhenius equation makes $E_a$ constant over its assumed range. This rate-based definition, and the need for care when identifying it with a threshold on a potential energy surface, are stated in the activation energy.

In a simple barrier picture, reaching a transition state costs energy relative to the starting reactants. The height of that barrier is different from the difference between the reactants and reaction products. An energetically favourable endpoint can still be reached slowly. An energy profile makes the distinction visible, while an Arrhenius plot tests temperature sensitivity from rate constants. A measured $E_a$ need not equal a barrier read from an arbitrarily chosen diagram.

For a computed example, assume a first-order reaction has $k(300\mathrm{K})=0.100{\mathrm{s}}^{-1}$, constant $E_a=80\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$, and a temperature-independent pre-exponential factor. The Arrhenius equation givesHence $k(310\mathrm{K})\approx 0.281398{\mathrm{s}}^{-1}$. Its half-life falls from $\ln 2/0.100\approx 6.93147\mathrm{s}$ to $\ln 2/0.281398\approx 2.46322\mathrm{s}$. These are chosen model parameters, not measurements of an identified chemical reaction. The exponent uses kelvin, not unconverted Celsius temperatures, and $80000\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}$ must match $R$ in $\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$.

A microscopic reaction energy barrier belongs to a specified path on a potential energy surface. An experimental activation energy belongs to a specified rate constant and its temperature dependence. A multi-step reaction mechanism can change which step controls the observed conversion, and a temperature-dependent pre-exponential factor contributes to the observed activation energy. Even the Gibbs energy of activation of transition state theory is a different quantity: it incorporates entropy as well as enthalpy and depends on the standard states. These distinctions prevent a fitted Arrhenius plot from being mistaken for a direct microscopic energy measurement.^[ref]

Catalysis illustrates the separation between a path and its endpoints. At the same temperature and with the same reactants, reaction products, and standard states, a regenerated catalyst changes the available reaction mechanism and rate constants, while the overall standard change in Gibbs free energy and the equilibrium constant remain the same:For a reversible first-order reaction pair $\mathrm{A}\rightleftharpoons \mathrm{B}$, consistency at chemical equilibrium requires $k_f/k_r=K$. A model that multiplies only $k_f$ by a catalytic speedup while holding $k_r$ fixed would change this ratio and misrepresent catalysis. This simple ratio applies to that specified rate law; a multi-step catalytic reaction mechanism needs its own rate equations. A potential-energy profile alone cannot determine $K$, because its endpoint difference is not generally the standard change in Gibbs free energy.^[ref]

In ammonia, the nitrogen chemical atom can lie on either side of the plane of the three hydrogen chemical atoms. Quantum tunnelling couples these configurations; the resulting pair of energy levels permits an electric dipole transition. This geometrical inversion is distinct from population inversion, which concerns how many molecules occupy each level.^[ref]

For an idealized pair of energy levels with equal quantum degeneracy, the Boltzmann distribution gives $N_u/N_l=\exp [-h{\nu }_0/(k_{\mathrm{B}}T)]<1$ at positive temperature $T$. A warm source therefore supplies upper-state molecules together with an unwanted lower-state population. Selecting the former can produce population inversion in the transmitted beam without exciting each transmitted molecule during selection. The energy per net downward transition isHere $h$ is the Planck constant. Taking the reported ${\nu }_0=23.87\mathrm{GHz}$ as a rounded input gives $ϵ\approx 1.582\times 10^{-23}\mathrm{J}$. This calculation is a conversion of the reported frequency, not a new measurement of the ammonia inversion transition.

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]}

For two individual microstates, the partition function cancels:At positive temperature, a higher-energy individual microstate is less probable, and equal-energy individual microstates are equally probable. For example, a gap $E_i-E_j=k_{B}T\ln 2$, using the natural logarithm, gives $p_i/p_j=1/2$. This comparison concerns individual microstates; a whole macrostate can contain several of them.

Changing the reference for energy leaves every probability unchanged. If all energies become $E_i^{\prime }=E_i+C$ for the same constant $C$, thenThe partition function itself changes, but its common factor cancels in the Boltzmann distribution. Thus choosing a ground state energy of zero simplifies calculations without changing populations. The expected value of energy shifts by $C$, as it must when its reference changes.

Keep the gap $\Delta$ fixed and vary positive temperature. With $x=k_{B}T/\Delta$, the Boltzmann factor is $q=e^{-1/x}$. In the four-microstate example, as $T\to 0^{+}$, $q\to 0$: the unique ground state has probability tending to one, and the other microstates have probabilities tending to zero. If several microstates instead shared the lowest energy, their state degeneracy would split that limiting population equally.

As $T\to \infty$ in this fixed finite model, all four individual Boltzmann factors tend to one, so all four individual microstates have probability tending to $1/4$. The three aggregate probabilities instead approachThe equal-population limit applies to this finite set of microstates. An infinite collection of microstates cannot all have the same positive normalized probability, so the statement cannot be extended to an unbounded energy spectrum. Real systems may also acquire additional accessible microstates as temperature rises.

The plotted crossing occurs at $x=1/\ln 2$, where $P(0)=P(\Delta )=4/9$. For larger positive temperatures, the middle macrostate is more probable than the lowest-energy macrostate, although the lowest-energy individual microstate remains more probable than either individual middle microstate. The expected value of energy rises from zero toward $\Delta$ in this finite model.

Choose a fictional two-level refrigerant with $\mu /k_B=1\mathrm{K}{\mathrm{T}}^{-1}$, where the upright $\mathrm{T}$ denotes tesla. This chosen moment is a convenient model parameter, not a fitted property of the historical salt. At bath temperature $T_b=1\mathrm{K}$, raise $B$ from $0.20\mathrm{T}$ to $1.00\mathrm{T}$. Then $x$ increases from $0.20$ to $1.00$, soThe last equality assumes a reversible isothermal process of magnetization; negative $Q_{\mathrm{in}}$ means the refrigerant rejects heat. The bath gains $0.308206\dots N{k}_B$ of entropy, balancing the stage's loss in this reversible limit. An appreciable temperature difference during real heat transfer adds entropy production.

Now open the heat switch and reduce $B$ from $1.00\mathrm{T}$ to $0.10\mathrm{T}$. Constant $x=1$ gives $T_f=0.10\mathrm{K}$. Throughout this ideal cooling leg, $p_{-}=0.880797\dots$, $p_{+}=0.119203\dots$, and $S_{\mathrm{mag}}=0.365334\dots N{k}_B$. Neither state becomes more probable. Field reduction shrinks the energy gap by ten while preserving the Boltzmann distribution's dimensionless ratio of gap to temperature.

There is a subtle energy check. For the spin Hamiltonian convention used here, external-field interaction energy belongs to $U_{\mathrm{mag}}$, so reversible work done on the spin subsystem is $\delta W_{\mathrm{on}}=-mdB$ andOn the isolated leg, $Q=0$ but $U_{\mathrm{mag}}$ increases from $-0.761594\dots N{k}_B\mathrm{K}$ to $-0.0761594\dots N{k}_B\mathrm{K}$, becoming less negative as the gap shrinks. Integrating $-mdB$ at constant $x=1$ supplies exactly $+0.685435\dots N{k}_B\mathrm{K}$. Cooling therefore need not decrease this chosen spin internal energy. An account of work exchanged with a magnet circuit must specify where field interaction energy is assigned before comparing signs. The occupation ratio and temperature prediction do not depend on that bookkeeping choice.

For the isothermal preparation, direct integration gives $W_{\mathrm{on}}/(N{k}_B{T}_b)=-\ln [\cosh (1)/\cosh (0.20)]=-0.413913\dots$. The endpoint internal energy difference is $-0.722119\dots N{k}_B{T}_b$, equal to $Q_{\mathrm{in}}+W_{\mathrm{on}}$. This independently checks both the rejected heat and the magnetic-work sign in the chosen convention.

For bound Kepler orbits, the Delaunay variables obey $0<G<L$, $-G<H<G$, and three periodic angles. A canonical transformation preserves phase space volume. At fixed $L$, integrating $H$ supplies the measure $2GdG$. Since $G=L\sqrt{1-e^2}$, the fraction with orbital eccentricity less than $e$ is $L^{-2}{\int }_{L\sqrt{1-e^2}}^{L}2GdG=e^2$. Thus the probability density function is $2e$, independently of any normalizable energy-only phase-space distribution function. A Boltzmann distribution is unnecessary. Even the time-dependent density $F(E(L))[1+\eta \cos (l-n(L)t)]$, $|\eta |<1$, has this same marginal because integrating the mean anomaly removes its oscillatory part; an eccentricity marginal alone does not imply equilibrium.

For positive diffusivity $D(x)$, a diffusion equation with potential-driven drift has diffusive flux $J=-D(c_x+\beta V^{\prime }c)$. On a finite interval with zero flux at both endpoints, the total amount is conserved. The substitution $c=e^{-\lambda t}{e}^{-\beta V}y$ converts the decay eigenvalue problem to the Sturm-Liouville problemThe Rayleigh quotient is $\lambda =\int p|y^{\prime }{|}^2/\int w|y{|}^2\ge 0$. Its zero eigenfunction is constant, giving the stationary solution $c_{\mathrm{eq}}=C{e}^{-\beta V}$. This is a spatial Boltzmann distribution; $D$ changes relaxation rates but cancels out of the zero-flux equilibrium profile. For a sufficiently regular positive $D$ and bounded $V$, the constant eigenfunction is the unique zero mode up to scale.

An electrolyte is a conducting medium in which mobile ions carry electric current, according to IUPAC Gold Book: electrolyte. In introductory chemistry, a substance that yields ions when dissolved in water is also called an electrolyte. Both usages point to the same mechanism: an electric field gives mobile cations and anions directed drift velocities.

Mobility matters as much as the presence of electric charge. In a simple ideal model of solid sodium chloride, ions are confined to their sites and do not carry a sustained transport electric current. Molten sodium chloride has mobile ions; dissolved sodium chloride gives mobile ${\mathrm{Na}}^{+}$ and ${\mathrm{Cl}}^{-}$ ions in water. Real solids can have some ionic transport through imperfections; the immobile-solid picture is an introductory limit, rather than a claim that all solids are perfect insulators.

A strong electrolyte yields essentially all ions from the dissolved substance, whereas a weak electrolyte yields only a fraction under the specified conditions. Dissolved sodium chloride is a strong electrolyte. Sucrose is a nonelectrolyte: it dissolves as neutral molecules, so dissolving something does not automatically create mobile ions. Water's own small population of hydronium and hydroxide remains possible. See 11.2 Electrolytes (§11.2: strong, weak, and nonelectrolytes) and Flowers et al., Chemistry 2e, §11.4, Table 11.3: sucrose and sodium chloride.

Acetic acid is a weak electrolyte and a weak acid. In water, some neutral acetic acid molecules transfer a proton to water, producing hydronium and the negatively charged partner:Both sides have net electric charge zero. At chemical equilibrium, neutral acetic acid coexists with these ions; the fraction converted depends on the conditions. Electrolyte strength describes this extent of ion production, while concentration describes how much substance is present per specified quantity of chemical solution. A dilute strong electrolyte can contain fewer ions per volume than a concentrated weak electrolyte. Neither the words strong and weak nor a conductivity comparison alone specifies concentration.

Now specify an electrolytic cell with a positive anode on the left, a negative cathode on the right, and an approximately uniform electric field pointing right in the bulk electrolyte. With a steady external source and neglecting diffusion and convection, positive cations drift right toward the negative cathode; negative anions drift left toward the positive anode. Yet both populations contribute a rightward conventional current! For number densities $n_{+},n_{-}$ and drift velocities ${\mathbf{v}}_{+},{\mathbf{v}}_{-}$, the transport current density isThe product of electric charge and velocity determines conventional current: for the anion, both the charge and its rightward velocity component are negative, making its rightward contribution positive. This is the same sign rule used for Electron conventional current. It follows directly from the existing charge carrier formula.

For an independent numerical example, stipulate $n_{+}=n_{-}=10^{25}{\mathrm{m}}^{-3}$, $v_{+,x}=2\times 10^{-4}\mathrm{m}{\mathrm{s}}^{-1}$, and $v_{-,x}=-3\times 10^{-4}\mathrm{m}{\mathrm{s}}^{-1}$. With the exact elementary charge $e=1.602176634\times 10^{-19}\mathrm{C}$, the contributions areEqual number densities of singly charged cations and anions give zero bulk charge density, while their motion gives nonzero current density. This is why a neutral electrolyte can conduct. These numbers are chosen toy inputs, not measurements of sodium chloride; maintaining sustained cell electric current also requires suitable electrode reactions. Blocking electrodes can instead lead to charge accumulation and decay of the electric current.

An electrolyte connects the microscopic motion of ions to observable electrical conductivity. A counterion-only model keeps just the mobile counterions needed to oppose a fixed electric charge; Poisson-Boltzmann equations use ideal Boltzmann distributions to model their equilibrium distributions, rather than the drift example above.

Freeze the magnitude of an order parameter and retain just one common orientation $\theta$ for $N$ rigidly aligned cells. Give this collective coordinate energy $E(\theta )=-Nh\cos \theta$ at positive temperature $T$. The Boltzmann distribution is the displayed normalized probability density function, where $k_B$ is the Boltzmann constant. Here $h$ is energy per cell, a different normalization from the source density in continuous vacuum selection by an external source. This is a deliberately restricted model of selection among already ordered orientations: it does not derive alignment or an ordering critical temperature from a spatial interaction.

Its mean orientation is an elementary ratio of integrals:At $h=0$ every angle has equal probability, and the numerator vanishes exactly. For small $x$, expansion of the exponential gives $m_N(h)=x/2+O(x^3)$ because the angular mean of ${\cos }^2\theta$ is $1/2$. For $x=1$ and $5$ the means are approximately $0.446390$ and $0.893383$. At fixed $h>0$, increasing $N$ concentrates the probability density function near $\theta =0$; expanding $\cos \theta =1-{\theta }^2/2+O({\theta }^4)$ gives angular variance approximately $1/x$. ConsequentlyThe uniform angular average remains a symmetric mixture if one sets $h=0$ first. The other order selects a definite direction. This is the continuous-orientation version of the symmetry-breaking order of limits; it illustrates selection without replacing the separate problem of stability against spatial fluctuations.

On a finite state space $F$ with finite objective function values, the fixed-$T$ target is a Boltzmann distribution:The partition function $Z_T$ normalizes the weights but cancels from the Metropolis–Hastings acceptance probability. A short detailed balance proof explains the target. With symmetric $q$ and $x\ne y$,Summing over $x$ gives ${\pi }_T{P}_T={\pi }_T$, so this is a stationary distribution. On a finite state space, a connected graph of positive symmetric proposals makes the accepted Markov chain an irreducible Markov chain. Convergence from any starting state additionally requires an aperiodic Markov chain; one positive self-transition suffices under irreducibility. Replacing $q$ by $(I+q)/2$ ensures this. For asymmetric $q$, require connectivity through bidirectional proposals: one-way proposals are rejected.^[ref]

In the five-state barrier example for simulated annealing, $T_{\ast }=1/\ln 2$ gives exact Boltzmann weights $2^{-E(x)}$. Multiplying by $16$ givesFor the uphill edge, detailed balance is visible numerically:Thus the global minimum has the largest individual stationary probability, $16/33$, but is far from certain at this positive $T_{\ast }$. At $T_{\ast }/2$ the exact stationary distribution is $(16,64,1,16,256)/353$.

For any finite $F$, let $M$ be the set of global minima and $E_{\min }$ the smallest objective function value. Cancelling $e^{-E_{\min }/T}$ from numerator and denominator shows that each nonminimal state's weight tends to zero as $T\downarrow 0$, while every state in $M$ has equal weight. Therefore ${\pi }_T$ tends to the uniform distribution on $M$. This is a statement about a family of equilibrium probability distributions; it does not prove that a cooling Markov chain follows that family.

Take three Ising spins $s_i\in \{-1,1\}$, with an antiferromagnetic interaction of energy scale $J$ on each edge of a triangle. An individually satisfied bond has $s_i{s}_j=-1$. But multiplying all three bond products gives $(s_1{s}_2)(s_2{s}_3)(s_3{s}_1)=1$, whereas three satisfied bonds would give $-1$. At least one bond must therefore remain unsatisfied. This is an exactly finite example of geometrical frustration.

There are six ground states with one unsatisfied bond and ground-state energy $-J$. The two all-aligned microstates have energy $3J$. No interaction is switched off in a ground state: the system minimizes their sum by accepting one local disappointment. Allowing continuous three-component spins instead of Ising spins changes this optimization problem and permits the $120$-degree arrangement described under geometrical frustration.

Writing $\beta =1/(k_{B}T)$ for inverse temperature, the exact partition function and mean internal energy areAt $T>0$ these expressions are smooth: a finite triangle has no thermodynamic phase transition. In zero magnetic field, spin inversion symmetry also gives $⟨s_i⟩=0$ at every site, even as the Boltzmann distribution concentrates on the six frustrated ground states.

For transverse coupling with equal mode amplitudes $\lambda$, Fermi's golden rule and the oscillator creation operator and annihilation operator givewhere the first index is the final spin and $g$ is the density of states per unit angular frequency. Decay creates a quantum and carries the factor $n_{\Omega }+1$; excitation annihilates one and carries $n_{\Omega }$. For nonuniform resonant occupations, $n_{\Omega }$ is their density-weighted mean.

With a thermal Bose-Einstein distribution, $n_{\Omega }/(n_{\Omega }+1)=e^{-\beta \hslash \Omega }$. The two-state master equation therefore has the Boltzmann distribution as its unique stationary populations when ${\Gamma }_0>0$. The population deviation relaxes at rate ${\Gamma }_0(2n_{\Omega }+1)$. If there is no resonant coupling, both rates vanish and stationarity does not determine the populations.

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]

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.

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]