Square lattice Ising model

Last updated July 08, 2024

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 for the special case that the external magnetic field H = 0.^[1] An analytical solution for the general case for $H\neq 0$ has yet to be found.

Defining the partition function

Consider a 2D Ising model on a square lattice $\Lambda$ with N sites and periodic boundary conditions in both the horizontal and vertical directions, which effectively reduces the topology of the model to a torus. Generally, the horizontal coupling $J$ and the vertical coupling $J^{*}$ are not equal. With $\textstyle \beta ={\frac {1}{kT}}$ and absolute temperature $T$ and the Boltzmann constant $k$ , the partition function

Z_{N}(K\equiv \beta J,L\equiv \beta J^{*})=\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).

Critical temperature

The critical temperature $T_{\text{c}}$ can be obtained from the Kramers–Wannier duality relation. Denoting the free energy per site as $F(K,L)$ , one has:

\beta F\left(K^{*},L^{*}\right)=\beta F\left(K,L\right)+{\frac {1}{2}}\log {\big [}\sinh \left(2K\right)\sinh \left(2L\right){\big ]}

where

\sinh \left(2K^{*}\right)\sinh \left(2L\right)=1

\sinh \left(2L^{*}\right)\sinh \left(2K\right)=1

Assuming that there is only one critical line in the (K, L) plane, the duality relation implies that this is given by:

\sinh \left(2K\right)\sinh \left(2L\right)=1

For the isotropic case $J=J^{*}$ , one finds the famous relation for the critical temperature $T_{c}$

{\frac {kT_{\text{c}}}{J}}={\frac {2}{\ln(1+{\sqrt {2}})}}\approx 2.26918531421

Dual lattice

Consider a configuration of spins $\{\sigma \}$ on the square lattice $\Lambda$ . Let r and s denote the number of unlike neighbours in the vertical and horizontal directions respectively. Then the summand in $Z_{N}$ corresponding to $\{\sigma \}$ is given by

e^{K(N-2s)+L(N-2r)}

Dual lattice Duallattice.png — Dual lattice

Construct a dual lattice $\Lambda _{D}$ as depicted in the diagram. For every configuration $\{\sigma \}$ , 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 $\Lambda$ 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.

Spin configuration on a dual lattice Duallattices.png — Spin configuration on a dual lattice

This reduces the partition function to

Z_{N}(K,L)=2e^{N(K+L)}\sum _{P\subset \Lambda _{D}}e^{-2Lr-2Ks}

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 $T\rightarrow 0,\ \ e^{-K},e^{-L}\rightarrow 0$ , so that

Z_{N}(K,L)=2e^{N(K+L)}\sum _{P\subset \Lambda _{D}}e^{-2Lr-2Ks}

defines a low temperature expansion of $Z_{N}(K,L)$ .

High-temperature expansion

Since $\sigma \sigma '=\pm 1$ one has

e^{K\sigma \sigma '}=\cosh K+\sinh K(\sigma \sigma ')=\cosh K(1+\tanh K(\sigma \sigma ')).

Therefore

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})

where $v=\tanh K$ and $w=\tanh L$ . Since there are N horizontal and vertical edges, there are a total of $2^{2N}$ 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 $v\sigma _{i}\sigma _{j}$ (or $w\sigma _{i}\sigma _{j})$ is chosen in the product. Summing over the configurations, using

\sum _{\sigma _{i}=\pm 1}\sigma _{i}^{n}={\begin{cases}0&{\mbox{for }}n{\mbox{ odd}}\\2&{\mbox{for }}n{\mbox{ even}}\end{cases}}

shows that only configurations with an even number of lines at each vertex (polygons) will contribute to the partition function, giving

Z_{N}(K,L)=2^{N}(\cosh K\cosh L)^{N}\sum _{P\subset \Lambda }v^{r}w^{s}

where the sum is over all polygons in the lattice. Since tanh K, tanh L $\rightarrow 0$ as $T\rightarrow \infty$ , this gives the high temperature expansion of $Z_{N}(K,L)$ .

The two expansions can be related using the Kramers–Wannier duality.

Exact solution

The free energy per site in the limit $N\to \infty$ is given as follows. Define the parameter $k$ as

k={\frac {1}{\sinh \left(2K\right)\sinh \left(2L\right)}}

The Helmholtz free energy per site $F$ can be expressed as

-\beta F={\frac {\log(2)}{2}}+{\frac {1}{2\pi }}\int _{0}^{\pi }\log \left[\cosh \left(2K\right)\cosh \left(2L\right)+{\frac {1}{k}}{\sqrt {1+k^{2}-2k\cos(2\theta )}}\right]d\theta

For the isotropic case $J=J^{*}$ , from the above expression one finds for the internal energy per site:

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-4k(1+k)^{-2}\sin ^{2}(\theta )}}}d\theta \right]

and the spontaneous magnetization is, for $T<T_{\text{c}}$ ,

M=\left[1-\sinh ^{-4}(2\beta J)\right]^{1/8}

and $M=0$ for $T\geq T_{\text{c}}$ .

Notes

↑ Onsager, Lars (1944-02-01). "Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition". Physical Review. 65 (3–4): 117–149. doi:10.1103/PhysRev.65.117.

Related Research Articles

In physics, the Lorentz transformations are a six-parameter family of linear transformations from a coordinate frame in spacetime to another frame that moves at a constant velocity relative to the former. The respective inverse transformation is then parameterized by the negative of this velocity. The transformations are named after the Dutch physicist Hendrik Lorentz.

<span class="mw-page-title-main">Lattice model (physics)</span>

In mathematical physics, a lattice model is a mathematical model of a physical system that is defined on a lattice, as opposed to a continuum, such as the continuum of space or spacetime. Lattice models originally occurred in the context of condensed matter physics, where the atoms 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 includes the presence of solitons. Techniques for solving these include the inverse scattering transform and the method of Lax pairs, the Yang–Baxter equation and quantum groups. 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 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 upper limit to the density of information, aka 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.

The Ising model, named after the physicists Ernst Ising and Wilhelm Lenz, is a mathematical model of ferromagnetism in statistical mechanics. The model consists of discrete variables that represent magnetic dipole moments of atomic "spins" that can be in one of two states. The spins are arranged in a graph, usually a lattice, allowing each spin to interact with its neighbors. Neighboring spins that agree have a lower energy than those that disagree; the system tends to the lowest energy but heat disturbs this tendency, thus creating the possibility of different structural phases. The model allows the identification of phase transitions as a simplified model of reality. The two-dimensional square-lattice Ising model is one of the simplest statistical models to show a phase transition.

In mathematical physics, n-dimensional de Sitter space is a maximally symmetric Lorentzian manifold with constant positive scalar curvature. It is the Lorentzian analogue of an n-sphere.

The classical XY model is a lattice model of statistical mechanics. In general, the XY model can be seen as a specialization of Stanley's n-vector model for $n = 2$ .

In statistical mechanics and mathematics, the Bethe lattice is an infinite connected cycle-free graph where all vertices have the same number of neighbors. The Bethe lattice was introduced into the physics literature by Hans Bethe in 1935. In such a graph, each node is connected to z neighbors; the number z is called either the coordination number or the degree, depending on the field.

In mathematics, the Weierstrass functions are special functions of a complex variable that are auxiliary to the Weierstrass elliptic function. They are named for Karl Weierstrass. The relation between the sigma, zeta, and $functions is analogous to that between the sine, cotangent, and squared cosecant functions: the logarithmic derivative of the sine is the cotangent, whose derivative is negative the squared cosecant.$

<span class="mw-page-title-main">Universal Transverse Mercator coordinate system</span> Map projection system

The Universal Transverse Mercator (UTM) is a map projection system for assigning coordinates to locations on the surface of the Earth. Like the traditional method of latitude and longitude, it is a horizontal position representation, which means it ignores altitude and treats the earth surface as a perfect ellipsoid. However, it differs from global latitude/longitude in that it divides earth into 60 zones and projects each to the plane as a basis for its coordinates. Specifying a location means specifying the zone and the x, y coordinate in that plane. The projection from spheroid to a UTM zone is some parameterization of the transverse Mercator projection. The parameters vary by nation or region or mapping system.

In physics and mathematics, the Gibbs measure, named after Josiah Willard Gibbs, is a probability measure frequently seen in many problems of probability theory and statistical mechanics. It is a generalization of the canonical ensemble to infinite systems. The canonical ensemble gives the probability of the system X being in state x as

The Kramers–Wannier duality is a symmetry in statistical physics. It relates the free energy of a two-dimensional square-lattice Ising model at a low temperature to that of another Ising model at a high temperature. It was discovered by Hendrik Kramers and Gregory Wannier in 1941. With the aid of this duality Kramers and Wannier found the exact location of the critical point for the Ising model on the square lattice.

The quantum Heisenberg model, developed by Werner Heisenberg, is a statistical mechanical model used in the study of critical points and phase transitions of magnetic systems, in which the spins of the magnetic systems are treated quantum mechanically. It is related to the prototypical Ising model, where at each site of a lattice, a spin $represents a microscopic magnetic dipole to which the magnetic moment is either up or down. Except the coupling between magnetic dipole moments, there is also a multipolar version of Heisenberg model called the multipolar exchange interaction.$

In condensed matter physics and crystallography, the static structure factor is a mathematical description of how a material scatters incident radiation. The structure factor is a critical tool in the interpretation of scattering patterns obtained in X-ray, electron and neutron diffraction experiments.

The spherical model is a model of ferromagnetism similar to the Ising model, which was solved in 1952 by T. H. Berlin and M. Kac. It has the remarkable property that for linear dimension d greater than four, the critical exponents that govern the behaviour of the system near the critical point are independent of d and the geometry of the system. It is one of the few models of ferromagnetism that can be solved exactly in the presence of an external field.

In complex analysis, a partial fraction expansion is a way of writing a meromorphic function $as an infinite sum of rational functions and polynomials. When is a rational function, this reduces to the usual method of partial fractions.$

<span class="mw-page-title-main">Genus of a multiplicative sequence</span> A ring homomorphism from the cobordism ring of manifolds to another ring

In mathematics, a genus of a multiplicative sequence is a ring homomorphism from the ring of smooth compact manifolds up to the equivalence of bounding a smooth manifold with boundary to another ring, usually the rational numbers, having the property that they are constructed from a sequence of polynomials in characteristic classes that arise as coefficients in formal power series with good multiplicative properties.

In mathematics, the Plancherel theorem for spherical functions is an important result in the representation theory of semisimple Lie groups, due in its final form to Harish-Chandra. It is a natural generalisation in non-commutative harmonic analysis of the Plancherel formula and Fourier inversion formula in the representation theory of the group of real numbers in classical harmonic analysis and has a similarly close interconnection with the theory of differential equations. It is the special case for zonal spherical functions of the general Plancherel theorem for semisimple Lie groups, also proved by Harish-Chandra. The Plancherel theorem gives the eigenfunction expansion of radial functions for the Laplacian operator on the associated symmetric space X; it also gives the direct integral decomposition into irreducible representations of the regular representation on $L 2 (X)$ . In the case of hyperbolic space, these expansions were known from prior results of Mehler, Weyl and Fock.

In statistical mechanics, the eight-vertex model is a generalisation of the ice-type (six-vertex) models; it was discussed by Sutherland, and Fan & Wu, and solved by Baxter in the zero-field case.

The Swendsen–Wang algorithm is the first non-local or cluster algorithm for Monte Carlo simulation for large systems near criticality. It has been introduced by Robert Swendsen and Jian-Sheng Wang in 1987 at Carnegie Mellon.

<span class="mw-page-title-main">Bending of plates</span> Deformation of slabs under load

Bending of plates, or plate bending, refers to the deflection of a plate perpendicular to the plane of the plate under the action of external forces and moments. The amount of deflection can be determined by solving the differential equations of an appropriate plate theory. The stresses in the plate can be calculated from these deflections. Once the stresses are known, failure theories can be used to determine whether a plate will fail under a given load.

In statistical mechanics, the Griffiths inequality, sometimes also called Griffiths–Kelly–Sherman inequality or GKS inequality, named after Robert B. Griffiths, is a correlation inequality for ferromagnetic spin systems. Informally, it says that in ferromagnetic spin systems, if the 'a-priori distribution' of the spin is invariant under spin flipping, the correlation of any monomial of the spins is non-negative; and the two point correlation of two monomial of the spins is non-negative.

References

Baxter, Rodney J. (1982), Exactly solved models in statistical mechanics (PDF), London: Academic Press Inc. [Harcourt Brace Jovanovich Publishers], ISBN 978-0-12-083180-7, MR 0690578
Baxter, Rodney J. (2016). "The bulk, surface and corner free energies of the square lattice Ising model". Journal of Physics A: Mathematical and Theoretical. 50 (1). IOP Publishing: 014001. arXiv: 1606.02029 . doi:10.1088/1751-8113/50/1/014001. ISSN 1751-8113. S2CID 2467419.
Kurt Binder (2001) [1994], "Ising model", Encyclopedia of Mathematics , EMS Press
BRUSH, STEPHEN G. (1967-10-01). "History of the Lenz-Ising Model". Reviews of Modern Physics. 39 (4). American Physical Society (APS): 883–893. Bibcode:1967RvMP...39..883B. doi:10.1103/revmodphys.39.883. ISSN 0034-6861.
Huang, Kerson (1987), Statistical mechanics (2nd edition), Wiley, ISBN 978-0471815181
Hucht, Alfred (2021). "The square lattice Ising model on the rectangle III: Hankel and Toeplitz determinants". Journal of Physics A: Mathematical and Theoretical. 54 (37). IOP Publishing: 375201. arXiv: 2103.10776 . Bibcode:2021JPhA...54K5201H. doi:10.1088/1751-8121/ac0983. ISSN 1751-8113. S2CID 232290629.
Ising, Ernst (1925), "Beitrag zur Theorie des Ferromagnetismus", Z. Phys., 31 (1): 253–258, Bibcode:1925ZPhy...31..253I, doi:10.1007/BF02980577, S2CID 122157319
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, Bibcode:1963JMP.....4..308M, doi:10.1063/1.1703955, ISSN 0022-2488, MR 0148406, archived from the original on 2013-01-12
Onsager, Lars (1944), "Crystal statistics. I. A two-dimensional model with an order-disorder transition", Phys. Rev., Series II, 65 (3–4): 117–149, Bibcode:1944PhRv...65..117O, doi:10.1103/PhysRev.65.117, MR 0010315
Onsager, Lars (1949), "Discussion", Supplemento al Nuovo Cimento, 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, Bibcode:1952PhRv...85..808Y, doi:10.1103/PhysRev.85.808, MR 0051740

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Onsager, Lars (1944-02-01). "Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition". Physical Review. 65 (3–4): 117–149. doi:10.1103/PhysRev.65.117.

[1]