We have already encountered (see The role of interaction range) the Hubbard-Stratonovich identity with the saddle point approximation to compute the partition function of a one-dimensional Ising model with long-range interactions. The same technique can also be used to formulate mean field theories for systems with short-ranged interactions, and it is one of the most useful ones.
We therefore start from the most general Hamiltonian of a two-dimensional Ising model with nearest-neighbour interactions:
whose partition function is:
The idea of Hubbard-Stratonovich mean field theory is the same as in The role of interaction range
, namely substituting the quadratic term
in the Hamiltonian with a linear one, introducing some auxiliary field over which we integrate.
In order to do so we must generalize the Hubbard-Stratonovich identity.
Let be a real symmetric matrix and an arbitrary vector. Then:
where for the sake of simplicity we have used Einstein's notation on repeated indices.
We omit the proof of this result.
and the partition function becomes:
where we have explicitly written that the trace is performed over the spin variables
Now that the terms in the exponential are linear in
the partition function is easier to compute; let us change variable defining:
Let us note that there is no trace left of the original degrees of freedom
in the partition function, and that it has the form of a functional integral
, since in the continuum limit (
and the distance between lattice sites tending to zero) the auxiliary fields become functions of the position,
Let us also note that until now we still have not done any approximation (we have just rewritten the partition function); this comes into play right now: since the partition function has the form of a functional integral we must find some approximate ways to compute it. The simplest of all possible approximations is the saddle point approximation
(see the appendix The saddle point approximation
), which essentially consists in approximating the integral with the maximum value of the integrand. In other words we approximate:
where the stationary solutions are those which satisfy:
In this case we will have:
multiplying both sides by
and summing over
, taking advantage of the fact that
and in the end renaming
Let us see the relation between the auxiliary fields
and the order parameter of the system
. We know that:
and from the saddle point approximation we have:
and plugging this into
We can note that this is a more general form of the self-consistency equation for the magnetization that we have found in Weiss mean field theory. In fact, if we set
(so that also
) and choose a nearest-neighbour interaction (i.e.
is equal to a constant
-th sites are nearest neighbours, otherwise is null), calling
the coordination number of the lattice we get exactly what we found in Weiss mean field theory for the Ising model
This means that in this approximation we deduce the same conclusions we have seen there that come from this self-consistency equation; in particular we will have that the temperature of the phase transition when
and that the critical exponents
are, respectively, 3 and 1.
We can thus already see that all the mean field theories are equivalent when applied to the same system (which is something reasonable since the approximation we make is always the same). In the following sections we will see many other mean field theories, applied also to fluids (see Mean field theories for fluids
), and we will always find the same values for the critical exponents. This is another expression of the fact that different systems, like an Ising model or a classical fluid, belong to the same universality class (even more, as we have shown in Ising model and fluids
they are actually equivalent) and thus behave similarly near a critical point.