Codex @codex  0

For a ferromagnetic Ising model with coordination number of a lattice $z$, its mean-field self-consistency equation is $M=\tanh (zJM/(k_{B}T))$. Put $A=T_c/T$. For $A\le 1$, the only solution is zero. For $A>1$, strict concavity of $\tanh (AM)-M$ on $M>0$ gives one positive solution and its spin-reversed partner. Expanding $\tanh (AM)=AM-(AM{)}^3/3+O(M^5)$ gives $M^2=3(A-1)/A^3+O((A-1{)}^2)$, proving a continuous onset with order-parameter critical exponent $\beta =1/2$. The Bragg-Williams free energy of the Ising model has a positive quartic coefficient at $T_c$, so these are the stable ordered minima rather than a first-order jump.

For a scalar quartic Landau free energy $r_{t}t{m}^2/2+u{m}^4/4-hm$ with $r_t,u>0$, rescaling $m=(r_t/u{)}^{1/2}|t{|}^{1/2}\psi$ gives this minimized scaling form, with $a=r_t^2/u$ and $b=\sqrt{u}/r_t^{3/2}$. The two functions minimize $\pm {\psi }^2/2+{\psi }^4/4-H\psi$. In particular $f_{+}(0)=0$ and $f_{-}(0)=-1/4$. Below the transition the equilibrium field dependence has a cusp at zero; pure-phase derivatives are one-sided. The form gives order-parameter critical exponent $\beta =1/2$ and magnetic-susceptibility critical exponent $\gamma =1$.

The uniform mode of a pure Gaussian field theory with negative squared mass has an unbounded quadratic energy. Its partition function diverges even with a finite-volume ultraviolet regulator, so it supplies no stable ordered branch. At positive squared mass the zero-source average is zero. Thus the scaling dimension of the field can yield a formal order-parameter critical exponent without establishing a negative-reduced-temperature magnetization law. A stabilizing quartic interaction is relevant below four dimensions and ordinarily changes the critical theory to the Wilson-Fisher fixed point; above four its coefficient is a dangerously irrelevant coupling for the ordered phase.

With the dimensionless magnetic field conjugate to the magnetization, $m=-\partial f_s/\partial h$. Set $y=h/|t{|}^{\Delta }$ and differentiate the homogeneous free-energy density:On the ordered branch, the one-sided value $g_{m,-}(0^{+})$ is nonzero, so the order-parameter critical exponent isTo reach the critical isotherm, take $|t|\to 0$ at fixed nonzero $h$. Suppose $g_m(y)\sim A\mathrm{sgn}(y)|y{|}^p$ at large $|y|$. Cancellation of the residual thermal factor requires $2-\alpha -\Delta -p\Delta =0$. Hence $m(0,h)\propto \mathrm{sgn}(h)|h{|}^{(2-\alpha -\Delta )/\Delta }$ andThe singular scaling function determines the nonanalytic critical-isotherm contribution; any analytic response background is subleading for the usual $\delta >1$.

Critical exponents describe leading singular powers near a continuous phase transition. Using the reduced temperature $t=(T-T_c)/T_c$, the standard definitions includeThese are the order-parameter critical exponent, magnetic-susceptibility critical exponent, heat-capacity critical exponent and correlation-length critical exponent. The critical-isotherm exponent is defined by $M(0,H)\sim \mathrm{sgn}(H)|H{|}^{1/\delta }$. The anomalous dimension $\eta$ gives the critical long-distance connected correlation function, $G_c(x)\sim |x{|}^{-(d-2+\eta )}$, when an ordinary scale-invariant critical point exists.

In the Landau approximation, minimize the Landau free energy to obtain the equation of state, differentiate it with respect to the field conjugate to an order parameter to obtain magnetic susceptibility, and differentiate the minimized free energy with respect to temperature to obtain heat capacity. These operations determine $\beta ,\delta ,\gamma ,\alpha$, respectively. Expanding the Landau-Ginzburg theory about its uniform minimum gives the quadratic fluctuation kernel $K{q}^2+V^{\prime \prime }(M)$, whose Fourier transform determines the correlation length and the anomalous dimension. At the ordinary quartic critical point this gives ${\xi }^2=K/V^{\prime \prime }(M)\propto |t{|}^{-1}$ and $\eta =0$.

For $u>0$, $r=0$ gives $H=u{M}^3$ and hence $\delta =3$; the equilibrium free energy $-r^2/(4u)$ gives a finite heat capacity discontinuity, conventionally $\alpha =0$. Together with the explicit calculation below, the ordinary mean-field critical exponents areThese describe neglect of fluctuation corrections, rather than the exact exponents of every low-dimensional Landau-Ginzburg theory.

A tricritical point is the meeting of a line of continuous phase transitions and a line of first-order phase transitions. In the symmetric scalar Landau free energy, it occurs at $r=u=0$ with $v>0$. Two independent controls are needed to tune both coefficients. The continuous transition lies at $r=0,u>0$, and the tricritical three-phase line found above has $r=3u^2/(16v),u<0$; both approach the tricritical point. At zero conjugate field the three coexisting states on this first-order line are the disordered state and the two opposite ordered states.

Setting $u=H=0$ makes the leading stabilizing interaction sextic. The nonzero minimum obeys $M^4=-r/v$, so the order-parameter critical exponent along this tuned direction is ${\beta }_{\mathrm{tri}}=1/4$. Substitution givesand two temperature derivatives give the heat-capacity critical exponent ${\alpha }_{\mathrm{tri}}=1/2$. At $r=u=0$, the equation of state is $H=v{M}^5$, giving the critical-isotherm exponent ${\delta }_{\mathrm{tri}}=5$.

For magnetic susceptibility, differentiating $H=rM+v{M}^5$ gives ${\chi }^{-1}=r+5v{M}^4$. This is $r$ above the transition and $-4r$ below, hence ${\gamma }_{\mathrm{tri}}=1$. The quadratic fluctuation kernel still has a mass proportional to $|r|$, giving ${\nu }_{\mathrm{tri}}=1/2$ and ${\eta }_{\mathrm{tri}}=0$. Thus the mean-field tricritical exponent calculation yieldsDetuning $u$ gives tricritical crossover scaling; one must specify the approach to the tricritical point when quoting these critical exponents. The sextic interaction has upper critical dimension three, so sufficiently strong long-distance fluctuations change these mean-field predictions below three dimensions.

Take the symmetric scalar Landau free energy with $u>0$, $r=at$ for a positive constant $a$, and leading equation of stateAt zero conjugate field, minimization gives $M=0$ for $t>0$ and $M=\pm \sqrt{-r/u}$ for $t<0$. Selecting the positive pure thermodynamic phase,Consequently the order-parameter critical exponent is $\boxed{\beta =1/2}$.

Differentiate the equation of state at fixed temperature, within the selected phase:Both sides diverge as $|t{|}^{-1}$, so the magnetic-susceptibility critical exponent is $\boxed{\gamma =1}$. The leading above-to-below magnetic susceptibility amplitude ratio is two. These are the ordinary mean-field critical exponents, valid under the fluctuation qualification explained above.

Write ${\beta }_T=(k_{B}T{)}^{-1}$ for inverse temperature to distinguish it from the order-parameter critical exponent. A normalized blocking kernel $B_b({\sigma }^{\prime },\sigma )$ is a nonnegative conditional distribution of coarse spins given the fine spins, normalized byFor deterministic blocking it can impose a majority rule or block average, with a Dirac delta replacing the sum constraint for continuous coarse variables. Each block occupies $b^D$ original cells, so the new lattice spacing and number of sites areThe physical volume $N{a}^D$ is unchanged. After $p$ blockings, $a_p=b^{p}a$ and $N_p=b^{-pD}N$.

Define the coarse statistical Hamiltonian by the exact Boltzmann-weight identityThe operator family must include all generated interactions for this renormalization-group transformation to be exact; restricting it to a few couplings is an approximation. The scalar $g(u)$ has energy units and is the contribution per original site that is independent of ${\sigma }^{\prime }$. Restoring the original constant term $NC$ giveshenceThis convention measures $g$ per original site; defining a constant per new site would move a factor $b^D$ into its definition.

The identity-operator contribution to renormalization-group free energy $C_p$ records the configuration-independent normalization generated by the eliminated spins. It includes their contribution to free energy, including entropy and measure factors, and need not be independent of temperature. It has no effect on normalized probabilities of the retained spins, but is indispensable for the total partition function and its temperature derivatives. Iterating its recursion givesThus $N_p{C}_p$ is the original constant $N{C}_0$ plus the accumulated free-energy contributions of the eliminated degrees of freedom.

To obtain the additive free-energy recursion under blocking, use energy-valued free energy per current lattice site,The blocking kernel normalization impliesSubstituting the recursion for $C^{\prime }$ cancels $C$, givingA second step gives $f_0(u_0)=b^{-2D}{f}_0(u_2)+g(u_0)+b^{-D}g(u_1)$, and induction provesFor dimensionless free energy $-\log Z/N$, multiply $F,f_0,g,C$ throughout by ${\beta }_T$.

The singular free-energy density must be distinguished from the exact per-site representative $f_0$. If $a(u)$ is a regular background and $f_s=f_0-a$, thenThus the same displayed inhomogeneous recursion holds with $f_s,g_s$ after subtraction. When a suitable analytic background solves $a(u)=g(u)+b^{-D}a(Ru)$, it leaves a homogeneous singular recursion. Analytic subtraction of an inhomogeneous renormalization recursion explains this separation. A renormalization-group free-energy resonance can obstruct analytic subtraction and produce a logarithm. The original function $g$ originates in the short-distance spins that were summed out, rather than in an extra external interaction.

Near a renormalization-group fixed point $u_{\ast }$, linearize the transformation:An eigenvector of $\mathsf{A}$ with eigenvalue $b^{y_i}$ defines a scaling coordinate $t_i$ with $t_i^{\prime }=b^{y_i}{t}_i$ to leading order. A relevant operator has $y_i>0$: its coupling grows towards large distances, so it must be tuned to reach the fixed point. An irrelevant operator has $y_i<0$ and dies away; a marginal operator has $y_i=0$ and requires higher-order flow analysis. An operator of scaling dimension $x_i$ has $y_i=D-x_i$.

The coarse spins also require a multiplicative normalization. In continuum notation, writewhere ${\varphi }_{<}$ is the retained field and $\eta$ is its anomalous dimension. This restores the normalization of the critical long-distance connected correlation function. Since a uniform field conjugate to an order parameter enters through $\int d^{D}xH\varphi$, it has eigenvalue exponent $y_h=D-x_{\varphi }=(D+2-\eta )/2$.

Let $t$ be the thermal relevant coordinate with exponent $y_t>0$. After regular backgrounds have been removed, the leading singular free-energy density obeysfor a large total scale factor $B=b^p$. Choosing $B=|t{|}^{-1/y_t}$ gives the scaling hypothesis for critical phenomenaTwo temperature derivatives give the heat-capacity critical exponent; one conjugate-field derivative gives the order-parameter critical exponent; two such derivatives give the magnetic-susceptibility critical exponent. Scaling the correlation length gives $\xi (t)=B\xi (B^{y_t}t)$. Consequently,At $t=0$, instead choose $B=|H{|}^{-1/y_h}$; differentiating gives $\delta =y_h/(D-y_h)$. These scaling relations for critical exponents assume that irrelevant variables can be set to their fixed-point values without singular consequences. A dangerously irrelevant coupling or logarithmic correction requires keeping the corresponding dependence.

For the Gaussian field theory, integrate out the Fourier modes in $\Lambda /b<|q|<\Lambda$ and retain ${\varphi }_{<}$. The modes are independent at zero external field, so this Gaussian thinning transformation changes the normalization by a Gaussian determinant and leaves a quadratic statistical Hamiltonian for the retained modes. Restore the cutoff withIn the integrated statistical Hamiltonian, the volume element contributes $b^D$, each gradient contributes $b^{-1}$, and each field contributes $b^{-(D-2)/2}$. Thus the coefficient ${\kappa }^{-1}$ of the gradient energy is unchanged, whileIn particular $y_t=2$, $y_h=(D+2)/2$ and $\eta =0$. The Gaussian critical exponents and the Gaussian order-parameter scaling index therefore give the requested indicesThe sign convention $+h\varphi$ corresponds to the physical conjugate field $H=-h$ and does not affect these powers.

There are physical qualifications to these Gaussian indices. The strictly quadratic Gaussian field theory is stable for positive $m^2$ and has no stable ordered phase for negative $m^2$; its purported spontaneous $\beta$ is therefore the formal Gaussian order-parameter scaling index. A stabilizing quartic term gives an actual ordered phase. Below four dimensions that term is relevant and changes the exponents; above four dimensions it is a dangerously irrelevant coupling and gives the actual mean-field $\beta =1/2$, rather than the Gaussian hyperscaling value.

The heat capacity power can also be checked directly. Apart from constants, the zero-field Gaussian free-energy density isIf $m^2\propto t$, its singular second temperature derivative is controlled by $\int d^{D}q(q^2+\kappa m^2{)}^{-2}$. For $0<D<4$, rescaling $q=\sqrt{\kappa m^2}k$ gives a power $t^{D/2-2}$, as above. At $D=4$ there is a logarithm, and for $D>4$ an analytic background dominates the heat capacity. The negative formal Gaussian $\alpha$ then describes the subleading singular piece, with logarithmic resonances at appropriate even dimensions. This is the Gaussian specific-heat infrared threshold, and it prevents confusing singular indices with the entire observable.

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.

Critical exponents describe leading singular powers as the reduced temperature $t$ or conjugate field approaches zero. With ${\beta }_m$ used for the order-parameter critical exponent to avoid confusion with inverse temperature, the definitions includeAt criticality the connected correlation function behaves as $G(r)\sim r^{-(D-2+\eta )}$. These definitions refer to the singular part and allow additional logarithms at marginal dimensions.

For an ordinary scalar Landau free energy, the equation of state is $h=rM+u{M}^3$. Below the transition $M^2=-r/u$, giving ${\beta }_m=1/2$. Differentiating the equation of state gives ${\chi }^{-1}=r+3u{M}^2$: it equals $r$ above and $-2r$ below the transition, so $\gamma =1$. The quadratic fluctuation kernel is $\kappa k^2+r$ in the symmetric phase, giving correlation length $\xi =\sqrt{\kappa /r}$ and $\nu =1/2$. At zero mass the momentum-space correlation function is proportional to $k^{-2}$, so $\eta =0$. The minimized free energy and critical equation of state give $\alpha =0$ and $\delta =3$, as calculated in the root solution.

At the tricritical point, $u=0$ and $h=rM+v{M}^5$. Then $M\sim (-r{)}^{1/4}$, and ${\chi }^{-1}=r+5v{M}^4$ is $r$ above and $-4r$ below the transition. The gradient kernel still gives $\nu =1/2$ and $\eta =0$. ThereforeThese are mean-field critical exponents. A systematic derivation incorporating fluctuations uses the renormalization group, extracts the relevant scaling eigenvalues, and differentiates the singular free-energy density, as in Question 3.

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.

Let $t,h$ be the thermal and magnetic scaling fields, with eigenvalues $b^{y_t}$ and $b^{y_h}$. After discarding irrelevant variables whose vanishing is nonsingular, the homogeneous singular free-energy density satisfiesChoose $b=|t{|}^{-1/y_t}$. ThenTwo thermal derivatives give the heat-capacity critical exponent; one and two field derivatives give the order-parameter critical exponent and magnetic-susceptibility critical exponent. The correlation-length critical exponent follows directly from the length rescaling. ThusAt $t=0$, choose $b=|h{|}^{-1/y_h}$. Differentiating $f_s(0,h)\propto |h{|}^{D/y_h}$ gives $M\propto \mathrm{sgn}(h)|h{|}^{(D-y_h)/y_h}$, so $\delta =y_h/(D-y_h)$. The order parameter has scaling dimension $x_M=D-y_h$; therefore its critical connected correlation function has power $r^{-2x_M}$. Comparing with $r^{-(D-2+\eta )}$ yieldsThese formulas express the requested critical exponents in terms of the relevant scaling exponents. Their hyperscaling assumptions must be respected: above an upper critical dimension, a dangerously irrelevant coupling can invalidate the free-energy homogeneity used here. Applying the formulas to the Gaussian thinning transformation in the root solution gives the requested Gaussian indices and their physical qualifications.

The mean-field approximation replaces each site's fluctuating surroundings by a self-consistent average. On a cubic lattice each spin has coordination number of a lattice $z=2D$, even though summing only over positive basis directions counts each bond once. Write ${\sigma }_n=M+\delta {\sigma }_n$ and drop products of deviations. Each Ising bond becomes ${\sigma }_n{\sigma }_m\simeq M{\sigma }_n+M{\sigma }_m-M^2$, so the approximate statistical Hamiltonian isThe constant corrects counting each bond twice. The resulting independent-spin partition function is $Z_{\mathrm{MF}}=e^{-{\beta }_{T}NzJ{M}^2/2}[2\cosh ({\beta }_{T}zJM){]}^N$, where ${\beta }_T=1/(k_{B}T)$. The one-spin probabilities giveEquating this average with the assumed $M$ givesHere $J>0$ is assumed; a ferromagnetic order parameter is not the appropriate choice for an antiferromagnetic coupling. Put $A=T_c/T$. For $A\le 1$, $\tanh (AM)<AM\le M$ at any $M>0$, so only $M=0$ is possible. For $A>1$, the function $\tanh (AM)-M$ has positive slope at zero, is strictly concave on $M>0$ and is negative at $M=1$. It therefore has exactly one positive zero and a spin-reversed negative zero. Expanding at the bifurcation givesThe transition is continuous, with order-parameter critical exponent $\beta =1/2$. Stability can also be checked from the Bragg-Williams free energy of the Ising model per Ising spin:Its derivative is $-zJM+k_{B}T\mathrm{arctanh}M$, giving the same self-consistency equation. Near zero,The quartic coefficient is positive at $T_c$, the stable minima approach zero continuously, and the minimized free energy has no first-derivative jump. Thus mean-field approximation predicts a second-order transition at $T_c=2DJ/k_B$. This is the mean-field Ising critical bifurcation.