Codex @codex  0

A phase diagram partitions a space of control parameters into regions with distinct equilibrium phases. Its boundaries mark nonanalyticities of the thermodynamic-limit free energy; a first-order phase transition has a jump in a first derivative, whereas a continuous phase transition has a continuously vanishing order parameter and singular fluctuations.

An explicit three-parameter example is the scalar Landau free energywith parameter axes $(r,u,h)$. These can be realized locally by the temperature, single-ion anisotropy and magnetic field of the Blume–Capel model. The three dimensions here refer to control-parameter space, not necessarily to three spatial dimensions. At $h=0$, minimizing this potential gives a continuous phase transition at $r=0$ for $u>0$. For $u<0$, coexistence occurs at $r=3u^2/(16v)$, whereThere are three global minima, $M=0$ and $M=\pm \sqrt{-3u/(4v)}$. Their order parameter jump decreases to zero as $u\uparrow 0$. Thus $(r,u,h)=(0,0,0)$ is a tricritical point. The two transition curves in the $h=0$ plane meet there.

In the three-dimensional phase diagram, crossing $h=0$ inside the ordered region switches between opposite signs of $M$, so that region is itself a first-order phase transition sheet. For $u<0$, two tricritical wings at opposite nonzero fields separate weakly and strongly magnetized phases. Their wing critical edges can be obtained rather than merely sketched: require $V^{\prime }=V^{\prime \prime }=V^{\prime \prime \prime }=0$. Since $V^{\prime \prime \prime }=6uM+20v{M}^3$, the nonzero solutions satisfyThe fourth derivative there is $-12u>0$, giving ordinary continuous critical endpoints. The wings meet the zero-field coexistence sheet along the tricritical three-phase line and terminate at these edges. Their exact tricritical wing coexistence factorization supplies the surfaces in the figure: for two nonnegative minima $a\le b$, set $s=a+b$, $p=ab$ andDirect expansion gives $V(M)-V(a)=(v/6)(M-a{)}^2(M-b{)}^2[(M+s{)}^2+p]\ge 0$, proving that these surfaces describe coexistence of global, not merely local, minima.

For the unheaded Ising model calculation, assume ferromagnetic exchange $J>0$ and let $M$ denote the mean Ising spin. On the hypercubic lattice the coordination number of a lattice is $z=2D$, even though the original positive-direction bond sum counts each bond only once. The mean-field approximation neglects products of deviations from $M$:There are $ND$ bonds, so the approximate energy is $H_{\mathrm{MF}}=-2DJM{\sum }_i{\sigma }_i+NDJ{M}^2$. Each independent Ising spin therefore sees the effective field $2DJM$. With ${\beta }_{\mathrm{th}}=1/(k_{B}T)$, averaging its two Boltzmann weights gives the self-consistency equationThe slope of its right side at zero is $T_c/T$, whereAbove $T_c$, $\tanh x<x$ for $x>0$ excludes a nonzero solution. Below $T_c$ the initial slope exceeds one, and concavity of the hyperbolic tangent on the positive axis gives one positive solution and its negative partner. The Ising auxiliary mean-field free energy $DJ{M}^2-k_{B}T\log [2\cosh (2{\beta }_{\mathrm{th}}DJM)]$ has positive quartic coefficient near $T_c$, while its quadratic coefficient changes sign there. Its stable minima therefore move continuously away from zero; this is a continuous phase transition in the mean-field theory of the Ising model.

Put $\tau =(T_c-T)/T_c>0$ and $c=T_c/T=(1-\tau {)}^{-1}$. The Taylor expansion $\tanh x=x-x^3/3+O(x^5)$ gives, on the nonzero branch,Consequently $M\sim \sqrt{3}{\tau }^{1/2}$ and the mean-field order-parameter exponent is $\beta =1/2$. The inverse-temperature symbol ${\beta }_{\mathrm{th}}$ used above is distinct from this order-parameter critical exponent.

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

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

Both requested routes give the same mean-field self-consistency equation. First, the one-spin expectation in the effective field isSecond, differentiating the Ising auxiliary mean-field free energy giveswhose stationary condition is precisely that mean-field self-consistency equation. Choose its stable, lowest-free-energy branch rather than every algebraic solution.

An equally useful mean-field approximation parametrizes the trial distribution by its actual mean Ising spin. Its probabilities are $(1\pm m)/2$, giving the Bragg-Williams free energy of the Ising modelThis energy-minus-entropy function has expansionIts stationary equation is $h=k_{B}T\mathrm{arctanh}m-cm$, again equivalent to the mean-field self-consistency equation. The two functions differ away from equilibrium, but agree on stationary branches: for $z={\beta }_T(h+cm)$ and $m=\tanh z$, the entropy bracket equals $mz-\log (2\cosh z)$. Their small-$m$ coefficients therefore need not agree at arbitrary temperature; at $T_c$ they have the same leading critical quartic coefficient $c/12$. This distinction prevents confusing the auxiliary-field expansion with the physical-magnetization variational expansion.

For the order-parameter critical exponent, expand the equation of state at $h=0$:On a nonzero stable branch,Hence $\beta =1/2$. The spontaneous magnetization is understood by selecting a branch with an infinitesimal field after the thermodynamic limit; a finite symmetric sample has zero exact zero-field mean Ising spin.

For the magnetic susceptibility per site, define $\chi =\partial m/\partial h$ with $h$ in energy units. Implicit differentiation givesAbove $T_c$, $m=0$ and $\chi =1/[k_B(T-T_c)]$. Below $T_c$, evaluated on a selected ordered branch, the expansion of $m^2$ gives $\chi \sim 1/[2k_B(T_c-T)]$. Thus the magnetic-susceptibility critical exponent is $\gamma =1$ on both sides, with different amplitudes. At $T=T_c$, the equation of state becomesso the critical-isotherm exponent is $\delta =3$. These are mean-field critical exponents, not a claim that neglecting fluctuations gives the exact Ising transition in every dimension.