Eight-vertex model explained

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

Description

As with the ice-type models, the eight-vertex model is a square lattice model, where each state is a configuration of arrows at a vertex. The allowed vertices have an even number of arrows pointing towards the vertex; these include the six inherited from the ice-type model (1-6), and sinks and sources (7, 8).

We consider a

N x N

lattice, with

N2

vertices and

2N2

edges. Imposing periodic boundary conditions requires that the states 7 and 8 occur equally often, as do states 5 and 6, and thus can be taken to have the same energy. For the zero-field case the same is true for the two other pairs of states. Each vertex

j

has an associated energy

\epsilonj

and Boltzmann weight
-\epsilonj
kT
w
j=e
, giving the partition function over the lattice as

Z=\sum\exp\left(-

\sumjnj\epsilonj
kT

\right)

where the summation is over all allowed configurations of vertices in the lattice. In this general form the partition function remains unsolved.

Solution in the zero-field case

The zero-field case of the model corresponds physically to the absence of external electric fields. Hence, the model remains unchanged under the reversal of all arrows; the states 1 and 2, and 3 and 4, consequently must occur as pairs. The vertices can be assigned arbitrary weights

\begin{align} w1=w2&=a\\ w3=w4&=b\\ w5=w6&=c\\ w7=w8&=d. \end{align}

The solution is based on the observation that rows in transfer matrices commute, for a certain parametrisation of these four Boltzmann weights. This came about as a modification of an alternate solution for the six-vertex model; it makes use of elliptic theta functions.

Commuting transfer matrices

The proof relies on the fact that when

\Delta'=\Delta

and

\Gamma'=\Gamma

, for quantities
\begin{align} \Delta&=a2+b2-c2-d2\\ \Gamma&=
2(ab+cd)
ab-cd
ab+cd

\end{align}

the transfer matrices

T

and

T'

(associated with the weights

a

,

b

,

c

,

d

and

a'

,

b'

,

c'

,

d'

) commute. Using the star-triangle relation, Baxter reformulated this condition as equivalent to a parametrisation of the weights given as

a:b:c:d=\operatorname{snh}(η-u):\operatorname{snh}(η+u):\operatorname{snh}():k\operatorname{snh}()\operatorname{snh}(η-u)\operatorname{snh}(η+u)

for fixed modulus

k

and

η

and variable

u

. Here snh is the hyperbolic analogue of sn, given by

\begin{align} \operatorname{snh}(u)&=-i\operatorname{sn}(iu)=i\operatorname{sn}(-iu)\\ where\operatorname{sn}(u)&=

H(u)
k1/2\Theta(u)

\end{align}

and

H(u)

and

\Theta(u)

are Theta functions of modulus

k

. The associated transfer matrix

T

thus is a function of

u

alone; for all

u

,

v

T(u)T(v)=T(v)T(u).

The matrix function

Q(u)

The other crucial part of the solution is the existence of a nonsingular matrix-valued function

Q

, such that for all complex

u

the matrices

Q(u),Q(u')

commute with each other and the transfer matrices, and satisfywhere

\begin{align} \zeta(u)&=[c-1H()\Theta(u)\Theta(u)]N\\ \phi(u)&=[\Theta(0)H(u)\Theta(u)]N. \end{align}

The existence and commutation relations of such a function are demonstrated by considering pair propagations through a vertex, and periodicity relations of the theta functions, in a similar way to the six-vertex model.

Explicit solution

The commutation of matrices in allow them to be diagonalised, and thus eigenvalues can be found. The partition function is calculated from the maximal eigenvalue, resulting in a free energy per site of

\begin{align} f=\epsilon5-2kT\sum

infty
n=1
\sinh2((\tau)n)(\cosh(nλ)-\cosh(n\alpha))
n\sinh(2n\tau)\cosh(nλ)

\end{align}

for
\begin{align} \tau&=\piK'\\ λ&=
2K
\piη\\ \alpha&=
iK
\piu
iK

\end{align}

where

K

and

K'

are the complete elliptic integrals of moduli

k

and

k'

.The eight vertex model was also solved in quasicrystals.

Equivalence with an Ising model

There is a natural correspondence between the eight-vertex model, and the Ising model with 2-spin and 4-spin nearest neighbour interactions. The states of this model are spins

\sigma=\pm1

on faces of a square lattice. The analogue of 'edges' in the eight-vertex model are products of spins on adjacent faces:

\begin{align} \alphaij&=\sigmaij\sigmai,j+1\\ \muij&=\sigmaij\sigmai+1,j. \end{align}

The most general form of the energy for this model is

\begin{align} \epsilon&=-\sumij(Jh\muij+Jv\alphaij+J\alphaij\muij+J'\alphai+1,j\muij+J''\alphaij\alphai+1,j) \end{align}

where

Jh

,

Jv

,

J

,

J'

describe the horizontal, vertical and two diagonal 2-spin interactions, and

J''

describes the 4-spin interaction between four faces at a vertex; the sum is over the whole lattice.

We denote horizontal and vertical spins (arrows on edges) in the eight-vertex model

\mu

,

\alpha

respectively, and define up and right as positive directions. The restriction on vertex states is that the product of four edges at a vertex is 1; this automatically holds for Ising 'edges'. Each

\sigma

configuration then corresponds to a unique

\mu

,

\alpha

configuration, whereas each

\mu

,

\alpha

configuration gives two choices of

\sigma

configurations.

Equating general forms of Boltzmann weights for each vertex

j

, the following relations between the

\epsilonj

and

Jh

,

Jv

,

J

,

J'

,

J''

define the correspondence between the lattice models:

\begin{align} \epsilon1&=-Jh-Jv-J-J'-J'',\epsilon2=Jh+Jv-J-J'-J''\\ \epsilon3&=-Jh+Jv+J+J'-J'',\epsilon4=Jh-Jv+J+J'-J''\\ \epsilon5&=\epsilon6=J-J'+J''\\ \epsilon7&=\epsilon8=-J+J'+J''. \end{align}

It follows that in the zero-field case of the eight-vertex model, the horizontal and vertical interactions in the corresponding Ising model vanish.

These relations gives the equivalence

ZI=2Z8V

between the partition functions of the eight-vertex model, and the 2,4-spin Ising model. Consequently a solution in either model would lead immediately to a solution in the other.

See also

Notes and References

  1. Sutherland . Bill . Two‐Dimensional Hydrogen Bonded Crystals without the Ice Rule . Journal of Mathematical Physics . AIP Publishing . 11 . 11 . 1970 . 0022-2488 . 10.1063/1.1665111 . 3183–3186. 1970JMP....11.3183S .
  2. Fan . Chungpeng . Wu . F. Y. . General Lattice Model of Phase Transitions . Physical Review B . American Physical Society (APS) . 2 . 3 . 1970-08-01 . 0556-2805 . 10.1103/physrevb.2.723 . 723–733. 1970PhRvB...2..723F .
  3. Baxter . R. J. . Eight-Vertex Model in Lattice Statistics . Physical Review Letters . American Physical Society (APS) . 26 . 14 . 1971-04-05 . 0031-9007 . 10.1103/physrevlett.26.832 . 832–833. 1971PhRvL..26..832B .