A number of lattice models can be described by the following data:
The Ising model is given by the usual cubic lattice graph G = ( Λ , E ) {\displaystyle G=(\Lambda ,E)} where Λ {\displaystyle \Lambda } is an infinite cubic lattice in R d {\displaystyle \mathbb {R} ^{d}} or a period n {\displaystyle n} cubic lattice in T d {\displaystyle T^{d}} , and E {\displaystyle E} is the edge set of nearest neighbours (the same letter is used for the energy functional but the different usages are distinguishable based on context). The spin-variable space is S = { + 1 , − 1 } = Z 2 {\displaystyle S=\{+1,-1\}=\mathbb {Z} _{2}} . The energy functional is
The spin-variable space can often be described as a coset. For example, for the Potts model we have S = Z n {\displaystyle S=\mathbb {Z} _{n}} . In the limit n → ∞ {\displaystyle n\rightarrow \infty } , we obtain the XY model which has S = S O ( 2 ) {\displaystyle S=SO(2)} . Generalising the XY model to higher dimensions gives the n {\displaystyle n} -vector model which has S = S n = S O ( n + 1 ) / S O ( n ) {\displaystyle S=S^{n}=SO(n+1)/SO(n)} .
We specialise to a lattice with a finite number of points, and a finite spin-variable space. This can be achieved by making the lattice periodic, with period n {\displaystyle n} in d {\displaystyle d} dimensions. Then the configuration space C {\displaystyle {\mathcal {C}}} is also finite. We can define the partition function
and there are no issues of convergence (like those which emerge in field theory) since the sum is finite. In theory, this sum can be computed to obtain an expression which is dependent only on the parameters { g i } {\displaystyle \{g_{i}\}} and β {\displaystyle \beta } . In practice, this is often difficult due to non-linear interactions between sites. Models with a closed-form expression for the partition function are known as exactly solvable.
Examples of exactly solvable models are the periodic 1D Ising model, and the periodic 2D Ising model with vanishing external magnetic field, H = 0 , {\displaystyle H=0,} but for dimension d > 2 {\displaystyle d>2} , the Ising model remains unsolved.
Due to the difficulty of deriving exact solutions, in order to obtain analytic results we often must resort to mean field theory. This mean field may be spatially varying, or global.
The configuration space C {\displaystyle {\mathcal {C}}} of functions σ {\displaystyle \sigma } is replaced by the convex hull of the spin space S {\displaystyle S} , when S {\displaystyle S} has a realisation in terms of a subset of R m {\displaystyle \mathbb {R} ^{m}} . We'll denote this by ⟨ C ⟩ {\displaystyle \langle {\mathcal {C}}\rangle } . This arises as in going to the mean value of the field, we have σ ↦ ⟨ σ ⟩ := 1 | Λ | ∑ v ∈ Λ σ ( v ) {\displaystyle \sigma \mapsto \langle \sigma \rangle :={\frac {1}{|\Lambda |}}\sum _{v\in \Lambda }\sigma (v)} .
As the number of lattice sites N = | Λ | → ∞ {\displaystyle N=|\Lambda |\rightarrow \infty } , the possible values of ⟨ σ ⟩ {\displaystyle \langle \sigma \rangle } fill out the convex hull of S {\displaystyle S} . By making a suitable approximation, the energy functional becomes a function of the mean field, that is, E ( σ ) ↦ E ( ⟨ σ ⟩ ) . {\displaystyle E(\sigma )\mapsto E(\langle \sigma \rangle ).} The partition function then becomes
As N → ∞ {\displaystyle N\rightarrow \infty } , that is, in the thermodynamic limit, the saddle point approximation tells us the integral is asymptotically dominated by the value at which f ( ⟨ σ ⟩ ) {\displaystyle f(\langle \sigma \rangle )} is minimised:
where ⟨ σ ⟩ 0 {\displaystyle \langle \sigma \rangle _{0}} is the argument minimising f {\displaystyle f} .
A simpler, but less mathematically rigorous approach which nevertheless sometimes gives correct results comes from linearising the theory about the mean field ⟨ σ ⟩ {\displaystyle \langle \sigma \rangle } . Writing configurations as σ ( v ) = ⟨ σ ⟩ + Δ σ ( v ) {\displaystyle \sigma (v)=\langle \sigma \rangle +\Delta \sigma (v)} , truncating terms of O ( Δ σ 2 ) {\displaystyle {\mathcal {O}}(\Delta \sigma ^{2})} then summing over configurations allows computation of the partition function.
Such an approach to the periodic Ising model in d {\displaystyle d} dimensions provides insight into phase transitions.
Suppose the continuum limit of the lattice Λ {\displaystyle \Lambda } is R d {\displaystyle \mathbb {R} ^{d}} . Instead of averaging over all of Λ {\displaystyle \Lambda } , we average over neighbourhoods of x ∈ R d {\displaystyle \mathbf {x} \in \mathbb {R} ^{d}} . This gives a spatially varying mean field ⟨ σ ⟩ : R d → ⟨ C ⟩ {\displaystyle \langle \sigma \rangle :\mathbb {R} ^{d}\rightarrow \langle {\mathcal {C}}\rangle } . We relabel ⟨ σ ⟩ {\displaystyle \langle \sigma \rangle } with ϕ {\displaystyle \phi } to bring the notation closer to field theory. This allows the partition function to be written as a path integral
where the free energy F [ ϕ ] {\displaystyle F[\phi ]} is a Wick rotated version of the action in quantum field theory.