Weiss mean field theory for the Ising model

Let us now study the mean field theory of an Ising model, called Weiss mean field theory. As usual, we start from the Hamiltonian of a nearest-neighbour interaction Ising model, so that the partition function of the system is[1]:

We define our order parameter as the single-spin magnetization:
and we introduce it in the partition function through the identity:
Note that the last term is , which exactly measures the fluctuations of the spin variables from their mean value. We now neglect it, so that in the partition function we can approximate:
If our system is isotropic then , so that we can write where is the coordination number of the lattice. This means that:
and the partition function reduces to:
Therefore the free energy of the system in the thermodynamic limit is:
and the magnetization is:
We have thus found the so called self-consistency equation for the magnetization:
which is almost the same as the one we have encountered in The role of interaction range. We now proceed exactly as before. Since we are interested in studying the possible occurrence of phase transitions and in particular the existence of spontaneous magnetization we focus on the case :
which is just like the self-consistency equation that we found in The role of interaction range, but taking into account also the coordination number of the lattice. In this case we will have that the phase transition occurs when , i.e. at the temperature:

Now, in order to better understand the physical properties of our system at different temperatures let us expand the free energy around , still in the case . Using the Taylor series for and we have , so:

Therefore, the behaviour of near depends strongly on the sign of the coefficient of ; in particular:

  • for , namely for the free energy has only one minimum for : this is the paramagnetic (disordered) phase
  • for , namely the free energy has two stable minima at:

This is the magnetized (ordered) phase.

Trend of for in the two cases considered

Note that in the computations we have just made we have never imposed a particular value for the dimensionality of the system. This means that the results of this approximation should be valid also for , but we know that in one dimension the Ising model does not exhibit a phase transition. This is an expression of the fact that in the one-dimensional case mean field theory is not a good approximation (again, the dimensionality of the system is still too low).

Let us finally note a fact. As we have previously stated mean field theories are characterized by the fact that the fluctuations of the order parameter are neglected, and we can see that the self-consistency equation can be indeed obtained by neglecting the fluctuations of the single spins. In other words we can find it if we consider all the spins "freezed" to the value ; in fact if we consider a single spin we will have:

(where because we are requiring that the mean values of all the spins are equal).
  1. The factor has been introduced for later convenience (it can be reabsorbed in the definition of ).