Physics:Square-lattice Ising model
In statistical mechanics, the two-dimensional square-lattice Ising model is a simple lattice model of interacting magnetic spins. The model is notable for having nontrivial interactions, yet having an analytical solution. The model was solved by Lars Onsager (Lars Onsager 1944) for the special case that the external magnetic field H = 0. An analytical solution for the general case for [math]\displaystyle{ H \neq 0 }[/math] has yet to be found.
Definition of the model
Consider the 2D Ising model on a square lattice [math]\displaystyle{ \Lambda }[/math] with N sites, with periodic boundary conditions in both the horizontal and vertical directions, which effectively reduces the topology of the model to a torus. In a general case, the horizontal coupling J is not equal to the coupling in the vertical direction, J*. With an equal number of rows and columns in the lattice, there will be N of each. In terms of
- [math]\displaystyle{ K = \beta J }[/math]
- [math]\displaystyle{ L = \beta J^{*} }[/math]
where [math]\displaystyle{ \beta = \frac{1}{kT} }[/math] where T is absolute temperature and k is Boltzmann's constant, the partition function [math]\displaystyle{ Z_N(K,L) }[/math] is given by
- [math]\displaystyle{ Z_N(K,L) = \sum_{\{\sigma\}} \exp \left( K \sum_{\langle ij \rangle_H} \sigma_i \sigma_j + L \sum_{\langle ij \rangle_V} \sigma_i \sigma_j \right). }[/math]
Critical temperature
The critical temperature [math]\displaystyle{ T_c }[/math] can be obtained from the Kramers–Wannier duality relation. Denoting the free energy per site as [math]\displaystyle{ F(K,L) }[/math], one has:
- [math]\displaystyle{ \beta F\left(K^{*}, L^{*}\right) = \beta F\left(K,L\right) + \frac{1}{2}\log\left[\sinh\left(2K\right)\sinh\left(2L\right)\right] }[/math]
where
- [math]\displaystyle{ \sinh\left(2K^{*}\right)\sinh\left(2L\right)=1 }[/math]
- [math]\displaystyle{ \sinh\left(2L^{*}\right)\sinh\left(2K\right)=1 }[/math]
Assuming there is only one critical line in the (K,L) plane, the duality relation implies that this is given by:
- [math]\displaystyle{ \sinh\left(2 K\right)\sinh\left(2 L\right)= 1 }[/math]
For the isotropic case [math]\displaystyle{ J = J^{*} }[/math], one finds the famous relation for the critical temperature [math]\displaystyle{ T_{c} }[/math]
- [math]\displaystyle{ \frac{k T_c}{J} = \frac{2}{\ln(1+\sqrt{2})} \approx 2.26918531421 }[/math]
Dual lattice
Consider a configuration of spins [math]\displaystyle{ \{ \sigma \} }[/math] on the square lattice [math]\displaystyle{ \Lambda }[/math]. Let r and s denote the number of unlike neighbours in the vertical and horizontal directions respectively. Then the summand in [math]\displaystyle{ Z_N }[/math] corresponding to [math]\displaystyle{ \{ \sigma \} }[/math] is given by
- [math]\displaystyle{ e^{K(N-2s) +L(N-2r)} }[/math]
Construct a dual lattice [math]\displaystyle{ \Lambda_D }[/math] as depicted in the diagram. For every configuration [math]\displaystyle{ \{ \sigma \} }[/math], a polygon is associated to the lattice by drawing a line on the edge of the dual lattice if the spins separated by the edge are unlike. Since by traversing a vertex of [math]\displaystyle{ \Lambda }[/math] the spins need to change an even number of times so that one arrives at the starting point with the same charge, every vertex of the dual lattice is connected to an even number of lines in the configuration, defining a polygon.
This reduces the partition function to
- [math]\displaystyle{ Z_N(K,L) = 2e^{N(K+L)} \sum_{P \subset \Lambda_D} e^{-2Lr-2Ks} }[/math]
summing over all polygons in the dual lattice, where r and s are the number of horizontal and vertical lines in the polygon, with the factor of 2 arising from the inversion of spin configuration.
Low-temperature expansion
At low temperatures, K, L approach infinity, so that as [math]\displaystyle{ T \rightarrow 0, \ \ e^{-K}, e^{-L} \rightarrow 0 }[/math], so that
- [math]\displaystyle{ Z_N(K,L) = 2 e^{N(K+L)} \sum_{ P \subset \Lambda_D} e^{-2Lr-2Ks} }[/math]
defines a low temperature expansion of [math]\displaystyle{ Z_N(K,L) }[/math].
High-temperature expansion
Since [math]\displaystyle{ \sigma \sigma' = \pm 1 }[/math] one has
- [math]\displaystyle{ e^{K \sigma \sigma'} = \cosh K + \sinh K(\sigma \sigma') = \cosh K(1+\tanh K(\sigma \sigma')). }[/math]
Therefore
- [math]\displaystyle{ Z_N(K,L) = (\cosh K \cosh L)^N \sum_{\{ \sigma \}} \prod_{\langle ij \rangle_H} (1+v \sigma_i \sigma_j) \prod_{\langle ij \rangle_V}(1+w\sigma_i \sigma_j) }[/math]
where [math]\displaystyle{ v =\tanh K }[/math] and [math]\displaystyle{ w = \tanh L }[/math]. Since there are N horizontal and vertical edges, there are a total of [math]\displaystyle{ 2^{2N} }[/math] terms in the expansion. Every term corresponds to a configuration of lines of the lattice, by associating a line connecting i and j if the term [math]\displaystyle{ v \sigma_i \sigma_j }[/math] (or [math]\displaystyle{ w \sigma_i \sigma_j) }[/math] is chosen in the product. Summing over the configurations, using
- [math]\displaystyle{ \sum_{\sigma_i = \pm 1} \sigma_i^n = \begin{cases} 0 & \mbox{for } n \mbox{ odd} \\ 2 & \mbox{for } n \mbox{ even} \end{cases} }[/math]
shows that only configurations with an even number of lines at each vertex (polygons) will contribute to the partition function, giving
- [math]\displaystyle{ Z_N(K,L) = 2^N(\cosh K \cosh L)^N \sum_{P \subset \Lambda} v^r w^s }[/math]
where the sum is over all polygons in the lattice. Since tanh K, tanh L [math]\displaystyle{ \rightarrow 0 }[/math] as [math]\displaystyle{ T \rightarrow \infty }[/math], this gives the high temperature expansion of [math]\displaystyle{ Z_N(K,L) }[/math].
The two expansions can be related using the Kramers–Wannier duality.
Exact solution
The free energy per site in the limit [math]\displaystyle{ N\to\infty }[/math] is given as follows. Define the parameter [math]\displaystyle{ k }[/math] as
- [math]\displaystyle{ k =\frac{1}{\sinh\left(2 K\right)\sinh\left(2 L\right)} }[/math]
The Helmholtz free energy per site [math]\displaystyle{ F }[/math] can be expressed as
- [math]\displaystyle{ -\beta F = \frac{\log(2)}{2} + \frac{1}{2\pi}\int_{0}^{\pi}\log\left[\cosh\left(2 K\right)\cosh\left(2 L\right)+\frac{1}{k}\sqrt{1+k^{2}-2k\cos(2\theta)}\right]d\theta }[/math]
For the isotropic case [math]\displaystyle{ J = J^{*} }[/math], from the above expression one finds for the internal energy per site:
- [math]\displaystyle{ U = - J \coth(2 \beta J) \left[ 1 + \frac{2}{\pi} (2 \tanh^2(2 \beta J) -1) \int_0^{\pi/2} \frac{1}{\sqrt{1 - 4 k (1+k)^{-2} \sin^2(\theta)}} d\theta \right] }[/math]
and the spontaneous magnetization is, for [math]\displaystyle{ T \lt T_c }[/math],
- [math]\displaystyle{ M = \left[ 1 - \sinh^{-4}(2 \beta J) \right]^{1/8} }[/math]
References
- Baxter, Rodney J. (1982), Exactly solved models in statistical mechanics, London: Academic Press Inc. [Harcourt Brace Jovanovich Publishers], ISBN 978-0-12-083180-7, https://physics.anu.edu.au/theophys/_files/Exactly.pdf
- Hazewinkel, Michiel, ed. (2001), "Ising model", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4, https://www.encyclopediaofmath.org/index.php?title=I/i120080
- Stephen G. Brush (1967), History of the Lenz-Ising Model. Reviews of Modern Physics (American Physical Society) vol. 39, pp 883–893. doi:10.1103/RevModPhys.39.883
- Huang, Kerson (1987), Statistical mechanics (2nd edition), Wiley, ISBN 978-0471815181
- Ising, E. (1925), "Beitrag zur Theorie des Ferromagnetismus", Z. Phys. 31 (1): 253–258, doi:10.1007/BF02980577, Bibcode: 1925ZPhy...31..253I
- Itzykson, Claude; Drouffe, Jean-Michel (1989), Théorie statistique des champs, Volume 1, Savoirs actuels (CNRS), EDP Sciences Editions, ISBN 978-2868833600
- Itzykson, Claude; Drouffe, Jean-Michel (1989), Statistical field theory, Volume 1: From Brownian motion to renormalization and lattice gauge theory, Cambridge University Press, ISBN 978-0521408059
- Barry M. McCoy and Tai Tsun Wu (1973), The Two-Dimensional Ising Model. Harvard University Press, Cambridge Massachusetts, ISBN:0-674-91440-6
- Montroll, Elliott W.; Potts, Renfrey B.; Ward, John C. (1963), "Correlations and spontaneous magnetization of the two-dimensional Ising model", Journal of Mathematical Physics 4 (2): 308–322, doi:10.1063/1.1703955, ISSN 0022-2488, Bibcode: 1963JMP.....4..308M, archived from the original on 2013-01-12, https://archive.is/20130112095848/http://link.aip.org/link/?JMAPAQ/4/308/1
- Onsager, Lars (1944), "Crystal statistics. I. A two-dimensional model with an order-disorder transition", Phys. Rev., Series II 65 (3–4): 117–149, doi:10.1103/PhysRev.65.117, Bibcode: 1944PhRv...65..117O
- Onsager, Lars (1949), "Discussion", Nuovo Cimento Supplement 6: 261
- John Palmer (2007), Planar Ising Correlations. Birkhäuser, Boston, ISBN:978-0-8176-4248-8.
- Yang, C. N. (1952), "The spontaneous magnetization of a two-dimensional Ising model", Physical Review, Series II 85 (5): 808–816, doi:10.1103/PhysRev.85.808, Bibcode: 1952PhRv...85..808Y