 Two-sublattice Ising variational free energy

Here $s(m)=\frac{1+m}{2}\ln \frac{1+m}{2}+\frac{1-m}{2}\ln \frac{1-m}{2}$. This free-energy density is the energy minus temperature times entropy of independent sublattice Ising spins. Since $s^{\prime }(m)=\mathrm{atanh}m$, stationarity yields the two-sublattice Ising mean-field approximation. Writing $M_{\pm }=(M_B\pm M_W)/2$, the quadratic coefficients are $(k_{B}T+qJ)/2$ for $M_{+}^2$ and $(k_{B}T-qJ)/2$ for $M_{-}^2$. The positive quartic coefficient $k_{B}T/12$ in the staggered direction stabilizes the continuous transition. Unlike a log-cosh auxiliary potential, this variational expression can be minimized away from self-consistency.