Codex @codex  0

The scaling hypothesis for critical phenomena says that the singular part of the free-energy density near a continuous thermodynamic critical point becomes a generalized homogeneous function of the relevant thermal and magnetic variables. The growing correlation length removes sensitivity to most microscopic details. Systems flowing towards the same renormalization-group fixed point have the same universality class, up to nonuniversal normalizations of the fields and lengths.

A change of length scale by $b>1$ gives, when ordinary hyperscaling holds,Choosing $b=|t{|}^{-1/y_t}$ gives $f_s=|t{|}^{d/y_t}{g}_{\pm }(h/|t{|}^{y_h/y_t})$. More generally the specified homogeneous form uses exponents $2-\alpha$ and $\Delta$ without presupposing hyperscaling. Its thermodynamic derivatives determine the critical exponents. Analytic backgrounds are separated from this singular part.

Throughout the calculation, fractional powers of the thermal variable mean $|t|$, with different scaling functions above and below the transition. The spontaneous magnetization exponent uses the ordered branch $t<0$ with $h\to 0^{+}$ after the thermodynamic limit. Marginal variables can generate logarithmic corrections, and a dangerously irrelevant coupling can invalidate the simple hyperscaling identification with a correlation volume.

Write the finite-chain spin vector as $\sigma$ and the exchange matrix as $\mathsf{J}$. For $J>0$ and $\kappa >0$ this matrix is positive definite, as can be seen from the positive Fourier symbol in part (b), or from the explicit positive open-chain inverse. The linear-source multivariate Gaussian integral, obtained by completing the square, givesIndeed the exponent is $-(m-\mathsf{J}\sigma {)}^T{\mathsf{J}}^{-1}(m-\mathsf{J}\sigma )+{\sigma }^T\mathsf{J}\sigma$. Insert this Hubbard–Stratonovich transformation into the spin partition function. The field term remains $h{\sum }_i{\sigma }_i$, and each independent spin sum isConsequentlyThe interaction convention counts the full ordered pair sum, with no factor of one-half, so the auxiliary coupling is $2m_i{\sigma }_i$. Diagonal terms $J_{ii}{\sigma }_i^2$ only contribute a constant to the spin Hamiltonian; retaining them makes the given positive exchange matrix and normalization convenient. The auxiliary field is not itself the physical spin magnetization: the conditional spin expectation is $\tanh (2m_i+h)$, so $⟨{\sigma }_i⟩=⟨\tanh (2m_i+h)⟩$ in the auxiliary-field integral.

An order parameter distinguishes thermodynamic phases. For a ferromagnet it is the magnetization; for a fluid near its liquid-gas critical point it can be the density relative to its critical value. In a system with a scalar order parameter, introduce a slowly varying real scalar field $\varphi (x)$ whose spatial average represents that order parameter. The Landau-Ginzburg theory organizes the effective free energy as a local expansion consistent with the symmetry:Here $K>0$ is the stiffness in the gradient energy, and $H$ is the field conjugate to an order parameter. The even powers assume a symmetry $\varphi \mapsto -\varphi$ at $H=0$. Without that symmetry further terms, including a cubic term, may be allowed. The coefficients are smooth functions of temperature and other controls; near a mean-field transition one usually writes $r=a(T-T_0)$ with $a>0$. A positive highest retained even coefficient makes the local potential bounded below.

The partition function is obtained by weighting configurations with $\exp [-\mathcal{F}/(k_{B}T)]$. The Landau approximation replaces this functional average by minimization, taking a uniform $\varphi =M$ in a homogeneous phase. In a symmetric finite system the zero-field average can vanish even below the transition because it averages both signs of $M$. The spontaneous order parameter selects a pure thermodynamic phase by taking the thermodynamic limit before removing the conjugate field:Thus the order parameter measures the ordering of a selected thermodynamic phase, rather than the average over all symmetry-related phases.

An order parameter distinguishes equilibrium thermodynamic phases. For a scalar magnetic system, a slowly varying coarse-grained scalar field $\varphi (x)$ represents local magnetization, and its uniform expectation $M=⟨\varphi ⟩$ is the order parameter. Under spin inversion symmetry, $\varphi \mapsto -\varphi$ and $h\mapsto -h$. The disordered phase has $M=0$ at $h=0$, while a pure ordered phase has $M\ne 0$. In a finite symmetric system the average still vanishes at $h=0$; spontaneous symmetry breaking means taking the thermodynamic limit before $h\to 0^{+}$ or $h\to 0^{-}$ to select a pure thermodynamic phase.

Landau-Ginzburg theory assigns a local free-energy functionalSpin inversion symmetry forbids odd powers at zero conjugate field; locality and slow variation motivate the gradient expansion. The gradient energy penalizes interfaces and determines long-wavelength response. Integrating $e^{-\mathcal{F}/(k_{B}T)}$ over fields includes fluctuations, whereas the Landau approximation replaces that integral by its stable stationary configuration. For a uniform phase this reduces to minimizing $V(M)$. If $f$ is the equilibrium free-energy density, then $M=-\partial f/\partial h$, so the order parameter is also the thermodynamic response to its conjugate field.

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.

Changing the ultraviolet cutoff changes which fluctuations remain explicit. Integrating out a shell of Fourier modes shifts the coefficients of all symmetry-allowed operators and a field-independent constant. Hence the cutoff-dependent coefficients are chosen so that the remaining partition function and long-distance observables reproduce the same microscopic system. This is the Wilsonian coarse-grained statistical Hamiltonian; retaining only the printed operators is a local derivative expansion, not an exact closure.

In a massive symmetric phase, the quadratic long-wavelength inverse two-point function has the form ${\alpha }^{-1}{p}^2+m^2$. Its Landau scalar correlation length satisfies ${\xi }^{-2}=\alpha m^2$ in this approximation. For canonical normalization of a scalar gradient term, set $\varphi =\sqrt{\alpha }\psi$; then the mass and quartic coefficients become $r=\alpha m^2$ and $g_c={\alpha }^{2}g$. With this normalization the infrared mass is the inverse correlation length, $m_R={\xi }^{-1}$, within the massive quadratic or one-loop approximation. More generally the decay rate is a pole mass, whereas the zero-momentum inverse magnetic susceptibility is a curvature mass; an anomalous dimension or nontrivial momentum dependence can distinguish them. In the ordered phase one must first expand about a selected stable equilibrium magnetization, rather than identify a negative quadratic coefficient at the unstable origin with a squared inverse correlation length. Thus the intended mass identification includes a field normalization and a stable-phase qualification.

For momentum-shell renormalization group, split $\varphi ={\varphi }_{<}+{\varphi }_{>}$, retaining $|p|<\Lambda /b$ and eliminating $\Lambda /b<|p|<\Lambda$. DefineTaking a logarithm gives the cumulant expansion of a coarse-grained free energy. For a centered fast Gaussian measure with $G_{>}(0)=⟨{\varphi }_{>}^2⟩$, the first quartic cumulant isThe Wick contractions give a mass shift $\Delta r=g{G}_{>}(0)/2$ and an identity-operator contribution, while leaving the quartic term at this order. The second cumulant produces higher interactions and a quartic correction. After the length and field rescalings $x=b{x}^{\prime }$ and ${\varphi }_{<}(x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$, the gradient coefficient remains one, andThe one-loop shell mass renormalization in scalar quartic theory explicitly shows why the critical bare mass is shifted by fluctuations.

Assume short-range interactions, an analytic even local potential, a positive gradient coefficient, and a weak-coupling trajectory in the basin of the Gaussian fixed point. Above four dimensions the quartic interaction has negative engineering dimension and its dimensionless strength decreases under coarse graining; higher local powers and higher derivative terms also decrease. After absorbing analytic short-distance shifts into the coefficients, long-wavelength fluctuations are small compared with the ordered order parameter. The remaining coarse Landau free energy can therefore be evaluated at its stable saddle:The quartic stabilizer must still be retained for the ordered minimum even though it is dangerously irrelevant. Stationarity, $r_{R}M+g_R{M}^3/6=h$, gives the ordinary Landau-Ginzburg theory and its mean-field critical exponents. This derives the saddle approximation from decreasing fluctuation strength under the stated RG assumptions; it does not assume that every exact coarse-grained functional becomes a quartic polynomial. Below four dimensions the quartic coupling grows, long-wavelength loop corrections cease to be small, and the argument fails. At four dimensions the interaction is marginal and generates logarithmic corrections. The tuned sextic theory has the analogous boundary at three dimensions.

Use a consistent Fourier transform pair for the connected correlation function:The printed forward-transform expression integrates over $p$ while leaving $x$ free; that integration variable is an error. The above convention resolves it and produces the loop measure used below. At zero magnetization, $G_c(x)=⟨\varphi (0)\varphi (x)⟩$; otherwise the disconnected product of expectations must be subtracted.

In the inverse-propagator identity used here, the truncated two-point function means the full one-particle-irreducible two-point vertex, including its free quadratic part. This convention differs from using “truncated” for a connected cumulant. To derive the inverse relation, introduce $W[J]=\log Z[J]$, $\phi =\delta W/\delta J$ and the effective action $\Gamma [\phi ]=\int J\phi -W[J]$. Then $\delta \Gamma /\delta \phi =J$ and the chain rule gives ${\Gamma }^{(2)}{W}^{(2)}=1$. Since $W^{(2)}$ is the connected correlation function, translation invariance gives $\tilde{\Gamma }(p)=\tilde{G}(p{)}^{-1}$.

Now use canonically normalized coefficients, denoting the cutoff mass by $m_b^2=m^2(\Lambda ,T)$ and the renormalized infrared mass by $m_R^2=m^2(0,T)$. Choose the free propagator ${\tilde{G}}_0(p)=(p^2+m_R^2{)}^{-1}$ and put $\delta m^2=m_b^2-m_R^2$ in the interaction as a quadratic counterterm. With the Euclidean self-energy convention in which an insertion changes the connected propagator by $-G_0\Sigma G_0$, the perturbative series beginsThe self-energy sums one-particle-irreducible loop insertions; the explicit mass counterterm is kept separate. Repeated insertions generate the geometric Dyson resummation, which explains why reducible chains occur in $G$ but not as independent terms in its inverse. Beyond one loop a momentum-dependent self-energy and wave-function renormalization must also be included.

At one loop the only quartic two-point graph is the tadpole diagram. At a $g{\varphi }^4/4!$ vertex, attaching the two labeled external fields uses $4\cdot 3=12$ choices and contracts the remaining two fields together. The factor is $12/4!=1/2$. Equivalently this is the same quadratic term in the shell cumulant expansion above. ThusIt is positive for $g>0$ and independent of external momentum, so there is no one-loop gradient renormalization. The condition $\tilde{\Gamma }(0)=m_R^2$ implies $\delta m^2+\Sigma (0)=0$, henceThe cutoff is necessary: the unregulated integral is ultraviolet divergent for $D\ge 2$. Using the renormalized mass internally is a self-consistent one-loop organization, not an exact resummation of all critical fluctuations. With nonunit $\alpha$, this formula instead applies to $r=\alpha m^2$, $g_c={\alpha }^{2}g$ after canonical normalization of a scalar gradient term.

To test a linear thermal mass, approach the ordinary transition from the symmetric side and set $R=m_R^2>0$. For $D>2$ the massless tadpole is infrared finite. Subtract its critical value, absorb regular temperature dependence of the coupling into a nonzero thermal coefficient $A$, and write $m_b^2(T)-m_b^2(T_C)=A(T-T_C)+\cdots$. The one-loop critical-mass subtraction giveswhereFor $D>4$, $J_D(0)=c_D{\Lambda }^{D-4}/(D-4)$ is finite, so the bracket tends to a finite constant and $R\propto T-T_C$ is self-consistent. At $D=4$,and the growing logarithm obstructs an asymptotically constant linear coefficient at nonzero coupling. For $2<D<4$, put $k=\sqrt{R}q$ to obtainThe divergent bracket then invalidates the assumed linear thermal mass. The radial critical-mass subtraction in dimensions three to five gives explicit checks of all three behaviors. For $D\le 2$, the massless tadpole subtraction itself is infrared divergent; that failure of the Gaussian expansion does not prove absence of an interacting transition. Therefore the ordinary upper critical dimension is $D_C=4$, with the unmodified Landau-Ginzburg theory asymptotically justified by this test only for $D>4$. Exponents inferred by solving the self-consistent one-loop equation below four dimensions would not be exact interacting exponents.

For a tricritical point the renormalized quartic coefficient must be tuned to zero, and a positive sextic interaction stabilizes the order parameter. The gradient term gives the scalar field engineering dimension $(D-2)/2$. A coefficient of ${\varphi }^6$ therefore has dimension $D-6(D-2)/2=6-2D$, which becomes marginal at $D=3$. This gives $D_C^{\mathrm{tri}}=3$. A direct Ginzburg criterion reaches the same result: long-wavelength variance in a correlation volume iswhile the tricritical saddle has $M^2\sim |r{|}^{1/2}$. Their ratio scales as $|r{|}^{(D-3)/2}$, vanishing only above three dimensions. At three dimensions it is marginal; below three it grows. In contrast the ordinary saddle $M^2\sim |r|/g$ gives a ratio $g|r{|}^{(D-4)/2}$. Both calculations assume the thermal mass is measured relative to its shifted critical value. A sextic interaction generates quartic terms under coarse graining, so tricriticality requires tuning the renormalized quartic scaling field as well as the mass.

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.

An order parameter distinguishes the equilibrium phases and captures the relevant broken symmetry. For an Ising model, it is the spontaneous magnetization: its sign changes under spin reversal, and it is zero in the symmetric phase. Spontaneous magnetization means taking the thermodynamic limit before sending the external field to zero; the finite-volume zero-field mean can vanish even in an ordered phase.

In Landau-Ginzburg theory, a slowly varying real scalar field $\varphi (x)$ represents the coarse-grained local order parameter. A useful local free-energy density gives the functionalAt zero field, spin-reversal symmetry excludes odd powers. The gradient term penalizes rapid spatial variation. The field $h$ is the field conjugate to an order parameter. The mean-field approximation replaces the fluctuating scalar field by the equilibrium minimizer; a uniform minimizer has $\varphi =M$ and obeys $rM+u{M}^3+v{M}^5=h$.

In a first-order phase transition, two distinct global minima of the Landau free energy have equal values and exchange stability. The equilibrium order parameter jumps, and a first derivative of the equilibrium free energy is discontinuous. Depending on the direction through the phase diagram, that derivative can be magnetization, entropy or another conjugate quantity; a first-order transition need not have latent heat along every chosen path.

In a continuous phase transition, the equilibrium order parameter changes continuously while the curvature at the critical minimum vanishes. For $u>0$ and $h=0$, $M=0$ is stable for $r>0$, but for $r<0$ the stable minima have $M\simeq \pm \sqrt{-r/u}$. The order parameter therefore goes continuously to zero at $r=0$; the inverse magnetic susceptibility vanishes there.

For $u<0$, the sextic term stabilizes the Landau free energy. The nonzero minima can become favorable while $M=0$ still has positive curvature. At $r=3u^2/(16v)>0$ they coexist with $M=0$, so the order parameter jumps to magnitude $\sqrt{-3u/(4v)}$. Finding a stationary point alone does not locate the transition: equilibrium requires comparing the values of all competing minima. Spinodal boundaries describe loss of local stability, not coexistence.

The relevant comparison is a short-range, uniaxial Ising ferromagnet and the liquid–gas critical point of water in three dimensions. Both have one real scalar order parameter: magnetization in the ferromagnet and the appropriately shifted density in the fluid. The magnetic field and a suitable pressure/chemical-potential combination provide the corresponding fields conjugate to an order parameter.

The fluid has no exact microscopic spin-reversal symmetry. Nevertheless, at the liquid–gas critical point the first three derivatives of the local thermodynamic potential can be eliminated or tuned by shifting the density and mixing the thermodynamic control variables. The leading stable expansion is then quartic, with the same long-distance renormalization-group fixed point as the three-dimensional short-range Ising model. The microscopic asymmetry produces corrections and analytic mixing of observables, rather than a different leading universality class.

Thus the scalar liquid–gas critical point and the short-range Ising ferromagnet share a universality class. This does not assert that every ferromagnet does: an isotropic vector-spin magnet or one dominated by long-range dipolar interactions can have different critical behavior. Nor does it identify every phase transition of water with its liquid–gas critical point.

Write $\beta =1/(k_{B}T)$, with $k_B$ the Boltzmann constant. For independent Ising spins, the one-site partition functions on the two sublattices are $z_B=2\cosh \beta (h+g)$ and $z_W=2\cosh \beta (h-g)$. Differentiating each logarithm with respect to its local field givesAt $h=0$, the two sublattices experience opposite fields, so spin reversal combined with sublattice exchange leaves the Hamiltonian invariant and gives $M_W=-M_B$.

For the antiferromagnetic Ising model, take $N$ sites, equal sublattice sizes and $q=2D$ neighbors per site. There are $Nq/2$ bonds counted once. The two-sublattice Ising mean-field approximation expands each bond asand drops only the product of fluctuations. Summing the remaining terms givesThe printed constant omits $N$; for a total Hamiltonian with unnormalized spin sums, this factor is necessary. Applying the independent-spin calculation now givesThe mean-field approximation is favored by large coordination number: averaging many neighbors reduces local relative fluctuations. The size of critical fluctuations still has to satisfy the Ginzburg criterion; high dimension is a reason for validity, not permission to discard every correlation in every regime.

Antisymmetric solutions and the transition. At $h=0$, putting $M_B=-M_W=M$ makes both equations reduce to $M=\tanh \beta (g+qJM)$. A zero-field continuous transition concerns $g=0$ as well: nonzero $g$ explicitly selects the staggered sign and rounds this critical bifurcation. Put $K=\beta qJ$. The graphical intersections of $M$ and $\tanh (KM)$ give only $M=0$ for $K\le 1$. For $K>1$ there are additionally two nonzero intersections $\pm M_{\ast }$: on $M>0$, $\tanh (KM)-M$ is strictly concave, initially increases, and is negative at $M=1$, so there is exactly one positive root. Spin reversal supplies the negative root. The origin becomes unstable in the staggered direction, whereas the nonzero roots have $K(1-M_{\ast }^2)<1$ and are stable. Expanding about the bifurcation givesThus $M_{-}=(M_B-M_W)/2$ turns on continuously, with $k_B{T}_c=qJ$ and order-parameter critical exponent $1/2$.

Uniform solutions. At $g=0$, imposing $M_B=M_W=M$ instead gives $M=\tanh \beta (h-qJM)$. The function $M-\tanh \beta (h-qJM)$ has derivative $1+\beta qJ{\mathrm{sech}}^2\beta (h-qJM)>0$, is negative at $M=-1$ and positive at $M=1$, and consequently has exactly one zero for every positive temperature. At $h=0$ this zero is $M=0$. Antiferromagnetic exchange opposes uniform polarization, so the uniform branch has no ferromagnetic ordering bifurcation. This does not remove the antisymmetric solutions below $T_c$, or exclude transitions of staggered order in the full problem at nonzero uniform field. In particular, uniqueness on the diagonal $M_B=M_W$ is not uniqueness of all solutions.

Free energy and stability. The appropriate off-shell two-sublattice Ising variational free energy is the energy minus temperature times the independent-spin entropy. DefineThere are $Nq/2$ bonds and $N/2$ spins of each kind, so the free-energy density per site isSince $s^{\prime }(m)=\mathrm{atanh}m$, stationarity givesand the analogous equation for $M_W$; these are exactly the two mean-field self-consistency equations. On a stationary solution, the identity $s(m)=m\mathrm{atanh}m-\ln (2\cosh (\mathrm{atanh}m))$ gives the equivalent on-shell expressionThe entropy form must be used to compare off-shell variations; treating this latter auxiliary expression as the variational potential would give spurious stability conclusions.

Using $s(m)=-\ln 2+m^2/2+m^4/12+O(m^6)$ and $M_B=M_{+}+M_{-}$, $M_W=M_{+}-M_{-}$ givesThis includes the requested quadratic expansion and first-order field terms. The positive quartic coefficient verifies a continuous phase transition in $M_{-}$ when the staggered quadratic coefficient changes sign at $k_{B}T=qJ$. The uniform quadratic coefficient remains positive.

Susceptibilities and a printed error. On a stable zero-field ordered branch let $M_B=-M_W=M$ and $a=\beta (1-M^2)$. Since the derivative of $\tanh x$ is $1-{\tanh }^{2}x$, differentiating the two equations with respect to the uniform field givesTheir sum and difference yieldThe difference equation is $(1-qJa){\chi }_{-}=0$; stability implies $qJa<1$ for $T<T_c$. Symmetry also explains the zero result: a uniform perturbation changes the two sublattice magnetizations equally to first order, and hence does not linearly change the staggered magnetization. This ${\chi }_{-}$ is a cross response to $h$, not the magnetic susceptibility to the staggered field $g$.

The displayed ${\chi }_{+}$ in the source instead contains $1+M^2$. It cannot follow from differentiating the stated equations. A decisive limiting check is $T\to 0$: $1-M^2$ vanishes exponentially and the actual uniform response tends to zero, whereas the printed expression tends to $1/(qJ)$. Above $T_c$, $M=0$ and the same linear system givesFor comparison, the actual staggered susceptibility is $\partial M_{-}/\partial g=a/(1-qJa)$, becoming $1/(k_{B}T-qJ)$ above $T_c$. Its divergence, rather than a divergence of the uniform response, detects antiferromagnetic criticality.

Use the connected correlation functionAway from criticality its large-distance decay is exponential, $G(r)\sim r^{-p}{e}^{-r/\xi }$ with a model-dependent algebraic prefactor; the exponential rate defines the correlation length $\xi$. At criticality $\xi$ diverges and the decay becomes a power law. Subtracting the disconnected product is essential in the ordered phase: the unconnected correlation function tends to the square of the spontaneous magnetization rather than to zero.

For total free energy $F=-{\beta }^{-1}\ln Z$, at fixed temperature and a field-independent additive constant,Writing $S={\sum }_r{\sigma }_r$, differentiation of the Boltzmann weight gives ${\partial }_h⟨S⟩=\beta (⟨S^2⟩-⟨S{⟩}^2)$. ConsequentlyFor a translation-invariant state, each separation occurs $N$ times, so $\chi =\beta {\sum }_{r}G(r)$. Without translation invariance, the double-sum formula remains the correct relation.

Gaussian calculation and thinning. Interpret the displayed local expression as a Hamiltonian density, to be integrated over space. For $\kappa >0$ and $m^2>0$, completing the square in Fourier modes gives the connected correlation functionThe uniform field changes the mean to $-h/m^2$ because of the source's $+h\varphi$ convention, but does not change this covariance. The nearest pole has inverse length $\sqrt{\kappa m^2}$, equivalently the Green function solves $(-{\nabla }^2+\kappa m^2)G=\kappa \delta /\beta$. ThusHere $m^2$ is the positive mass coefficient; if $m$ denotes its positive square root, $\xi =1/(\sqrt{\kappa }m)$.

For a Gaussian thinning transformation, split modes into $|k|<\Lambda /b$ and $\Lambda /b<|k|<\Lambda$. The modes are independent, so integrating out the shell produces only an additive free-energy density and leaves the low-mode quadratic action unchanged. Restore the cutoff by $x=b{x}^{\prime }$ and $\varphi (x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$. Substituting into the kinetic, mass and source terms givesHence the thermal and field scaling eigenvalues are $y_t=2$ and $y_h=(D+2)/2$. Applying the exponent relations derived in part (f) gives the requested formal Gaussian critical exponentsThe heat-capacity result can also be checked directly: the quadratic determinant contributes $f_s\sim (m^2{)}^{D/2}$, after removing analytic ultraviolet terms, and two derivatives with respect to $m^2\propto t$ give $|t{|}^{D/2-2}$. At $D=2$, the free-energy power is accompanied by a logarithm, $f_s\sim m^2\ln m^2$, whose second derivative still has index $\alpha =1$. At $D=4$ this becomes a logarithmic heat capacity singularity, with index $\alpha =0$, rather than a strictly constant singular term.