Codex @codex  0

Let $A(M)=V(M)+hM$ be the zero-field constrained order-parameter free energy. At coexistence in a conjugate field $h_c$, two global minima $M_1<M_2$ obeyThus the same line of slope $h_c$ is tangent to $A$ at both endpoints: the common-tangent construction for phase coexistence. Equivalently,This is the Maxwell equal-area construction on the equation-of-state curve $h=A^{\prime }(M)$: the signed areas on either side of the horizontal coexistence line cancel. Equal slope alone finds stationary phases; the integral additionally enforces equal free energy.

If the mean order parameter is constrained between the endpoints, the system can separate into two macroscopic phases. With $\lambda =(M-M_1)/(M_2-M_1)$, their bulk free-energy density is $(1-\lambda )A(M_1)+\lambda A(M_2)$. The interface contributes subextensively in the thermodynamic limit. The resulting constrained order-parameter free energy is the convex envelope of $A$, with a straight coexistence segment. This explains why a nonconvex uniform Landau free energy must be convexified when describing stable constrained equilibrium, and why a spinodal point is not the Maxwell construction boundary.

Write ${\beta }_T=(k_{B}T{)}^{-1}$ to distinguish inverse temperature from an order-parameter critical exponent. For the Ising model, define the full two-point correlation function $G(r)=⟨{\sigma }_0{\sigma }_r⟩$ and its connected correlation function $G_c(r)=G(r)-M^2$ in a translation-invariant phase. Away from a thermodynamic critical point, the large-distance envelope decays exponentially, possibly multiplied by a power of distance:This defines the correlation length in physical distance units when the limit exists. At a continuous phase transition, $\xi \to \infty$ and the critical connected correlation function typically has algebraic decay $G_c(r)\sim |r{|}^{-(D-2+\eta )}$, where $\eta$ is the anomalous dimension. Below the transition the full $G(r)$ tends to $M^2$, so the decay definition must use $G_c$.

Let $S={\sum }_r{\sigma }_r$ and $F=-{\beta }_T^{-1}\log Z$ be the extensive Helmholtz free energy. Holding the field-independent constant $C$ and other physical couplings fixed,Because the Hamiltonian contains $-hS$, ${\partial }_h\log Z={\beta }_T⟨S⟩$ and ${\partial }_h^2\log Z={\beta }_T^2(⟨S^2⟩-⟨S{⟩}^2)$. Translation invariance then proves the correlation-function susceptibility sum rule:In the ordered zero-field phase these statements refer to a selected pure thermodynamic phase; averaging equally over opposite magnetizations introduces the separate spin-mixture contribution to zero-field susceptibility.

A normalized blocking kernel $K_b({\sigma }^{\prime },\sigma )\ge 0$ satisfies ${\sum }_{{\sigma }^{\prime }}{K}_b({\sigma }^{\prime },\sigma )=1$. For example, assign one block spin to each block of $b^D$ microscopic spins using majority sign, with equal probabilities for a tie. Define the blocked Hamiltonian, retaining the entire generated operator space, byThe field-independent contribution $N{g}_E$ belongs to the identity operator. With $N^{\prime }=N{b}^{-D}$ andsumming over ${\sigma }^{\prime }$ proves $Z(u,C,N)=Z(u^{\prime },C^{\prime },N^{\prime })$ exactly. Iterating gives $N_p=N{b}^{-pD}$ and the same partition function with $(u_p,C_p,N_p)$. Observables at large scales are reproduced by appropriately blocked observables and sources. Truncating the generated interactions makes this an approximation; dropping the additive constant already destroys the equality of free energies even when normalized spin expectations remain correct.

The extensive Helmholtz free energy is unchanged: $F(u_0,C_0,N)=F(u_p,C_p,N_p)$. In contrast, its value per site changes because the number of sites changes. Define the reduced zero-constant free-energy density $f_0(u)={\beta }_{T}F(u,0,N)/N$ in the thermodynamic limit. The preceding equality yieldsand iteration proves the additive free-energy recursion under blocking:Equivalently $C_p=b^{pD}[C_0+{\sum }_{j=0}^{p-1}{b}^{-jD}{g}_E(u_j)]$. Thus the inhomogeneous source records the identity-operator contribution to renormalization-group free energy from eliminated short-distance degrees of freedom.

To express this recursion for the singular part, write $f_0=a+f_s$, where $a$ is the regular analytic background. Then the same displayed recursion holds with $f_s$ andThis is the analytic subtraction of an inhomogeneous renormalization recursion. A purely analytic source can often be absorbed into $a$, after which it does not control the leading singular powers. It is not legitimate simply to discard an analytic source in the full free energy, where it can dominate the regular background. Moreover an analytic monomial $c{t}^m{h}^n$ contributesIf $m{l}_t+n{l}_h=D$, the sum is $p$ and hence generates a logarithm when $p\sim -\log |t|/(l_t\log b)$: a renormalization-group free-energy resonance. Such a term cannot be discarded when describing that logarithmic singularity. Neglect of the inhomogeneous term in the leading power-law scaling therefore assumes that its regular part has been subtracted, that any remaining source only changes finite scaling amplitudes, and that there are no relevant resonances or marginal logarithms. Omitting an irrelevant operator further requires absence of a dangerously irrelevant coupling, whose vanishing limit can make a scaling function singular.

A renormalization-group fixed point obeys $R_b{u}_{\ast }=u_{\ast }$. Linearization gives $\delta u^{\prime }=A\delta u$; in eigenvector coordinates $w_i^{\prime }=b^{l_i}{w}_i$. Positive $l_i$ label relevant operators, negative $l_i$ label irrelevant operators, and zero $l_i$ labels a marginal operator requiring nonlinear analysis. The critical surface is the stable manifold obtained by tuning all relevant scaling fields to zero. A repulsive renormalization-group trajectory leaves the fixed point along a relevant direction as the observation length increases. Locally a typical zero-field section has $dw/d\ell =-\omega w$, $dt/d\ell =l_{t}t$, $\ell =\log b$, so the critical surface $t=0$ attracts along $w$, while trajectories with either sign of $t$ depart from it. A magnetic scaling field adds another relevant direction.

Under the stated homogeneous-scaling assumptions and after fixing irrelevant fields at their limiting values, two relevant scaling fields obey $t^{\prime }=b^{l_t}t$, $h^{\prime }=b^{l_h}h$. These fields are analytic coordinates near the renormalization-group fixed point, proportional to the physical reduced temperature and conjugate field to leading order. The singular free-energy density consequently obeysHere $F_s$ denotes the singular free energy per original site or unit volume, with regular dimensional factors absorbed; the extensive singular part is $N{F}_s$. Iterate and choose $b^p=|t{|}^{-1/l_t}$ to obtain the scaling hypothesis for critical phenomena:The sign labels the two sides of the transition: $f_{+}$ describes $t>0$, the disordered side, while $f_{-}$ describes $t<0$, the ordered side. The two functions need not coincide; $f_{-}$ has the zero-field one-sided derivative appropriate to spontaneous magnetization. The numbers $l_t,l_h$ are scaling eigenvalues in logarithmic form: the corresponding discrete RG multipliers are $b^{l_t},b^{l_h}$.

The correlation length obeys $\xi (t,h)=b\xi (b^{l_t}t,b^{l_h}h)$. The same scale choice gives $\nu =1/l_t$. Differentiating $F_s$ once and twice with respect to the conjugate field, and twice with respect to temperature, givesAt $t=0$, instead choose the scale $b^p=|h{|}^{-1/l_h}$. Then $F_s(0,h)\sim |h{|}^{D/l_h}$ and $M\sim \mathrm{sgn}(h)|h{|}^{D/l_h-1}$, so $\delta =l_h/(D-l_h)$. Thus differentiation and scale matching establishThe first is the Widom scaling relation and the second the hyperscaling relation. The hyperscaling relation requires precisely the regular homogeneous scaling used here. Above the ordinary upper critical dimension, the quartic coupling is dangerously irrelevant in the ordered phase and this derivation cannot discard it: the mean-field critical exponents $\nu =1/2$, $\alpha =0$ do not obey $\alpha =2-D\nu$ for $D>4$. At marginal dimensions logarithmic factors likewise require care beyond the displayed pure powers.

An order parameter distinguishes phases macroscopically. In an Ising model it is the local coarse-grained magnetization; the uniform mean is $M=V^{-1}\int \varphi (x)d^{D}x$. In a selected ordered phase $M\ne 0$, whereas the disordered phase has $M=0$ at zero conjugate field. The ordered signs are related by spin inversion symmetry. A finite zero-field system may average both signs to zero, so spontaneous symmetry breaking is defined using the thermodynamic limit before removing a symmetry-selecting field. For other systems the order parameter can be a vector, a density difference, or a tensor rather than a single magnetization.

Landau-Ginzburg theory makes the order parameter a spatial field and expands a local free energy consistent with the symmetries. For a short-range scalar system with spin inversion symmetry at $h=0$, a useful expression isThe linear source is the field conjugate to an order parameter. Odd terms at zero field are forbidden by spin inversion symmetry. The coefficients depend smoothly on control parameters away from singular fluctuation effects. With $u>0$ the quartic term suffices locally; if $u<0$, take $v>0$ for stability. The Landau approximation minimizes this functional; statistical field theory also integrates fluctuations about those minima.

The Maxwell construction replaces a nonconvex homogeneous free-energy density by its coexistence envelope. If $U_0(M)$ excludes the source, the coexisting values $M_1,M_2$ satisfyThus one straight line is tangent at both endpoints. Equivalently,the equal-area form of the common-tangent construction. At fixed average $\overline{M}$ between the endpoints, the lever rule gives $\overline{M}=p{M}_1+(1-p)M_2$. Such a mixture has the tangent-line bulk energy, below the energy of an unstable homogeneous state.

The two phases occupy spatial domains separated by domain walls. The gradient term gives a positive wall energy, but a macroscopic interface has area of order $L^{D-1}$ rather than bulk volume $L^D$, so its cost per volume vanishes in the thermodynamic limit. For the symmetric quartic example a scalar quartic domain wall obeys $\kappa {\varphi }^{\prime \prime }=r\varphi +u{\varphi }^3$ and interpolates between $\pm M_0$ as $\varphi =M_0\tanh (x/\ell )$, with $M_0^2=|r|/u$ and $\ell =\sqrt{2\kappa /|r|}$. Multiplication of the differential equation by ${\varphi }^{\prime }$ gives the first integral $\kappa ({\varphi }^{\prime }{)}^2/2=u({\varphi }^2-M_0^2{)}^2/4$, explaining the finite positive interface cost. A uniform unconstrained system need not form both domains: phase separation is particularly relevant when the average order parameter is fixed or boundaries enforce competing phases.

Hysteresis occurs when the system follows a metastable local minimum after another minimum has become the equilibrium one. For $r<0,u>0$, the equation of state is $h=rM+u{M}^3$. The hysteresis spinodals of a scalar quartic potential follow from $r+3u{M}^2=0$:Following local minima until these limits gives different switching fields on increasing and decreasing $h$. Equilibrium instead switches at $h=0$. Actual switching can occur before a mean-field spinodal through nucleation, and depends on time scales; hysteresis is not a property of a globally minimized equilibrium potential. The second panel below distinguishes the local branches and their limiting switches.

The partition function sums the Boltzmann weights of all statistical configurations. For discrete states $s$ of energy $E_s$, $Z={\sum }_s{e}^{-{\beta }_{\mathrm{th}}{E}_s}$; for a regulated scalar field it is a functional integral, $Z=\int \mathcal{D}\varphi e^{-{\beta }_{\mathrm{th}}H[\varphi ]}$. It normalizes probabilities and generates thermodynamic averages. Its Helmholtz free energy isA quantum version uses $Z=\mathrm{Tr}{e}^{-{\beta }_{\mathrm{th}}\hat{H}}$. In the following identification, $F(M)$ is a constrained free-energy density, so its units match those of a Hamiltonian density rather than a total extensive energy.

Fix a microscopic cutoff ${\Lambda }_0$. Split the field into retained modes ${\varphi }_{<}$ with momenta below $\Lambda$ and eliminated modes ${\varphi }_{>}$ between $\Lambda$ and ${\Lambda }_0$. Define the Wilsonian coarse-grained statistical Hamiltonian exactly byIts coefficients and additive constant depend on $\Lambda$ so that integrating the retained modes gives precisely the original $Z$. With sources and observables carried along, the same construction preserves long-wavelength responses. Thus changing this intermediate ultraviolet cutoff is a change of description, not a change of the physical system. It is advantageous because one can integrate out fluctuations irrelevant to the observable length scale, organize a local gradient expansion, and use renormalization-group flow toward simpler universal fixed-point descriptions. Exact invariance retains all generated operators; truncating them is an approximation.

The quadratic gradient coefficient generally becomes $K_{\Lambda }$, not a constant. Field renormalization sets a chosen kinetic normalization, for example ${\varphi }_R=\sqrt{K_{\Lambda }}\varphi$ when this makes the gradient coefficient one in the energy convention. Sources and observables must be transformed inversely so that their physical meaning is preserved. Under rescaled coordinates $x^{\prime }=x/b$, the critical field scaling isIts anomalous part is the anomalous dimension; it gives $G(r)\sim r^{-2x_{\varphi }}$ and the conjugate-field exponent $y_h=D-x_{\varphi }$. Different conventions for the wave-function factor require corresponding changes in its logarithmic derivative. A fixed physical magnetization must not be confused with a cutoff-dependent dimensionless renormalized field.

To justify the zero-cutoff identification of constrained free energy, work first in a finite periodic volume $V$. Define the physical spatial mean $M=V^{-1}\int \varphi (x)d^{D}x$ and constrain that mean in the microscopic integral:Retain the constant mode $M$ while eliminating every nonzero mode. When $\Lambda$ is below the smallest nonzero box momentum, the remaining energy is $H_{\Lambda }[M]=V{U}_{\Lambda }(M)$. Consequently $e^{-{\beta }_{\mathrm{th}}V{U}_{\Lambda }(M)}=Z_M$ with a consistently fixed, possibly field-independent measure normalization. This provesFor a uniform field, $U_{\Lambda }(M)$ is the Hamiltonian density $\mathcal{H}(\Lambda ,M)$ in the question. The statement assumes that the physical zero mode is retained, all other modes are integrated exactly, field-rescaling factors are undone, and density and additive-normalization conventions agree. Taking the cutoff limit in the finite box before its thermodynamic limit makes this derivation unambiguous. With an applied uniform source, $Z(h)=\int dM{e}^{-{\beta }_{\mathrm{th}}V[F_c(M)-hM]}$ up to the same measure convention; its large-volume saddle minimizes $F_c(M)-hM$.

The Landau approximation further assumes that this potential can be approximated by a regular local polynomial on the homogeneous branches, and that the remaining fluctuations are negligible for the observables sought. It is not an assertion that an exact zero-cutoff free energy is the bare polynomial. In a coexistence region the exact constrained density can be convexified by domains, while a nonconvex Landau branch potential describes local phases and their metastability.

For the breakdown criterion use a short-range scalar theory with positive quadratic gradient stiffness. At the Gaussian fixed point, invariance of $\int (\nabla \varphi {)}^2{d}^{D}x$ gives field dimension $(D-2)/2$. A coupling $g_{2n}$ multiplying ${\varphi }^{2n}$ therefore has engineering scaling exponentIt is irrelevant only above its upper critical dimension. For the ordinary critical point the leading stabilizing interaction is quartic, while at a tricritical point the renormalized quartic term is tuned away and the leading interaction is sextic. ThereforeOne can also see the failure through the Ginzburg criterion. Long-wavelength fluctuations averaged over a correlation volume obey $⟨(\delta M{)}^2{⟩}_{\xi }\sim {\xi }^{2-D}$, up to temperature and stiffness factors. The quartic saddle has $M^2\sim {\xi }^{-2}$, giving a fluctuation/mean-square ratio proportional to ${\xi }^{4-D}$. The tricritical sextic saddle has $M^2\sim {\xi }^{-1}$, giving a ratio proportional to ${\xi }^{3-D}$. Both ratios grow without bound below their respective critical dimension, invalidating the fluctuation-free Landau exponents sufficiently close to a transition.

At equality the ratios are marginal, not power-law divergent. The loop expansion resolves this boundary: the quartic vertex correction contains $\int d^{D}q/(q^2+{\xi }^{-2}{)}^2$, logarithmic at $D=4$; the sextic vertex correction contains a two-loop three-propagator integral whose overall degree is $2D-6$, logarithmic at $D=3$. Marginal running therefore generally adds logarithmic corrections at these dimensions. Below them, where the corresponding transition exists, interactions control the long-distance critical behavior. Above them, leading thermodynamic powers are mean-field, although stabilizing dangerously irrelevant couplings explain why naive hyperscaling need not hold.

An order parameter distinguishes phases by a macroscopic quantity that changes at the transition. For an Ising model ferromagnet it is the magnetization per Ising spin, $M=N^{-1}{\sum }_n⟨{\sigma }_n⟩$. At zero magnetic field the statistical Hamiltonian is invariant under spin inversion symmetry. The disordered phase preserves this symmetry, whereas a selected ordered phase has $M\ne 0$ and a spin-reversed partner with $-M$. At finite volume with symmetric boundary conditions the expected value is zero even below the transition; spontaneous magnetization means taking the thermodynamic limit before sending a symmetry-breaking field to zero.

Landau-Ginzburg theory introduces a coarse-grained local order parameter $\varphi (x)$ and a symmetry-constrained free-energy functional. For a short-range system with scalar spin inversion symmetry, a useful expansion isThe coefficients vary smoothly with microscopic control parameters at the level of the Landau approximation; normally $r=A(T-T_c)$ with $A>0$. The gradient term penalizes spatial variation. The magnetic field is conjugate to the order parameter, and odd powers at $h=0$ are excluded by spin inversion symmetry. With $u>0$, the sextic term may be omitted locally; with $u<0$, a stabilizing positive $v$ is essential. The equilibrium uniform order parameter minimizes the local potential $V(M)$, givingThis equation of state follows from minimization, rather than being imposed as an independent assumption. Vector or tensor order parameters and their allowed invariants describe other symmetry classes.

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.

An order parameter distinguishes phases by an equilibrium expectation value. For an Ising model, the local field $M(x)$ is the coarse-grained magnetization, and the uniform value is $M=N^{-1}⟨{\sum }_r{\sigma }_r⟩$. The spin inversion symmetry maps $M$ to $-M$, so at zero field the Landau free energy contains even powers. Its disordered minimum is $M=0$; ordered pure thermodynamic phases have $M=\pm M_0$. Their sign is selected by taking the thermodynamic limit before sending $h$ to zero. A symmetric mixture can have zero average without being a single disordered phase. The term $-hM$ identifies $h$ as the field conjugate to an order parameter.

A limiting argument starts with a family of microscopic models $T_N$, a specified macroscopic quantity $M_N$ and a topology or operational criterion for convergence as $N\to \infty$. It is not enough to say that a large system has many constituents: one must identify what is held fixed, what observables are compared and whether predictions converge uniformly. A thermodynamic limit, for example, normally holds density and intensive parameters fixed while size increases. Other limits remove a resolution cutoff or enlarge an observable algebra.

A limiting theory can be both deduced from microscopic models and qualitatively novel. The compatibility depends on distinguishing an exact property at the mathematical limit from the approximate collective behaviour of actual finite systems. The conceptual distinction is developed below, followed by three complementary examples rather than a claim that every form of emergence has the same mechanism.

For a ferromagnet the scalar order parameter is the spontaneous spin magnetization. It distinguishes a symmetry-invariant phase with zero mean from one of two phases with a nonzero mean of opposite signs. At zero field, a finite system with spin inversion symmetry has exactly zero mean. The spontaneous value is defined by taking the thermodynamic limit first and then sending the symmetry-breaking conjugate field to zero from one side.

A Landau-Ginzburg theory promotes the coarse magnetization to a scalar field $M(x)$ and usesAt $h=0$, spin inversion symmetry permits only even local powers. The gradient term penalizes spatial variations. The equilibrium statistical weight is $e^{-\mathcal{F}/(k_{B}T)}$; the Landau approximation chooses a stable minimum of $\mathcal{F}$, whereas the full functional integral retains fluctuations around it. In a homogeneous phase this gives the equation of state $h=rM+u{M}^3+v{M}^5$. The field $h$ is conjugate to the order parameter, with $M=-\partial f_{\mathrm{eq}}/\partial h$.

The microcanonical ensemble describes an isolated system at fixed energy, volume and particle number. The canonical ensemble describes fixed particle number and volume with heat exchange at a prescribed temperature. The grand canonical ensemble additionally allows particles to exchange with a reservoir, prescribing temperature, volume and chemical potential. Their respective equilibrium probabilities are uniform on the accessible energy shell, proportional to $e^{-\beta E_n}$, and proportional to $e^{-\beta (E_n-\mu N_n)}$. For stable short-range systems in the thermodynamic limit, concentration of relative energy and particle-number fluctuations and concavity of entropy make their bulk predictions equivalent; phase coexistence or nonadditive interactions require care.

For the microcanonical case restrict the allowed states to the fixed-energy shell of size $\Omega$, imposing $p_n=0$ outside it. Maximize the Gibbs entropy subject to $\sum p_n=1$. A Lagrange multiplier gives $-k_B(\log p_n+1)+\lambda =0$ on the shell, soFor the canonical case impose both normalization and $\sum p_n{E}_n=U$. Stationarity of $-k_B\sum p_n\log p_n-\lambda (\sum p_n-1)-k_B\beta (\sum p_n{E}_n-U)$ gives $p_n\propto e^{-\beta E_n}$. ThusThe last identification follows from $\partial S/\partial U=1/T$. Entropy is strictly concave in the probabilities, so these stationary distributions are maxima.

For the independent three-level particles, put $z=ϵ/(k_{B}T)$ and $Z_1=1+2\cosh z$. The full partition function is $Z_1^N$ because the particles occupy labeled lattice sites. Differentiation givesThe heat capacity also equals $N\mathrm{Var}(E_1)/(k_B{T}^2)$, confirming positivity. As $T\to +\infty$, the three occupations become equal:As $T\downarrow 0$, all particles occupy their unique lowest state: