Square lattice Ising model

Short description: Model in statistical mechanics

In statistical mechanics, the two-dimensional square lattice Ising model is a simple lattice model of interacting magnetic spins, an example of the class of Ising models. 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≠0 has yet to be found.

## Contents

* 1 Defining the partition function
* 2 Critical temperature
* 3 Dual lattice
* 4 Low-temperature expansion
* 5 High-temperature expansion
* 6 Exact solution
* 6.1 Transfer matrix
* 6.2 T in terms of Pauli matrices
* 6.3 Spin flip creation and annihilation operators
* 7 Notes
* 8 References

## Defining the partition function

Consider a 2D Ising model on a square lattice Λ 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 β=1kT and absolute temperature T and the Boltzmann constant k, the partition function

ZN(K≡βJ,L≡βJ*)=∑{σ}exp⁡(K∑⟨ij⟩Hσiσj+L∑⟨ij⟩Vσiσj).

## Critical temperature

The critical temperature Tc can be obtained from the Kramers–Wannier duality relation. Denoting the free energy per site as F(K,L), one has:

βF(K*,L*)=βF(K,L)+12log⁡[sinh⁡(2K)sinh⁡(2L)]

where

sinh⁡(2K*)sinh⁡(2L)=1
sinh⁡(2L*)sinh⁡(2K)=1

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

sinh⁡(2K)sinh⁡(2L)=1

For the isotropic case J=J*, one finds the famous relation for the critical temperature Tc

kTcJ=2ln⁡(1+2)≈2.26918531421

## Dual lattice

Consider a configuration of spins {σ} on the square lattice Λ. Let r and s denote the number of unlike neighbours in the vertical and horizontal directions respectively. Then the summand in ZN corresponding to {σ} is given by

eK(N−2s)+L(N−2r)
Dual lattice

Construct a dual lattice ΛD as depicted in the diagram. For every configuration {σ}, 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 Λ 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

This reduces the partition function to

ZN(K,L)=2eN(K+L)∑P⊂ΛDe−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->0, e−K,e−L->0, so that

ZN(K,L)=2eN(K+L)∑P⊂ΛDe−2Lr−2Ks

defines a low temperature expansion of ZN(K,L).

## High-temperature expansion

Since σσ′=±1 one has

eKσσ′=cosh⁡K+sinh⁡K(σσ′)=cosh⁡K(1+tanh⁡K(σσ′)).

Therefore

ZN(K,L)=(cosh⁡Kcosh⁡L)N∑{σ}∏⟨ij⟩H(1+vσiσj)∏⟨ij⟩V(1+wσiσj)

where v=tanh⁡K and w=tanh⁡L. Since there are N horizontal and vertical edges, there are a total of 22N 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σiσj (or wσiσj) is chosen in the product. Summing over the configurations, using

∑σi=±1σin={0for n odd2for n even

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

ZN(K,L)=2N(cosh⁡Kcosh⁡L)N∑P⊂Λvrws

where the sum is over all polygons in the lattice. Since tanh K, tanh L ->0 as T->∞, this gives the high temperature expansion of ZN(K,L).

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

## Exact solution

The free energy per site in the limit N->∞ is given as follows. Define the parameter k as

k=1sinh⁡(2K)sinh⁡(2L)

The Helmholtz free energy per site F can be expressed as

−βF=log⁡(2)2+12π∫0πlog⁡[cosh⁡(2K)cosh⁡(2L)+1k1+k2−2kcos⁡(2θ)]dθ

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

U=−Jcoth⁡(2βJ)[1+2π(2tanh2(2βJ)−1)∫0π/211−4k(1+k)−2sin2(θ)dθ]

and the spontaneous magnetization is, for T<Tc,

M=[1−sinh−4(2βJ)]1/8

and M=0 for T≥Tc.

### Transfer matrix

Start with an analogy with quantum mechanics. The Ising model on a long periodic lattice has a partition function

∑{S}exp⁡(∑ijSi,j(Si,j+1+Si+1,j)).

Think of the i direction as space, and the j direction as time. This is an independent sum over all the values that the spins can take at each time slice. This is a type of path integral, it is the sum over all spin histories.

A path integral can be rewritten as a Hamiltonian evolution. The Hamiltonian steps through time by performing a unitary rotation between time t and time t \+ Δt: U=eiHΔt

The product of the U matrices, one after the other, is the total time evolution operator, which is the path integral we started with.

UN=(eiHΔt)N=∫DXeiL

where N is the number of time slices. The sum over all paths is given by a product of matrices, each matrix element is the transition probability from one slice to the next.

Similarly, one can divide the sum over all partition function configurations into slices, where each slice is the one-dimensional configuration at time 1. This defines the transfer matrix: TC1C2.

The configuration in each slice is a one-dimensional collection of spins. At each time slice, T has matrix elements between two configurations of spins, one in the immediate future and one in the immediate past. These two configurations are C1 and C2, and they are all one-dimensional spin configurations. We can think of the vector space that T acts on as all complex linear combinations of these. Using quantum mechanical notation: |A⟩=∑SA(S)|S⟩

where each basis vector |S⟩ is a spin configuration of a one-dimensional Ising model.

Like the Hamiltonian, the transfer matrix acts on all linear combinations of states. The partition function is a matrix function of T, which is defined by the sum over all histories which come back to the original configuration after N steps: Z=tr(TN).

Since this is a matrix equation, it can be evaluated in any basis. So if we can diagonalize the matrix T, we can find Z.

### T in terms of Pauli matrices

The contribution to the partition function for each past/future pair of configurations on a slice is the sum of two terms. There is the number of spin flips in the past slice and there is the number of spin flips between the past and future slice. Define an operator on configurations which flips the spin at site i:

σix.

In the usual Ising basis, acting on any linear combination of past configurations, it produces the same linear combination but with the spin at position i of each basis vector flipped.

Define a second operator which multiplies the basis vector by +1 and −1 according to the spin at position i:

σiz.

T can be written in terms of these:

∑iAσix+Bσizσi+1z

where A and B are constants which are to be determined so as to reproduce the partition function. The interpretation is that the statistical configuration at this slice contributes according to both the number of spin flips in the slice, and whether or not the spin at position i has flipped.

### Spin flip creation and annihilation operators

Just as in the one-dimensional case, we will shift attention from the spins to the spin-flips. The σz term in T counts the number of spin flips, which we can write in terms of spin-flip creation and annihilation operators:

∑Cψi†ψi.

The first term flips a spin, so depending on the basis state it either:

1. moves a spin-flip one unit to the right
2. moves a spin-flip one unit to the left
3. produces two spin-flips on neighboring sites
4. destroys two spin-flips on neighboring sites.

Writing this out in terms of creation and annihilation operators: σix=Dψ†iψi+1+D*ψ†iψi−1+Cψiψi+1+C*ψ†iψ†i+1.

Ignore the constant coefficients, and focus attention on the form. They are all quadratic. Since the coefficients are constant, this means that the T matrix can be diagonalized by Fourier transforms.

Carrying out the diagonalization produces the Onsager free energy.

## Notes

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. Bibcode: 1944PhRv...65..117O. https://link.aps.org/doi/10.1103/PhysRev.65.117.

## 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
* Baxter, Rodney J. (2016). "The bulk, surface and corner free energies of the square lattice Ising model". Journal of Physics A: Mathematical and Theoretical (IOP Publishing) 50 (1). doi:10.1088/1751-8113/50/1/014001. ISSN 1751-8113.
* 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
* BRUSH, STEPHEN G. (1967-10-01). "History of the Lenz-Ising Model". Reviews of Modern Physics (American Physical Society (APS)) 39 (4): 883–893. doi:10.1103/revmodphys.39.883. ISSN 0034-6861. Bibcode: 1967RvMP...39..883B.
* Huang, Kerson (1987), Statistical mechanics (2nd edition), Wiley, ISBN 978-0-471-81518-1
* Hucht, Alfred (2021). "The square lattice Ising model on the rectangle III: Hankel and Toeplitz determinants". Journal of Physics A: Mathematical and Theoretical (IOP Publishing) 54 (37): 375201. doi:10.1088/1751-8121/ac0983. ISSN 1751-8113. Bibcode: 2021JPhA...54K5201H.
* Ising, Ernst (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-2-86883-360-0
* 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-0-521-40805-9
* 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, 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", 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, doi:10.1103/PhysRev.85.808, Bibcode: 1952PhRv...85..808Y

0.00

(0 votes)

Original source: https://en.wikipedia.org/wiki/Square lattice Ising model. Read more

|

Retrieved from "https://handwiki.org/wiki/index.php?title=Physics:Square_lattice_Ising_model&oldid=4184515"