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

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.

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$.

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.

If the molar latent heat $L_m$ is approximately constant on the chosen temperature interval, the ideal-vapor Clausius-Clapeyron approximation integrates toFor $L_m>0$, higher absolute temperature gives higher saturation pressure. A plot of $\ln (p/p_1)$ against inverse temperature $1/T$ is a straight line of slope $-L_m/R$. One reference pair $(T_1,p_1)$ fixes its intercept. On a per-unit-mass basis replace $L_m/R$ by $L/R_s$, with $R_s=R/M$ and $M$ the molar mass; using $L$ in $\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$ with the molar $R$ would give incompatible units.

If only the constant-latent heat approximation is relaxed, the same ideal-vapor model instead givesConsequently an apparent straight-line fit across a broad temperature range need not represent a single material constant. The simple exponential also approaches a finite value as $T\to \infty$; this is a property of the stipulated formula, not a prediction that liquid–gas phase coexistence extends beyond a thermodynamic critical point.

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]

Within the ideal-vapor Clausius-Clapeyron approximation, a local slope of the natural logarithm of saturation pressure against inverse temperature estimates the molar latent heat:In the stipulated vapor-pressure calculation, the slope is $-4811.16\mathrm{K}$, which recovers $40.0\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$. For variable $L_m(T)$, a two-point estimate gives an average weighted by $T^{-2}$, rather than the latent heat at an arbitrary midpoint. Differentiating noisy pressure measurements amplifies their errors. Outside the ideal-vapor regime, use independently determined molar volumes in $L_m=T\Delta v_{m}dp/dT$; a straight-looking graph alone cannot establish that the gas is ideal.

A macrostate groups microstates sharing prescribed coarse data. The grouping depends on what is retained: here retain only total energy. For the two distinguishable two-level systems, the macrostates $E=0,\Delta ,2\Delta$ have state degeneracies $1,2,1$. Although two units are not macroscopic in size, they make this loss of information visible: knowing $E=\Delta$ does not tell us which unit is excited.

At inverse temperature $\beta \Delta =\ln 3$, each unit's excited-state probability is $p=1/4$. Their independence gives microstate probabilities $9/16,3/16,3/16,1/16$ in the order above. Adding within each macrostate givesThe total expected value of energy is therefore $\Delta (6/16)+2\Delta (1/16)=\Delta /2$, matching two units each with expected value of energy $\Delta /4$. Equal state degeneracies would give the wrong answer: macrostate probabilities combine counting with the Boltzmann factor. At very high positive temperature, all four microstates become equally likely, but these three macrostate probabilities approach $1/4,1/2,1/4$.

With no pairing, the Hamiltonian is diagonal in the fermionic occupation numbers. The Pauli exclusion principle allows each momentum mode to be either empty or occupied once. Without a fixed fermion number constraint, occupying a mode lowers the energy precisely when its single-particle energy is negative. Therefore, with $|0⟩$ the Fock vacuum, a ground state isChoose a fixed ordering of the fermionic creation operators; changing it alters only the overall phase. Any zero-energy modes may be occupied or empty without changing the ground-state energy.

To identify the occupied modes, complete the square:Here $\beta$ is the dispersion parameter, not inverse temperature. If ${\beta }^2<4\mu$, every mode has positive energy and the unique ground state is the Fock vacuum, with $E_{\mathrm{GS}}=0$. If ${\beta }^2=4\mu$, there are still no negative-energy modes, but $p=\pm \sqrt{\beta /2}$ have zero energy if they lie on the allowed momentum grid. Each such mode doubles the ground-state degeneracy.

For ${\beta }^2>4\mu$, defineBoth are positive, and the occupied Fermi sea consists of the two intervals $p_{-}<|p|<p_{+}$. With $p=2\pi n/L$, the requested finite-size ground-state energy is explicitlyThe factor two counts opposite momenta. Endpoints are excluded from this sum because they have zero energy; when an endpoint is allowed, its optional occupation adds ground-state degeneracy. A narrow negative-energy interval can contain no allowed finite-size momenta, in which case the sum is empty. The zero-momentum mode is always empty since $ϵ(0)=\alpha \mu >0$.

Write $K=\beta J$, where $\beta =1/(k_{B}T)$ is the inverse temperature, $T$ is temperature and $k_B$ is the Boltzmann constant. In the canonical ensemble, each spin configuration has probability $e^{-\beta H}/Z_N$. The nearest-neighbor Boltzmann factors multiply, so the spin-chain transfer matrix, in the order $+1,-1$, isMultiplying its entries around the ring and summing over every Ising spin yields the partition function:The trace implements the periodic boundary conditions; replacing it by an endpoint sum would describe a different partition function. The normalized eigenvectors and corresponding eigenvalues areThe negative eigenvalue is essential for the antiferromagnetic Ising model. Introduce their ratioThe exact finite-$N$ partition function and Helmholtz free energy areFor finite $K>0$, $|t|<1$, so $Z_N>0$ for both even and odd $N$. The last term retains the finite-ring effect of geometric frustration.

To insert an Ising spin into this trace, use the diagonal matrix $\mathsf{S}=\mathrm{diag}(1,-1)$. In the eigenbasis of $\mathsf{T}$, it exchanges $|+⟩$ and $|-⟩$, and hence has vanishing diagonal entries. The one-point correlation function is thereforeThe same conclusion follows directly from spin inversion symmetry: reversing every Ising spin preserves the energy and the Boltzmann factor while reversing $S_1$.

For the two-point correlation function, the two insertions divide the ring into paths of $n$ and $N-n$ bonds. ThusIn particular, the exact finite-ring antiferromagnetic Ising correlation isBecause the one-point correlation functions vanish, this is also the connected correlation function. It gives $C_0=1$ and $C_n=C_{N-n}$, as required by $S_1^2=1$ and exchanging the two paths around the ring. At $K=0$, distinct Ising spins are independent and $C_n=0$ for $1\le n\le N-1$; $C_0=1$ is understood separately if a power $0^0$ occurs in the displayed expression.

Make the Wick rotation $t^{\prime }=-i\tau$, and write $T_E=\beta \hslash$ for the final imaginary time, where $\beta$ is inverse temperature and $\hslash$ is the reduced Planck constant. Then $d{t}^{\prime }=-id\tau$ and $dx/d{t}^{\prime }=idx/d\tau$. Thus the continued action isBoth the kinetic energy and the potential energy enter this Euclidean action with positive signs relative to their real-time definitions. Also $e^{-iH(-i\beta \hslash )/\hslash }=e^{-\beta H}$, so the Euclidean path integral isIts boundary conditions fix the two endpoint positions. They need not coincide: identifying and integrating over them would instead give the thermal trace.

The normalization of $\mathcal{D}x$ can be fixed without continuing a formal measure. Apply the Trotter product formula to $H={\hat{p}}^2/(2m)+V(\hat{x})$, divide $T_E$ into $M$ intervals of size $ϵ=T_E/M$, and integrate the intermediate positions. The resulting time-sliced configuration-space path integral isThis also explains the upper limit $\beta \hslash$ rather than $\beta$ when the reduced Planck constant has not been set to one.

Count oscillator polarizations at fixed center-of-mass momentum. In a critical closed bosonic string, each chiral sector has 24 transverse oscillator species. A species at frequency $m$ can be occupied $0,1,2,\dots$ times, so its contribution is $(1-q^m{)}^{-1}$. The full chiral generating function is thereforeThe two sectors are independent, while closed-string level matching restricts their oscillator numbers to the same $N$. Thus the number of states at the closed-string mass level is $g_N=d_N^2$, not $d_{2N}$.

Put $q=e^{-\beta }$, with $\beta >0$. The Dedekind eta function obeysThe transformed argument $2\pi i/\beta$ has large imaginary part, so its product approaches one and $\eta (2\pi i/\beta )\sim e^{-{\pi }^2/(6\beta )}$. HenceThis establishes the singular exponential as well as its power prefactor. The corrections to the transformed eta product are exponentially small in $1/\beta$.

The Cauchy coefficient formula, followed by $q=e^{-\beta }$, gives a steepest-descent integral with leading formLet $n=N-1$. The saddle of $S(\beta )=n\beta +4{\pi }^2/\beta$ isAlong the vertical steepest-descent direction $\beta ={\beta }_{\ast }+it$, the quadratic part is $-S^{\prime \prime }({\beta }_{\ast })t^2/2$. Thus the saddle-point approximation gives the Gaussian prefactorThe point $q=1$ supplies the largest exponential contribution. At another fixed root of unity of order $k>1$, the corresponding modular singularity gives a smaller exponential $e^{4\pi \sqrt{n}/k}$, so it does not change this leading asymptotic. Squaring the chiral result gives the large-level degeneracy of a closed bosonic stringReplacing $N-1$ by $N$ gives the equivalent usual leading asymptotic. The power $-27/2$ matters: the two-sector degeneracy is not just an unspecified exponential.

Since $M=2\sqrt{(N-1)/{\alpha }^{\prime }}$, the exponential is $e^{{\beta }_{H}M}$ withHere Boltzmann's constant is one. Coarse-graining the levels into a density per unit rest mass introduces the Jacobian $dN/dM={\alpha }^{\prime }M/2$. ConsequentlyThis is distinct from the degeneracy per level, whose power in $M$ is $-27$. The internal-state contribution to a canonical partition function behaves at large mass asIt converges for $\beta >{\beta }_H$ and diverges for $\beta <{\beta }_H$. At $\beta ={\beta }_H$, this particular internal rest-mass integral converges because of the power $M^{-26}$; one should not discard the prefactor and claim divergence there without specifying the full thermodynamic ensemble. For example, including continuous momenta in 25 spatial dimensions multiplies the high-mass integrand by a factor proportional to $M^{25/2}$, still leaving convergence of this ideal one-string integral at the endpoint, although sufficiently high derivatives are singular. Multi-string effects, volume and interactions affect the detailed limiting thermodynamics.

The leading microcanonical entropy grows as $S(E)\sim {\beta }_{H}E$, with logarithmic corrections, so increasingly large energies are stored in increasingly excited, long strings rather than in an arbitrarily hot ordinary gas. The Hagedorn temperature is the boundary of canonical convergence from above in inverse temperature; extrapolating the free-string canonical description to $T>T_H$ fails. This counting does not settle the interacting phase above that scale. Also, the bosonic ground-state tachyon is a separate vacuum instability; the high-level asymptotic is meaningful as a formal spectrum calculation without asserting a stable bosonic thermal vacuum.

Use ${\beta }_T=1/(k_{B}T)$ for inverse temperature, to distinguish it from the order-parameter critical exponent $\beta$. In a translation-invariant pure thermodynamic phase, the connected two-point correlation function isFor finite correlation length, at $|r|\gg \xi$ it decays exponentially, possibly with an algebraic prefactor; ${\xi }^{-1}$ is the asymptotic exponential decay rate. Below an ordered transition the pure phase is selected by taking the thermodynamic limit before $h\to 0^{+}$ or $h\to 0^{-}$. A zero-field symmetric mixture of the two ordered phases has a nondecaying connected term and does not describe the usual pure-phase magnetic susceptibility.

Let $\mathcal{F}=-k_{B}T\log Z$ denote extensive physical free energy and $S={\sum }_r{\sigma }_r$. Since the field contributes $-hS$ to the statistical Hamiltonian, differentiation givesExpanding the variance and using translation invariance proves the correlation-function susceptibility sum rule:The sum includes $r=0$. If the source is instead defined as the dimensionless quantity ${\beta }_{T}h$, the corresponding derivative omits the displayed inverse-temperature factor.

A real-space renormalization group step of length factor $b>1$ uses a local blocking kernel $B_b({\sigma }^{\prime },\sigma )\ge 0$ satisfying ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )=1$. For example, a deterministic majority-spin rule assigns one coarse Ising spin to each block; tied blocks can be resolved with equal probabilities. Define the blocked statistical Hamiltonian byThe spacing becomes $a^{\prime }=ba$. All generated interactions, including multispin terms and the field-independent constant, must be retained for an exact transformation. Summing over ${\sigma }^{\prime }$ uses normalization to prove $Z(u^{\prime },C^{\prime },N^{\prime })=Z(u,C,N)$. Iteration gives $N_p=b^{-pD}N$ and $Z(u_p,C_p,N_p)=Z(u_0,C_0,N)$.

The coarse probability law is the exact marginal law of these block variables. Indeed, for any coarse observable $A$, its expected value is the microscopic expected value of ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )A({\sigma }^{\prime })$. A blocking rule that retains the long-wavelength order-parameter mode, with source and observable transformations tracked, therefore retains long-distance response and correlation functions. Their field normalization may change, while their physical length scales and critical singularities do not. Partition-function equality alone would not justify this assertion for an arbitrary kernel that discards all slow variables. Truncating the generated operator family is likewise an approximation.

For the free-energy recursion use dimensionless free energy per site,in the thermodynamic limit. Exact partition-function invariance gives $F(u_0,C_0)=b^{-pD}F(u_p,C_p)$. If one step sends $C$ to $C^{\prime }=b^{D}C+c(u)$, the coefficient $c(u)$ records the identity-operator contribution from eliminated fluctuations and normalization. ConsequentlyRepeated substitution derives the additive free-energy recursion under blocking:In physical energy units, divide this convention by ${\beta }_T$. The additive term is essential for absolute free energies even when it does not affect the leading critical powers.

A renormalization-group fixed point satisfies $R_b(u_{\ast })=u_{\ast }$. Linearize in scaling coordinates $v_i$ so that $v_i^{\prime }=b^{y_i}{v}_i+O(v^2)$. A relevant operator has $y_i>0$, an irrelevant operator has $y_i<0$, and a marginal operator has $y_i=0$ with behaviour decided by nonlinear terms. Thus the eigenvalues of the discrete linear map are $b^{y_i}$; the sign criterion applies to the generator exponents $y_i$, not to the discrete multipliers. The critical surface is the stable manifold on which every relevant scaling field has been tuned away. Irrelevant coordinates flow to the renormalization-group fixed point along it. A repulsive renormalization-group trajectory leaves the renormalization-group fixed point along an unstable relevant direction, ultimately entering a noncritical phase or another scaling regime.

To separate the singular free-energy density from its regular background, suppose the eliminated-mode contribution $g$ is analytic near the renormalization-group fixed point. Look for an analytic background $a(u)$ satisfying $a(u)=b^{-D}a(R_{b}u)+g(u)$. Then $F_s=f-a$ satisfies $F_s(u)=b^{-D}{F}_s(R_{b}u)$. In linear thermal and field coordinates, an analytic monomial $g_{mn}{t}^m{h}^n$ has background coefficientprovided the denominator is nonzero. This explicitly explains when analytic subtraction of an inhomogeneous renormalization recursion works. The printed inhomogeneous equation may describe the free-energy function before this subtraction, or a singular part defined only up to analytic terms. The inhomogeneous sum need not be small compared with the full free energy: an analytic constant can dominate a vanishing singular contribution. It can be omitted in deriving leading singular powers after regular subtraction, provided no renormalization-group free-energy resonance, marginal logarithm or dangerously irrelevant coupling obstructs the homogeneous critical asymptotics. At a renormalization-group free-energy resonance $m{\lambda }_t+n{\lambda }_h=D$, additive logarithms can remain, so simply dropping the term would miss them.

Suppose there are exactly two relevant scaling fields, thermal $t$ and magnetic $h$, with positive exponents ${\lambda }_t,{\lambda }_h$, and the remaining fields approach a regular limit. Physical reduced temperature and magnetic field agree with these coordinates to leading order, up to nonuniversal metric factors. The homogeneous scaling hypothesis for critical phenomena isas an asymptotic scaling equation. Choose $\ell =|t{|}^{-1/{\lambda }_t}$ to make the thermal argument $+1$ or $-1$. Defining $f_{\pm }(z)=F_s(\pm 1,z)$ givesThe two functions describe respectively $t>0$ and $t<0$: approaching the transition from disordered and ordered sides need not give the same function. In particular the ordered zero-field function has one-sided derivatives associated with spontaneous magnetization. Metric factors multiply this formula in different microscopic normalizations.

The correlation length in Bravais lattice units obeys $\xi (t,h)=\ell \xi ({\ell }^{{\lambda }_t}t,{\ell }^{{\lambda }_h}h)$. At $h=0$, the same choice of $\ell$ gives $\nu =1/{\lambda }_t$. Two thermal derivatives of $F_s(t,0)$ give $2-\alpha =D/{\lambda }_t$. One and two field derivatives give, respectively, $\beta =(D-{\lambda }_h)/{\lambda }_t$ and $\gamma =(2{\lambda }_h-D)/{\lambda }_t$, assuming the relevant branch amplitudes are nonzero. At $t=0$, choose $\ell =|h{|}^{-1/{\lambda }_h}$; then $M\propto \mathrm{sgn}(h)|h{|}^{(D-{\lambda }_h)/{\lambda }_h}$, so $\delta ={\lambda }_h/(D-{\lambda }_h)$. HenceSimple algebra now gives the Rushbrooke scaling relation $\alpha +2\beta +\gamma =2$ and the hyperscaling relation $\alpha =2-D\nu$. Also $\beta \delta ={\lambda }_h/{\lambda }_t=\beta +\gamma$, giving the Widom scaling relation. These statements refer to leading singular contributions with the homogeneous two-field scaling assumptions. Above the upper critical dimension, the quartic coupling is dangerously irrelevant because it is still needed to stabilize the ordered phase; ordinary mean-field critical exponents then satisfy Rushbrooke and Widom but generally not hyperscaling. Marginal logarithms also require qualifications at the upper critical dimension.

For the connected-correlation scaling form $G(r)=|r{|}^{-(D-2+\eta )}{f}_G(|r|/\xi )$, the function $f_G(x)$ tends to a finite nonzero value as $x\to 0$ and decays exponentially at large $x$, up to powers of $x$, in a massive short-range phase. The scaling ansatz alone does not fix that prefactor. If the usual massive Gaussian field correlation tail $G(r)\propto r^{-(D-1)/2}{e}^{-r/\xi }$ applies, matching the powers gives $f_G(x)\sim B{x}^{(D-3)/2+\eta }{e}^{-x}$, with a correlation-length-dependent amplitude in $G$ fixed by matching the critical scaling form.

In physical coordinates the lattice sum has site-density factor $a^{-D}$. Its long-distance singular part therefore giveswhere ${\Omega }_{D-1}$ is the area of the unit sphere in $D$ dimensions. For $\eta <2$, finite $f_G(0)$, an integrable exponential tail and a nonzero integral, the last factor tends to a finite nonzero constant as $\xi \to \infty$. Microscopic distances contribute a regular background. Thus ${\chi }_s\propto {\xi }^{2-\eta }\propto |t{|}^{-(2-\eta )\nu }$, proving the Fisher scaling relationEquivalently the field has scaling dimension $(D-2+\eta )/2$, so its conjugate field has ${\lambda }_h=(D+2-\eta )/2$, which gives the same identity from the exponent formulas above.

A partition function sums the Boltzmann weights of all states: $Z={\sum }_s{e}^{-{\beta }_{T}H(s)}$ for discrete states, or the corresponding integral for continuous fields. It normalizes equilibrium probabilities. Its physical Helmholtz free energy is $\mathcal{F}=-k_{B}T\log Z$, and its free-energy density is $\mathcal{F}/V$. With a momentum ultraviolet cutoff ${\Lambda }_0$, a scalar-field example isThe normalization of the functional measure and field-independent energy terms must be tracked when discussing absolute free energy.

Choose a smaller cutoff $\Lambda$ and split the field into retained modes ${\varphi }_{<}$ with $|q|<\Lambda$ and eliminated modes ${\varphi }_{>}$ with $\Lambda <|q|<{\Lambda }_0$. Define the Wilsonian coarse-grained statistical Hamiltonian bySubstitution into the partition function proves $Z=\int \mathcal{D}{\varphi }_{<}{e}^{-{\beta }_T{H}_{\Lambda }[{\varphi }_{<}]}$ for every intermediate cutoff. A source retained in this definition similarly preserves long-wavelength connected correlation functions and response. Thus the coefficients in a Hamiltonian density must run with $\Lambda$: the dependence induced by eliminated modes compensates the dependence of the remaining field integral. Holding the original bare coefficients fixed while changing the cutoff would generally change the system. Exact integration generates all interactions allowed by symmetry, including higher powers, derivative terms and an additive constant; a local derivative expansion or a truncation to a few couplings is an extra approximation.

Cutoff independence makes arbitrary intermediate coarse-graining scales unphysical. It allows calculations at the scale of the observable, avoids repeatedly resolving irrelevant short-distance details and exposes the renormalization-group fixed points responsible for universality. A microscopic cutoff, such as a Bravais lattice spacing, can still label a physical microscopic model; taking it to infinity as a continuum limit requires tuning appropriate bare couplings. This is distinct from lowering an intermediate Wilsonian cutoff in a fixed model.

To interpret the zero-cutoff limit, define $M=V^{-1}\int \varphi (x)d^{D}x$, the spatially uniform mode of the order parameter. Keep this mode fixed and integrate the others, rather than integrating over $M$ too. Write $\varphi =M+\eta$, with $\int \eta =0$, and define the constrained order-parameter free energy byIn a finite periodic box, once $\Lambda$ is below the smallest nonzero momentum, all nonconstant modes have been eliminated. The remaining effective energy is $H_{\Lambda }[M]=V{U}_{\Lambda }(M)$, so its defining integral is exactly $Z_M$. Therefore, with consistent normalization,When $\mathcal{H}(\Lambda ,M)$ denotes this uniform effective Hamiltonian density, this is the stated identification of $F(M)$ with its zero-cutoff limit. It requires a free-energy density constrained at fixed $M$, not the unconstrained extensive free energy; all integrated fluctuations and additive constants are included, the zero mode in field theory is retained, and a density for the uniform configuration is well defined. The finite-volume argument avoids an interchange of the thermodynamic and cutoff limits without justification. On a homogeneous pure-phase branch one may use a local potential description, while a thermodynamic limit allowing phase separation can convexify the constrained potential. A bare nonconvex Landau polynomial is accordingly a homogeneous-branch approximation, not a claim that the exact equilibrium Legendre potential stays nonconvex.

In a field $h$, the remaining integral is $Z(h)=\int dM{e}^{-{\beta }_{T}V[F_c(M)-hM]}$, up to the consistently chosen zero-mode measure. The thermodynamic saddle then gives $f(h)={\inf }_M\{F_c(M)-hM\}$ and selects the equilibrium magnetization. Landau theory assumes that the relevant uniform potential can be approximated by a regular analytic expansion in $M$, with symmetry-constrained coefficients, and minimizes it. At tree level the integrated fluctuations are neglected, recovering the familiar bare polynomial. With fluctuations retained, the coefficients run and can become singular near criticality; the zero-cutoff identification of constrained free energy alone does not justify mean-field critical exponents.

For an ordinary short-range scalar thermodynamic critical point, the leading interaction is quartic. Absorb inverse temperature and the gradient coefficient into field normalization and writeScale $x=b{x}^{\prime }$ and $\varphi (x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$ to leave the gradient term invariant. The quartic term changes to $g{b}^{4-D}\int {\varphi }^{\prime 4}/4!$, so its engineering scaling exponent is $4-D$. It is irrelevant at the Gaussian fixed point for $D>4$, relevant for $D<4$, and marginal for $D=4$. HenceThe loop expansion makes the breakdown concrete. A quartic-vertex correction is proportional to $g^2{I}_2(r)$, whereFor $0<D<4$, rescaling $q=\sqrt{r}k$ gives $I_2(r)\sim C_D{r}^{(D-4)/2}$, diverging as $r\to 0$. At $D=4$, radial integration giveswhich is logarithmically divergent. Above four dimensions its infrared limit is finite; the leading cutoff-dependent constant can be absorbed into a coupling renormalization. The effective expansion parameter is thus of order $g{\xi }^{4-D}$ below four, or $g\log (\Lambda \xi )$ at four, since $\xi \sim r^{-1/2}$. It eventually becomes large for any fixed bare interaction in the lower-dimensional critical limit, defeating the uncorrected mean-field expansion. The Ginzburg criterion gives the same threshold by comparing correlation-volume fluctuations with the squared mean order parameter.

Below four dimensions an interacting renormalization-group fixed point can replace the Gaussian fixed point and change critical exponents, provided a continuous transition exists. At four dimensions the interaction is marginal and produces logarithmic corrections: the leading power exponents can remain their mean-field values, so “breakdown” there means that pure mean-field power laws omit those corrections, rather than necessarily that every power exponent changes. Above four the stabilizing quartic interaction is dangerously irrelevant for the ordered phase, explaining why mean-field thermodynamic exponents need not obey naive hyperscaling with the spatial dimension.

Use inverse temperature ${\beta }_T=(k_{B}T{)}^{-1}$, reserving $\beta$ without a subscript for the critical exponent. Write $c=qJ$; each undirected bond occurs once, so the number of bonds is $Nq/2$. In the mean-field theory of the Ising model, write each Ising spin as ${\sigma }_i=m+\delta {\sigma }_i$ and neglect the product of fluctuations:The resulting independent-spin statistical Hamiltonian isThe constant corrects the double counting of interaction energy. The Ising spin sums now factorize and can all be evaluated:This is the approximate partition function at an assumed mean field; equilibrium fixes $m$ self-consistently. The corresponding Ising auxiliary mean-field free energy isIts equilibrium value gives the mean-field Helmholtz free energy in the imposed field. Away from a stationary point its parameter $m$ is an assumed field variable, not necessarily the actual mean Ising spin of that independent-spin distribution.

For the full small-$m$ expansion at fixed $h$, set $t_h=\tanh ({\beta }_{T}h)$ and $s_h={\mathrm{sech}}^2({\beta }_{T}h)$. Differentiating $\log \cosh z$ gives successive derivatives $t_h$, $s_h$, $-2t_h{s}_h$, and $2s_h(2-3s_h)$ at $z={\beta }_{T}h$. ThereforeAt zero field, spin inversion symmetry eliminates odd powers and this simplifies toThe quadratic coefficient changes sign and the quartic coefficient is positive atThus the zero-field mean-field prediction is a continuous, continuous phase transition: the stable zero spin magnetization develops two symmetry-related nonzero minima continuously below $T_c$.