Mean-field theory of the Ising model ($m=\tanh [{\beta }_T(h+qJm)]$)

The mean-field approximation neglects products of Ising spin deviations, replacing each Ising spin's neighbours by the common spin magnetization. With coordination number of a lattice $q$, the effective field is $h+qJm$. The self-consistency equation equates the assumed mean Ising spin with its independent-spin expectation in that effective field. It predicts $k_B{T}_c=qJ$ and the usual quartic mean-field critical exponents, while neglecting spatially correlated fluctuations.

A product trial distribution with actual mean Ising spin $m$ givesIts entropy is the sum of independent-spin entropies, with $0\log 0=0$. Differentiation gives $k_{B}T\mathrm{arctanh}m=qJm+h$, the self-consistency equation. The variational mean parameter is distinct from the assumed field parameter of the Ising auxiliary mean-field free energy until the equilibrium condition is imposed.

The approximate partition function obtained by decoupling each Ising bond has this free energy. The quadratic constant compensates double counting. Its stationary condition is the self-consistency equation; away from stationarity $m$ parametrizes the assumed effective field and need not equal that distribution's mean Ising spin. It agrees at stationary points with the Bragg-Williams free energy of the Ising model, although their off-equilibrium expansions differ.