Codex @codex  0

At a continuous phase transition the equilibrium order parameter tends continuously to its critical value, there is no latent heat, and the correlation length diverges. For $u>0$ at zero field, minimizing $V(M)=r{M}^2/2+u{M}^4/4$ gives $M=0$ when $r\ge 0$ and $M=\pm \sqrt{-r/u}$ when $r<0$. These minima join continuously at $r=0$. The minimized singular free-energy density is zero above and $-r^2/(4u)$ below, so its first thermal derivative is continuous, while its second derivative has a finite jump in mean-field approximation.

A first-order phase transition instead has competing minima whose free energies cross while their order parameters remain distinct. The first derivatives of the equilibrium free energy generally jump; along a temperature-driven crossing this gives latent heat. For the symmetric sextic potential with $u<0,v>0$, let $y=M^2>0$. At phase coexistence, the stationary points and equality of the two free energies requireEliminating $r$ gives $-u{y}^2/4-v{y}^3/3=0$, henceAt coexistence $V(M)=v{M}^2(M^2+3u/(4v){)}^2/6$, proving that all three degenerate minima are global. The equilibrium order parameter jumps from zero to one of the two nonzero values even though $r_{\mathrm{coex}}>0$. The disordered minimum loses local stability at $r=0$, and the nonzero stationary points first appear at $r=u^2/(4v)$. These spinodal points delimit metastability and are not the equilibrium transition line. Cubic terms, when symmetry permits them, provide another common route to a first-order transition.

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.

Universality classes identify systems with the same long-distance renormalization-group fixed point. Spatial dimension, order parameter symmetry and number of components, and the range of interactions help determine the equilibrium class. Microscopic lattice details and many local interactions become irrelevant operators under coarse-graining. Universal quantities include critical exponents, scaling functions after fixing metric normalizations, and suitable amplitude ratios. The critical temperature, bare couplings, absolute correlation-length and susceptibility amplitudes, analytic free-energy density backgrounds, and microscopic length scales are generally not universal. Dynamic universality can additionally depend on conservation laws, so it should not be inferred solely from a static scalar potential.

Set $h=0$ in the even sextic Landau free energy. For $u>0$, near $r=0$ the nonzero stationary branch satisfies $M^2=-r/u+O(r^2)$ for $r<0$. The minimum value is $V_{\min }=-r^2/(4u)+O(r^3)$ below the transition and zero above. The order parameter is continuous, and if $r$ crosses zero linearly with temperature, the first temperature derivative of the minimized free-energy density is continuous too. There is no latent heat; a second derivative has a finite jump at this level.

For $u<0$, solve the two conditions $V^{\prime }(M)=0$ and $V(M)=V(0)$ with $M\ne 0$. Setting $q=M^2$ gives $r+uq+v{q}^2=0$ and $r/2+uq/4+v{q}^2/6=0$. Eliminating $r$ givesBoth the disordered and ordered minima have positive curvature there, but their values coincide. Crossing this phase coexistence boundary changes the global minimum discontinuously. The order parameter jumps; for $r=r(T)$ at fixed $u,v$, ${\partial }_T{V}_{\min }=r^{\prime }(T)M^2/2$ jumps as well, generically yielding a nonzero latent heat. Local extrema can persist past coexistence as metastable states, which explains hysteresis in a dynamically delayed transition. The spinodal points, where local curvature vanishes, do not locate the equilibrium transition.

The scaling hypothesis for critical phenomena says that the singular equilibrium free-energy density is a generalized homogeneous function of the thermal and field variables, after analytic backgrounds are removed. To derive the ordinary mean-field form, put $r=r_{t}t$, $u>0$ and neglect the higher powers near the transition. Rescale $m=\sqrt{r_t/u}|t{|}^{1/2}\psi$. Minimizing the quartic potential becomesThis is the mean-field scalar free-energy scaling form with $a=r_t^2/u$, $b=\sqrt{u}/r_t^{3/2}$. The two functions correspond to $t>0$ and $t<0$; the latter has two pure-phase branches and a cusp at zero field. The question's $>$ and $<$ subscripts distinguish these temperatures above and below $T_c$.

Differentiating with respect to $h$ gives $m=-{\partial }_h{A}_s\propto |t{|}^{1/2}$ on an ordered pure branch, and differentiating again gives $\chi =-{\partial }_h^2{A}_s\propto |t{|}^{-1}$. Thus the displayed scaling form independently yields $\beta =1/2$ and $\gamma =1$.

To allow anomalous powers, replace it byThen $\beta =2-\alpha -\Delta$ and $\gamma =2\Delta -(2-\alpha )$. Taking $t\to 0$ at fixed small $h$ gives $A_s\propto |h{|}^{(2-\alpha )/\Delta }$, so $1/\delta =(2-\alpha -\Delta )/\Delta =\beta /\Delta$. Eliminating $\Delta$ proves the Rushbrooke scaling relation and Widom scaling relation:These power-law relations assume the simple scaling form; marginal corrections or dangerously irrelevant couplings can require additional scaling variables.

Use the connected correlation function $G(r)=⟨{\sigma }_0{\sigma }_r⟩-⟨{\sigma }_0⟩⟨{\sigma }_r⟩$. Away from criticality its large-distance envelope decays exponentially on the correlation length $\xi$, often with an algebraic prefactor. At criticality $\xi$ diverges and a short-range scalar model instead has the power law $G(r)\sim r^{-(D-2+\eta )}$, with anomalous dimension $\eta$.

A blocking kernel $K_b({\sigma }^{\prime },\sigma )$ assigns coarse spins to each block of side $b$ lattice spacings, with ${\sum }_{{\sigma }^{\prime }}{K}_b({\sigma }^{\prime },\sigma )=1$. Define the coarse Hamiltonian byThe operator basis must be sufficiently complete to include interactions generated by blocking. Summing over coarse spins preserves the partition function, while the physical lattice spacing and site count change to $a^{\prime }=ba$, $N^{\prime }=b^{-D}N$. After $p$ steps they are $a_p=b^{p}a$ and $N_p=b^{-pD}N$. Exact blocking preserves the long-distance predictions; truncating the generated interaction space introduces an approximation.

Define the dimensionless free energy per site $F(u,C)=-N^{-1}\log Z(u,C,N)=f(u)+\beta C$. Preservation of $Z$ gives $F(u_0,C_0)=b^{-pD}F(u_p,C_p)$. For one step, removing the explicitly additive constant givesThe function $g$ records the additive free-energy density generated by eliminated degrees of freedom; it is usually regular near a finite blocking step's critical point. Iteration provesThis is the additive free-energy recursion under blocking. After removing regular backgrounds, its singular part obeys homogeneous scaling.

A renormalization-group fixed point satisfies $R_b(u^{\ast })=u^{\ast }$. Linearizing the map gives scaling coordinates $v_i^{\prime }=b^{{\lambda }_i}{v}_i$. A relevant operator has ${\lambda }_i>0$, an irrelevant operator has ${\lambda }_i<0$, and a marginal operator has ${\lambda }_i=0$ and needs nonlinear analysis. The critical surface consists of points flowing into the critical fixed point after all relevant directions are tuned. An outgoing or repulsive trajectory lies in the unstable manifold and describes perturbations that grow under coarse-graining.

With relevant thermal and magnetic coordinates $t,h$, the singular free-energy density obeys $F_s(t,h)=b^{-D}{F}_s(b^{{\lambda }_t}t,b^{{\lambda }_h}h)$. Choose $b=|t{|}^{-1/{\lambda }_t}$ to obtainHere ${\lambda }_t,{\lambda }_h$ are the positive scaling eigenvalues, and $f_{\pm }$ refer to the signs of $t$. Since $\xi (t,h)=b\xi (b^{{\lambda }_t}t,b^{{\lambda }_h}h)$, the correlation-length exponent is $\nu =1/{\lambda }_t$. Two temperature derivatives of $F_s(t,0)$ give $\alpha =2-D/{\lambda }_t$, establishing the hyperscaling relationThe magnetization exponent similarly is $\beta =(D-{\lambda }_h)/{\lambda }_t$. These expressions assume no additional dangerous scaling variables.

For the Gaussian field theory, decompose the field into low and high Fourier modes, integrate out $\Lambda /b<|p|<\Lambda$, and restore the cutoff using $x=b{x}^{\prime }$ and $\varphi (x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$. The gradient coefficient remains fixed, the mass coefficient transforms as $m^{\prime 2}=b^2{m}^2$, and the uniform source as $h^{\prime }=b^{(D+2)/2}h$. Hence ${\lambda }_t=2$, ${\lambda }_h=(D+2)/2$, and the Gaussian momentum-shell scaling givesThe sign convention $+h\varphi$ reverses the magnetization's sign but not these exponents. For $2<D<4$, these are Gaussian scaling predictions, not the exponents of a generic interacting scalar critical point: the quartic interaction is relevant there. The free quadratic Hamiltonian also has no stable ordered phase at negative $m^2$, so the quoted $\beta$ is its formal order-parameter scaling index, not a spontaneously magnetized state of the unstabilized Gaussian integral. The relevance distinction is also explained in Tong's renormalization-group lecture notes.

In Landau-Ginzburg theory, a scalar order parameter $M$ distinguishes macroscopic phases: for an Ising model it is the magnetization per site, while near a liquid–gas endpoint it can be a shifted density. An ordered pure thermodynamic phase has a nonzero value even when the conjugate field is removed. For a symmetry-related pair, take the thermodynamic limit before $h\to 0^{\pm }$; a finite-volume symmetric ensemble instead has zero average.

Coarse-grain to a slowly varying real field and use a local, symmetry-constrained free-energy density. A useful convention isFor spin inversion symmetry, odd powers other than the externally imposed source are absent. Take $v>0$ when the quartic coefficient can be negative; for an ordinary transition $u>0$ already stabilizes the quartic truncation. The gradient energy penalizes spatial variations. In the Landau approximation, minimize this functional, usually first over constant fields. The corresponding statistical field theory also integrates fluctuations with weight $e^{-\mathcal{F}/(k_{B}T)}$. The order parameter identifies the phase; its spatial field describes both equilibrium order and its fluctuations.

A critical exponent specifies a leading singular power in a response or order parameter as the thermodynamic critical point is approached. Analytic backgrounds and the selected thermodynamic branch must be distinguished from that singular part. Write $r=at$, $a>0$, and first take $u>0$. The equation of state is $h=rM+u{M}^3$.

At zero field on the ordered branch, $M=(-r/u{)}^{1/2}$, giving the order-parameter critical exponent $\beta =1/2$. Differentiating the equation of state gives $\chi =(r+3u{M}^2{)}^{-1}$: above the transition $\chi =1/r$, below it $\chi =1/(2|r|)$, so the magnetic-susceptibility critical exponent is $\gamma =1$. At $r=0$, $M=\mathrm{sgn}(h)(|h|/u{)}^{1/3}$, giving the critical-isotherm exponent $\delta =3$.

The minimized singular free-energy density is $f_s=0$ above and $f_s=-r^2/(4u)$ below. Two temperature derivatives produce a finite specific-heat jump, corresponding to the heat-capacity critical exponent $\alpha =0$. The inverse quadratic fluctuation kernel is $K{p}^2+r+3u{M}^2$, so $\xi =\sqrt{K/(r+3u{M}^2)}$ gives $\nu =1/2$; the critical propagator proportional to $p^{-2}$ gives $\eta =0$. ThusThese mean-field critical exponents follow from minimization and linear response. More generally, linearize the renormalization group and differentiate its singular free-energy scaling form to obtain the actual fluctuation-controlled exponents.

At the tricritical point, tune $u=0$, hold $v>0$, and let $r=at$ cross zero. The equation of state becomes $h=rM+v{M}^5$. At $h=0,r<0$, $M^4=-r/v$, while the inverse susceptibility on the ordered branch is $r+5v{M}^4=4|r|$. On the disordered side it is $r$. At $r=0$, $M=\mathrm{sgn}(h)(|h|/v{)}^{1/5}$. Therefore $\beta =1/4$, $\gamma =1$ and $\delta =5$.

The minimized singular free-energy density isTwo temperature derivatives give a singular heat capacity proportional to $|t{|}^{-1/2}$ on the ordered side. Consequently the mean-field tricritical exponent calculation givesThe displayed scaling relations for critical exponents hold numerically:These are Landau theory values with both relevant coefficients tuned: along a generic path with fixed $u>0$, the asymptotic transition is ordinary rather than tricritical. The sextic upper critical dimension is three, so fluctuations can change the powers below three dimensions, and logarithmic corrections are possible at three. The vanishing disordered mean-field singular amplitude does not invalidate the ordered-side exponent.

Use a selected translationally invariant pure thermodynamic phase and define the connected lattice correlation functionAway from criticality its long-distance tail has an exponential factor $e^{-r/\xi }$, possibly multiplied by a power. Equivalently ${\xi }^{-1}=-{\lim }_{r\to \infty }{r}^{-1}\log |G(r)|$ when this limit exists; different directions can have different lengths on an anisotropic lattice. At a thermodynamic critical point the correlation length diverges and the tail becomes algebraic. Differentiate $M=N^{-1}⟨{\sum }_n{\sigma }_n⟩$ with respect to the physical field appearing in the energy Statistical Hamiltonian. With ${\beta }_{\mathrm{th}}=1/(k_{B}T)$, the correlation-function susceptibility sum rule isIn physical continuum coordinates the sum becomes $a^{-D}\int d^{D}rG(r)$. Differentiating with respect to the dimensionless source ${\beta }_{\mathrm{th}}h$ instead removes the explicit inverse-temperature factor. Below the transition, take a one-sided pure-phase limit rather than a macroscopic symmetric mixture.

A normalized blocking kernel $B_b({\sigma }^{\prime },\sigma )$ assigns a probability or a delta constraint to a coarse configuration for each microscopic one, with ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )=1$. Define the blocked Boltzmann weight bySumming over ${\sigma }^{\prime }$ proves exact equality of the partition functions. Repeated real-space renormalization group steps compose the kernels and integrate out shorter-scale variables. To retain equality one must keep all generated interactions and the field-independent constant; a finite truncation is an approximation. Long-distance observables are mapped through the corresponding coarse observables, not identified blindly with microscopic spins. In physical units,A subsequent coordinate rescaling restores the reference lattice spacing while expressing couplings in cutoff units.

Let $F(u,C)=-({\beta }_{\mathrm{th}}N{)}^{-1}\log Z=C+f(u)$ be free energy per microscopic site. Invariance of $Z$ givesThe constant transforms as $C_1=b^D{C}_0+c(u_0)$. Therefore, with $g(u)=b^{-D}c(u)$, the identity-operator contribution to renormalization-group free energy givesThe second equality follows by substitution and induction. The inhomogeneous term records eliminated-mode entropy, determinants, and normalization factors, all multiplying the identity operator; it is not an interaction with the retained spins. Strictly, this $f$ can contain regular background terms. After an analytic subtraction of an inhomogeneous renormalization recursion, the genuinely singular part satisfies the homogeneous scaling equation, except for possible logarithmic resonances.

A renormalization-group fixed point satisfies $R_b(u_{\ast })=u_{\ast }$. Linearize in scaling coordinates $s_i$: $s_i^{\prime }=b^{{\lambda }_i}{s}_i$. A relevant operator has ${\lambda }_i>0$, an irrelevant operator has ${\lambda }_i<0$, and a marginal operator has ${\lambda }_i=0$ and needs nonlinear analysis. These ${\lambda }_i$ are logarithmic scaling eigenvalues; the discrete Jacobian eigenvalues are $b^{{\lambda }_i}$, not ${\lambda }_i$. The critical surface is the stable manifold obtained by tuning all relevant coordinates to zero. A repulsive renormalization-group trajectory leaves the fixed point along a relevant direction as the coarse-graining length grows. Stable and unstable manifolds become curved away from linear order.

The left section has $h=0$, one irrelevant coordinate $w$ and a relevant thermal coordinate $t$. Its critical surface is $t=0$ in that section. With two relevant variables, the full critical surface has codimension two: the second section shows repulsion in the $(t,h)$ plane at $w=0$.

Assume an isolated fixed point with two relevant scaling fields, analytic changes from physical controls to $(t,h)$, short-range isotropic scaling, no dangerously irrelevant variable in the thermodynamic scaling function, and no marginal or resonant logarithm. Let ${\lambda }_t,{\lambda }_h>0$ denote their scaling eigenvalues. The homogeneous singular free-energy density obeysChoose a stopping scale $L=b^p\simeq |t{|}^{-1/{\lambda }_t}$ so the renormalized thermal field is order one. Irrelevant fields vanish at this scale. With $s=D/{\lambda }_t$ and $\Delta ={\lambda }_h/{\lambda }_t$,Continuous scale interpolation absorbs the immaterial integer stopping-scale convention. Nonuniversal thermal, field and free-energy metric factors can be included explicitly; the normalized functions depend only on the universality class. The $+$ and $-$ labels refer to $t>0$ and $t<0$, which flow to different sides of the critical surface.

For clarity about the additive terms, let $a(t,h)$ solve $a(t,h)=b^{-D}a(b^{{\lambda }_t}t,b^{{\lambda }_h}h)+g(t,h)$. If $g=\sum g_{mn}{t}^m{h}^n$, its analytic coefficients areprovided the denominators do not vanish. Subtracting $a$ makes the recursion homogeneous. A vanishing denominator is a resonance and generically produces logarithms, which the simple power-law hypothesis excludes. In a convention where the regular field-dependent background has been subtracted, write ${\Phi }_{\pm }=f_{\pm }+C_{\pm }$ and $I_{\pm }(t)=C_0+a(t,0)$. This yields the requested formThe $C_{\pm }$ contain branch-dependent matching contributions; only $f_{\pm }+C_{\pm }$ is physically determined. The $I_{\pm }$ are nonsingular functions analytic through the thermodynamic critical point and represent the same regular Taylor background on the two sides, though either branch may be used to describe them. For the full free energy at nonzero field, the regular background is generally $I(t,h)$, not just $I(t)$. Adding an analytic identity-operator contribution proportional to $h^2$ provides a direct counterexample to an exact unrestricted $I(t)$ formula; it changes no singular exponents. The printed expression is justified for the leading singular scaling form with that analytic field background removed. Confluent corrections from irrelevant fields are also omitted in this leading form.

Differentiate the singular scaling form, with one-sided derivatives in the ordered phase. It givesThe length transforms as $\xi (t,h)=L\xi (L^{{\lambda }_t}t,L^{{\lambda }_h}h)$, so $\nu =1/{\lambda }_t$. Thus the Rushbrooke scaling relation and hyperscaling relation follow by direct substitution:The same reasoning also gives $\gamma =\beta (\delta -1)$. Hyperscaling relation relies on the assumptions above and can fail above an upper critical dimension because of a stabilizing dangerously irrelevant coupling.

For the universal specific-heat amplitude ratio, write $F_s=A_f|a_{t}t{|}^s{\Phi }_{\pm }(0)$ at $h=0$. Two temperature derivatives multiply both branches by the same dimensional and thermal metric factors. Away from logarithmic resonances and exceptional zero amplitudes, $A_{+}/A_{-}={\Phi }_{+}(0)/{\Phi }_{-}(0)$ with the same heat-capacity sign convention on both sides. HenceThe functions and their matched branch constants are tied to the fixed-point theory; the constants are not independent arbitrary experimental parameters. If $\alpha <0$, this means the amplitudes of the subleading singular contribution after subtracting the analytic specific-heat background. A finite background ratio itself need not be universal.

The scaling hypothesis for critical phenomena says that the singular part of the equilibrium free energy is a generalized homogeneous function of its thermal and field-like controls. Smooth backgrounds must first be removed; the hypothesis does not claim that the entire free energy, including arbitrary regular terms, has a pure scaling form.

For the ordinary scalar quartic LG theory, set $r=r_{t}t$ with $r_t>0$ and $u>0$, and rescale the uniform order parameter asThe order-parameter-dependent free-energy density becomesDefine $f_{\pm }(H)$ as the minimum of the braces, with the sign of the quadratic term respectively positive or negative. ConsequentlyFor an extensive $A_s$, $a$ additionally contains the system volume. The printed double-inequality subscript labels the two temperature branches: the upper-temperature function is $f_{+}$ and the lower-temperature function is $f_{-}$. It does not classify positive and negative magnetic fields. In particular $f_{+}(0)=0$ and $f_{-}(0)=-1/4$; below the transition the field dependence has a cusp at zero, so derivatives are taken on a selected branch.

For the following derivatives, take $A_s$ to be a density, so $m$ is the magnetization density; for total free energy the derivative gives total magnetization instead. Differentiate this mean-field scalar free-energy scaling form. The magnetization is $m=-\partial A_s/\partial h=-ab|t{|}^{1/2}{f}_{\pm }^{\prime }(H)$, giving $m_0\propto (-t{)}^{1/2}$ and $\beta =1/2$. The magnetic susceptibility is $\chi =-{\partial }_h^2{A}_s=-a{b}^2|t{|}^{-1}{f}_{\pm }^{\prime \prime }(H)$, giving $\gamma =1$. At zero field the nonzero curvature amplitudes are $-f_{+}^{\prime \prime }(0)=1$ and $-f_{-}^{\prime \prime }(0^{+})=1/2$.