text:
# Introduction to solvable lattice models in statistical and mathematical physics 11footnote 1 In ”Classical and Quantum Integrable Systems: Theory and Applications” (ISBN 07503 09598), Institute of Physics Publishing, (2003) http://bookmarkphysics.iop.org/ (Chap. 5 , pp. 113-151.)

## 1 Introduction

We introduce the six-vertex model defined on a two-dimensional square lattice. We describe the model in detail, since it gives an important prototype of many solvable lattice models defined on two-dimensional lattices . The transfer matrix of the six-vertex model generalizes the XXZ quantum spin chain which plays a central role among integrable quantum spin chains . The eight-vertex model, which generalizes the six-vertex model directly, may be considered as the most important exactly solvable model in statistical mechanics . Moreover, many mathematical theories such as the algebraic Bethe ansatz  and quantum groups  are closely related to the six-vertex model. Starting from the six-vertex model, one may have a wide viewpoint on various physical and mathematical topics related to solvable models. There are quite a large number of topics related to exactly solvable models in physics and mathematics .

We explain in §2 some features of the six-vertex model defined on a square lattice. We introduce the Boltzmann weights and the transfer matrix for the six-vertex model. We review a method for diagonalizing the transfer matrix, which is called the coordinate Bethe ansatz, and we show the expressions of the free energy per site in the ferroelectric, the anti-ferroelectric, and the disordered phases, respectively . The disordered phase is gapless, while the ferroelectric and the anti-ferroelectric phases have gaps . We derive a critical singularity appearing at the phase transition from the anti-ferroelectric to the disordered phases. We review the calculation of the singular part of the free energy through the analytic continuation, as shown in Ref. . The critical singularity is very weak and has the essential singularity similar to the Kosterlitz-Thouless transition. We have thus derived the KT-like singularity through exact calculation. After reviewing the finite-size analysis of conformal invariance , we discuss that the massless phase of the six-vertex model is related to the conformal field theory of $`c=1`$ which has $`U(1)`$ symmetry. The $`c=1`$ CFT has the critical line where critical exponents change continuously with respect to some parameter of the model . There are quite a few papers on the finite-size corrections of the integrable models . (For a review, see .) The critical line is also characteristic to the Tomonaga-Luttinger liquid . The existence of a critical line was first discovered by R.J. Baxter through the exact solution of the eight-vertex model .

In §3, we review various integrable models in statistical mechanics . We briefly introduce the Ising model , the Potts models , and the chiral Potts model  and then the eight-vertex model  and IRF models . In §4, we solve explicitly the Yang-Baxter equations for the six-vertex model. We introduce the algebraic Bethe ansatz . Here we show that the Yang-Baxter equation of the algebraic Bethe ansatz can be expressed by graphs. In §5, we discuss some mathematical theories associated with integrable models such as the braid group  and the quantum groups .

There have been novel developments in the mathematical physics associated with the six-vertex model . The integrable vertex models associated with various Lie algebras have been obtained, which generalize the six-vertex model . The crystal basis of the quantum groups is derived from the mathematical analysis of the corner transfer matrix which is fundamental for calculating the one-point functions of the vertex models and IRF models . Through the $`q`$-vertex operators, correlation functions of the XXZ spin chain or the six vertex model are obtained . The dynamical Yang-Baxter equation  and the elliptic quantum groups  have also been extensively discussed. In fact, we can derive the $`R`$ matrix of the eight-vertex model systematically from the elliptic quantum group through the twists . Furthermore, the correlation functions of the XXZ model calculated with the $`q`$-vertex operators have been rederived for large but finite chains through the algebraic Bethe ansatz with the Drinfeld twists . Here we note that the $`q`$-vertex operator can be defined only on the infinite chain, while the algebraic Bethe ansatz with the Drinfeld twists can be applied to any finite chain. By taking the thermodynamic limit, it has been shown that the two approaches indeed give the same results. These papers indeed illustrate nontrivial physical applications of the Drinfeld twits. It is recently found that the symmetry of the six-vertex model is enhanced at some particular coupling constants: the transfer matrix commutes with the generators of the $`sl_2`$ loop algebra for the six-vertex model at roots of unity .

Let us discuss some physical motivations for the six-vertex model. The exact solution of the six-vertex model was originally introduced for studying the statistical mechanics of ferro-electrics such as the residual entropy and the ferroelectric transitions . However, it seems that the physical motivation of the six-vertex model for ferro-electricity has decreased. On the other hand, there are many different physical applications of of the six-vertex model. Here we consider a few examples, the domain wall theory , crystal growth  and the thermodynamics of the XXZ spin chain through the quantum transfer matrix . The crystal growth on surfaces has been discussed by applying exact solutions of the six-vertex model  and some extensions . The free energy of the six vertex model gives the equilibrium crystal shapes . The finite-temperature thermodynamics of the XXZ spin chain has been studied extensively through the quantum transfer matrix, which is a version of the inhomogeneous six-vertex transfer matrix . The quantum transfer matrix is obtained by regarding one direction of the square lattice as the imaginary time or the inverse temperature . There have been considerable efforts to evaluate thermal quantities analytically or numerically. Several functional equations on the eigenvalues of the transfer matrix have been devised . Finally, we remark that a universal relation between the dispersion curve and the ground-state correlation length in quantum spin chains are discussed by using the exact solutions of the vertex models .

In §2, we employ mainly the notation of Baxter’s textbook  except for the transfer matrix. In §4, however, we briefly show how the notation in statistical mechanics is related to the notation of the algebraic Bethe ansatz or the Quantum Inverse Scattering Method. The graphical illustration should be useful. Finally in §5 we discuss many connections of the six-vertex model to several mathematical developments such as the quantum groups.

## 2 Solvable vertex models

### 2.1 The six-vertex model

#### 2.1.1 Ice rule

Let us consider a square lattice as a model of two-dimensional ferroelectric crystal. Molecules are placed on the vertices of the lattice. Arrows are placed on the edges of the lattice, which correspond to directions of dipole moments of hydrogen bonds. As a crystal of hydrogen bonding, we may consider ice, i.e., the crystal of water molecules. In this review, however, we simplify the molecular background of the model (For instance, see ). We assume that the dipole moments defined on edges take only two values $`\pm 1`$.

At a vertex in the lattice, there are four edges. There are possible 16 configurations of the four edges around the vertex since each of the edges takes two values $`\pm 1`$.

Let symbols $`\alpha ,\beta ,\gamma ,\delta `$ denote the values of dipole moments around the vertex. Due to charge neutrality, they should satisfy the following condition

$$\alpha +\beta =\gamma +\delta $$
(2.1)
For an illustration, let us consider the case when $`\alpha ,\beta ,\gamma `$ and $`\delta `$ are given by +1. Then, $`\alpha `$ and $`\beta `$ give + 2 to the vertex, while $`\gamma `$ and $`\delta `$ deprive + 2 from it, so that the net charge around the vertex is kept neutral: $`\alpha +\beta \gamma \delta =0`$.

There are only 6 configurations satisfying the condition. The other configurations that do not satisfy the condition are not allowed in thermal equilibrium. Thus, we call the model the six-vertex model. The condition (2.1) is sometimes called ice rule, since ice as a crystal consists of water molecules connected by hydrogen bonding.

We denote by 1 and 2 the values of polarization 1 and $`1`$, respectively. The symbols 1 and 2 are useful for matrix notation. Let $`p`$ denote the notation of $`\pm 1`$ and $`k`$ 1 and 2. Then, they are related by the relation: $`k=1+(1p)/2`$.

#### 2.1.2 Boltzmann weights

It is a key idea in exactly solvable models that we define the model by the Boltzmann weights not by the energies of configurations. Let us introduce the Boltzmann weights for configurations around a vertex. For a vertex configuration $`\alpha ,\beta ,\gamma ,\delta `$, we denote by $`ϵ(\alpha ,\beta |\gamma ,\delta )`$ the energy at the vertex. Then, the Boltzmann weight for a temperature $`T`$ is given by

$$w(\alpha ,\beta |\gamma ,\delta )=\mathrm{exp}(ϵ(\alpha ,\beta |\gamma ,\delta )/k_BT)$$
(2.2)
Under the ice rule, there are only six configurations allowed round a vertex. Here, it is assumed that the energy of a configuration violating the ice rule should be infinite. We denote by $`ϵ_j`$ the energy of the $`j`$th vertex configuration shown in Fig. 2.2.

Under no external field, the Boltzmann weights must be invariant when reversing all the polarizations simultaneously. Thus, we have $`ϵ_1=ϵ_2`$, $`ϵ_3=ϵ_4`$ and $`ϵ_5=ϵ_6`$, when there is no external field. We denote the Boltzmann weights as follows.

$`w(1,1|1,1)`$ $`=`$ $`w(2,2|2,2)=w_1=a`$
$`w(1,2|2,1)`$ $`=`$ $`w(2,1|1,2)=w_2=b`$
$`w(1,2|1,2)`$ $`=`$ $`w(2,1|2,1)=w_3=c`$ (2.3)
The Boltzmann weights of the zero-field six-vertex model have essentially only two parameters. For instance, we may choose $`a/c`$ and $`b/c`$. Note that the probability for the vertex configuration of $`a`$ is given by $`a/(a+b+c)`$, which does not change by replacing $`a,b`$ and $`c`$ with $`\rho a`$, $`\rho b`$ and $`\rho c`$.

Let us consider $`\pi /2`$ rotation of the square lattice. If we rotate vertex configuration (1) of Fig. 2.2 by the angle $`\pi /2`$ in the counterclockwise direction, then, it becomes vertex configuration (4). Under the $`\pi /2`$ rotation, the weight $`a`$ is exchanged with the weight $`b`$, while the weight $`c`$ does not change.

### 2.2 Partition function and the transfer matrix

Let us discuss the partition function of the system. We now set the boundary conditions. Here, we consider the periodic boundary conditions for the two-dimensional lattice. We take a product of the Boltzmann weights over all the vertices of the lattice, and sum up the product over all the allowed configurations of arrows on the lattice.

$$Z=\underset{config}{}\underset{j:vertices}{}w(a_j,b_j|c_j,d_j)$$
(2.4)
The partition function of the square lattice can be formulated as the trace of the products of the transfer matrices. Let us define the transfer matrix $`\tau `$ of the six-vertex model. The matrix elements of the transfer matrix $`\tau `$ acting on $`N`$ lattice sites are given by

$$\tau _{b_1,\mathrm{},b_N}^{a_1,\mathrm{},a_N}=\underset{c_1,\mathrm{},c_N}{}w(c_1,b_1|a_1,c_2)w(c_2,b_2|a_2,c_3)\mathrm{}w(c_N,b_N|a_N,c_1)$$
(2.5)
Under the periodic boundary conditions, the partition function $`Z_{NN^{^{}}}`$ of $`N\times N^{^{}}`$ lattice is given by the trace of the $`N^{^{}}`$th power of the transfer matrices:

$$Z_{NN^{^{}}}=\mathrm{Tr}\left(\tau ^N^{^{}}\right)=\underset{a_1,\mathrm{},a_N}{}\left(\tau ^N^{^{}}\right)_{a_1,\mathrm{},a_N}^{a_1,\mathrm{},a_N}=\mathrm{\Lambda }_1^N^{^{}}+\mathrm{\Lambda }_2^N^{^{}}+\mathrm{}+\mathrm{\Lambda }_{2^N}^N^{^{}}$$
(2.6)
Here $`\mathrm{\Lambda }_j`$ denotes the eigenvalue of the transfer matrix $`\tau `$.

The free energy per site $`f`$ is given by

$$f=k_BT\mathrm{log}Z_{NN^{^{}}}/(NN^{^{}}).$$
(2.7)
In the thermodynamic limit: $`N,N^{^{}}\mathrm{}`$, the free energy per site is given by the largest eigenvalue $`\mathrm{\Lambda }_{max}`$ of the transfer matrix $`\tau `$.

### 2.3 Diagonalization of the transfer matrix

#### 2.3.1 The Yang-Baxter relations for six-vertex model

Let us consider three sets of the Boltzmann weights: $`(w_1,w_2,w_3)=(a,b,c)`$, $`(a^{^{}},b^{^{}},c^{^{}})`$, and $`(a^{^{\prime \prime }},b^{^{\prime \prime }},c^{^{\prime \prime }})`$. We denote by $`\tau ^{^{}}`$ and $`\tau ^{^{\prime \prime }}`$ the transfer matrices constructed from the sets of the Boltzmann weights $`(a^{^{}},b^{^{}},c^{^{}})`$ and $`(a^{^{\prime \prime }},b^{^{\prime \prime }},c^{^{\prime \prime }})`$, respectively. If the three sets of the Boltzmann weights satisfy the Yang-Baxter equation

$`{\displaystyle \underset{\alpha ,\beta ,\gamma }{}}w(\alpha ,\gamma |a_1,a_2)w^{^{}}(\beta ,b_3|\gamma ,a_3)w^{^{\prime \prime }}(b_1,b_2|\alpha ,\beta )`$ (2.8)
$`=`$ $`{\displaystyle \underset{\alpha ,\beta ,\gamma }{}}w^{^{\prime \prime }}(\beta ,\alpha |a_2,a_3)w^{^{}}(b_1,\gamma |a_1,\beta )w(b_2,b_3|\gamma ,\alpha ),`$
then the transfer matrices $`\tau ^{^{}}`$ and $`\tau ^{^{\prime \prime }}`$ commute. The derivation of the commutation relation is given in Appendix A. We note that a graphical presentation of the Yang-Baxter equation (2.8) will be shown in Fig. 4.1.

Let us define the parameter $`\mathrm{\Delta }`$ as follows

$$\mathrm{\Delta }=\frac{a^2+b^2c^2}{2ab}.$$
(2.9)
For the zero-field six-vertex model, we can show that if the two sets of the Boltzmann weights have the same value of the parameter $`\mathrm{\Delta }`$, then their transfer matrices commute. We shall explicitly discuss in §4 that it is indeed derived from the Yang-Baxter equations (2.8).

#### 2.3.2 The coordinate Bethe ansatz

Let us consider the matrix element $`\tau _{b_1,\mathrm{}b_N}^{a_1,\mathrm{},a_N}`$ of the transfer matrix $`\tau `$. Due to the ice rule, we may express the suffix $`a_1,\mathrm{},a_N`$ by the positions of the value 2, as follows. Suppose that there are $`n`$ suffices given by the value 2 among the $`N`$ suffices $`a_1,\mathrm{},a_N`$. The $`n`$ suffices are expressed as $`a_{x_1},a_{x_2},\mathrm{},a_{x_n}`$ where $`x_j`$’s are in increasing order: $`x_1<x_2<\mathrm{}<x_n`$ . Then, the entry $`a_1,\mathrm{},a_N`$ is equivalent to the set of the $`x_j`$’s: $`x_1,\mathrm{},x_n`$. For an illustration, let us consider the case $`N=5`$ and $`n=3`$. Then, $`(x_1,x_2,x_3)=(1,3,4)`$ corresponds to $`(a_1,a_2,a_3,a_4,a_5)=(2,1,2,2,1)`$. Thus, the matrix element $`\tau _{b_1,\mathrm{},b_N}^{a_1\mathrm{}a_N}`$ can be denoted briefly as $`\tau _{y_1,\mathrm{},y_n}^{x_1,\mathrm{},x_n}`$.

Let us now discuss how to solve the secular equation: $`\tau g=\mathrm{\Lambda }g`$. Here, the transfer matrix $`\tau `$ is a $`2^N\times 2^N`$ matrix, $`g`$ is a $`2^N`$-dimensional eigenvector with eigenvalue $`\mathrm{\Lambda }`$. In terms of matrix elements, the secular equation is written as

$$\underset{y_1,\mathrm{},y_n}{}\tau _{y_1,\mathrm{},y_n}^{x_1,\mathrm{},x_n}g(y_1,\mathrm{},y_n)=\mathrm{\Lambda }g(x_1,\mathrm{},x_n)$$
(2.10)
Here, $`g(x_1,\mathrm{},x_n)`$ denotes the matrix element of the vector $`g`$ for the entry $`(x_1,\mathrm{},x_n)`$. In the coordinate Bethe ansatz, we assume the following form for the matrix element of the possible eigenvector $`g`$

$$g(x_1,\mathrm{},x_n)=\underset{PS_n}{}A_P\mathrm{exp}\left(k_{P1}x_1+\mathrm{}+k_{Pn}x_n\right)$$
(2.11)
Here, $`S_n`$ denotes the symmetric group of order $`n`$, and $`P`$ is a permutation of $`n`$ letters, $`1,2,\mathrm{},n`$, where $`P`$ maps $`j`$ into $`Pj`$. The expression (2.11) is called the Bethe ansatz wave-function. If the vector $`g`$ whose elements are of the form (2.11) is an eigenvector of the transfer matrix, then we call it a Bethe ansatz eigenvector.

For general $`n`$, the vector (2.11) becomes an eigenvector of the transfer matrix, if the wave-numbers $`k_j`$’s satisfy the Bethe ansatz equations. They are given by the following

$$\mathrm{exp}(iNk_j)=(1)^{n1}\underset{\mathrm{}=1}{\overset{n}{}}\mathrm{exp}(i\mathrm{\Theta }(k_j,k_{\mathrm{}})),\mathrm{for}j=1,\mathrm{},n,$$
(2.12)
where $`\mathrm{\Theta }(p,q)`$ is defined by

$$\mathrm{exp}(i\mathrm{\Theta }(p,q))=\frac{12\mathrm{\Delta }e^{ip}+e^{i(p+q)}}{12\mathrm{\Delta }e^{iq}+e^{i(p+q)}}$$
(2.13)
For the solutions $`k_j`$’s to the Bethe ansatz equations, the eigenvalue $`\mathrm{\Lambda }`$ of the transfer matrix is given by

$$\mathrm{\Lambda }(k_1,\mathrm{},k_n)=a^NL(z_1)\mathrm{}L(z_n)+b^NM(z_1)\mathrm{}M(z_n),$$
(2.14)
where $`z_j=\mathrm{exp}(ik_j)`$ for $`j=1,\mathrm{},n`$ and the functions $`L(z)`$ and $`M(z)`$ are defined by

$$L(z)=\frac{ab+(c^2b^2)z}{a(abz)},M(z)=\frac{a^2c^2abz}{b(abz)}.$$
(2.15)
When we discuss the spectrum of an integrable model through the coordinate Bethe ansatz, we often assume that all the eigenvectors of the transfer matrix are characterized by the Bethe ansatz wave-function (2.11). However, it is not certain whether the assumption is valid or not. Thus, we have to check it by other methods. In fact, there are several numerical studies on the validity of the completeness of the Bethe ansatz for some integrable models.

On the other hand, there is no doubt on the mathematical structure of the Bethe ansatz wave-function. We can derive the expression (2.11) by the algebraic Bethe ansatz through the ‘two-site’ model . (There is an instructive notice in .) It was shown that the matrix elements of the product of $`B`$ operators acting on the vacuum are given by the Bethe ansatz wave-function (2.11) with $`k_j`$’s being generic.

The Yang-Baxter relation leads to not only the integrability of the six-vertex model but also the systematic construction of the eigenvectors. In fact, we shall see in §4 that the algebraic Bethe ansatz is solely based on the Yang-Baxter equation.

#### 2.3.3 An example of the eigenvector

For an illustration, let us consider the eigenvector $`g`$ for the case of $`n=1`$: $`g(x)=A\mathrm{exp}(ikx)`$. Through a direct calculation, we have

$`{\displaystyle \underset{y=1}{\overset{N}{}}}\tau _y^xg(y)`$ $`=`$ $`{\displaystyle \underset{y=1}{\overset{x1}{}}}a^{xy1}b^{Nx+y1}c^2g(y)+{\displaystyle \underset{y=x+1}{\overset{N}{}}}a^{N+xy1}b^{yx1}c^2g(y)`$ (2.16)
$`=`$ $`\left(a^NL(z)+b^NM(z)\right)g(x)+{\displaystyle \frac{a^{N1}b^{Nx}c^2z}{abz}}(1z^N)`$
Therefore, $`g(x)`$ becomes an eigenvector if the Bethe ansatz equation $`z^N=1`$ is satisfied.

### 2.4 The free energy of the six-vertex model

#### 2.4.1 Three phases of the six-vertex model

There are three phases for the zero-field six-vertex model. They are given by the regions of the parameter $`\mathrm{\Delta }`$: ferroelectric phase when $`\mathrm{\Delta }>1`$; anti-ferroelectric phase when $`\mathrm{\Delta }<1`$; disordered phase when $`1<\mathrm{\Delta }<1`$. It is found that the disordered phase ($`1<\mathrm{\Delta }<1`$) is gapless (massless), while the ferroelectric phase ($`\mathrm{\Delta }>1`$) and the anti-ferroelectric phase ($`\mathrm{\Delta }<1`$) are gapful (massive).

Let us consider the phase diagram of the six vertex model shown in Fig. 2.4. We recall that the ratios $`a/c`$ and $`b/c`$ determine the model. Here we note that the number of independent parameters is given by two, since the overall normalization factor is arbitrary. The regimes I and II give the ferroelectric phase, the regimes III and IV are anti-ferroelectric and disordered, respectively. In terms of the Boltzmann weights, the regime I is given by $`a>b+c`$, the regime II by $`b>a+c`$, and the regime IV by $`a+b<c`$ . Then, the regime III is given by $`a<b+c`$, $`b<c+a`$ and $`c<a+b`$.

Let us derive the three phases through an intuitive argument. For the ferroelectric regime of $`\mathrm{\Delta }>1`$, we have $`a^2+b^2c^2>2ab`$ , which leads to the inequality: $`|ab|>c`$. When $`a>b`$, we have $`a>b+c`$. Thus, the configuration where all the Boltzmann weights are given by the weight $`a`$ should be the largest contribution to the partition function $`Z`$. In fact, when $`n=0`$ in eq. (2.14), we have $`\mathrm{\Lambda }=a^N+b^N`$. When $`a>b`$, the free energy per site, $`f`$, is given by $`ϵ_1`$: $`f=k_BT\mathrm{log}a`$. When $`b>a`$, we have $`f=ϵ_3`$ , similarly.

For the anti-ferroelectric regime of $`\mathrm{\Delta }<1`$, we have the inequality $`a+b<c`$. Thus, the vertex configurations for $`c`$ should be more favorable than those of $`a`$ or $`b`$. In fact, it is shown that the transfer matrix has the largest eigenvalue when $`n=N/2`$. Furthermore, if we send $`|\mathrm{\Delta }|`$ to infinity, all the vertex configurations should be given by those of $`c`$. The phase is thus called anti-ferroelectric.

We may explain the reason why it is called anti-ferroelectric. Let us consider the configuration (5) of Fig. 2.2. We see that the arrow coming from the left goes upward, while the arrow coming from the right goes downward. If all the vertex configurations on the square lattice are given by (5), then the lines coming from South West to North East and the lines coming from North East to South West occupy the lattice alternatively. This gives the anti-ferroelectric order.

#### 2.4.2 Parametrization of the Boltzmann weights

It is nontrivial how to parametrize the Boltzmann weights. Recall that there are two independent parameters for the six-vertex model. Thus, if we consider $`\mathrm{\Delta }`$ as a parameter, there is only another one. As we shall see later, it is related to the spectral parameter.

We recall that the phases of the zero-field six-vertex model are classified into the following: $`\mathrm{\Delta }<1`$, $`1<\mathrm{\Delta }<1`$ and $`1<\mathrm{\Delta }`$.

(1) Anti-ferroelectric phase

For $`\mathrm{\Delta }<1`$, we define a real parameter $`\lambda `$ by

$$\mathrm{\Delta }=\mathrm{cosh}\lambda (0<\lambda )$$
(2.17)
and we parametrize the Boltzmann weights as follows

$$a=\rho \mathrm{sinh}(\frac{\lambda v}{2}),b=\rho \mathrm{sinh}(\frac{\lambda +v}{2}),c=\rho \mathrm{sinh}\lambda ,(\lambda <v<\lambda )$$
(2.18)
where $`\rho `$ is the normalization factor. Let us define the rapidity $`\alpha `$ for the wavenumber $`k`$ as

$$\mathrm{exp}(ik)=\frac{e^\lambda e^{i\alpha }}{e^{\lambda i\alpha }1}=\frac{\mathrm{sin}\frac{1}{2}(\alpha i\lambda )}{\mathrm{sin}\frac{1}{2}(\alpha +i\lambda )}$$
(2.19)
By replacing the wave numbers $`p`$ and $`q`$ with the rapidities $`\alpha `$ and $`\beta `$ in eq. (2.13), the phase factor $`\mathrm{\Theta }(p,q)`$ is written as follows

$$\mathrm{exp}(i\mathrm{\Theta }(p,q))=\frac{e^{2\lambda }e^{i(\alpha \beta )}}{e^{2\lambda i(\alpha \beta )}1}=\frac{\mathrm{sin}\frac{1}{2}((\alpha \beta )2i\lambda )}{\mathrm{sin}\frac{1}{2}((\alpha \beta )+2i\lambda )}$$
(2.20)
(2) Disordered phase

For $`1<\mathrm{\Delta }<1`$, we define a positive real parameter $`\mu `$ by

$$\mathrm{\Delta }=\mathrm{cos}\mu (0<\mu <\pi )$$
(2.21)
and we parameterize the Boltzmann weights as follows

$$a=\rho \mathrm{sin}(\frac{\mu w}{2}),b=\rho \mathrm{sin}(\frac{\mu +w}{2}),c=\rho \mathrm{sin}\mu ,(\mu <w<\mu ).$$
(2.22)
Here $`\rho `$ is the normalization factor. We define the rapidity $`\alpha `$ for the wavenumber $`k`$ by

$$\mathrm{exp}(ik)=\frac{e^{i\mu }e^\alpha }{e^{i\mu +\alpha }1}=\frac{\mathrm{sinh}\frac{1}{2}(\alpha i\mu )}{\mathrm{sinh}\frac{1}{2}(\alpha +i\mu )}$$
(2.23)
In terms of rapidities $`\alpha `$ and $`\beta `$, the phase factor $`\mathrm{\Theta }(p,q)`$ is expressed as

$$\mathrm{exp}(i\mathrm{\Theta }(p,q))=\frac{e^{2i\mu }e^{\alpha \beta }}{e^{2i\mu +\alpha \beta }1}=\frac{\mathrm{sinh}\frac{1}{2}((\alpha \beta )2i\mu )}{\mathrm{sinh}\frac{1}{2}((\alpha \beta )+2i\mu )}$$
(2.24)
(3) Ferroelectric phase

For $`\mathrm{\Delta }>1`$, we define the real parameter $`\lambda `$ by

$$\mathrm{\Delta }=\mathrm{cosh}\lambda (0<\lambda )$$
(2.25)
and we may parameterize the Boltzmann weights as follows: When $`a>b`$ (the regime I),

$$a=\rho \mathrm{sinh}(\frac{\lambda v}{2}),b=\rho \mathrm{sinh}(\frac{\lambda +v}{2}),c=\rho \mathrm{sinh}(\lambda ),(v<\lambda ).$$
(2.26)
When $`a<b`$ (the regime II),

$$a=\rho \mathrm{sinh}(\frac{\lambda v}{2}),b=\rho \mathrm{sinh}(\frac{\lambda +v}{2}),c=\rho \mathrm{sinh}(\lambda ),(v>\lambda ).$$
(2.27)
Here we recall that $`\rho `$ is the normalization factor. The wavenumber $`k`$ is related to the rapidity $`\alpha `$ by

$$\mathrm{exp}(ik)=\frac{e^\lambda e^{i\alpha }}{e^{\lambda i\alpha }1}=\frac{\mathrm{sin}\frac{1}{2}(\alpha i\lambda )}{\mathrm{sin}\frac{1}{2}(\alpha +i\lambda )}$$
(2.28)
The phase factor $`\mathrm{\Theta }(p,q)`$ is written as

$$\mathrm{exp}(i\mathrm{\Theta }(p,q))=\frac{e^{2\lambda }e^{i(\alpha \beta )}}{e^{2\lambda i(\alpha \beta )}1}=\frac{\mathrm{sin}\frac{1}{2}((\alpha \beta )2i\lambda )}{\mathrm{sin}\frac{1}{2}((\alpha \beta )+2i\lambda )}$$
(2.29)
Let us consider curved lines given by changing the spectral parameter $`w`$ or $`v`$ continuously. All the curves pass through the points $`(a/c,b/c)`$=$`(1,0)`$ and $`(0,1)`$. Except for the two points, however, any points on the horizontal axis: $`b/c=0`$ or the vertical axis: $`a/c=0`$ are never reached by the above parametrizations with finite values.

#### 2.4.3 Expressions of the free energy

(1) Anti-ferroelectric phase

When $`\mathrm{\Delta }<1`$, the system is in the anti-ferroelectric phase. The free energy per site $`f`$ is given by

$`f`$ $`=`$ $`k_BT\mathrm{log}ak_BT\left({\displaystyle \frac{\lambda +v}{2}}+{\displaystyle \underset{m=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{e^{m\lambda }\mathrm{sinh}m(\lambda +v)}{m\mathrm{cosh}m\lambda }}\right)`$ (2.30)
$`=`$ $`k_BT\mathrm{log}bk_BT\left({\displaystyle \frac{\lambda v}{2}}+{\displaystyle \underset{m=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{e^{m\lambda }\mathrm{sinh}m(\lambda v)}{m\mathrm{cosh}m\lambda }}\right)`$ (2.31)
for $`\lambda <v<\lambda `$.

(2) Disordered phase

When $`1<\mathrm{\Delta }<1`$, we have

$`f`$ $`=`$ $`k_BT\mathrm{log}ak_BT{\displaystyle _{\mathrm{}}^{\mathrm{}}}{\displaystyle \frac{\mathrm{sinh}(\mu +w)x\mathrm{sinh}(\pi \mu )x}{2x\mathrm{sinh}\pi x\mathrm{cosh}\mu x}}𝑑x`$ (2.32)
$`=`$ $`k_BT\mathrm{log}bk_BT{\displaystyle _{\mathrm{}}^{\mathrm{}}}{\displaystyle \frac{\mathrm{sinh}(\mu w)x\mathrm{sinh}(\pi \mu )x}{2x\mathrm{sinh}\pi x\mathrm{cosh}\mu x}}𝑑x`$ (2.33)
for $`\mu <w<\mu `$.

(3) Ferroelectric phase

When $`a>b`$ (the regime I), we have $`f=k_BT\mathrm{log}a`$, while when $`b>a`$ (the regime II), we have $`f=k_BT\mathrm{log}b`$.

#### 2.4.4 Correlation length

Correlation length has been calculated for the six-vertex model . In fact, it is calculated for the eight-vertex model. In the ferroelectric and the anti-ferroelectric regimes of the six-vertex model, the correlation length is finite. It becomes very large near their phase boundaries to the disordered phase. In the disordered phase of the six-vertex model, the correlation length diverges throughout the regime. Thus, the disordered phase is critical.

### 2.5 Critical singularity in the anti-ferro regime near the phase boundary

We now discuss a singular behavior of the free energy shown at the phase transition from the anti-ferroelectric to the disordered phases . Here we note that the former one is gapful, while the latter is gapless. Recall that the phase boundary between the regimes III and IV is given by $`a+b=c`$. When $`T>T_c`$, the system should be in the disordered regime ($`1<\mathrm{\Delta }<1)`$, while when $`T<T_c`$ it is in the anti-ferroelectric regime ($`\mathrm{\Delta }<1)`$. In the lower temperature, the system should be ordered.

Let us calculate the analytic continuation of the high-temperature free energy (2.32) into the low-temperature phase, so that we can single out the singularity of the free energy near $`T_c`$ and for $`T<T_c`$. First, we reformulate the integral of eq. (2.32) as follows.

$$_{\mathrm{}}^{\mathrm{}}\frac{\mathrm{sinh}(\mu +w)x\mathrm{sinh}(\pi \mu )x}{2x\mathrm{sinh}\pi x\mathrm{cosh}\mu x}𝑑x=𝒫_{\mathrm{}}^{\mathrm{}}\frac{\mathrm{sinh}(\mu +w)x\mathrm{exp}(\pi \mu )x}{2x\mathrm{sinh}\pi x\mathrm{cosh}\mu x}𝑑x$$
(2.34)
Here $`𝒫`$ denote the principal value integral. Second, we set $`\lambda `$ to be a very small positive real number, and take a value $`v`$ satisfying $`\lambda <v<\lambda `$. We consider the path in the complex $`\mu `$-plane:

$$\mu =\lambda \mathrm{exp}(i\theta ),\mathrm{for}0\theta \pi /2.$$
(2.35)
Along the path, we calculate the analytic continuation of the high-temperature free energy (2.32). Here, we also consider the path of $`w`$: $`w=v\mathrm{exp}(i\theta )`$ for $`0\theta \pi /2`$. Then, we have

$`𝒫{\displaystyle _{\mathrm{}}^{\mathrm{}}}{\displaystyle \frac{i\mathrm{sin}(\lambda +v)x\mathrm{exp}(\pi +i\lambda )x}{2x\mathrm{sinh}\pi x\mathrm{cos}\lambda x}}𝑑x`$ (2.36)
$`=`$ $`{\displaystyle \frac{\lambda +v}{2}}+{\displaystyle \underset{m=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{e^{m\lambda }\mathrm{sinh}m(\lambda +v)}{m\mathrm{cosh}m\lambda }}i{\displaystyle \underset{m=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(1)^m\mathrm{cos}\left((m1/2)\pi v/\lambda \right)e^{(m1/2)\pi ^2/\lambda }}{(m1/2)\mathrm{sinh}\left((m1/2)\pi ^2/\lambda \right)}}`$
The real part of the analytic continuation corresponds to the expression of the anti-ferroelectric free energy. Therefore, we obtain the singular part of the free energy

$$f_{sing}=ik_BT\underset{m=1}{\overset{\mathrm{}}{}}\frac{(1)^m\mathrm{cos}\left((m1/2)\pi v/\lambda \right)\mathrm{exp}\left((m1/2)\pi ^2/\lambda \right)}{(m1/2)\mathrm{sinh}\left((m1/2)\pi ^2/\lambda \right)}$$
(2.37)
We define the reduced temperature $`t`$ by

$$t=(a+bc)/c.$$
(2.38)
In the low-temperature phase ($`T<T_c`$) and near $`T_c`$, $`t`$ is given by

$$t\frac{1}{8}(\lambda ^2v^2)$$
(2.39)
Approximately, $`t`$ is given by $`t\lambda ^2/8`$. Near $`T_c`$, we have

$$f_{sing}4ik_BTe^{\pi ^2/\lambda }\mathrm{cos}(\frac{\pi v}{2\lambda })$$
(2.40)
Thus, we have

$$f_{sing}\mathrm{exp}\left(\frac{\mathrm{constant}}{\sqrt{t}}\right)$$
(2.41)
Near $`T_c`$, the free energy has an essential singularity.

The singularity of the free energy is very close to that of the Kosterlitz-Thouless transition. In fact, calculating exactly, we can show that the correlation length $`\xi `$ diverges at $`T_c`$ as $`\xi \mathrm{exp}\left(\mathrm{constant}/\sqrt{t}\right)`$, when $`T`$ approaches $`T_c`$ in the anti-ferroelectric phase.

### 2.6 XXZ spin chain and the transfer matrix

The logarithmic derivative of the transfer matrix of the six-vertex model gives the Hamiltonian of the XXZ spin chain.

$$\frac{d}{dv}\mathrm{log}\tau |_{v=\lambda }=\tau ^1\frac{d}{dv}\tau H_{XXZ}+constant$$
(2.42)
where $`H_{XXZ}`$ is given by

$$H_{XXZ}=J\underset{j=1}{\overset{L}{}}\left(\sigma _j^X\sigma _{j+1}^X+\sigma _j^Y\sigma _{j+1}^Y+\mathrm{\Delta }\sigma _j^Z\sigma _{j+1}^Z\right).$$
(2.43)
Intuitively, we may express it by $`\tau _{6V}(v)\mathrm{exp}(vH_{XXZ})`$.

The XXZ spin chain and the six-vertex transfer matrix have the same eigenvectors in common thanks to eq. (2.42). Taking logarithm of the Bethe ansatz equations, we have

$$Nk_j=2\pi I_j\underset{\mathrm{}=1}{\overset{M}{}}\mathrm{\Theta }(k_j,k_{\mathrm{}}),\mathrm{for}j=1,\mathrm{},M.$$
(2.44)
where $`M`$ is the number of down spins. (We assume $`2MN`$.) Here $`I_j`$ is an integer if $`M`$ is odd and half an integer if $`M`$ is even.

The ground state of the XXZ spin chain for $`\mathrm{\Delta }<1`$ was obtained by Yang and Yang . The ground state is specified by the integers $`I_j=j(M+1)/2`$ for $`j=1,\mathrm{},M`$. When $`N`$ is very large, the distribution of $`k_j`$’s becomes continuous. The number of $`k_j`$’s between $`k`$ and $`k+dk`$ can be approximated by $`N\rho (k)dk`$. Thus, we have the integral equation of $`\rho (k)`$

$$2\pi \rho (k)=1+_Q^Q\frac{\mathrm{\Theta }(k,k^{^{}})}{k}\rho (k^{^{}})𝑑k^{^{}}$$
(2.45)
where $`Q`$ is determined by the normalization condition

$$_Q^Q\rho (k)𝑑k=M/N.$$
(2.46)
The integral equation (2.45) can be solved by changing the variable $`k`$ to rapidity $`\alpha `$, and then by taking the Fourier transform for the half-filling case $`M/N=1/2`$. When $`M/N`$ is close to $`1/2`$, the integral equation can be solved by the Wiener-Hopf method .

### 2.7 Low-lying excited spectrum of the transfer matrix and conformal field theory

In this section, we assume that the low excited spectra of the transfer matrix of gapless models should be characterized by the conformal invariance, if the system size is large enough. The assumption is not rigorous, however, there are many studies which confirm it numerically for integrable models. We review the finite-size corrections for the XXZ spin chain and the six-vertex model.

#### 2.7.1 Finite-size corrections

Let us consider a conformally invariant field theory defined in the two-dimensional Euclidean space with coordinates $`r_1`$ and $`r_2`$. The energy momentum tensor $`T_{\mu \nu }`$ for $`\mu ,\nu =1,2`$ should be symmetric and traceless due to the conformal symmetry. Introducing the complex coordinates $`z=r_1+ir_2`$, $`\overline{z}=r_1ir_2`$, we define the chiral operator : $`T=(T_{11}T_{22}2iT_{12})/4`$, and the anti-chiral operator: $`\overline{T}=(T_{11}T_{22}+2iT_{12})/4`$. The operator $`T`$ (or $`\overline{T}`$) depends only on the variable $`z`$ (or $`\overline{z}`$).

The energy momentum tensor has the operator product expansion

$$T(z_1)T(z_2)=\frac{c/2}{(z_1z_2)^4}+\frac{2T(z_2)}{(z_1z_2)^2}+\frac{T(z_2)}{z_1z_2}+\mathrm{}$$
(2.47)
Here $`c`$ is called the central charge. We define the operators $`L_n`$ by the expansion: $`T(z)=_{n=\mathrm{}}^{\mathrm{}}L_nz^{n2}`$. The operator product expansion (2.47) corresponds to the Virasoro algebra

$$[L_n,L_m]=(mn)L_{m+n}+\frac{c}{12}(m^3m)\delta _{m+n,0}$$
(2.48)
Under a conformal transformation $`zw`$, the energy momentum tensor is transformed as

$$T(z)=\left(\frac{dw}{dz}\right)^2\stackrel{~}{T}(w)+\frac{c}{12}\{w,z\}$$
(2.49)
where the symbol $`\{w,z\}`$ denotes the Schwarzian derivative: $`(d^3w/dz^3)/(dw/dz)`$ \- $`(3/2)(d^2w/dz^2)^2/(dw/dz)^2`$.

Let us consider the conformal mapping from the $`z`$-plane to a cylinder of circumference $`L`$:

$$zw=\frac{L}{2\pi }\mathrm{log}z$$
(2.50)
Here $`w=\tau ix`$ with imaginary time $`\tau =it`$. The Hamiltonian $`\widehat{H}`$ on the cylinder is given by the space integral of the (1,1) component of the energy momentum tensor $`(T_{cyl})_{\mu \nu }`$

$`\widehat{H}`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}{\displaystyle _0^L}𝑑x\left(T_{cyl}(w)+\overline{T}_{cyl}(w)\right)`$ (2.51)
$`=`$ $`{\displaystyle \frac{2\pi }{L}}(L_0+\overline{L}_0){\displaystyle \frac{\pi c}{6L}}`$
Here we have used (2.49). For the momentum operator on the cylinder, we have

$`\widehat{P}`$ $`=`$ $`{\displaystyle \frac{1}{2\pi }}{\displaystyle _0^L}𝑑x\left(T_{cyl}(w)\overline{T}_{cyl}(w)\right)`$ (2.52)
$`=`$ $`{\displaystyle \frac{2\pi }{L}}(L_0\overline{L}_0)`$
Let us now discuss the application of the formulas (2.51) and (2.52) to the quantum spin chains. We assume that the low excited energies should be gapless and conformally invariant. In other words, we assume that the excitations near the ground state have a linear dispersion relation. Let $`v`$ denote the velocity of the linear dispersion. Then, for the ground-state energy $`E_0`$, we have

$$E_0=Le_{\mathrm{}}\frac{\pi vc}{6L}$$
(2.53)
and for the excited energy $`E_{ex}`$ and the momentum $`P_{ex}`$ we have

$`E_{ex}E_0`$ $`=`$ $`{\displaystyle \frac{2\pi v}{L}}(h+\overline{h}+N+\overline{N})`$
$`P_{ex}P_0`$ $`=`$ $`{\displaystyle \frac{2\pi }{L}}(h\overline{h}+N\overline{N})`$ (2.54)
where $`h`$ and $`\overline{h}`$ are the conformal weights related to the zero modes of the field, and the eigenvalues of $`N`$ and $`\overline{N}`$ are given by non-negative integers.

There is another viewpoint on the finite-size scaling. Let us consider the $`t`$ axis as the space axis for an infinitely long quantum spin chain, and $`x`$ axis as the imaginary time axis. Here we assume that $`L=v\beta =v/T`$. Thus, our system now becomes the quantum spin chain in the finite temperature $`T`$. Replacing $`E_0`$ with $`v\beta f`$, where $`f`$ denotes the finite-temperature free energy of the spin chain, we have $`f=\pi cT^2/6v`$. Thus, we may calculate the specific heat $`C`$ by the formula $`C=T^2f/T^2`$, and we have

$$C=\frac{\pi c}{3v}T$$
(2.55)
#### 2.7.2 The free Boson: CFT with $`c=1`$

Let us consider a free Bose field $`\phi (x,t)`$ defined on a cylinder of circumference $`L`$. The Lagrangian is given by

$$=\frac{1}{2}g𝑑x\left\{\left(_t\phi \right)^2\left(_x\phi \right)^2\right\}$$
(2.56)
We define the mode $`\phi _n`$ by the Fourier expansion $`\phi (x,t)=\phi _n(t)\mathrm{exp}(2\pi inx/L)`$ . From the canonical quantization, we have the conjugate momentum $`\pi _n=gL\dot{\phi }_n`$, and the commutation relation $`[\phi _n,\pi _m]=i\delta _{nm}`$ . With the operators $`a_n`$ and $`\overline{a}_n`$ for $`n0`$ satisfying

$$[a_n,a_m]=n\delta _{n+m}[a_n,\overline{a}_m]=0[\overline{a}_n,\overline{a}_m]=n\delta _{n+m},$$
(2.57)
the Fourier mode is expressed as $`\phi _n=i\left(a_n\overline{a}_n\right)/(n\sqrt{4\pi g})`$ , for $`n0`$. The Hamiltonian is given by

$$=\frac{1}{2gL}\pi _0^2+\frac{2\pi }{2L}\underset{n0}{}(a_na_n+\overline{a}_n\overline{a}_n)$$
(2.58)
Hereafter, we assume $`g=1/4\pi `$. The convention is consistent with the conformally invariant partition functions.

We now discuss the compactification of the Boson with radius $`R`$. Suppose that the field operator $`\phi `$ takes its value only on the circle of radius $`R`$. In other words, we may identify $`\phi `$ with $`\phi +2\pi R`$. Then, the eigenvalue of the momentum $`\pi _0`$ conjugate to $`\phi _0`$ is given by $`n/R`$ for an integer $`n`$. Here we recall that the wavenumber of a one-dimensional system of size $`L`$ is given by $`2\pi n/L`$ ($`n𝐙`$), and also that the range of $`\phi _0`$ is given by $`2\pi R`$, which corresponds to $`L`$. Furthermore, we may assign on the operator $`\phi `$ the boundary condition for an integer $`m`$

$$\phi (x+L,t)=\phi (x,t)+2\pi mR$$
(2.59)
Here, the integer $`m`$ is called the winding number. The mode expansion of $`\phi `$ is given by

$$\phi (x,t)=\phi _0+\frac{4\pi }{L}\pi _0t+\frac{2\pi Rm}{L}x+i\underset{n0}{}\frac{1}{n}\left(a_ne^{2\pi in(xt)/L}\overline{a}_ne^{2\pi in(x+t)/L}\right)$$
(2.60)
In terms of the coordinates $`z=\mathrm{exp}(2\pi (\tau ix)/L)`$ and $`\overline{z}=\mathrm{exp}(2\pi (\tau +ix)/L)`$, $`\phi (x,t)`$ is given by the sum of the holomorphic and anti-holomorphic parts: $`\phi (z,\overline{z})=\varphi (z)+\overline{\varphi }(\overline{z})`$ . Here, they are given by

$$\varphi (z)=\frac{\phi _0}{2}ia_0\mathrm{log}(z)+i\underset{k0}{}\frac{1}{k}a_kz^k\overline{\varphi }(\overline{z})=\frac{\phi _0}{2}i\overline{a}_0\mathrm{log}(\overline{z})+i\underset{k0}{}\frac{1}{k}\overline{a}_k\overline{z}^k$$
(2.61)
with $`a_0=n/R+mR/2`$ and $`\overline{a}_0=n/RmR/2`$ . Here we can show that the operator $`J(z)=i\varphi (z)/z`$ is the $`U(1)`$ current operator.

Making use of Noether’s theorem, we have

$$T(z)=\frac{1}{2}:(\frac{\varphi (z)}{z})^2:,\overline{T}(\overline{z})=\frac{1}{2}:(\overline{}\overline{\varphi }(\overline{z}))^2:$$
(2.62)
Here $`::`$ denotes a proper normal ordering. Then, through the Laurent expansion of powers of $`z`$, we have

$$L_0=\frac{1}{2}a_0^2+\underset{n=1}{\overset{\mathrm{}}{}}a_na_n$$
(2.63)
Thus, the conformal weights $`h_{nm}`$ and $`\overline{h}_{nm}`$ are given by

$$h_{n,m}=\frac{1}{2}\left(\frac{n}{R}+\frac{1}{2}mR\right)^2,\overline{h}_{n,m}=\frac{1}{2}\left(\frac{n}{R}\frac{1}{2}mR\right)^2$$
(2.64)
#### 2.7.3 The XXZ spin chain and CFT with $`c=1`$

We discuss the finite-size corrections to the XXZ spin chain. The finite-size corrections to the ground state energy is calculated in  based on the method with the Euler-MacLaurin formula . (For a review, see .) The result is

$`E_{ex}`$ $`=`$ $`Le_{\mathrm{}}{\displaystyle \frac{\pi v}{6L}}+{\displaystyle \frac{2\pi v}{L}}\left((\mathrm{\Delta }D)^2\xi ^2+{\displaystyle \frac{(\mathrm{\Delta }M)^2}{4\xi ^2}}+N+\overline{N}\right)`$ (2.65)
$`P_{ex}P_0`$ $`=`$ $`2k_F\mathrm{\Delta }D+{\displaystyle \frac{2\pi }{L}}(\mathrm{\Delta }D\mathrm{\Delta }M+N\overline{N})`$ (2.66)
Here the term $`e_{\mathrm{}}`$ denotes the ground-state energy per site. The Fermi velocity is obtained by the derivative of the dressed energy  with respect to the rapidity at the Fermi level.

The central charge $`c`$ is given by 1. $`\mathrm{\Delta }D`$ and $`\mathrm{\Delta }M`$ are integers. $`\mathrm{\Delta }M`$ denotes the change in the number of down spins, and $`\mathrm{\Delta }D`$ the number of particles jumping over the Fermi sea through the backscattering. We note that the difference of the conformal weights (2.64) is given by $`h_{nm}\overline{h}_{nm}=nm`$. Thus, $`\mathrm{\Delta }M`$ and $`\mathrm{\Delta }D`$ correspond to $`n`$ and $`m`$ of the $`c=1`$ CFT, respectively. $`N`$ and $`\overline{N}`$ are derived from the particle-hole excitations near the Fermi surface. The Fermi wavenumber $`k_F`$ is given by $`k_F=\pi M/L`$ where $`M`$ is the number of down spins. If the dispersion is linear, $`k_F`$ is consistent with the number of particles $`M`$.

The parameter $`\xi `$ is given by the dressed charge, which is defined by an integral equation. We note that the sum of the conformal weights (2.64) is given by $`h_{n,m}+\overline{h}_{n,m}`$ = $`n^2/R^2+m^2R^2/4`$. Thus, the dressed charge $`\xi `$ corresponds to the radius $`R`$ of the $`c=1`$ CFT, $`\xi =R/2`$. Under zero magnetic field, the dressed charge $`\xi `$ or the radius $`R`$ is given by

$$R=\left(\frac{2\pi }{\pi \mu }\right)^{1/2}(0\mu \pi ).$$
(2.67)
## 3 Various integrable models on two-dimensional lattices

### 3.1 Ising model and Potts model

#### 3.1.1 Ising model

Let us consider the Ising model defined on a square lattice. Each lattice site has spin variable which takes the two values $`\pm 1`$. We denote by $`\sigma _j`$ the spin variable of lattice site $`j`$. The Hamiltonian of the Ising model is given by

$$=J\underset{i,j}{}\sigma _i\sigma _j$$
(3.1)
where the symbol $`i,j`$ denote that sites $`i`$ and $`j`$ are nearest neighbors, and we take the sum over all the pairs of adjacent sites on the lattice.

There have been many papers written on the two-dimensional Ising model . However, it should be remarked that correlation functions are calculated exactly for the two-dimensional Ising model in the scaling limit .

We remark that exact solutions are discussed for the Ising model defined on various two-dimensional lattices such as the Kagome lattice (For a review, see Ref.  ).

#### 3.1.2 Self-dual Potts model

The Potts model generalizes the Ising model into a $`p`$-state model with $`p>2`$ . Let us consider the three-state Potts model defined on a square lattice . Each lattice site has spin variable which takes three values 1, 2, 3. We denote by $`\sigma _j`$ the spin variable of lattice site $`j`$. The Hamiltonian of the Potts model is given by

$$=J\underset{i,j}{}\delta (\sigma _i,\sigma _j)$$
(3.2)
Here the symbol $`\delta (a,b)`$ denote the Kronecker delta.

In general, the Potts model is not solvable. However, at its criticality, it is equivalent to some variant of the six-vertex model, and is solvable. Thus, the self-dual Potts model is integrable .

#### 3.1.3 Ashkin-Teller model

With each site $`i`$ we associate two spins: $`s_i`$ and $`t_i`$. They take values $`\pm 1`$. The Hamiltonian of the model is given by

$$=\underset{i,j}{}\left[K_2(s_is_j+t_it_j)+K_4s_is_jt_it_j\right]$$
(3.3)
It is known that the Ashkin-Teller model and the eight-vertex model are in the same universality class . The universality class is described by the $`c=1`$ CFT with the twisted boson .

### 3.2 Chiral Potts model

#### 3.2.1 General case

There is another version of the Potts model . We may assign the chiral symmetry, or $`𝐙/p𝐙`$ symmetry on the Potts model, which we call the chiral Potts model.

We now explain the most general chiral Potts model defined on a square lattice . Let $`a`$ and $`b`$ denote the spin variables defined on two nearest-neighboring sites. The interaction energy between the spins depends on the difference $`n=ab`$ (mod $`N`$) as

$$(n)=\underset{j=1}{\overset{N1}{}}E_j\omega ^{jn}\omega =e^{2\pi /N}$$
(3.4)
We note that $`(N+n)=(n)`$. The parameter $`E_j`$ can be written as

$$\frac{E_j}{k_BT}=K_j\omega ^{\mathrm{\Delta }_j},\mathrm{for}j=1,\mathrm{}[[\frac{N}{2}]]$$
(3.5)
where $`K_j`$ and $`\mathrm{\Delta }_j`$ constitute $`N1`$ independent variables, and the symbol $`[[]]`$ denotes the Gaussian symbol. For a real number $`x`$, $`[[x]]`$ denotes the biggest integer not larger than $`x`$. We also assume that $`E_{Nj}`$ is complex conjugate to $`E_j`$: $`E_{Nj}=E_j^{}`$. When $`N`$ is odd, we have

$$\frac{(n)}{k_BT}=\underset{j=1}{\overset{[[(N1)/2]]}{}}2K_j\mathrm{cos}\left(\frac{2\pi }{N}(jn+\mathrm{\Delta }_j)\right)+K_{N/2}(1)^n\frac{((1)^N+1)}{2}$$
(3.6)
#### 3.2.2 Integrable chiral Potts model

Let us discuss the integrable restriction of the most general chiral Potts model defined on the square lattice. It has horizontal and vertical couplings. Suppose that spin variables $`a`$ and $`b`$ are located on two neighboring sites connected by a horizontal line. When the line goes rightward from $`a`$ to $`b`$, the horizontal coupling has the energy $`_{pq}(ab)`$ and the Boltzmann weight $`W_{pq}(n)`$. Here $`p`$ and $`q`$ are ‘rapidity’ parameters. For spin variables $`c`$ and $`d`$ located on two neighboring sites connected by a vertical line, the vertical coupling has the energy $`\overline{}_{pq}(cd)`$ and the Boltzmann weight $`\overline{W}_{pq}(cd)`$ when the vertical line goes upward from $`c`$ to $`d`$.

The model is called solvable if the Boltzmann weights $`W`$’s and $`\overline{W}`$’s satisfy the star-triangle equation

$$\underset{d=1}{\overset{N}{}}\overline{W}_{qr}(bd)W_{pr}(ad)\overline{W}_{pq}(dc)=R_{pqr}W_{pq}(ab)\overline{W}_{pr})(bc)W_{qr}(ac)$$
(3.7)
The solution to the above equation is given by

$`W_{pq}(n)`$ $`=`$ $`W_{pq}(0){\displaystyle \underset{j=1}{\overset{n}{}}}\left({\displaystyle \frac{\mu _p}{\mu _q}}{\displaystyle \frac{y_qx_p\omega ^j}{y_px_q\omega ^j}}\right)`$
$`\overline{W}_{pq}(n)`$ $`=`$ $`\overline{W}_{pq}(0){\displaystyle \underset{j=1}{\overset{n}{}}}\left(\mu _p\mu _q{\displaystyle \frac{\omega x_px_q\omega ^j}{y_qy_p\omega ^j}}\right)`$ (3.8)
with a constant $`R_{pqr}`$ depending on the three rapidity variables. The constraint that the Boltzmann weight should have periodicity modulo $`N`$: $`W_{pq}(n+N)=W_{pq}(n)`$ gives for all rapidity pairs $`p`$ and $`q`$

$$\left(\frac{\mu _p}{\mu _q}\right)^N=\frac{y_p^Nx_q^N}{y_q^Nx_p^N},\left(\mu _p\mu _q\right)^N=\frac{y_q^Ny_p^N}{x_p^Nx_q^N}$$
(3.9)
We can define $`k`$ and $`k^{^{}}`$ such that

$`\mu _p`$ $`=`$ $`{\displaystyle \frac{k^{^{}}}{1kx_p^N}}={\displaystyle \frac{1ky_p^N}{k^{^{}}}}`$ (3.11)
$`x_p^N+y_p^N=k(1+x_p^Ny_p^N)`$
where $`k^2+(k^{^{}})^2=1`$. Thus, the rapidities are placed on a curve of genus $`g>1`$ of Fermat type.

We make some comments. Multiplying the two equations of (3.9) and noting that $`p`$ and $`q`$ are independent, we can show eq. (3.11) and then eq. (3.11). The star-triangle equation is proven illustratively in the appendix of Ref. . We note that the vertex-type formulation of the chiral Potts model is related to the tetrahedron equation . The integrable chiral Potts model generalizes the self-dual $`Z_N`$ model given by Fateev and Zamolodchikov . An elliptic extension of the self-dual $`Z_N`$ model is introduced in Ref. , and the Yang-Baxter equation is proven for the model in Ref. . Some nontrivial connections among higher-rank chiral Potts models, elliptic IRF models and Belavin’s $`Z_N`$ symmetric model are explicitly discussed in Ref. .

### 3.3 The eight-vertex model

Let us explain the eight-vertex model which was solved by Baxter in 1972 . The Boltzmann weights $`w(\alpha ,\beta |\gamma ,\delta ;u)`$ of are nonzero if the charge is conserved modulo 2: $`\alpha +\beta =\gamma +\delta `$ (mod) 2. We have the six configurations around a vertex shown in Fig. 2.2 and the two in Fig. 3.1 . We assume that there is no external field. The weights therefore become symmetric: $`ϵ_1=ϵ_2`$, $`ϵ_3=ϵ_4`$, $`ϵ_5=ϵ_6`$, and $`ϵ_7=ϵ_8`$.

$`w(1,1|1,1)`$ $`=`$ $`w(2,2|2,2)=w_1=a_{8V}`$
$`w(1,2|2,1)`$ $`=`$ $`w(2,1|1,2)=w_2=b_{8V}`$
$`w(1,2|1,2)`$ $`=`$ $`w(2,1|2,1)=w_3=c_{8V}`$
$`w(1,1|2,2)`$ $`=`$ $`w(2,2|1,1)=w_4=d_{8V}`$ (3.12)
We give a parametrization of the Boltzmann weights. We define the theta function

$$\theta (z;\tau )=2p^{1/4}\mathrm{sin}\pi z\underset{n=1}{\overset{\mathrm{}}{}}(1p^{2n})(1p^{2n}\mathrm{exp}(2\pi iz))(1p^{2n}\mathrm{exp}(2\pi iz)),$$
(3.13)
where the nome $`p`$ is related to the parameter $`\tau `$ by $`p=\mathrm{exp}(\pi i\tau )`$ with $`\mathrm{Im}\tau >0`$ . We also define the theta functions $`\theta _0(z)`$ and $`\theta _1(z)`$ satisfying $`\theta _\alpha (z+1)=(1)^\alpha \theta _\alpha (z)`$ and $`\theta _\alpha (z+\tau )=ie^{\pi i(z+\tau /2)}\theta _{1\alpha }(z)`$ for $`\alpha =0,1`$ , where we define $`\theta _1(z)`$ by $`\theta _1(z;\tau )=\theta (z;2\tau )`$. The Boltzmann weights $`a_{8V}(z),b_{8V}(z),c_{8V}(z)`$ and $`d_{8V}(z)`$ are expressed as

$`a_{8V}(z)`$ $`=`$ $`{\displaystyle \frac{\theta _0(z)\theta _0(2\eta )}{\theta _0(z2\eta )\theta _0(0)}},b_{8V}(z)={\displaystyle \frac{\theta _1(z)\theta _0(2\eta )}{\theta _1(z2\eta )\theta _0(0)}},`$
$`c_{8V}(z)`$ $`=`$ $`{\displaystyle \frac{\theta _0(z)\theta _1(2\eta )}{\theta _1(z2\eta )\theta _0(0)}},d_{8V}(z)={\displaystyle \frac{\theta _1(z)\theta _1(2\eta )}{\theta _0(z2\eta )\theta _0(0)}}.`$ (3.14)
### 3.4 IRF models

#### 3.4.1 Unrestricted 8V SOS model

We introduce unrestricted Solid-on-Solid models. They are also called Interaction-Round-a-Face models or IRF models, briefly .

To each site $`i`$ of a two-dimensional square lattice, a spin $`a_i`$ is associated. Let $`i,j,k`$, and $`\mathrm{}`$ be the lattice sites surrounding a face (or a square), where $`i,j,k,\mathrm{}`$ are placed counterclockwise from the southwest corner. We assume that an elementary configuration is given by that of the four spin variables around the face, and the probability of having $`a_i,a_j,a_k,a_{\mathrm{}}`$ is denoted by the Boltzmann weight $`w(a_i,a_j,a_k,a_{\mathrm{}};z)`$. Here the variable $`z`$ is called the spectral parameter.

For unrestricted 8V SOS model, $`a_i`$ can take any integer. When sites $`i`$ and $`j`$ are nearest neighboring, then the states $`a_i`$ and $`a_j`$ are said to be admissible if and only if $`|a_ia_j|=1`$ . The Yang-Baxter equations are given by

$`{\displaystyle \underset{g}{}}w(a,b,g,f;zw)w(f,g,d,e;z)w(g,b,c,d;w)`$ (3.15)
$`=`$ $`{\displaystyle \underset{g}{}}w(f,a,g,e;w)w(a,b,c,g;z)w(g,c,d,e;zw),`$
where the summation of the variable $`g`$ is taken over all the admissible states.

The Boltzmann weights are given by

$`w(d+1,d+2,d+1,d;z,w_0)`$ $`=`$ $`w(d,d1,d2,d1;z)={\displaystyle \frac{\theta (2\eta z)}{\theta (2\eta )}}`$
$`w(d1,d,d+1,d;z,w_0)`$ $`=`$ $`w(d+1,d,d1,d;z,w_0)`$
$`=`$ $`{\displaystyle \frac{\theta (z)}{\theta (2\eta )}}{\displaystyle \frac{\sqrt{\theta (2\eta (d+1)+w_0)\theta (2\eta (d1)+w_0)}}{\theta (2\eta d+w_0)}}`$
$`w(d+1,d,d+1,d;z,w_0)`$ $`=`$ $`{\displaystyle \frac{\theta (z+2\eta d+w_0)}{\theta (2\eta d+w_0)}}`$
$`w(d1,d,d1,d;z,w_0)`$ $`=`$ $`{\displaystyle \frac{\theta (z2\eta dw_0)}{\theta (2\eta d+w_0)}}`$ (3.16)
#### 3.4.2 RSOS models

Let us explain restricted Solid-on-Solid models (RSOS models) . Let $`s`$ denote the number of elements in $`S`$. Consider a $`s\times s`$ matrix $`C`$ satisfying the following conditions :

(i) $`C_{ab}=C_{ba}=0`$ or 1

(ii) $`C_{aa}=0`$

(iii) For each $`aS`$, there should exist $`bS`$ such that $`C_{ab}=1`$

For such choice of $`C`$, we impose a restriction that two states $`a`$ and $`b`$ can occupy the neighboring lattice sites if and only if $`C_{ab}=1`$. We call such a pair of the states $`(a,b)`$ admissible. For the case of unrestricted models, the infinite matrix $`C`$ satisfies the conditions (i), (ii), and (iii) with an infinite set $`S`$.

For an illustration, let us consider the restricted eight-vertex Solid-on-Solid model (the restricted 8V SOS model), which we also call the ABF model . For the $`N`$-state case, we have $`S=\{1,2,\mathrm{},N\}`$. The nonzero matrix elements of $`C`$ are given by $`C_{j,j+1}=C_{j+1,j}=1`$ for $`j=1,2,\mathrm{},N1`$; other matrix elements such as $`C_{1,N}`$ and $`C_{N,1}`$ are given by zero. Setting $`w_0=0`$, then we have the Boltzmann weights of the ABF model. The Boltzmann weights satisfy the Yang-Baxter relations (3.15) with the finite set: $`S=\{1,2,\mathrm{},N\}`$.

Let us explain CSOS models, another type of RSOS models . Here we assume that $`2N\eta =m_1`$, where integer $`m_1`$ has no common divisor with $`N`$. If we set $`w_00`$, then we have the Boltzmann weights of the cyclic SOS model (CSOS model). We can show that the Boltzmann weights satisfy the Yang-Baxter relations with the finite set: $`S=\{1,2,\mathrm{},N\}`$ and the cyclic admissible conditions: $`C_{1,N}=C_{N,1}=1`$.

We remark that the connection between the 6j symbols and the Boltzmann weights of IRF models is first discussed in Ref. .

#### 3.4.3 Fusion IRF models and ABCD IRF models

The 8VSOS model was generalized into the fusion IRF models.  IRF models associated with $`A_n^{(1)}`$ Lie algebra are constructed.  The IRF models associated with $`B^{(1)}`$ $`C^{(1)}`$ $`D^{(1)}`$ typed Lie algebra are also obtained .

In the IRF models, the one-point function, which is the magnetization per site, can be calculated by the corner transfer matrix method invented by Baxter .

#### 3.4.4 Gauge transformations

It is sometimes convenient to employ a gauge transformation

$$w(a,b,c,d;z)w(a,b,c,d;z)\frac{g_c}{g_a}.$$
(3.17)
The transformed Boltzmann weights also satisfy the Yang-Baxter relations (3.15). For instance, we may set $`g_a=\mathrm{exp}(\pi ia/2)\sqrt{\theta (2\eta a+w_0)}`$ ($`a𝐙`$).

## 4 Yang-Baxter equation and the algebraic Bethe Ansatz

### 4.1 Solutions to the Yang-Baxter equation

#### 4.1.1 Derivation of a solution for the six-vertex model

Let us solve the Yang-Baxter equation for the six-vertex model. We recall that it is given by

$`{\displaystyle \underset{\alpha ,\beta ,\gamma }{}}w(\alpha ,\gamma |a_1,a_2)w^{^{}}(\beta ,b_3|\gamma ,a_3)w^{^{\prime \prime }}(b_1,b_2|\alpha ,\beta )`$ (4.1)
$`=`$ $`{\displaystyle \underset{\alpha ,\beta ,\gamma }{}}w^{^{\prime \prime }}(\beta ,\alpha |a_2,a_3)w^{^{}}(b_1,\gamma |a_1,\beta )w(b_2,b_3|\gamma ,\alpha ).`$
The Yang-Baxter equation is illustrated in Fig. 4.1.

There are $`2^3\times 2^3=64`$ cases for the entries $`(a_1,a_2,a_3)`$ and $`(b_1,b_2,b_3)`$. Due to the ice rule, however, the Yang-Baxter equation is trivial unless $`a_1+a_2+a_3=b_1+b_2+b_3`$. We have thus only 20 entries: $`\left({}_{3}{}^{}C_{0}^{}\right)^2+\left({}_{3}{}^{}C_{1}^{}\right)^2+\left({}_{3}{}^{}C_{2}^{}\right)^2+\left({}_{3}{}^{}C_{3}^{}\right)^2=1+9+9+1=20`$.

Let us consider symmetries of the Boltzmann weights given in (2.3). If we exchange $`1`$ and $`2`$, the Boltzmann weights do not change. Thus, we have reduced 20 cases into 10 cases.

$`(a_1,a_2,a_3;b_1,b_2,b_3)=(1,1,1;1,1,1),`$
$`(1,1,2;1,1,2),(1,1,2;1,2,1),(1,1,2;2,1,1),`$
$`(1,2,1;1,1,2),(1,2,1;1,2,1),(1,2,1;2,1,1),`$
$`(2,1,1;1,1,2),(2,1,1;1,2,1),(2,1,1;2,1,1)`$ (4.2)
We now recall that the Boltzmann weights (2.3) have the symmetry

$$w(\alpha ,\beta |\gamma ,\delta )=w(\gamma ,\delta |\alpha ,\beta )=w(\beta ,\alpha |\delta ,\gamma )$$
(4.3)
Combining these symmetries we can show that the Yang-Baxter equation for the two entries (1) and (2) are equivalent: (1) $`(a_1,a_2,a_3;b_1,b_2,b_3)`$ ; (2) $`(b_3,b_2,b_1;a_3,a_2,a_1)`$. Precisely, the L.H.S. (or R.H.S) of the Yang-Baxter equation of the case $`(1)`$ corresponds to the R.H.S. (or L.H.S) of eq. (4.1). Thus, we have the three cases as follows.

$$(1,1,2;1,1,2),(1,1,2;1,2,1),(1,2,1;1,1,2)$$
(4.4)
For the three cases, the Yang-Baxter equations are given by

$`ac^{^{}}a^{^{\prime \prime }}`$ $`=`$ $`bc^{^{}}b^{^{\prime \prime }}+ca^{^{}}c^{^{\prime \prime }}`$
$`ab^{^{}}c^{^{\prime \prime }}`$ $`=`$ $`ba^{^{}}c^{^{\prime \prime }}+cc^{^{}}b^{^{\prime \prime }}`$
$`cb^{^{}}a^{^{\prime \prime }}`$ $`=`$ $`ca^{^{}}b^{^{\prime \prime }}+bc^{^{}}c^{^{\prime \prime }}`$ (4.5)
A nontrivial solution $`(a^{^{\prime \prime }},b^{^{\prime \prime }},c^{^{\prime \prime }})`$ exists only if the determinant vanishes

$$\left|\begin{array}{ccc}ac^{^{}}& bc^{^{}}& ca^{^{}}\\ 0& cc^{^{}}& ba^{^{}}ab^{^{}}\\ cb^{^{}}& ca^{^{}}& bc^{^{}}\end{array}\right|=abca^{^{}}b^{^{}}c^{^{}}\left(\frac{(a^{^{}})^2+(b^{^{}})^2(c^{^{}})^2}{a^{^{}}b^{^{}}}\frac{a^2+b^2c^2}{ab}\right)$$
(4.6)
We define the parameter $`\mathrm{\Delta }`$ by

$$\mathrm{\Delta }=\frac{a^2+b^2c^2}{2ab}$$
(4.7)
The condition that the determinant vanishes is given by

$$\mathrm{\Delta }=\mathrm{\Delta }^{^{}}$$
(4.8)
Thus, the transfer matrices $`\tau `$ and $`\tau ^{^{}}`$ commute, if the two sets of weights $`(a,b,c)`$ and $`(a^{^{}},b^{^{}},c^{^{}})`$ have the same $`\mathrm{\Delta }`$.

In terms of the spectral parameter $`u`$, we can parametrize the three Boltzmann weights $`a`$, $`b`$ and $`c`$. We express the weight as $`w(\alpha ,\beta |\gamma ,\delta ;u)`$. Let $`u`$ and $`v`$ be arbitrary. We denote $`w(\alpha ,\beta |\gamma ,\delta )`$, $`w^{^{}}(\alpha ,\beta |\gamma ,\delta )`$ , and $`w^{^{\prime \prime }}(\alpha ,\beta |\gamma ,\delta )`$ as $`w(\alpha ,\beta |\gamma ,\delta ;u)`$, $`w(\alpha ,\beta |\gamma ,\delta ;u+v)`$, and $`w(\alpha ,\beta |\gamma ,\delta ;v)`$, respectively. The Yang-Baxter equations are depicted in Fig. 4.2. As a solution, we may have $`(a,b,c)=(\rho \mathrm{sinh}(u+2\eta ),=\rho \mathrm{sinh}u,=\rho \mathrm{sinh}2\eta `$ . Here we set $`\mathrm{\Delta }=\mathrm{cosh}2\eta `$ . The transfer matrices $`\tau (u)`$ and $`\tau (v)`$ commute: $`\tau (u)\tau (v)=\tau (v)\tau (u)`$.

#### 4.1.2 Gauge transformations for vertex models

Let us suppose that $`w`$’s satisfy the Yang-Baxter equation. Then we can show that transformed weights $`\stackrel{~}{w}`$’s defined by

$$\stackrel{~}{w}(\alpha ,\beta |\gamma ,\delta ;u)=(ϵ)^{\alpha +\gamma }\mathrm{exp}(\kappa (\alpha +\gamma \beta \delta )u)w(\alpha ,\beta |\gamma ,\delta ;u)$$
(4.9)
also satisfy the Yang-Baxter equations . Here $`ϵ=\pm 1`$, and the number $`\kappa `$ is arbitrary.

The gauge transformation is important in the derivation of the Jones polynomial from the symmetric Boltzmann weights of the six-vertex model under zero field . It is also quite useful when we discuss the relation of the six-vertex model to the quantum group, as we shall see in §5 (see also the appendix of Ref. ).

### 4.2 Algebraic Bethe ansatz

#### 4.2.1 $`R`$ matrix and the $`L`$ operator

Let us diagonalize the transfer matrix by the method of the algebraic Bethe ansatz.

Let us introduce the notation of the matrix tensor product. We define the direct product $`AB`$ of matrices $`A`$ and $`B`$. Let $`A_k^j`$ denote the matrix element for the entry of column $`j`$ and row $`k`$ of the matrix $`A`$. Then, the matrix element of column $`(j_1,j_2)`$ and $`(k_1,k_2)`$ is defined by

$$\left(AB\right)_{k_1,k_2}^{j_1,j_2}=A_{k_1}^{j_1}B_{k_2}^{j_2}$$
(4.10)
We now define the R matrix of the XXZ spin chain. The element $`R_{cd}^{ab}`$ corresponds to the entry of column $`(a,b)`$ and row $`(c,d)`$.

$$R(z)=\left(\begin{array}{cccc}R(z)_{11}^{11}& R(z)_{12}^{11}& R(z)_{21}^{11}& R(z)_{22}^{11}\\ R(z)_{11}^{12}& R(z)_{12}^{12}& R(z)_{21}^{12}& R(z)_{22}^{12}\\ R(z)_{11}^{21}& R(z)_{12}^{21}& R(z)_{21}^{21}& R(z)_{22}^{21}\\ R(z)_{11}^{22}& R(z)_{12}^{22}& R(z)_{21}^{22}& R(z)_{22}^{22}\end{array}\right)=\left(\begin{array}{cccc}a(z)& 0& 0& 0\\ 0& c(z)& b(z)& 0\\ 0& b(z)& c(z)& 0\\ 0& 0& 0& a(z)\end{array}\right)$$
(4.11)
Here $`a(z)`$, $`b(z)`$ and $`c(z)`$ are given by

$$a(z)=\mathrm{sinh}(z+2\eta ),b(z)=\mathrm{sinh}z,c(z)=\mathrm{sinh}2\eta .$$
(4.12)
Here, the functions $`a(z)`$, $`b(z)`$ and $`c(z)`$ are equivalent to the Boltzmann weights $`a`$, $`b`$ and $`c`$ in §2

We now introduce $`L`$ operators. We write the matrix element the $`L`$ operator with entry $`(j,k)`$ as $`\left(L_n(z)\right)_{jk}`$ or $`L_n(z)_k^j`$ . The $`L`$ operator for the XXZ spin chain is given by

$$L_n(z)=\left(\begin{array}{cc}L_n(z)_1^1& L_n(z)_2^1\\ L_n(z)_1^2& L_n(z)_2^2\end{array}\right)=\left(\begin{array}{cc}\mathrm{sinh}\left(zI_n+\eta \sigma _n^z\right)& \mathrm{sinh}2\eta \sigma _n^{}\\ \mathrm{sinh}2\eta \sigma _n^+& \mathrm{sinh}\left(zI_n\eta \sigma _n^z\right)\end{array}\right)$$
(4.13)
Here $`I_n`$ and $`\sigma _n^a`$ ($`n=1,\mathrm{},L`$) are acting on the $`n`$ th vector space $`V_n`$. The $`L`$ operator is an operator-valued matrix which acts on the auxiliary vector space $`V_0`$. The symbols $`\sigma ^\pm `$ denote $`\sigma ^+=E_{12}`$ and $`\sigma ^{}=E_{21}`$, and $`\sigma ^x,\sigma ^y,\sigma ^z`$ are the Pauli matrices

In terms of the $`R`$ matrix and $`L`$ operators, the Yang-Baxter equation is expressed as

$$R(zt)\left(L_n(z)L_n(t)\right)=\left(L_n(t)L_n(z)\right)R(zt)$$
(4.14)
Here the tensor symbol in $`L_n(z)L_n(t)`$ denotes the tensor product of auxiliary spaces.

The Yang-Baxter equation (4.14) gives the relation between the two products of $`4\times 4`$ matrices. For an illustration, we consider the L.H.S. of (4.14).

$`\left[R(zt)L_n(z)L_n(t)\right]_{b_1,b_2}^{a_1,a_2}`$ $`=`$ $`{\displaystyle \underset{c_1,c_2}{}}R(zt)_{c_1,c_2}^{a_1,a_2}\left(L_n(z)L_n(t)\right)_{b_1,b_2}^{c_1,c_2}`$ (4.15)
$`=`$ $`{\displaystyle \underset{c_1,c_2}{}}R(zt)_{c_1,c_2}^{a_1,a_2}L_n(z)_{b_1}^{c_1}\times L_n(t)_{b_2}^{c_2}`$
Here the symbol $`\times `$ denotes the product of matrices acting on the $`n`$ th space $`V_n`$.

Expressing the operator products in the $`n`$th space $`V_n`$, the L.H.S. of (4.14) is written as follows

$`\left[R(zt)L_n(z)L_n(t)\right]_{b_1,b_2}^{a_1,a_2}|_{\alpha _n,\beta _n}`$ $`=`$ $`{\displaystyle \underset{c_1,c_2}{}}R(zt)_{c_1,c_2}^{a_1,a_2}\left(L_n(z)L_n(t)\right)_{b_1,b_2}^{c_1,c_2}|_{\alpha _n,\beta _n}`$ (4.16)
$`=`$ $`{\displaystyle \underset{c_1,c_2}{}}{\displaystyle \underset{\gamma _n}{}}R(zt)_{c_1,c_2}^{a_1,a_2}L_n(z)_{b_1}^{c_1}|_{\alpha _n,\gamma _n}L_n(t)_{b_2}^{c_2}|_{\gamma _n,\beta _n}`$
#### 4.2.2 Monodromy matrix and the construction of the eigenvector

We define the monodromy matrix by the product of the $`L`$ operators (see also Fig. A.1)

$$T(z)=L_N(z)\mathrm{}L_2(z)L_1(z).$$
(4.17)
The transfer matrix $`\tau _{6V}(z)`$ of the six vertex model is given by the trace of $`T(z)`$

$$\tau _{6V}(z)=\mathrm{tr}T(z)=A(z)+D(z),\mathrm{where}T(z)=\left(\begin{array}{cc}A(z)& B(z)\\ C(z)& D(z)\end{array}\right)$$
(4.18)
The Yang-Baxter equation (4.14) leads to the commutation relation: $`R(zt)(T(z)T(t))=(T(t)T(z))R(zt)`$, from which we have many relations among the operators $`A`$, $`B`$, $`C`$, and $`D`$. For instance, we have $`B(z)B(t)=B(t)B(z)`$. Furthermore, we have

$`A(z)B(t)`$ $`=`$ $`{\displaystyle \frac{a(tz)}{b(tz)}}B(t)A(z){\displaystyle \frac{c(tz)}{b(tz)}}B(z)A(t)`$ (4.19)
$`D(z)B(t)`$ $`=`$ $`{\displaystyle \frac{a(zt)}{b(zt)}}B(t)D(z){\displaystyle \frac{c(zt)}{b(zt)}}B(z)D(t)`$ (4.20)
We define the ‘vacuum’ by

$$|0=\stackrel{N}{\stackrel{}{|_1|_2\mathrm{}|_N}}$$
(4.21)
Multiplying $`A(z)`$ and $`D(z)`$ on the vacuum, we have

$$A(z)|0=a(z\eta )^N|0,D(z)|0=b(z\eta )^N|0.$$
(4.22)
Let us consider the vector generated by the product of $`B`$ operators.

$$|M=B(t_1)\mathrm{}B(t_M)|0$$
(4.23)
Then, through the commutation relations such as (4.19) and (4.20) we can show  that the vector $`|M`$ gives an eigenvector of the transfer matrix $`\tau _{6V}(z)`$ if rapidities $`t_1,t_2,\mathrm{}`$, $`t_M`$ satisfy the set of equations

$$\left(\frac{a(t_j\eta )}{b(t_j\eta )}\right)^N=\frac{c(t_kt_j)}{b(t_kt_j)}\frac{b(t_jt_k)}{c(t_jt_k)}\underset{k=1;kj}{\overset{M}{}}\left(\frac{a(t_jt_k)}{b(t_jt_k)}\frac{b(t_kt_j)}{a(t_kt_j)}\right),\text{ for}j=1,\mathrm{},M.$$
(4.24)
They are the Bethe ansatz equations (2.12) with different parametrization. For a set of solutions, $`t_1,t_2,\mathrm{},t_M`$ to eqs. (4.24), the eigenvalue of the transfer matrix $`\tau _{6V}(z)`$ is given by

$`\mathrm{\Lambda }(z;t_1,t_2,\mathrm{},t_M)`$ $`=`$ $`a(z\eta )^N{\displaystyle \underset{j=1}{\overset{M}{}}}{\displaystyle \frac{a(t_jz)}{b(t_jz)}}+b(z\eta )^N{\displaystyle \underset{j=1}{\overset{M}{}}}{\displaystyle \frac{a(zt_j)}{b(zt_j)}}`$
$`=`$ $`\mathrm{sinh}^N(z+\eta ){\displaystyle \underset{j=1}{\overset{M}{}}}{\displaystyle \frac{\mathrm{sinh}(t_jz+2\eta )}{\mathrm{sinh}(t_jz)}}+\mathrm{sinh}^N(z\eta ){\displaystyle \underset{j=1}{\overset{M}{}}}{\displaystyle \frac{\mathrm{sinh}(t_jz2\eta )}{\mathrm{sinh}(t_jz)}}.`$
#### 4.2.3 Connection to the coordinate Bethe ansatz result

Let us compare the result of the coordinate Bethe ansatz in §2. We consider the disordered phase: $`1<\mathrm{\Delta }<1`$. First, we change the variables $`w`$ and $`\alpha `$ defined in §2 into $`u`$ and $`\zeta `$ by $`w=2u\mu `$ and $`\alpha =2\zeta i\mu `$, respectively. Thus, from the expressions (2.22) we have $`(a,b,c)`$ = $`(\mathrm{sin}(\mu u),\mathrm{sin}u,\mathrm{sin}\mu )`$ = $`i(\mathrm{sinh}(i\mu +iu),\mathrm{sinh}(iu),\mathrm{sinh}(i\mu ))`$ and

$$L^{Baxter}(z_j)=\frac{ab+(c^2b^2)z_j}{a(abz_j)}=\frac{\mathrm{sinh}((\alpha _jiw2i\mu )/2)}{\mathrm{sinh}((\alpha _jiw)/2)}=\frac{\mathrm{sinh}((\zeta _jiu)+i\mu )}{\mathrm{sinh}(\zeta _jiu)}$$
(4.26)
Here we recall that the symbol $`L^{Baxter}(z)`$ has been defined in eq. (2.15) in order to denote the eigenvalue of the transfer matrix of the six-vertex model. The expression (4.26) can be derived from the formula (LABEL:evABA) of the algebraic Bethe ansatz as follows. First, we take the gauge transformation: $`b(z)b(z)`$, which corresponds to the case $`ϵ=1`$ and $`\kappa =0`$ in eq. (4.9). Then, we replace the variables $`z`$, $`t_j`$’s, and $`2\eta `$ in (LABEL:evABA) by $`z\eta iu`$, $`t_j\eta \zeta _j`$, and $`2\eta i\mu `$, respectively. We have the following:

$$\frac{a(t_jz)}{b(t_jz)}\frac{a(t_jz)}{b(t_jz)}=\frac{\mathrm{sinh}(t_jz+2\eta )}{\mathrm{sinh}(t_jz)}\frac{\mathrm{sinh}((\zeta _jiu)+i\mu )}{\mathrm{sinh}(\zeta _jiu)}$$
(4.27)
and $`\mathrm{sinh}^N(z+\eta )\mathrm{sinh}^N(iu+i\mu )`$ . Thus, the formula (LABEL:evABA) reproduces the expression (2.14) for the eigenvalues of the transfer matrix except for the normalization factor $`i^N`$.

## 5 Mathematical structures of integrable lattice models

### 5.1 Braid group

#### 5.1.1 The Yang-Baxter equation in operator formalism

The Yang-Baxter equation in §2 gives a sufficient condition for the existence of commuting transfer matrices. However, there are other viewpoints on the Yang-Baxter equation.

Let $`E_k^j`$ denote the matrix given by

$$\left(E_k^j\right)_b^a=\delta _{a,j}\delta _{b,k}a,b=1,2.$$
(5.1)
We define operators $`X_j(u)`$’s by

$$X_j(z)=\underset{a,b,c,d}{}w(a,b|c,d;z)\stackrel{(j1)}{\stackrel{}{I\mathrm{}I}}E_a^cE_b^d\stackrel{(Nj1)}{\stackrel{}{I\mathrm{}I}}\text{ for}j=1,\mathrm{},N1.$$
(5.2)
Here, the symbol $`w(a,b|c,d;z)`$ corresponds to $`R_{a,b}^{c,d}(z)`$, and the essential part of $`X_j(z)`$ is given by

$$X(z)=\underset{a,b,c,d}{}w(a,b|c,d;z)E_a^cE_b^d$$
(5.3)
which is equivalent to $`R(z)`$. In terms of the $`X_j(z)`$’s, the Yang-Baxter equation is expressed by the following

$$X_j(z)X_{j+1}(z+t)X_j(t)=X_{j+1}(t)X_j(z+t)X_{j+1}(z),\mathrm{for}j=1,\mathrm{},N1.$$
(5.4)
#### 5.1.2 The braid group

The braid group $`B_N`$ for $`N`$ strings is an infinite group which is generated by the generators $`b_1,\mathrm{},b_{N1}`$ satisfying the defining relations

$`b_jb_{j+1}b_j`$ $`=`$ $`b_{j+1}b_jb_{j+1}`$
$`b_ib_j`$ $`=`$ $`b_jb_i\mathrm{for}|ij|>1`$ (5.5)
Let us assume that the limit: $`lim_z\mathrm{}X_j(z)`$ exists. Then, eq. (5.4) becomes

$$X_j(\mathrm{})X_{j+1}(\mathrm{})X_j(\mathrm{})=X_{j+1}(\mathrm{})X_j(\mathrm{})X_{j+1}(\mathrm{}),\mathrm{for}j=1,\mathrm{},N1.$$
(5.6)
This is nothing but the defining relations of the braid group. Thus, the Boltzmann weights of solvable models expressed in terms of the spectral parameter lead to representations of the braid group.

We now show that from a given exactly solvable model, one can derive two different representations of the braid group . Here, the gauge transformation (4.9) plyas a central role. This technical point is quite fundamental when we make connections of exactly solvable models to quantum groups and the Temperley-Lieb algebra (for instance, see the appendix of Ref. ).

We first consider the Boltzmann weights of the zero-field six-vertex model given by eqs. (4.12). Taking the infinite limit to them, we have the following

$$\underset{z\mathrm{}}{lim}X(z)/\mathrm{sinh}(z+2\eta )=\left(\begin{array}{cccc}1& 0& 0& 0\\ 0& 0& \mathrm{exp}(2\eta )& 0\\ 0& \mathrm{exp}(2\eta )& 0& 0\\ 0& 0& 0& 1\end{array}\right)$$
(5.7)
The representation of the braid group (5.7) leads to a link polynomial equivalent to the linking number.

Let us now apply the gauge transformation (4.9) with $`ϵ=1`$ and $`\kappa =1/2`$ to the weights given by eqs. (4.12)  . Then, from the transformed Boltzmann weights we have

$$\underset{z\mathrm{}}{lim}\stackrel{~}{X}(z)/\mathrm{sinh}(z+2\eta )=\left(\begin{array}{cccc}1& 0& 0& 0\\ 0& 0& \mathrm{exp}(2\eta )& 0\\ 0& \mathrm{exp}(2\eta )& 1\mathrm{exp}(4\eta )& 0\\ 0& 0& 0& 1\end{array}\right)$$
(5.8)
This matrix representation of the braid group leads to the Jones polynomial with $`q=\mathrm{exp}(2\eta )`$ .

We remark that based on the representations of the braid group which are derived from the Boltzmann weights of exactly solvable models, we can construct various invariants of knots and links (see, for reviews, Refs. ).

### 5.2 Quantum groups (Hopf algebras)

From the quantum groups (the Hopf algebras) we can systematically construct representations of the braid group such as derived from the six-vertex model. Almost all the solutions of the Yang-Baxter equations can be constructed in some framework of quantum groups. Furthermore, the connection of solvable models to the quantum groups is useful for investigating nontrivial properties of integrable models.

Let us introduce the quantum group $`U_q(sl_2)`$ , which is a $`q`$-analog of the universal enveloping algebra of $`sl_2`$. Generators $`X^\pm ,H`$ satisfy

$$KX^\pm K^1=q^{\pm 2}X^\pm ,[X^+,X^{}]=\frac{KK^1}{qq^1}$$
(5.9)
We may express $`K`$ as $`K=q^H`$. Taking the limit of $`q`$ to unity, the relations are reduced into the commutation relations of $`sl_2`$.

The tensor product is defined by the following

$`\mathrm{\Delta }(K)`$ $`=`$ $`KK`$
$`\mathrm{\Delta }(X^+)`$ $`=`$ $`X^+I+KX^+`$
$`\mathrm{\Delta }(X^{})`$ $`=`$ $`X^{}K^1+IX^{}`$ (5.10)
The operation $`\mathrm{\Delta }()`$ is called the comultiplication. In the quantum group, the comultiplication does not commute with the exchange operator $`\sigma `$ defined by $`\sigma ab=ba`$. However, there is an operator $`R`$ which satisfy the following

$$R\mathrm{\Delta }(x)=\sigma \mathrm{\Delta }(x)R,\mathrm{for}xU_q(sl_2)$$
(5.11)
Thus, the tensor product $`V_1V_2`$ can be related to $`V_2V_1`$ through the R matrix of the quantum group.

In the $`U_q(sl_2)`$, the R matrix can be constructed in the operator formalism. If operators $`X^\pm `$ and $`H`$ satisfy the defining relations of the $`U_q(sl_2)`$, then the operator $``$ defined by the following satisfy the intertwining relation (5.11)

$$=q^{HH/2}\mathrm{exp}_q\left((qq^1)K^1X^+X^{}K\right)$$
(5.12)
where $`\mathrm{exp}_qx`$ denotes the infinite series

$$\mathrm{exp}_qx=\underset{n=0}{\overset{\mathrm{}}{}}q^{n(n1)/2}\frac{x^n}{[n]!}$$
(5.13)
The operator $``$ is called the universal $`R`$ matrix.

The representation (5.8) of the braid group corresponds to the representation of the universal $`R`$ matrix on the tensor product of two fundamental representations.

We note that the universal $`R`$ matrix can be constructed canonically through Drinfeld’s quantum double construction (for instance, see ). This is similar to the Sugawara construction which derives the energy momentum tensor from the current operator.

## Appendix A Commuting transfer matrices and the Yang-Baxter equations

We show that if given two sets of the Boltzmann weights of the six vertex model satisfy the Yang-Baxter relation, then their transfer matrices commute. We consider three sets of the Boltzmann weights: $`(w_1,w_2,w_3)=(a,b,c)`$, $`(a^{^{}},b^{^{}},c^{^{}})`$, and $`(a^{^{\prime \prime }},b^{^{\prime \prime }},c^{^{\prime \prime }})`$. Let us denote by $`\tau ^{^{}}`$ and $`\tau ^{^{\prime \prime }}`$ the transfer matrices constructed from the sets of the Boltzmann weights $`(a^{^{}},b^{^{}},c^{^{}})`$ and $`(a^{^{\prime \prime }},b^{^{\prime \prime }},c^{^{\prime \prime }})`$, respectively. Then, we can show that if the three sets of the Boltzmann weights satisfy the Yang-Baxter equations given in §2.3, then the transfer matrices $`\tau ^{^{}}`$ and $`\tau ^{^{\prime \prime }}`$ commute.

Let us now explicitly discuss the commutation relation. We first introduce the monodromy matrix. It is an $`N`$ ranked tensor, whose $`(\alpha ,\beta )`$ elements are defined as follows

$$\left(T_{\alpha ,\beta }\right)_{b_1,\mathrm{},b_N}^{a_1,\mathrm{},a_N}=\underset{c_2,\mathrm{},c_N}{}w(\beta ,b_1|a_1,c_2)w(c_2,b_2|a_2,c_3)\mathrm{}w(c_N,b_N|a_N,\alpha )$$
(A.1)
The transfer matrix is given by the trace of the monodromy matrix

$$\tau =\mathrm{tr}\left(T\right)=\underset{\alpha =1,2}{}T_{\alpha ,\alpha }$$
(A.2)
In terms of the matrix elements, we have

$$\left(\tau \right)_{b_1,\mathrm{},b_N}^{a_1,\mathrm{},a_N}=\underset{\alpha =1,2}{}\left(T_{\alpha ,\alpha }\right)_{b_1,\mathrm{},b_N}^{a_1,\mathrm{},a_N}$$
(A.3)
Let us denote by $`T^{^{}}`$ and $`T^{^{\prime \prime }}`$ the monodromy matrices for the sets of the Boltzmann weights $`(a^{^{}},b^{^{}},c^{^{}})`$ and $`(a^{^{\prime \prime }},b^{^{\prime \prime }},c^{^{\prime \prime }})`$, respectively. We consider the product of the matrix elements of the two monodromy matrices: $`T_{\alpha ,\beta }^{^{}}T_{\gamma ,\delta }^{^{\prime \prime }}`$. The entry of $`(a_1,\mathrm{},a_N)`$ and $`(b_1,\mathrm{},b_N)`$ of the product $`T_{\alpha ,\beta }^{^{}}T_{\gamma ,\delta }^{^{\prime \prime }}`$ is given by

$`\left(T_{\gamma _1,\gamma _{N+1}}^{^{}}T_{\delta _1,\delta _{N+1}}^{^{\prime \prime }}\right)_{b_1,\mathrm{},b_N}^{a_1,\mathrm{},a_N}={\displaystyle \underset{e_1,\mathrm{},e_N}{}}\left(T_{\gamma _1,\gamma _{N+1}}^{^{}}\right)_{e_1,\mathrm{},e_N}^{a_1,\mathrm{},a_N}\left(T_{\delta _1,\delta _{N+1}}^{^{\prime \prime }}\right)_{b_1,\mathrm{},b_N}^{e_1,\mathrm{},e_N}`$ (A.4)
$`=`$ $`{\displaystyle \underset{e_1,\mathrm{},e_N}{}}{\displaystyle \underset{c_2,\mathrm{},c_N}{}}w^{^{}}(\gamma _1,e_1|a_1,c_2)w^{^{}}(c_2,e_2|a_2,c_3)\mathrm{}w^{^{}}(c_N,e_N|a_N,\gamma _{N+1})`$
$`\times {\displaystyle \underset{d_2,\mathrm{},d_N}{}}w^{^{\prime \prime }}(\delta _1,b_1|e_1,d_2)w^{^{\prime \prime }}(d_2,b_2|e_2,d_3)\mathrm{}w^{^{\prime \prime }}(d_N,b_N|e_N,\delta _{N+1})`$
$`=`$ $`{\displaystyle \underset{c_2,\mathrm{},c_N}{}}{\displaystyle \underset{d_2,\mathrm{},d_N}{}}{\displaystyle \underset{e_1,\mathrm{},e_N}{}}\left(w^{^{}}(\gamma _1,e_1|a_1,c_2)w^{^{\prime \prime }}(\delta _1,b_1|e_1,d_2)\right)`$
$`\left(w^{^{}}(c_2,e_2|a_2,c_3)w^{^{\prime \prime }}(d_2,b_2|e_2,d_3)\right)\mathrm{}\left(w^{^{}}(c_N,e_N|a_N,\gamma _{N+1})w^{^{\prime \prime }}(d_N,b_N|e_N,\delta _{N+1})\right)`$
$`=`$ $`{\displaystyle \underset{c_2,\mathrm{},c_N}{}}{\displaystyle \underset{d_2,\mathrm{},d_N}{}}S(a_1,b_1)_{\delta _1,d_2}^{\gamma _1,c_2}S(a_2,b_2)_{d_2,d_3}^{c_2,c_3}\mathrm{}S(a_N,b_N)_{d_N,\delta _{N+1}}^{c_N,\gamma _{N+1}}`$
Here $`S(a_j,b_j)_{d_j,d_{j+1}}^{c_j,c_{j+1}}`$ has been defined by

$$S(a_j,b_j)_{d_j,d_{j+1}}^{c_j,c_{j+1}}=\underset{e_j}{}w^{^{}}(c_j,e_j|a_j,c_{j+1})w^{^{\prime \prime }}(d_j,b_j|e_j,d_{j+1})$$
(A.5)
We define the matrix element $`M_{d_0,d_1}^{c_0,c_1}`$ as follows

$$M_{d_0,d_1}^{c_0,c_1}=w(d_0,d_1|c_0,c_1)$$
(A.6)
Here we assume that $`M_{d_0,d_1}^{c_0,c_1}`$ denotes the matrix element for column $`(c_0,d_0)`$ and row $`(c_1,d_1)`$ of the matrix $`M`$. Multiplying the matrix $`M`$ to the product $`T_{\alpha ,\beta }^{^{}}T_{\gamma ,\delta }^{^{\prime \prime }}`$ and applying the Yang-Baxter relation $`N`$ times, we can derive the following

$$\underset{c_1,d_1}{}M_{\delta _0,d_1}^{\gamma _0,c_1}T_{c_1,\gamma _{N+1}}^{^{}}T_{d_1,\delta _{N+1}}^{^{\prime \prime }}=\underset{c_N,d_N}{}T_{c_1,c_{N+1}}^{^{\prime \prime }}T_{d_1,d_{N+1}}^{^{}}M_{d_N,\delta _{N+1}}^{c_N,\gamma _{N+1}}$$
(A.7)
Let us briefly discuss the derivation of relation (A.7). It is depicted in Fig. A.2. In the first equality of Fig. 2.6, we have applied the Yang-Baxter relation formulated as follows

$$\underset{c_1,d_1}{}M_{\delta _0,d_1}^{\gamma _0,c_1}S(a_1,b_1)_{d_1,d_2}^{c_1,c_2}=\underset{c_1,d_1}{}S^{^{}}(a_1,b_1)_{\delta _0,d_1}^{\gamma _0,c_1}M_{d_1,d_2}^{c_1,c_2}$$
(A.8)
Here the symbol $`S^{^{}}(a_j,b_j)_{d_j,d_{j+1}}^{c_j,c_{j+1}}`$ has been defined by

$$S^{^{}}(a_j,b_j)_{d_j,d_{j+1}}^{c_j,c_{j+1}}=\underset{e_j}{}w^{^{\prime \prime }}(c_j,e_j|a_j,c_{j+1})w^{^{}}(d_j,b_j|e_j,d_{j+1})$$
(A.9)
We also note that the LHS of (A.7) corresponds to the sum

$$\underset{c_1,c_2,\mathrm{},c_N}{}\underset{d_1,d_2,\mathrm{},d_N}{}M_{\delta _0,d_1}^{\gamma _0,c_1}S(a_1,b_1)_{d_1,d_2}^{c_1,c_2}S(a_2,b_2)_{d_2,d_3}^{c_2,c_3}\mathrm{}S(a_N,b_N)_{d_N,\delta _{N+1}}^{c_N,\gamma _{N+1}}$$
(A.10)
Let us consider the inverse of the matrix $`M`$

$$\left(M^1M\right)_{d_1,d_2}^{c_1,c_2}=\left(MM^1\right)_{d_1,d_2}^{c_1,c_2}=\delta _{c_1,c_2}\delta _{d_1,d_2}$$
(A.11)
Multiplying the inverse $`M^1`$ to both hand sides of (A.7) we have

$$MT^{^{}}T^{^{\prime \prime }}M^1=T^{^{\prime \prime }}T^{^{}}$$
(A.12)
Noting $`tr(MT^{^{}}T^{^{\prime \prime }}M^1)=tr(T^{^{}}T^{^{\prime \prime }})`$, we obtain the commutation relation of the transfer matrices

$$\tau ^{^{}}\tau ^{^{\prime \prime }}=\tau ^{^{\prime \prime }}\tau ^{^{}}$$
(A.13)