4.4: Ising model - Weiss molecular-field theory
( \newcommand{\kernel}{\mathrm{null}\,}\)
The Landau mean-field theory is phenomenological in the sense that even within the range of its validity, it tells us nothing about the value of the critical temperature Tc and other parameters (in Equation (4.3.6), the coefficients a, b, and c), so that they have to be found from a particular “microscopic” model of the system under analysis. In this course, we would have time to discuss only the Ising model (4.2.3) for various dimensionalities d.
F=−(NJd)η2−Nhη.
This energy is plotted in Figure 4.4.1a as a function of η, for several values of h.
The plots show that at h=0, the system may be in either of two stable states, with η=±1, corresponding to two different spin directions (i.e. two different directions of magnetization), with equal energy.35 (Formally, the state with η=0 is also stationary, because at this point ∂F/∂η=0, but it is unstable, because for the ferromagnetic interaction, J>0, the second derivative ∂2F/∂η2 is always negative.)
As the external field is increased, it tilts the potential profile, and finally at the critical field,
h=hc≡2Jd,
So, this simplest mean-field theory (???) does give a (crude) description of the ferromagnetic ordering. However, this theory grossly overestimates the stability of these states with respect to thermal fluctuations. Indeed, in this theory, there is no thermally-induced randomness at all, until T becomes comparable with the height of the energy barrier separating two stable states,
ΔF≡F(η=0)−F(η=±1)=NJd,
which is proportional to the number of particles. At N→∞, this value diverges, and in this sense, the critical temperature is infinite, while numerical experiments and more refined theories of the Ising model show that actually its ferromagnetic phase is suppressed at T>Tc∼Jd – see below.
The accuracy of this theory may be dramatically improved by even an approximate account for thermally-induced randomness. In this approach (suggested in 1907 by Pierre-Ernest Weiss), called the molecular-field theory,38 random deviations of individual spin values from the lattice average,
˜sk≡sk−η, with η≡⟨sk⟩,
are allowed, but considered small, |˜sk<<η. This assumption allows us, after plugging the resulting expression sk=η+˜sk to the first term on the right-hand side of Equation (4.2.3),
Em=−J∑{k,k}(η+˜sk)(η+˜sk′)−h∑ksk≡−J∑{k,k′}[η2+η(˜sk+˜sk′)+˜sk˜sk′]−h∑ksk,
ignore the last term in the square brackets. Making the replacement (???) in the terms proportional to ˜sk, we may rewrite the result as
Em≈E′m≡(NJd)η2−hef∑ksk,
where hef is defined as the sum
hef≡h+(2Jd)η.
This sum may be interpreted as the effective external field, which takes into account (besides the genuine external field h) the effect that would be exerted on spin sk by its 2d next neighbors if they all had non-fluctuating (but possibly continuous) spin values sk′=η. Such addition to the external field,
Weiss molecular field:
hmol≡hef−h=(2Jd)η,
is called the molecular field – giving its name to the Weiss theory.
From the point of view of statistical physics, at fixed parameters of the system (including the order parameter η), the first term on the right-hand side of Equation (???) is merely a constant energy offset, and hef is just another constant, so that
E′m= const +∑kεk, with εk=−hefsk≡{−hef, for sk=+1+hef, for sk=−1.
Such separability of the energy means that in the molecular-field approximation the fluctuations of different spins are independent of each other, and their statistics may be examined individually, using the energy spectrum εk. But this is exactly the two-level system that was the subject of Problems 2.2- 2.4. Actually, its statistics is so simple that it is easier to redo this fundamental problem starting from scratch, rather than to use the results of those exercises (which would require changing notation).
Indeed, according to the Gibbs distribution (2.4.7)-(2.4.8), the equilibrium probabilities of the states sk=±1 may be found as
W±=1Ze±hef/T, with Z=exp{+hefT}+exp{−hefT}≡2coshhefT.
From here, we may readily calculate F=–TlnZ and all other thermodynamic variables, but let us immediately use Equation (???) to calculate the statistical average of sj, i.e. the order parameter:
η≡⟨sj⟩=(+1)W++(−1)W−=e+hef/Te−hef/T2cosh(hef/T)≡tanhhefT.
Now comes the punch line of the Weiss' approach: plugging this result back into Equation (???), we may write the condition of self-consistency of the molecular-field theory:
Self-consistency equation:
hef−h=2JdtanhhefT.
This is a transcendental equation, which evades an explicit analytical solution, but whose properties may be readily analyzed by plotting both its sides as functions of the same argument, so that the stationary state(s) of the system corresponds to the intersection point(s) of these plots.
First of all, let us explore the field-free case (h=0), when hef=hmol≡2dJη, so that Equation (???) is reduced to
η=tanh(2JdTη),
giving one of the patterns sketched in Figure 4.4.2, depending on the dimensionless parameter 2Jd/T.
If this parameter is small, the right-hand side of Equation (???) grows slowly with η (see the red line in Figure 4.4.2), and there is only one intersection point with the left-hand side plot, at η=0. This means that the spin system has no spontaneous magnetization; this is the so-called paramagnetic phase. However, if the parameter 2Jd/T exceeds 1, i.e. if T is decreased below the following critical value,
Critical ("Curie") temperature:
Tc=2Jd,
the right-hand side of Equation (???) grows, at small η, faster than its left-hand side, so that their plots intersect it in 3 points: η=0 and η=±η0 – see the blue line in Figure 4.4.2. It is almost evident that the former stationary point is unstable, while the two latter points are stable. (This fact may be readily verified by using Equation (???) to calculate F. Now the condition ∂F/∂η|h=0=0 returns us to Equation (???), while calculating the second derivative, for T<Tc we get ∂2F/∂η2>0 at η=±η0, and ∂2F/∂η2<0 at η=0). Thus, below Tc the system is in the ferromagnetic phase, with one of two possible directions of the average spontaneous magnetization, so that the critical (Curie39) temperature, given by Equation (???), marks the transition between the paramagnetic and ferromagnetic phases. (Since the stable minimum value of the free energy F is a continuous function of temperature at T=Tc, this phase transition is continuous.)
Now let us repeat this graphics analysis to examine how each of these phases responds to an external magnetic field h≠0. According to Equation (???), the effect of h is just a horizontal shift of the straight-line plot of its left-hand side – see Figure 4.4.3. (Note a different, here more convenient, normalization of both axes.)

In the paramagnetic case (Figure 4.4.3a) the resulting dependence hef(h) is evidently continuous, but the coupling effect (J>0) makes it steeper than it would be without spin interaction. This effect may be quantified by the calculation of the low-field susceptibility defined by Equation (4.2.9). To calculate it, let us notice that for small h, and hence small hef, the function tanh in Equation (???) is approximately equal to its argument so that Equation (???) is reduced to
hef−h=2JdThef, for |2JdThef|<<1.
Solving this equation for hef, and then using Equation (???), we get
hef=h1−2Jd/T≡h1−Tc/T.
Recalling Equation (???), we can rewrite this result for the order parameter:
η=hef−hTc=hT−Tc,
so that the low-field susceptibility
Curie-Weiss law:
χ≡∂η∂h|h=0=1T−Tc, for T>Tc.
This is the famous Curie-Weiss law, which shows that the susceptibility diverges at the approach to the Curie temperature Tc.
In the ferromagnetic case, the graphical solution (Figure 4.4.3b) of Equation (???) gives a qualitatively different result. A field increase leads, depending on the spontaneous magnetization, either to the further saturation of hmol (with the order parameter η gradually approaching 1), or, if the initial η was negative, to a jump to positive η at some critical (coercive) field hc. In contrast with the crude approximation (???), at T>0 the coercive field is smaller than that given by Equation (???), and the magnetization saturation is gradual, in a good (semi-qualitative) accordance with experiment.
To summarize, the Weiss molecular-field theory gives an approximate but realistic description of the ferromagnetic and paramagnetic phases in the Ising model, and a very simple prediction (???) of the temperature of the phase transition between them, for an arbitrary dimensionality d of the cubic lattice. It also enables calculation of other parameters of Landau's mean-field theory for this model – an easy exercise left for the reader. Moreover, the molecular-field approach allows one to obtain analytical (if approximate) results for other models of phase transitions – see, e.g., Problem 18.