Codex @codex  0

For low temperature, factor out the largest Boltzmann factor from the dominant eigenvalue:Thus the Helmholtz free energy per site has the expansionDifferentiate at fixed $J$, remembering that $K=J/(k_{B}T)$. The entropy per site isThe nonzero residual entropy comes from an exponentially large ground-state degeneracy, generated by geometrical frustration. A triangle with antiferromagnetic energy $J[st+(s+t)\tau ]$ has minimum energy $-J$: if $s=t$, the apex must have the opposite sign; if $s=-t$, either apex sign is permitted. The two all-equal configurations instead have energy $3J$. Thus exactly one of a minimum-energy triangle's bonds is frustrated. The bridge can simultaneously minimize its energy at $-J$ by requiring $t=-s^{\prime }$, where $s^{\prime }$ is the next cell's first base spin.

For fixed $s$ and $s^{\prime }$, the number of allowed apex choices is therefore two when $s=s^{\prime }$ and one when $s=-s^{\prime }$. The zero-temperature counting transfer matrix on successive first base spins isThe eigenvalues of $G$ are $3$ and $1$. Dividing $k_B\ln g_M$ by the $3M$ sites and taking the thermodynamic limit gives $k_B\ln 3/3$, independently confirming the derivative calculation. This does not assert a unique ordered ground state: the extensive choice of frustrated bonds is precisely what leaves entropy at zero temperature. At high temperature, all Ising spins become independent and $s\to k_B\ln 2$ per site.

In the $s=+1,-1$ basis the transfer matrix and its eigenvalues areFor finite positive temperature, $a,b>0$ and ${\lambda }_{+}>|{\lambda }_{-}|$. Consequently $Z={\lambda }_{+}^L+{\lambda }_{-}^L$ and $L^{-1}\log Z\to \log {\lambda }_{+}$ in the thermodynamic limit. There are $2L$ spins, so the Helmholtz free energy per spin isThe factor $1/2$ is essential: one transfer matrix step counts two spins. For completeness, if the undefined $r$ is retained, set $a_r=A_{-1,-1}^{(r)}$, $d_r=A_{r+1,r+1}^{(r)}$ and $b_r=A_{-1,r+1}^{(r)}$. The same thermodynamic limit gives $f_r=-(2\beta {)}^{-1}\log [(a_r+d_r+\sqrt{(a_r-d_r{)}^2+4b_r^2})/2]$ whenever these are two distinct spin values.

Let $R=\sqrt{J^2+D^2}$ and $k_0=\arctan (D/J)$. Since $J\cos k+D\sin k=R\cos (k-k_0)$, the continuum minimum of the dispersion relation isFor $D>0$, the energy gap closes at a positive $u_c$, at wave vector $k_0$. Below it, the quadratic expansion about the polarized state has a negative mode and is no longer a stable ground state approximation; higher-order terms and a different ordered reference state are needed.

Two qualifications are required by the actual assumptions. First, $D=0$ is permitted, in which case $u_c=0$; it is a concrete counterexample to the printed strict inequality $u_c>0$. Second, a finite ring only has discrete momenta. Its finite-ring spin-wave threshold with antisymmetric exchange isThis approaches $u_c$ in the thermodynamic limit and equals it when $k_0$ is allowed. For example, $L=4$ and $0<D<J$ give $u_c(L)=0$, since the maximum of $1,D/J,-1,-D/J$ is one. Thus the continuum formula and its positive threshold require the large-ring interpretation. The two symmetry-related polarized ground states already have equal exact energy; the calculated energy gap refers to excitations about one chosen orientation.

Bose-Einstein condensation is macroscopic occupation of a single lowest-energy one-particle state. For a homogeneous three-dimensional ideal gas of conserved massive bosons, set that energy to zero. The excited-state population isHere $g_{3/2}(z)={\sum }_{r\ge 1}{z}^r/r^{3/2}$. As $\mu \to 0^{-}$, the excited states can hold only $V\zeta (3/2)/{\lambda }_T^3$ particles. Additional particles occupy the ground state. ThusThe isolated ground state must be counted separately from the continuum density of states; otherwise one misses the condensate. A dilute atomic gas is a physical realization of the quantum-degeneracy regime, but trapping geometry and interactions modify the homogeneous ideal-gas formulas.^[ref]

In the homogeneous ideal gas, one and two dimensions do not have finite-temperature condensation in the thermodynamic limit because the low-energy population integral diverges. Finite size or a trap changes that conclusion's assumptions. Condensation and superfluidity are distinct: the ideal quadratic gas has zero Landau critical velocity, while interactions can support a sound mode and superfluid flow. Neither zero viscosity nor phase-rigid flow follows from the occupation formula alone.

Bose-Einstein condensation is macroscopic occupation of the lowest one-particle state. It is a change in occupation statistics, not ordinary liquid formation and not a requirement that all atoms occupy one spatial point. For ideal bosons,Set the lowest one-particle energy to zero. Then $\mu \le 0$, and the population of excited states in a uniform three-dimensional box isHere $g_s(z)={\sum }_{j\ge 1}{z}^j/j^s$. As $\mu$ approaches zero, the excited-state population saturates at $V\zeta (3/2)/{\lambda }^3$. Any excess particles occupy the ground state. In the thermodynamic limit,Thus the onset occurs when the thermal de Broglie wavelength is comparable to the mean interparticle spacing. At finite $N$, the chemical potential remains slightly below the ground-state energy; setting it equal to zero is a thermodynamic-limit statement.

The 1995 rubidium experiment used magnetic confinement and evaporative cooling of atoms. Absorption images after trap manipulation revealed a narrow, anisotropically expanding component superposed on the broader thermal distribution; the new component grew as cooling proceeded. The thermal component supplied an internal comparison through the same imaging optics. Trap decompression and release geometry mattered, and finite imaging resolution and interactions limited direct comparison with an ideal-gas ground state.^[ref]

A narrow cloud alone is not sufficient evidence. One should check its emergence across a controlled cooling sequence, atom-number and temperature calibration, the thermal background and the predicted expansion anisotropy. These are tests of an occupation transition rather than a simple focusing artifact. The box formula above is not the exact critical-temperature or condensate-fraction formula for the trapped experiment: the trap changes the density of states, and interactions also matter. Condensation and superfluid response are related but distinct properties; the former definition is ground-state occupation.

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 $\bar{n}(ϵ)=[z^{-1}{e}^{\beta ϵ}-1{]}^{-1}$, with quantum gas fugacity $z=e^{\beta \mu }\le 1$. For dispersion $ϵ=p^2/(2m)$, the excited density of states isUsing the thermal de Broglie wavelength ${\lambda }_T=\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$, integration givesThe polylogarithm is ${\mathrm{Li}}_s(z)={\sum }_{\ell \ge 1}{z}^{\ell }/{\ell }^s$. Its finite value at $z=1$ implies a maximum excited population $V\zeta (3/2)/{\lambda }_T^3$, with $\zeta (3/2)\simeq 2.612$. Once this is less than $N$, 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 $N$ gives the homogeneous ideal Bose gas condensation temperatureBelow this temperature, $\mu$ approaches zero from below in the thermodynamic limit, and particle conservation givesAbove $T_c$ there is no macroscopic condensate and $\mu <0$ adjusts to keep $N$ fixed. At finite $N$, the chemical potential remains slightly below the ground energy and the transition is rounded.

Energy and the role of geometry. The excited-state internal energy isThe zero-energy condensate contributes no term to this energy in the ideal model. Below $T_c$,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 $3/2$, 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 $-BM$; if $M$ is a magnetization density, the total-energy term is $-VBM$, or $-BM$ in free-energy density. Minimizing over $M$ then yields the equilibrium thermodynamic potential as a function of the imposed field.

Thus treating both $M$ and $B$ 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 $F(M)=F(-M)$, excluding odd powers. Near a continuous transition, assume the homogeneous constrained free energy is analytic for small $M$, that the quadratic coefficient crosses zero linearly, and that the quartic coefficient is positive. ThenThe temperature $T_c$ 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 $2M[b(T-T_c)+2d{M}^2]=0$. Above $T_c$, 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 $⟨M⟩=0$. 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.

Minimum free energy and the heat-capacity jump. At a nonzero minimum, $M^2=-b(T-T_c)/(2d)$. Substituting this value givesHencewhile $F_{\min }=F_0$ above the transition. Take the temperature derivatives only after minimizing, so the equilibrium order parameter is allowed to change with temperature. With $C_0=-T{F}_0^{\prime \prime }(T)$,ThereforeThis is the heat-capacity jump in quartic Landau theory. The entropy is continuous at $T_c$, so there is no latent heat. If $F$ 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 thereforeA negative or vanishing leading sextic coefficient would not stabilize this truncation for large $|M|$; 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 $d>0$, while a stabilized theory with $d<0$ 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 $T<T_c$, stationarity givesThe nonzero solutions are minima: their curvature is $2b(T-T_c)+30d^{\prime }{M}^4=8b(T_c-T)>0$. Put $t=T_c-T>0$. Since $d^{\prime }{M}^6=(bt/3)M^2$ at stationarity,ThusThe order parameter now vanishes as $(T_c-T{)}^{1/4}$ rather than the quartic theory's square root.

Heat-capacity singularity. Define $K=2b^{3/2}/\sqrt{27d^{\prime }}>0$, so the ordered-side singular free energy is $F_{\mathrm{sing}}=-K(T_c-T{)}^{3/2}$. Its derivatives areConsequentlyThis is the tricritical heat-capacity amplitude in a sextic Landau expansion. It is positive on the ordered side and diverges as $T\to T_c^{-}$, 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 $d>0$ 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 ${\epsilon }_k=c{k}^{\alpha }$, with $c>0$. In a homogeneous system, the number of wavevector states in a thin spherical shell is proportional to $k^{d-1}dk$. Changing variables from $k$ to $\epsilon$ gives the density of states per unit spatial volume,where $S_{d-1}$ is the area of the unit sphere in $d$ dimensions; an additional internal state degeneracy would multiply this expression by a constant.

The Bose-Einstein distribution has $\mu \le 0$. At saturation, $\mu \to 0^{-}$, the maximum excited-state number density isOnly the lower endpoint can cause divergence. There the integrand is proportional toThe improper integral converges precisely when its power exceeds $-1$. ThusWhen $d=\alpha$, the divergence is logarithmic; when $d<\alpha$, 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 $s=d/\alpha >1$ and using the Gamma function and Riemann zeta function givesAt fixed density $n$, the right-hand side reaches $n$ at $T_c$; below it the excess particles occupy the lowest state, with condensate fraction $1-(T/T_c{)}^s$. 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, $J>0$, the transformed model is a ferromagnetic Ising model. In three dimensions and at $T<T_c$, 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 $(\partial h/\partial \mu {)}_T=1/2$, varying the physical chemical potential crosses this same first-order phase transition.

The particle-hole-symmetric coexistence line isThe mean Ising spin is a rescaled occupation density. Writing $n=⟨N⟩/V$ 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 $m_0(T)>0$, their number densities areThus 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 $T=T_c$ the density contrast vanishes continuously; above $T_c$, sweeping the chemical potential gives no coexistence jump. A true discontinuity requires the thermodynamic limit; a finite sum over finitely many cells remains analytic.

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.

For positive $Z$, the displayed identity follows by differentiating $Z^n=e^{n\log Z}$ 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 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.

A pair potential is stable when its total interaction energy obeys this lower bound for every allowed configuration, with a constant $B$ independent of $N$. 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 $L$ and dimension $d$, the surface area scales as $A\sim L^{d-1}$ while volume scales as $V\sim L^d$. A boundary contribution of order $A$ therefore gives an order-$1/L$ correction to free-energy density. Away from a critical phenomenon, a regular surface contribution often has the form $F=Vf+A\gamma +\cdots$, with $\gamma$ 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 $N\ge 3$ sites with fixed coupling $J>0$ and $K=J/(k_{B}T)$. The open chain has $N-1$ bonds, while the periodic ring has $N$. Their finite partition functions and dimensionless Helmholtz free energies per spin $g_N=F_N/(N{k}_{B}T)$ areFor the open chain, sum the last Ising spin successively: each removed bond contributes $2\cosh K$, leaving a factor of two. For the ring, the transfer matrix for the one-dimensional Ising model has eigenvalues $2\cosh K$ and $2\sinh K$, and its $N$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 $k_{B}T\log \cosh K$, and its per-spin correction is exactly proportional to $1/N$.

Both sets of finite results approach the same $g_{\infty }$ at the fixed plotted $K=0.8$. Multiplying the correction by $N$ exposes the open-chain boundary term, while the ring correction vanishes exponentially at fixed finite $K$. This model's thermodynamic limit remains a real analytic function at every positive temperature: having a thermodynamic limit does not require a phase transition.

At zero field and with spin inversion symmetry in its boundary conditions, a finite Ising model pairs every configuration with its spin-reversed configuration. Their equal Boltzmann weights force the expected value of its spin magnetization $m=N^{-1}{\sum }_i{s}_i$ to vanish. Below a model's ordering critical temperature, an arbitrarily small positive field can instead select the positive spontaneous magnetization if size is taken to infinity first:The distinction is the order of limits and the stated symmetric boundary conditions; plus boundary conditions already break spin inversion symmetry. A zero-field symmetric mixture of the two ordered states can retain zero mean spin magnetization even in the thermodynamic limit. Also, $⟨|m|⟩>0$ in a finite system is not by itself proof of spontaneous symmetry breaking.^[ref]

A standard sufficient tempering condition is this tail bound beyond a fixed radius, with $ϵ>0$ in spatial dimension $d$. It makes the distant interaction tail integrable; finite-range pair potentials satisfy it. Together with a stable pair potential and suitable growing containers, such assumptions support the conventional thermodynamic limit. A Curie–Weiss model with its size-dependent coupling uses a different prescription.

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

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

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

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

Two local Hamiltonian operators describe the same gapped phase when an allowed continuous path $H(s)$, $0\le s\le 1$, connects them while retaining a nonzero bulk spectral gap. The path must keep locality bounds under control. For a family of finite systems with an isolated low-energy sector, a standard sufficient condition is a lower boundindependent of system size $L$ and path parameter $s$. A positive gap in each small finite system alone does not establish a gapped thermodynamic limit.

For sufficiently regular paths of local interactions, quasi-adiabatic continuation transports the low-energy subspace by a quasi-local evolution. Lieb-Robinson bounds control its spatial tails. The resulting equivalence concerns locality-preserving transport, not arbitrary unitary operators on a finite Hilbert space. ^[ref]

The permitted paths define what is being classified. In a symmetry-protected topological phase, they must preserve the specified symmetry. Allowing symmetry-breaking terms can join phases that were distinct under that restriction; the Haldane phase is an important example of why the symmetry condition matters.

For an ideal Bose gas in a two-dimensional isotropic harmonic trap, the leading density of states measured above its ground state is $D(E)=E/(\hslash \omega {)}^2$. At saturation of the Bose-Einstein distribution, $N_{\mathrm{ex}}=(k_{B}T/(\hslash \omega ){)}^2{\int }_0^{\infty }u/(e^u-1)du=({\pi }^2/6)(k_{B}T/(\hslash \omega ){)}^2$. The finite excited-state capacity gives the displayed critical temperature and condensate fraction $N_0/N=1-(T/T_c{)}^2$ below it. The appropriate trapped-gas thermodynamic limit fixes $N{\omega }^2$ as $N\to \infty$ and $\omega \to 0$; a finite trap exhibits finite-size rounding. The number density is inhomogeneous, and the physical chemical potential approaches the oscillator ground state energy, not an unshifted zero by convention.