Editing Lattice model (physics)

{{Short description|A physical model defined on a lattice}}
{{refimprove|date=February 2016}}
{{other uses|Lattice model (disambiguation)}}

[[File:LatticemodelAB.png|thumb|A three-dimensional lattice filled with two molecules A and B, here shown as black and white spheres. Lattices such as this are used - for example - in the [[Flory–Huggins solution theory]]]]
In [[mathematical physics]], a '''lattice model''' is a [[mathematical model]] of a physical system that is defined on a [[lattice (group)|lattice]], as opposed to a [[continuum (theory)|continuum]], such as the continuum of [[space]] or [[spacetime]]. Lattice models originally occurred in the context of [[condensed matter physics]], where the [[atom]]s of a [[crystal]] automatically form a lattice. Currently, lattice models are quite popular in [[theoretical physics]], for many reasons. Some models are [[exactly solvable]], and thus offer insight into physics beyond what can be learned from [[perturbation theory]]. Lattice models are also ideal for study by the methods of [[computational physics]], as the discretization of any continuum model automatically turns it into a lattice model. The exact solution to many of these models (when they are solvable) includes the presence of [[soliton]]s. Techniques for solving these include the [[inverse scattering transform]] and the method of [[Lax pair]]s, the [[Yang–Baxter equation]] and [[quantum group]]s. The solution of these models has given insights into the nature of [[phase transitions]], [[magnetization]] and [[scaling behaviour]], as well as insights into the nature of [[quantum field theory]]. Physical lattice models frequently occur as an approximation to a continuum theory, either to give an [[ultraviolet cutoff]] to the theory to prevent divergences or to perform [[numerical analysis|numerical computations]]. An example of a continuum theory that is widely studied by lattice models is the [[QCD lattice model]], a discretization of [[quantum chromodynamics]]. However, [[digital physics]] considers nature fundamentally discrete at the Planck scale, which imposes [[Talk:Bekenstein bound#Cellular automata|upper limit to the density of information]], aka [[Holographic principle#Limit on information density|Holographic principle]]. More generally, [[lattice gauge theory]] and [[lattice field theory]] are areas of study. Lattice models are also used to simulate the structure and dynamics of polymers.

==Mathematical description==
A number of lattice models can be described by the following data:
* A [[lattice (group)|lattice]] <math>\Lambda</math>, often taken to be a lattice in <math>d</math>-dimensional Euclidean space <math>\mathbb{R}^d</math> or the <math>d</math>-dimensional torus if the lattice is periodic. Concretely, <math>\Lambda</math> is often the [[integer lattice|cubic lattice]]. If two points on the lattice are considered 'nearest neighbours', then they can be connected by an edge, turning the lattice into a [[lattice graph]]. The vertices of <math>\Lambda</math> are sometimes referred to as sites.
* A spin-variable space <math>S</math>. The configuration space <math>\mathcal{C}</math> of possible system states is then the space of functions <math>\sigma: \Lambda \rightarrow S</math>. For some models, we might instead consider instead the space of functions <math>\sigma: E \rightarrow S</math> where <math>E</math> is the edge set of the graph defined above.
* An energy functional <math>E:\mathcal{C}\rightarrow\mathbb{R}</math>, which might depend on a set of additional parameters or 'coupling constants' <math>\{g_i\}</math>.


===Examples===
The [[Ising model]] is given by the usual cubic lattice graph <math>G = (\Lambda, E)</math> where <math>\Lambda</math> is an infinite cubic lattice in <math>\mathbb{R}^d</math> or a period <math>n</math> cubic lattice in <math>T^d</math>, and <math>E</math> 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 <math>S = \{+1,-1\} = \mathbb{Z}_2</math>. The energy functional is 
:<math>E(\sigma) = -H\sum_{v\in\Lambda}\sigma(v) - J\sum_{\{v_1,v_2\}\in E}\sigma(v_1)\sigma(v_2).</math>

The spin-variable space can often be described as a [[coset]]. For example, for the Potts model we have <math>S = \mathbb{Z}_n</math>. In the limit <math>n\rightarrow \infty</math>, we obtain the XY model which has <math>S = SO(2)</math>. Generalising the XY model to higher dimensions gives the <math>n</math>-vector model which has <math>S = S^n = SO(n+1)/SO(n)</math>.

===Solvable models===
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 <math>n</math> in <math>d</math> dimensions. Then the configuration space <math>\mathcal{C}</math> is also finite. We can define the [[partition function (statistical mechanics)|partition function]]
:<math>Z = \sum_{\sigma \in \mathcal{C}}\exp(-\beta E(\sigma))</math>
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 <math>\{g_i\}</math> and <math>\beta</math>. 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, <math>H = 0,</math> but for dimension <math>d>2</math>, the Ising model remains unsolved.

===Mean field theory===
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. 
====Global mean field====
The configuration space <math>\mathcal{C}</math> of functions <math>\sigma</math> is replaced by the [[convex hull]] of the spin space <math>S</math>, when <math>S</math> has a realisation in terms of a subset of <math>\mathbb{R}^m</math>. We'll denote this by <math>\langle\mathcal{C}\rangle</math>. This arises as in going to the mean value of the field, we have <math>\sigma \mapsto \langle \sigma \rangle := \frac{1}{|\Lambda|}\sum_{v\in\Lambda}\sigma(v)</math>. 

As the number of lattice sites <math>N = |\Lambda|\rightarrow \infty</math>, the possible values of <math>\langle \sigma \rangle</math> fill out the convex hull of <math>S</math>. By making a suitable approximation, the energy functional becomes a function of the mean field, that is, <math>E(\sigma)\mapsto E(\langle \sigma \rangle).</math> The partition function then becomes
:<math>Z = \int_{\langle\mathcal{C}\rangle}d\langle\sigma\rangle e^{-\beta E(\langle\sigma\rangle)}\Omega(\langle\sigma\rangle) =: \int_{\langle\mathcal{C}\rangle}d\langle\sigma\rangle e^{-N\beta f(\langle\sigma\rangle)}.</math>
As <math>N\rightarrow \infty</math>, that is, in the [[thermodynamic limit]], the [[saddle point approximation]] tells us the integral is asymptotically dominated by the value at which <math>f(\langle\sigma\rangle)</math> is minimised:
:<math>Z \sim e^{-N\beta f(\langle\sigma\rangle_0)}</math>
where <math>\langle\sigma\rangle_0</math> is the argument minimising <math>f</math>.

A simpler, but less mathematically rigorous approach which nevertheless sometimes gives correct results comes from linearising the theory about the mean field <math>\langle\sigma\rangle</math>. Writing configurations as <math>\sigma(v)=\langle\sigma\rangle + \Delta\sigma(v)</math>, truncating terms of <math>\mathcal{O}(\Delta\sigma^2)</math> then summing over configurations allows computation of the partition function.

Such an approach to the periodic Ising model in <math>d</math> dimensions provides insight into [[phase transitions]].
====Spatially varying mean field====
Suppose the [[continuum limit]] of the lattice <math>\Lambda</math> is <math>\mathbb{R}^d</math>. Instead of averaging over all of <math>\Lambda</math>, we average over neighbourhoods of <math>\mathbf{x}\in\mathbb{R}^d</math>. This gives a spatially varying mean field <math>\langle\sigma\rangle:\mathbb{R}^d\rightarrow \langle\mathcal{C}\rangle</math>. We relabel <math>\langle\sigma\rangle</math> with <math>\phi</math> to bring the notation closer to field theory. This allows the partition function to be written as a [[path integral formulation|path integral]]
:<math>Z = \int \mathcal{D}\phi e^{-\beta F[\phi]}</math>
where the free energy <math>F[\phi]</math> is a [[Wick rotated]] version of the action in [[quantum field theory]].

==Examples==
===Condensed matter physics===
*[[Ising model]]
*[[ANNNI model]]
*[[Potts model]]
*[[Chiral Potts model]]
*[[XY model]]
*[[Classical Heisenberg model]]
*[[n-vector model]]
*[[Vertex model]]
*[[Toda lattice]]
*[[cellular automata]]

===Polymer physics===
* [[Bond fluctuation model]]
* [[2nd model]]

===High energy physics===
* [[QCD lattice model]]

==See also==
* [[Crystal structure]]
* [[Continuum limit]]
* [[QCD matter]]
* [[Lattice gas]]

==References==
*{{Citation | last1=Baxter | first1=Rodney J. | authorlink1=Rodney_Baxter | title=Exactly solved models in statistical mechanics | url=https://physics.anu.edu.au/theophys/_files/Exactly.pdf | publisher=Academic Press Inc. [Harcourt Brace Jovanovich Publishers] | location=London | isbn=978-0-12-083180-7 | mr=690578  | year=1982}}

{{DEFAULTSORT:Lattice Model (Physics)}}
[[Category:Lattice models| ]]