As we have alreadt stated, the Ising model can also be used to describe systems different from magnets. The most important example is the correspondence that can be established between an Ising model and a fluid through a lattice gas model. We will now show how this equivalence can be defined; in order to do so, we will proceed in two steps: we will first show how an Ising model is equivalent to a lattice gas, and then show (qualitatively) that this lattice gas model is equivalent to the classical model for a fluid.
Before doing so, we briefly review the formalism used in classical statistical mechanics in order to describe fluids (see The canonical ensemble at work, for example).
Let us then consider a system of particles subjected to some generic potentials, so that its Hamiltonian can be written as:
-body potential, which generally depends only on the distances between the particles if it involves two or more of them.
Working in the grand canonical ensemble, the grand partition function will be:
is the Planck constant, and
the dimensionality of the system
). It is convenient to separate the contributes due to the kinetic and configurational terms:
we mean all the possible potentials acting on the system. However (see, again, The canonical ensemble at work
) we have:
is the thermal wavelength
The grand free energy is thus
and as usual the thermodynamic limit of
is defined as:
with the constraint:
Ising model and lattice gas[edit | edit source]
The basic idea of the lattice gas model is to describe a fluid where the particles are located on the sites of a lattice instead of continuously occupying any position in space; it is a sort of "discretization" of the classical description of fluids. The correspondence with the Ising model is established relating the local density of a fluid with the local magnetization density of an Ising model.
Let us therefore consider a -dimensional lattice with coordination number . Each site of the lattice can be occupied by a particle, so if we call the occupation number of a site we will have either or , and the total number of particles will be:
is the number of sites. In analogy with the continuum case we can guess a Hamiltonian of the form:
where the interaction potential
is symmetric, namely
We have written only the configurational term, because we know that the kinetic part always contributes to the grand partition function with the thermal wavelength.
This way, in the grand canonical ensemble we have:
Considering now an Ising model defined on the same lattice, we can establish a correspondence with the lattice gas model defining:
is the spin variable of the Ising model defined on the same lattice: in fact, doing so we have
What we now want to show is that substituting
we obtain a Hamiltonian that leads to the grand partition function of the Ising model, so that the equivalence between the two models is made explicit
From now on we will neglect any potential that couples more than two particles; therefore, the first term of expression of
while the second term:
(where in the second summation has been obtained relabelling indices and using the symmetry of
). If we now suppose that
is a nearest-neighbour interaction, i.e.:
If we now set
(i.e. there are no external fields acting on the fluid) we have:
For the sake of simplicity, we set
We thus have that the grand partition function of the lattice gas can be written as:
is the grand partition function of the Ising model.
Therefore, we see that the grand partition functions of both the Ising model and the lattice gas are equal, except for the irrelevant constant rescaling factor
This way we can see explicitly that the two models are perfectly equivalent.
Lattice gas and continuous fluids[edit | edit source]
We now want to show (although not really rigorously) that the lattice gas model can be derived from the "classical" model of a fluid.
Consider the configurational sum:
We can approximate it dividing our system in
cells of linear dimension
, such that the probability to find more than one particle inside one cell is negligible; in other words, if our particles interact via a potential that has a hard core component we can take
of the order of the hard core radius. This way the integral in
can be rewritten as:
labels the particles while
labels the cells.
Since we have focused our attention towards the cells, we can substitute the interaction
between particles with an interaction between occupied cells:
As we have done before we consider only two-particle interactions and suppose that the relative potential is short-ranged, namely
acts only between nearest neighbours and in that case is equal to a constant value
where of course
is the occupation number of the
Now, for each specified configuration of the occupation numbers
there can be
possible configurations for the positions
: in fact, the configuration
of the occupied cells only tells us where
we can find a particle, but not which
specific particle is in a given cell; furthermore, the system is left unchanged if we interchange the particles. We can thus write:
means that the sum on
must be subject to the constraint that the total number of particle be fixed and equal to
where "c.f." stands for "continuum fluid". The presence of the sum
strongly simplifies the calculations: in fact, if
is a generic function, then:
where the last sum in unconstrained. Therefore we can write
is the grand partition function of the lattice gas model, while:
is the "classical" chemical potential. However, since the chemical potential has the meaning of an energy, any shift in
by a constant quantity
is physically irrelevant, so
Therefore, the grand partition functions of the continuum fluid and the lattice gas model are the same: we thus see that these two models are equivalent.
- ↑ This of course is valid for classical mechanics, but the correspondence we will establish holds also in quantum statistical mechanics.
- ↑ The only big difference with The canonical ensemble at work is that we are considering generic instead of equal to 3.
- ↑ A small remark: for finite and "reasonable" systems the grand free energy is not singular even if it involves an infinite sum over . The reason is that generally (in the "reasonable" cases we have just mentioned) the interaction potentials have a hard-core component that prevents the particles from overlapping: therefore, a finite system will be able to contain only a finite number of particles, so that the sum has in reality an upper limit and is not infinite.
- ↑ Of course, if two systems have the same partition function, their thermodynamics will coincide, so they are at all effects equivalent.
- ↑ From these definition we see that in reality the precise values of the coefficients , , that we have encountered are absolutely irrelevant.
- ↑ Remember that our system is at fixed temperature, so is constant.
- ↑ Equivalently, we can note that , so that we also have : the grand partition functions of the fluid and the lattice gas differ for a constant rescaling factor, similarly to what we have seen before.