For a simple homogeneous thermodynamic system, internal energy is extensive: . Euler theorem for homogeneous functions and give , hence
Also , so at fixed ,
At coexistence gas and liquid have the same temperature, pressure and chemical potential. Integrating along the intervening isotherm gives . Integration by parts, with endpoints at the same coexistence pressure , gives the Maxwell construction
The horizontal line at replaces the mechanically unstable loop of the van der Waals isotherm. The extensivity argument presupposes the ordinary thermodynamic limit and neglects macroscopic interface contributions.
In the microcanonical ensemble, the microcanonical energy-shell entropy is the Boltzmann entropy
For a macroscopic system, changing the small energy resolution by a fixed factor changes only by a nonextensive additive term, while . Its relative effect therefore vanishes in the thermodynamic limit.
The thermodynamic temperature, pressure, and chemical potential are defined in the entropy representation by
Thus
or, equivalently, the first law of thermodynamics in fundamental form is
Expanding the trace in the spin basis gives
so the transfer matrix for the one-dimensional Ising model reproduces the partition function.
In the ordered basis ,
Its eigenvalues are
Since at every positive temperature, the thermodynamic limit gives the Helmholtz free energy per spin
This is an analytic function of and for , so the one-dimensional short-range Ising model has no finite-temperature phase transition.
Differentiate the bound with respect to one , using and :
Thus an interior minimizing Gaussian variational kernel for a gradient quartic interaction has
Use the half-space mode density given in the question. In the thermodynamic limit, the coupled self-consistency equations are
The admissible solution must have . The coefficient is unchanged because the interaction's derivative contributes only a constant and a term.
The formal continuum integrals require the coarse-graining ultraviolet cutoff in three dimensions: the first radial integrand tends to a constant at large , so taking would leave divergent. The low-wavevector fluctuation argument is independent of this short-distance regularization.
Use the stable parameters . The bulk stationary points satisfy . For , has negative curvature, while the two minima are
Both have and . To establish minimization at fixed total composition, write
For mean composition with , a mixture of the two minima with volume fraction in the positive phase attains this lower bound in the bulk. Macroscopic slabs have only a subextensive interfacial excess free energy, so the bound is attained per unit volume in the thermodynamic limit. Thus constrained minimization gives phase coexistence in the mean-field approximation; at the endpoints only one phase is needed.
The two equilibrium conditions also follow directly by constrained optimization. Introduce a Lagrange multiplier for in the mixture's bulk free-energy density. Varying and gives ; varying gives . These equate the chemical potentials and pressures, respectively.
Equivalently, their bulk chemical potentials are both zero and their pressures are equal:
These are precisely the conditions of the common-tangent construction for phase coexistence: the horizontal supporting tangent touches both minima. Equal chemical potentials prevent net composition exchange, and equal pressure gives mechanical equilibrium at a flat interface.
Let the pair potential be , so that in the dilute limit only one interacting pair need be retained at a time. The leading Mayer f-function expansion is
For each pair, translation of its center coordinate supplies and the other particles supply . Neglecting boundary corrections,
The Taylor expansion and in the thermodynamic limit give
Using and the ideal gas free energy,
With number density , this is , so is precisely the second virial coefficient.
The conserved number operator is
Substitution of into the Landau-Ginzburg theory gives
The imaginary term is a total derivative, so
Thus is the number density and is the canonical momentum conjugate to . The canonical commutation relation is . Consequently, for in volume ,
Equivalently, the averaged phase obeys . This number-phase conjugacy means that a state of sharp has no sharp phase, whereas a phase-selected state exhibiting spontaneous symmetry breaking must superpose different number sectors. Such sectors become effectively degenerate in the thermodynamic limit.
Translational invariance makes the pointwise variance independent of . Applying Parseval identity to the stated Fourier convention gives
In the thermodynamic limit, the reciprocal-lattice sum becomes
At each point is a centered Gaussian random variable, so the supplied absolute third moment gives
Temperature as an equilibrium relation. Classical thermodynamics begins with the zeroth law of thermodynamics: thermal equilibrium is transitive. A thermometer can therefore assign a common empirical temperature to systems in mutual thermal equilibrium. Different thermometer properties initially give different monotone empirical scales; a resistance reading by itself does not define the absolute scale.
The Second law of thermodynamics provides a substance-independent thermodynamic scale. A reversible engine operating between reservoirs has a heat ratio depending only on their temperatures. One chooses absolute temperature so that , with the positive heat magnitudes. Its reversible efficiency is . Equivalently, reversible heat defines entropy through . For an equilibrium simple compressible system,
Classical thermodynamics uses this relation without needing a microscopic account of entropy. A dilute-gas thermometer realizes an empirical scale approaching the thermodynamic scale after corrections for nonideal behavior.
Statistical meaning of the same derivative. In the microcanonical description, an equilibrium energy shell containing accessible states has . Its statistical temperature is consequently
To see why this is an equilibrium temperature, let two weakly coupled systems exchange a small energy while their total energy stays fixed. The combined multiplicity is . Maximizing its logarithm gives
which is equality of the two inverse temperatures. Stability requires the total entropy maximum, not only stationarity. The overwhelmingly dominant macrostate in a large system therefore realizes the thermodynamic equilibrium condition.
A large reservoir gives a small subsystem the canonical probabilities , with . In that ensemble, and agree with macroscopic thermodynamic relations in their usual equilibrium regime. The equipartition identity for a classical dilute monatomic gas, , is a useful operational connection between microscopic motion and temperature, but is not a universal definition: quantum modes need not obey equipartition, and a system away from equilibrium need not have one temperature at all.
For finite systems one must specify the entropy convention and energy resolution; surface and fluctuation corrections can distinguish definitions that agree in the ordinary thermodynamic limit. Bounded-energy population-inverted systems can also have a negative entropy derivative and hence a negative statistical temperature. That possibility lies outside an ordinary positive-temperature ideal-gas thermometer; it does not mean temperature is merely an arbitrary average kinetic energy.
Obtaining the energy–temperature conversion. The Boltzmann constant connects entropy measured in with the dimensionless logarithm of a state count. Historically it could be determined from the ideal gas law : measure the molar gas constant and independently determine the Avogadro constant , then obtain . The independent particle-number determination matters; using alone would not fix energy per particle.
An electrical method uses Johnson-Nyquist noise. A resistor of resistance at known temperature has, in the classical range , a one-sided open-circuit voltage-noise spectrum . Over a calibrated flat bandwidth ,
A practical measurement uses an equilibrium resistor, a calibrated low-noise amplifier and a known transfer function, subtracts or independently characterizes amplifier noise, and replaces a nominal bandwidth by the actual effective noise bandwidth. Varying temperature and fitting the mean-square voltage is more informative than treating one amplified voltage sample as the thermal variance. Other primary methods relate temperature to sound speed in a dilute gas or to calibrated particle diffusion and mobility.
There is a historically relevant unit-definition qualification: by the date of this exam, the revised SI had already taken effect on 20 May 2019. The kelvin is now defined by fixing
The same experiments now realize or compare thermodynamic temperatures using this fixed constant, whereas historically they measured the conversion constant on an independently defined temperature scale. The physical distinction between equilibrium temperature and microscopic state counting is unchanged.[ref]
At zero field, write the right-hand side of the self-consistency equation as , with . The origin is always a solution, and the slope there is
The critical crossing occurs when this equals one, giving
This is a continuous phase transition in the mean-field approximation. For ,
so the positive branch is a strictly concave function, while is odd and bounded by one. For , the slope at zero is at most one and strict concavity puts below for every , leaving only the unmagnetized solution. For , it initially lies above but has ; hence there is one positive intersection and, by spin inversion symmetry, one negative intersection.
The curvature of the Helmholtz free energy in the mean-field approximation is per site. It is negative at below and positive at the nonzero intersections, identifying the ordered minima and the unstable zero solution.
The critical temperature and branches here belong to the mean-field approximation. In the exact finite system, spin inversion symmetry forces zero average magnetization at zero field; selecting a broken-symmetry phase requires the thermodynamic limit with a vanishing field of specified sign.
Excited-state saturation. Consider first a uniform three-dimensional noninteracting gas of identical spinless bosons in a large box. Set the one-particle ground energy to zero. The Bose-Einstein distribution is , with quantum gas fugacity . For dispersion , the excited density of states is
Using the thermal de Broglie wavelength , integration gives
The polylogarithm is . Its finite value at implies a maximum excited population , with . Once this is less than , the excess particles occupy the ground state macroscopically. This is Bose-Einstein condensation, rather than occupation of a state by merely one or a few atoms.
Transition and condensate fraction. Setting the saturation number equal to gives the homogeneous ideal Bose gas condensation temperature
Below this temperature, approaches zero from below in the thermodynamic limit, and particle conservation gives
Above there is no macroscopic condensate and adjusts to keep fixed. At finite , the chemical potential remains slightly below the ground energy and the transition is rounded.
Figure 1.
Original uniform ideal-gas condensate fraction and schematic accumulation in the lowest state below the critical temperature; the exponent is for a three-dimensional box
. Details.
Energy and the role of geometry. The excited-state internal energy is
The zero-energy condensate contributes no term to this energy in the ideal model. Below ,
Condensation depends on low-energy state counting: for a uniform quadratic-dispersion gas, a finite-temperature bulk condensate requires spatial dimension greater than two. A trap changes the density of states and therefore the temperature law; a three-dimensional harmonic trap gives a cubic, rather than , condensate-fraction exponent. Finally, an ideal-gas condensate should not automatically be called a superfluid: its free quadratic spectrum has zero critical speed by the Landau criterion. Interactions matter for the robust collective-flow properties of real dilute condensates.
Order parameter and the free-energy landscape. In Landau theory, take a scalar magnetization component at fixed volume and composition. It is an order parameter: its nonzero spontaneous value distinguishes an ordered phase with ferromagnetism from the symmetric high-temperature phase. One first constructs a constrained order-parameter free energy at specified magnetization by summing over microscopic states compatible with it. An external magnetic field is an independent control, so it can also appear in that constrained function. For a total magnetic moment the conjugate term is ; if is a magnetization density, the total-energy term is , or in free-energy density. Minimizing over then yields the equilibrium thermodynamic potential as a function of the imposed field.
Thus treating both and as arguments of a Landau function does not treat them as independent natural variables of the already-minimized equilibrium free energy. The first is the fluctuating or constrained order parameter, and the second its conjugate field.
At zero field, reversal symmetry requires , excluding odd powers. Near a continuous transition, assume the homogeneous constrained free energy is analytic for small , that the quadratic coefficient crosses zero linearly, and that the quartic coefficient is positive. Then
The temperature is where the zero-field quadratic coefficient vanishes and spontaneous ordering begins in this approximation. Higher powers are negligible sufficiently close to it; spatial fluctuations are neglected in this mean-field treatment.
Stationarity gives . Above , the only minimum is at zero. Below it, zero becomes a maximum and two minima appear:
The two signs describe spontaneous symmetry breaking: they are the average magnetizations in the symmetry-broken pure phases. A strictly finite system at exactly zero field samples the two wells symmetrically and has ensemble average . A spontaneous pure-phase value is defined by selecting a phase, for example taking the thermodynamic limit before removing a small field. This qualification distinguishes a minimum's location from the average over a symmetric two-well distribution.
Figure 1.
Quartic Landau free energy above and below the transition, showing the disordered and symmetry-broken minima
. Details.
Minimum free energy and the heat-capacity jump. At a nonzero minimum, . Substituting this value gives
Hence
while above the transition. Take the temperature derivatives only after minimizing, so the equilibrium order parameter is allowed to change with temperature. With ,
Therefore
This is the heat-capacity jump in quartic Landau theory. The entropy is continuous at , so there is no latent heat. If denotes free-energy density rather than total free energy, these expressions denote heat capacity per unit volume with the same normalization.
Stabilizing a tricritical expansion. At a tricritical point, the quartic coefficient is tuned to zero at the same point where the quadratic coefficient changes sign. Reversal symmetry still excludes odd powers. The next available stabilizing term is therefore
A negative or vanishing leading sextic coefficient would not stabilize this truncation for large ; a higher positive term would then be needed and the present tricritical form would be insufficient. The ordinary quartic theory has a continuous transition for , while a stabilized theory with can have a first-order transition. Their meeting at zero field is what makes the simultaneous tuning tricritical.
The tricritical minimum. At zero field and , stationarity gives
The nonzero solutions are minima: their curvature is . Put . Since at stationarity,
Thus
The order parameter now vanishes as rather than the quartic theory's square root.
Heat-capacity singularity. Define , so the ordered-side singular free energy is . Its derivatives are
Consequently
This is the tricritical heat-capacity amplitude in a sextic Landau expansion. It is positive on the ordered side and diverges as , rather than remaining a finite jump. This value of the heat-capacity critical exponent is the mean-field result along the path keeping the quartic term zero; a detuning to ultimately restores ordinary quartic behavior. Fluctuation corrections to the mean-field theory are not included in the stipulated polynomial minimization.
Measure single-particle energy from the lowest state and write , with . In a homogeneous system, the number of wavevector states in a thin spherical shell is proportional to . Changing variables from to gives the density of states per unit spatial volume,
where is the area of the unit sphere in dimensions; an additional internal state degeneracy would multiply this expression by a constant.
The Bose-Einstein distribution has . At saturation, , the maximum excited-state number density is
Only the lower endpoint can cause divergence. There the integrand is proportional to
The improper integral converges precisely when its power exceeds . Thus
When , the divergence is logarithmic; when , it is a power divergence. In either case the excited states can accommodate any fixed finite density at every positive temperature, so there is no bulk Bose-Einstein condensation at nonzero temperature.
For completeness, putting and using the Gamma function and Riemann zeta function gives
At fixed density , the right-hand side reaches at ; below it the excess particles occupy the lowest state, with condensate fraction . This is Bose condensation with power-law dispersion in the homogeneous thermodynamic limit. A finite vessel or a trapping potential changes the low-energy state count and must not be silently substituted for that setting.
A phase transition is a qualitative change of thermodynamic phase as control parameters such as temperature, pressure or chemical potential are varied. A genuine equilibrium singularity is a nonanalyticity of the appropriate thermodynamic potential in the thermodynamic limit; finite systems generally show rounded crossovers instead.
An order parameter is a macroscopic quantity that distinguishes the phases. It often vanishes in a symmetric disordered phase and becomes nonzero after spontaneous symmetry breaking, although density-driven transitions need not have that simple symmetry interpretation. It can be a scalar, vector or tensor. At a first-order phase transition an appropriate first derivative of the free energy is discontinuous, and the order parameter commonly jumps; at a continuous phase transition the order parameter can approach zero continuously while a susceptibility or a higher derivative becomes singular.
Three examples are:
These examples also illustrate why an order parameter must reflect the physical distinction between phases rather than merely be any quantity that changes with temperature.
For attractive nearest-neighbour interactions, , the transformed model is a ferromagnetic Ising model. In three dimensions and at , its infinite-volume spontaneous magnetization has two branches. As the effective field crosses zero, the equilibrium mean spin changes discontinuously from a negative value to a positive value. Since , varying the physical chemical potential crosses this same first-order phase transition.
The particle-hole-symmetric coexistence line is
The mean Ising spin is a rescaled occupation density. Writing for the physical number density,
The negative-spin phase is a dilute gas and the positive-spin phase is a dense, liquid-like phase. If the spontaneous spin magnitude at coexistence is , their number densities are
Thus the physical order parameter is a density or occupation imbalance, not an actual magnetic moment of these particles. In a closed system with an intermediate mean density, the same transition appears as phase coexistence between the two densities. At the density contrast vanishes continuously; above , sweeping the chemical potential gives no coexistence jump. A true discontinuity requires the thermodynamic limit; a finite sum over finitely many cells remains analytic.
Pure thermodynamic phase Created 2026-10-06 Updated 2026-10-14
A pure thermodynamic phase is an extremal equilibrium Gibbs state rather than a probabilistic mixture of distinct macroscopic phases. On an ordered Ising branch, taking the thermodynamic limit before removing a selecting conjugate field gives a nonzero magnetization. A symmetric mixture can instead have zero mean but a nondecaying two-point contribution. The connected correlation function and correlation length describing ordinary critical fluctuations should be evaluated in a chosen pure phase when this distinction matters.
Replica trick 2026-10-18
For positive , the displayed identity follows by differentiating at zero.
Integer powers of a partition function represent independent copies with the same disorder. Integer moments alone do not justify analytic continuation or interchanging disorder averages and thermodynamic limits.[ref]
Spontaneous symmetry breaking Created 2026-09-24 Updated 2026-10-18
Spontaneous symmetry breaking occurs when the action and equations have a symmetry but a chosen ground state is invariant only under a proper subgroup. Acting with the broken symmetry generates a degenerate vacuum manifold. An order parameter distinguishes the selected states: the equations do not prefer its direction, yet a particular state does. This lets one use the symmetry of the laws to predict excitations around a state that has less symmetry, notably the Goldstone bosons of a broken continuous symmetry.
A finite system with symmetry-preserving boundary conditions can have a symmetric ground state; a symmetric statistical mixture of broken states also has vanishing mean order parameter. A nonzero mean therefore requires a stated selection procedure. The symmetry-breaking order of limits takes the thermodynamic limit in a small applied source before removing that source. The finite-volume angular source selection example makes the distinction explicit for a continuous symmetry.
Stable pair potential 2026-10-14
A pair potential is stable when its total interaction energy obeys this lower bound for every allowed configuration, with a constant independent of . It prevents collapse with an arbitrarily negative energy per particle. Stability alone does not control the long-distance contribution needed for the ordinary thermodynamic limit.
For finite-range interactions and bounded wall couplings, only sites within an interaction range of the wall have altered local environments. In a regular container of linear size and dimension , the surface area scales as while volume scales as . A boundary contribution of order therefore gives an order- correction to free-energy density. Away from a critical phenomenon, a regular surface contribution often has the form , with a wall Helmholtz free energy per area. Critical finite-size rounding, edges and long-range interactions require separate care; vanishing relative boundary size does not guarantee a simple surface expansion in every regime.
An exact check uses the zero-field nearest-neighbour Ising model on a chain of sites with fixed coupling and . The open chain has bonds, while the periodic ring has . Their finite partition functions and dimensionless Helmholtz free energies per spin are
For the open chain, sum the last Ising spin successively: each removed bond contributes , leaving a factor of two. For the ring, the transfer matrix for the one-dimensional Ising model has eigenvalues and , and its th power's trace gives the displayed sum. The two ends are a constant-size boundary in one dimension, so the open-chain total correction is , and its per-spin correction is exactly proportional to .
Figure 1.
Exact boundary corrections for a finite Ising chain and ring
. Details.
Both sets of finite results approach the same at the fixed plotted . Multiplying the correction by exposes the open-chain boundary term, while the ring correction vanishes exponentially at fixed finite . This model's thermodynamic limit remains a real analytic function at every positive temperature: having a thermodynamic limit does not require a phase transition.