arxiv: v3 [math-ph] 17 Aug 2016

Transcription

1 arxiv: v3 [math-ph] 17 Aug 2016 Exact soution of the trigonometric SU(3) spin chain with generic off-diagona boundary refections Guang-Liang Li a, Junpeng Cao b,c,d, Kun Hao e,f, Fakai Wen e,f, Wen-Li Yang e,f,g1 and Kangjie Shi e,f a Department of Appied Physics, Xian Jiaotong University, Xian , China b Institute of Physics, Chinese Academy of Sciences, Beijing , China c Schoo of Physica Sciences, University of Chinese Academy of Sciences, Beijing, China d Coaborative Innovation Center of Quantum Matter, Beijing, China e Institute of Modern Physics, Northwest University, Xian , China f Shaanxi Key Laboratory for Theoretica Physics Frontiers, Northwest University, Xian , China g Beijing Center for Mathematics and Information Interdiscipinary Sciences, Beijing, , China Abstract The nested off-diagona Bethe ansatz is generaized to study the quantum spin chain associated with the SU q (3) R-matrix and generic integrabe non-diagona boundary conditions. By using the fusion technique, certain cosed operator identities among the fused transfer matrices at the inhomogeneous points are derived. The corresponding asymptotic behaviors of the transfer matrices and their vaues at some specia points are given in detai. Based on the functiona anaysis, a nested inhomogeneous T Q reations and Bethe ansatz equations of the system are obtained. These resuts can be naturay generaized to cases reated to the SU q (n) agebra. PACS: Pq, Ik, Pm Keywords: Spin chain; refection equation; Bethe ansatz; T Q reation 1 Corresponding author: wyang@nwu.edu.cn

2 1 Introduction Exact soution is a very important issue in studies of statistica mechanics, condensed matter physics, quantum fied theory and mathematica physics [1, 2] since those resuts can provide important benchmarks for understanding physica effects in a variety of systems. The coordinate Bethe ansatz and the agebraic Bethe ansatz are two powerfu methods to obtain the exact soution of the integrabe systems [3, 4, 5, 6, 7]. With these methods, many interesting exacty sovabe modes, such as the one-dimensiona Hubbard mode, supersymmetric t J mode, Heisenberg spin chain and the δ-potentia quantum gas mode, were exacty soved. For integrabe systems with U(1) symmetry, it is easy to find a reference state and these conventiona Bethe ansatz can be appied to. Indeed, most of the previous studies focus on periodic or diagona open boundary conditions without breaking the U(1) symmetry. However, there exists another kind of integrabe systems which does not have the U(1) symmetry, such as the integrabe systems with generic off-diagona boundary refections. Because the reference state of this kind of integrabe system is absent, the conventiona Bethe ansatz methods are faied. On the other hand, many interesting phenomena arise in this kind of systems, such as the topoogica eementary excitations in the spin-1/2 torus [8], spira phase in the Heisenberg mode with unparae boundary magnetic fied [9] and stochastic process in non-equiibrium statistica mechanics [10, 11, 12]. Motivated by these important appications, many interesting methods such as the q-onsager agebra [13, 14, 15], the modified agebraic Bethe ansatz [16, 17, 18, 19] and the Skyanin s separation of variabes (SoV) [20, 21, 22, 23, 24] were aso appied to some integrabe modes without U(1) symmetry. Other interesting progress can be found in [25, 26, 27, 28, 29]. Recenty, a new approach, i.e., the off-diagona Bethe ansatz (ODBA) [8] was proposed to obtain exact soutions of generic integrabe modes either with or without U(1) symmetry. Severa ong-standing probems were then soved [30, 31, 32, 33, 34, 35, 36] via this method. For comprehensive introduction to this method we refer the readers to [37]. In order to study the high rank integrabe modes, the nested version of ODBA has been proposed for the isotropic (or rationa) modes [33]. In this paper, we study the anisotropic rank-2 spin mode with generic integrabe boundary conditions. Here the R-matrix is the trigonometric one associated with the SU q (3) agebra and the boundary refection matrices are the most generic refection matrices which have non-vanishing off-diagona eements. Because the off- 2

3 diagona eements of the refection matrices break the U(1) symmetry, the exact soution of the system has been missing even its integrabiity was known for many years ago. By using the fusion technique and nested ODBA, we successfuy obtain the cosed operator identities, the vaues at the specia points and the asymptotic behaviors. Based on them, we construct the nested inhomogeneous T Q reation and obtain the eigenvaue of the transfer matrix thus the energy spectrum of the system. These resuts can be generaized to mutipe components spin chains reated to more higher rank agebra cases. The paper is organized as foows. Section 2 is the genera description of the mode. The SU q (3) R-matrixand corresponding generic integra non-diagona boundary refection matrices are introduced. In Section 3, by using the fusion technique, we derive the cosed operator identities for the fused transfer matrices and the quantum determinant. The asymptotic behaviors of the fused transfer matrix and their vaues at specia points are aso obtained. In section 4, we ist some necessary functiona reations which are used to determine the eigenvaues. Section 5 is devoted to the construction of the nested inhomogeneous T Q reation and the Bethe ansatz equations. In section 6, we summarize our resuts and give some discussions. Some resuts reated to the other types of the genera off-diagona boundary refections are given in Appendix. 2 The mode Throughout, V denotes a three-dimensiona inear space and et { i,i = 1,2,3} be an orthonormabasis ofit. Wesha adoptthestandardnotations. ForanymatrixA End(V), A j is an embedding operator in the tensor space V V, which acts as A on the j-th space and as identity on the other factor spaces. For B End(V V), B ij is an embedding operator of B in the tensor space, which acts as identity on the factor spaces except for the i-th and j-th ones. The R-matrix R(u) End(V V) used in this paper is the trigonometric one associated withthesu q (3)agebra, whichwasfirstproposedbyperkandshutz[38]andfurtherstudied 3

4 in [39, 40, 41, 42, 43], a(u) R 12 (u) = b(u) d(u) b(u) d(u) c(u) b(u) a(u) b(u) d(u) c(u) b(u) c(u) b(u) a(u), (2.1) where the matrix eements are a(u) = sinh(u+η), b(u) = sinh(u), (2.2) c(u) = e u sinhη, d(u) = e u sinhη. (2.3) The R-matrix satisfies the quantum Yang-Baxter equation (QYBE) R 12 (u 1 u 2 )R 13 (u 1 u 3 )R 23 (u 2 u 3 ) = R 23 (u 2 u 3 )R 13 (u 1 u 3 )R 12 (u 1 u 2 ), (2.4) and possesses the foowing properties, Initia condition : R 12 (0) = sinhηp 12, (2.5) Unitarity reation : R 12 (u)r 21 ( u) = ρ 1 (u) id, (2.6) Crossing Unitarity reation : R t 1 12 (u)m 1 R t 1 21 ( u 3η)M 1 1 = ρ 2 (u) id, (2.7) PT-symmetry : R 21 (u) = R t 1t 2 12 (u), (2.8) Periodicity : R 12 (u+iπ) = R 12 (u). (2.9) Here R 21 (u) = P 12 R 12 (u)p 12 with P 12 being the usua permutation operator and t i denotes transposition in the i-th space. The functions ρ 1 (u), ρ 2 (u) and the crossing matrix M are given by ρ 1 (u) = sinh(u η)sinh(u+η), (2.10) ρ 2 (u) = sinh(u)sinh(u+3η), (2.11) e 4η M = e 2η. (2.12) 1 4

5 It is easy to check that the R-matrix (2.1) aso has the foowing properties M 1 M 2 R 12 (u)m 1 1 M 1 2 = R 12 (u). (2.13) Let usintroduce therefection matrixk (u)anditsduaonek + (u). Theformersatisfies the refection equation (RE) and the atter satisfies the dua RE R 12 (u 1 u 2 )K 1 (u 1 )R 21 (u 1 +u 2 )K 2 (u 2 ) = K 2 (u 2 )R 12 (u 1 +u 2 )K 1 (u 1 )R 21 (u 1 u 2 ), (2.14) R 12 (u 2 u 1 )K + 1 (u 1 )M 1 1 R 21 ( u 1 u 2 3η)M 1 K + 2 (u 2 ) = K + 2 (u 2)M 1 2 R 12( u 1 u 2 3η)M 2 K + 1 (u 1)R 21 (u 2 u 1 ). (2.15) In this paper we consider the generic non-diagona K-matrices K (u) found in [44, 45, 46]. There are three kinds of refecting K-matrix: e u sinh(ζ u)+ce 2u sinh(2u) 0 0 (I) : K (u) = 0 e u sinh(ζ u) c 1 sinh(2u),(2.16) 0 c 2 sinh(2u) e u sinh(ζ +u) with the constraint c 2 = c 1 c 2 +ce ζ. Thus the four boundary parameters c, c 1 c 2 and ζ are not independent with each other. e u sinh(ζ u) 0 c 1 sinh(2u) (II) : K (u) = 0 e u sinh(ζ u)+csinh(2u) 0, (2.17) c 2 sinh(2u) 0 e u sinh(ζ +u) with the constraint c 2 = c 1 c 2 +ce ζ. (III) : K (u) = e u sinh(ζ u) c 1 sinh(2u) 0 c 2 sinh(2u) e u sinh(ζ +u) e u sinh(ζ +u)+ce 2u sinh(2u) (2.18), 5

6 with the constraint c 2 = c 1 c 2 ce ζ. The dua non-diagona refection matrix K + (u) is given by K + (u) = MK ( u 3η/2) (ζ,c,c1,c 2 ) (ζ,c,c 1,c 2 ). (2.19) In order to construct the mode s Hamitonian of the system, we first introduce the rowto-row (or one-row) monodromy matrices T 0 (u) and ˆT 0 (u) T 0 (u) = R 0N (u θ N )R 0N 1 (u θ N 1 ) R 01 (u θ 1 ), (2.20) ˆT 0 (u) = R 10 (u+θ 1 )R 20 (u+θ 2 ) R N0 (u+θ N ), (2.21) where {θ j,j = 1,,N} are the inhomogeneous parameters and N is the number of sites. The one-row monodromy matrices are the 3 3 matrices in the auxiary space 0 and their eements act on the quantum space V N. For the system with open boundaries, we need to define the doube-row monodromy matrix T 0 (u) Then the transfer matrix of the system is constructed as [7] T 0 (u) = T 0 (u)k 0 (u)ˆt 0 (u). (2.22) t(u) = tr 0 {K + 0 (u)t 0(u)}. (2.23) FromtheQYBE(2.4),RE(2.14)andduaRE(2.15),onecanprovethatthetransfermatrices with different spectra parameters commute with each other, [t(u), t(v)] = 0. Therefore, t(u) serves as the generating functiona of a the conserved quantities of the system. The mode Hamitonian can be constructed by taking the derivative of the ogarithm of the transfer matrix of the system 3 Fusion H = sinhη nt(u) u=0,{θj }=0. (2.24) u Foowing [33], we appy the fusion technique [47, 48, 49] to study the present mode. The fusion procedure wi ead to the desired operator identities to determine the spectrum of the 6

7 transfer matrix t(u) given by (2.23). For this purpose, et us introduce the foowing vectors in the tensor space V V simiary as [36] Φ (1) 12 = 1 2eη coshη ( 1,2 eη 2,1 ), Φ (2) 12 = 1 2eη coshη ( 1,3 eη 3,1 ), (3.1) Φ (3) 12 = 1 2eη coshη ( 2,3 eη 3,2 ), in the tensor space V V and Φ 123 = 1 2e 3η (2coshη +cosh3η) ( 1,2,3 eη 1,3,2 +e 2η 3,1,2 e η 2,1,3 +e 2η 2,3,1 e 3η 3,2,1 ), (3.2) in the tensor space V V V. The associated projectors 2 are P 12 = Φ(1) 12 Φ(1) 12 + Φ(2) 12 Φ(2) 12 + Φ(3) 12 Φ(3) 12, (3.3) P 123 = Φ 123 Φ 123. (3.4) Direct cacuation shows that the R-matrix given by (2.1) at some degenerate points are proportiona to the projectors, R 12 ( η) = P 12 S 12, R 12 ( η)r 13 ( 2η)R 23 ( η) = P 123 S 123, (3.5) where the diagona matrices S 12 and S 123 are given by S 12 = sinh2η Diag[1,e η,e η,e η,1,e η,e η,e η,1], (3.6) S 123 = 2sinh2ηsinh 2 η(2coshη+cosh3η) Diag[1,1,1,1,1,e 3η,1,e η,1,1,1,e η, 1,1,1,e η,1,1,1,e η,1,e 3η,1,1,1,1,1]. (3.7) The fused transfer matrices are defined as t m (u) = tr 12 m {K + <12 m> (u)t <12 m>(u)k <12 m> (u)ˆt <12 m> (u)}, m = 1,2,3, (3.8) 2 We note that in contrast to most of rationa modes, here P12 P 21. Therefore, the orders of sub-indices in (3.9)-(3.14) are crucia (c.f. [33]). 7

8 where K + 12 m (u) = K+ <2 m> (u η)m 1 2 R 1m( 2u+(m 1)η 3η) R 1m 1 ( 2u+(m 2)η 3η) R 12 ( 2u+η 3η)M 2 K 1 + (u), (3.9) K + <12 m> (u) = P 12 m K+ 12 m (u)p mm 1 1, (3.10) K 12 m(u) = K 1 (u)r 21 (2u η) R 31 (2u 2η) R m1 (2u (m 1)η)K <2 m>(u η), (3.11) K <12 m> (u) = P mm 1 1 K 12 m (u)p 12 m, (3.12) T <12 m> (u) = P mm 1 1 T 1(u)T 2 (u η) T m (u (m 1)η)P mm 1 1, (3.13) ˆT <12 m> (u) = P 12 mˆt 1 (u)ˆt 2 (u η) ˆT m (u (m 1)η)P 12 m, (3.14) and the notation t 1 (u) = t(u) is used. By repeatedy using the QYBE (2.4), the RE (2.14), the dua RE (2.15) and the definition (3.8), one can prove that a these fused transfer matrices are commutative with each other [t m (u),t k (v)] = 0, m,k = 1,2,3. (3.15) Thus they have the common eigenstates. Furthermore, we find that the transfer matrix given by (3.8) satisfies the foowing operator production identities m t(±θ j )t m (±θ j η) = t m+1 (±θ j ) ρ 1 2 (±2θ j kη), j = 1,...,N, m = 1,2,(3.16) k=1 t 2 (±θ j +η) = 0, j = 1,,N. (3.17) We note that the fused transfer matrix t 3 (u) equas to its quantum determinant mutipying the unity matrix. Thus the operators production identities (3.16) are cosed. The expicit form of the fused transfer matrix t 3 (u) reads t 3 (u) = q (u) id = q {T(u)} q {ˆT(u)} q {K + (u)} q {K (u)} id, (3.18) where q {T(u)}, q {ˆT(u)}, q {K + (u)} and q {K (u)} are the quantum determinants of the matrices T(u), ˆT(u), K + (u) and K (u), respectivey. The quantum determinants of the one-row monodromy matrices are q {T(u)} = N sinh(u θ +η)sinh(u θ η)sinh(u θ 2η), (3.19) =1 8

9 N q {ˆT(u)} = sinh(u+θ +η)sinh(u+θ η)sinh(u+θ 2η). (3.20) =1 The quantum determinant of the refecting matrix (I) given by (2.16) is q {K (u)} = (e u η sinh(ζ u+η)+ce 2u 2η sinh(2u 2η)) (sinh(ζ u)sinh(ζ +u) c 1 c 2 sinh(2u)sinh(2u)) sinh(2u 2η)sinh(2u 3η)sinh(2u 4η) = (e u η sinh(ζ u+η)+ce 2u 2η sinh(2u 2η)) (e u sinh(ζ u)+ce 2u sinh(2u)) (e u sinh(ζ +u) ce 2u sinh(2u)) sinh(2u 2η)sinh(2u 3η)sinh(2u 4η). (3.21) The quantum determinant of the refecting matrix (II) given by (2.17) is q {K (u)} = (e u η sinh(ζ u+η)+csinh(2u 2η)) (sinh(ζ u)sinh(ζ +u) c 1 c 2 sinh(2u)sinh(2u)) sinh(2u 2η)sinh(2u 3η)sinh(2u 4η) = (e u η sinh(ζ u+η)+csinh(2u 2η)) (e u sinh(ζ u)+csinh(2u)) (e u sinh(ζ +u) csinh(2u)) sinh(2u 2η)sinh(2u 3η)sinh(2u 4η). (3.22) The quantum determinant of the refecting matrix (III) given by (2.18) is q {K (u)} = (e u+η sinh(ζ +u η)+ce 2u+2η sinh(2u 2η)) (sinh(ζ u)sinh(ζ +u) c 1 c 2 sinh(2u)sinh(2u)) sinh(2u 2η)sinh(2u 3η)sinh(2u 4η) = (e u+η sinh(ζ +u η)+ce 2u+2η sinh(2u 2η)) (e u sinh(ζ u) ce 2u sinh(2u)) (e u sinh(ζ +u)+ce 2u sinh(2u)) sinh(2u 2η)sinh(2u 3η)sinh(2u 4η). (3.23) 9

10 The quantum determinant of the dua refecting matrices K + (u) can be obtained by the mapping q {K + (u)} = e 6η q {K ( u+ η 2 )} (c,ζ c,ζ ). Then the equation (3.18) can be proved easiy based on the facts T <123> (u) = q {T(u)}P 321, (3.24) ˆT <123> (u) = q {ˆT(u)}P 123, (3.25) K <123> (u) = q{k (u)} Φ 321 Φ 123, (3.26) K + <123>(u) = q {K + (u)} Φ 123 Φ 321. (3.27) Form the definition of fused transfer matrices (3.8), the corresponding asymptotic behaviors can be cacuated directy. Obviousy, different refection parameters wi give different asymptotic behaviors. Without osing the generaity, we consider the case corresponding to the refection matrices K ± (u) given by (2.16) and (2.19) and the detais for the resuts for the other cases wi be presented in Appendix. Then the asymptotic behaviors read t 1 (u) u + = 1 [ e η e 2ηQ(1) 4 N+1e(2N+4)u+3η c + (c 1 c 2 +c 1 c 2e 2η )e Nη e ηq(1)] +, (3.28) t 1 (u) u = 1 4 N+1e (2N+4)u 3η [ c c 1c 2 c eη e 2ηQ(1) + (c 1 c 2 +c 1 c 2e 2η )e Nη e ηq(1)] +, (3.29) t 2 (u) u + = 1 [ c c 1c 2 4 2N+3e4(N+3)u+4η c c eη e 2ηQ(1) + (c 1 c 2 +c 1 c 2e 2η )e Nη e ηq(1)] +, (3.30) t 2 (u) u = 1 [ c c 1c 2 e η e 2ηQ(1) 4 2N+3e 4(N+3)u 2η c c + (c 1 c 2 +c 1 c 2e 2η )e Nη e ηq(1)] +, (3.31) where the operator Q (1) is Q (1) = N =1 E 11, E 11 = (3.32) 10

11 In the derivation, the reations c(c e ζ ) = c 1 c 2 and c (c e ζ ) = c 1c 2 are used. It is remarked that the non-diagona K-matrices (given by (2.16) and (2.19)) ony break two of the origina three U(1)-symmetries for the diagona K-matrices or periodica case, and that the system sti has a remaining U(1) -symmetry which is generated by the operator Q (1). The fused transfer matrices t m (u) have other usefu properties. For exampe, their vaues at some specia points can be cacuated directy by using the properties of the R-matrix and the refection matrices K ±. We ist them in the foowing t m (u+iπ) = t m (u), m = 1,2, (3.33) N t(0) = ρ 1 ( θ )tr{k + (0)}K (0) id, (3.34) =1 t( iπ 2 ) = N ( 1)N ρ 1 ( θ + iπ 2 )tr{k+ ( iπ 2 )}K ( iπ ) id, (3.35) 2 =1 t( 3η N 2 ) = ρ 2 ( θ 3η 2 )tr{k ( 3η 2 )M}M 1 K + ( 3η ) id, (3.36) 2 =1 t( 3η 2 + iπ 2 ) = N ( 1)N ρ 2 ( θ 3η 2 iπ 2 )tr{k ( 3η 2 + iπ 2 )M} =1 M 1 K + ( 3η 2 + iπ ) id, (3.37) 2 t 2 ( η { } 2 ) = tr 12 P12 K+ 2 ( η 2 )M 1 2 R 12( 3η)M 2 K 1 + (η 2 )R 12(0)P12 ρ K (η 2 ) t 2 ( η 2 + iπ 2 ) = tr 12 N ρ 1 ( η 2 θ )ρ 1 ( η 2 θ ) id, (3.38) =1 { P12 K+ 2 ( η 2 + iπ 2 )M 1 2 R 12( 3η)M 2 K 1 + (η 2 + iπ } 2 )R 12(0)P12 ρ K (η 2 + iπ N 2 ) ρ 1 ( η 2 θ + iπ 2 )ρ 1( η 2 θ + iπ ) id, (3.39) 2 =1 t 2 ( η) = tr 12 { P 12 R 12 (0)K 1 ( η)r 21( 3η)K 2 ( 2η)M 1M 2 P 12 t 2 ( η + iπ 2 ) = tr 12 } ρ + K (η 2 ) N ρ 2 ( θ 2η)ρ 2 ( θ η) id, (3.40) =1 { P12 R 12(0)K1 ( η + iπ 2 )R 21( 3η)K2 ( 2η + iπ } 2 )M 1M 2 P12 11

12 t 2 (0) = b( η)b( 2η)K (0) ρ + K (η 2 + iπ N 2 ) ρ 2 ( θ 2η + iπ 2 )ρ 2( θ η + iπ ) id, (3.41) 2 =1 N ρ 1 ( θ )tr{k + (0)}t( η), (3.42) =1 t 2 ( iπ 2 ) = b( η)b( 2η)K ( iπ N 2 )( 1)N ρ 1 ( iπ 2 θ )tr{k + ( iπ 2 )}t( η + iπ ), (3.43) 2 =1 t 2 ( η 2 ) = b( η)b( 2η)M 1 K + ( 3η N 2 ) ρ 2 ( θ 3η 2 )tr{k ( 3 2 η)m}t( η 2 ),(3.44) =1 t 2 ( η 2 + iπ 2 ) = b( η)b( 2η)M 1 K + ( 3η 2 + iπ N 2 )( 1)N ρ 2 ( θ 3η 2 + iπ 2 ) tr{k ( 3η 2 + iπ 2 )M}t( η 2 + iπ ), (3.45) 2 t 2 (η) = t 2 (η + iπ 2 ) = t 2( 3η 2 ) = t 2( 3η 2 + iπ ) = 0, (3.46) 2 where the notations ρ K (u) and ρ+ K (u) are defined as ρ K (u) = K (u)k ( u), ρ + K (u) = ρ K (u) c,ζ,c 1,c 2 c,ζ,c 1,c 2. (3.47) In the derivation, we have used the reations K ± (u+πi) = K ± (u), R(u+πi) = R(u), (3.48) T 0 (u)ˆt 0 ( u) = =1 N ρ 1 (u θ ) id, (3.49) =1 T t 0 0 (u)m 0ˆTt 0 0 ( u 3η)M 1 0 = 4 Functiona reations N ρ 2 (u θ ) id. (3.50) Because the fused transfer matrices t m (u) commute with each other, they have the common eigenstates. Let Ψ be a common eigenstate of t m (u), which dose not depend upon u, with the eigenvaues Λ m (u), =1 t m (u) Ψ = Λ m (u) Ψ, m = 1,2,3. Again, we use the notation Λ 1 (u) = Λ(u), which represents the eigenvaue of transfer matrix t(u) given by (2.23). The Λ(u), as an entire function of u, is a trigonometric poynomia 12

13 of degree 2N + 4, which can be competey determined by 2N + 5 conditions. The Λ 2 (u), as an entire function of u, is a trigonometric poynomia of degree 4N + 12, which can be competey determined by 2N +13 conditions 3. From the operator production identities (3.16) and (3.33), we have Λ m (u+iπ) = Λ m (u), m = 1,2, (4.1) m Λ(±θ j )Λ m (±θ j η) = Λ m+1 (±θ j ) ρ 1 2 (±2θ j kη), j = 1,...,N, m = 1,2, (4.2) k=1 Λ 2 (±θ j +η) = 0, j = 1,...,N, (4.3) Λ 3 (u) = q (u). (4.4) The vaues of Λ(u) at the specia points 0, iπ 2, 3η 2, 3η 2 + iπ 2, (4.5) shoud be the same as those given by (3.34)-(3.37) of the transfer matrix t(u). At the same time, the vaues of Λ 2 (u) at the specia points 0, iπ 2, η 2, η 2 + iπ 2, η, η + iπ 2, η 2, η 2 + iπ 2, η, η + iπ 2, 3η 2, 3η 2 + iπ 2, (4.6) shoud be the same as those given by (3.38)-(3.46) of the fused transfer matrix t 2 (u). The asymptotic behaviors of Λ m (u) can be obtained by acting the operators in (3.28)- (3.31) on the corresponding eigenstates. The asymptotic behaviors (3.28)-(3.31) aows us to decompose the whoe Hibert space H into N +1 subspaces, i.e., H = N M=0 H(M) according to the action of the operator Q (1) given by (3.32): Q (1) H (M) = M H (M), M = 0,1,,N. (4.7) The commutativity of the transfer matrices and the operator Q (1) impies that each of the subspace is invariant under t m (u). Hence thewhoe set of eigenvaues of the transfer matrices can be decomposed into N +1 series. Acting the operators in (3.28)-(3.31) on any subspace H (M), we obtain the asymptotic behaviors of the corresponding Λ m (u) [ ] Λ 1 (u) u + = e(2n+4)u+3η e 2Mη+η +(c 4 N+1 c 1 c 2 +c 1 c 2e 2η )e (N M)η 3 It is noted that the reations (3.17) give the other 2N conditions. +, (4.8) 13

14 [ Λ 1 (u) u = e (2N+4)u 3η c c ] 1c 2 4 N+1 c e 2Mη+η +(c 1 c 2 +c 1c 2 e 2η )e (M N)η +, (4.9) Λ 2 (u) u + = 1 [ c c 1c 2 4 2N+3e4(N+3)u+4η c c e 2Mη+η +(c 1 c 2 +c 1 c 2e 2η )e (M N)η] +, (4.10) Λ 2 (u) u = 1 [ c c 1c 2 e 2Mη+η 4 2N+3e 4(N+3)u 2η c c +(c 1 c 2 +c 1c 2 e 2η )e (N M)η] +. (4.11) Therefore, the functiona reations (4.1)-(4.4), the vaues at specia points (4.5)-(4.6) and the asymptotic behaviors 4 (4.8)-(4.11) can provide us sufficient conditions to competey determine the corresponding eigenvaues Λ m (u). 5 Nested inhomogeneous T Q reation Now we construct the eigenvaues Λ m (u) of the fused transfer matrices t m (u). For simpicity, we define some functions b 0 (u) = N sinh(u θ j )sinh(u+θ j ), a 0 (u) = b 0 (u+η), (5.1) j=1 L k Q (k) (u) = =1 sinh(u λ (k) )sinh(u+λ (k) +kη), k = 1,2, (5.2) where L 1 and L 2 are non-negative integers. Due to the survived U(1) conserved charge Q (1) in the system, the number of one kind of Bethe roots can be chosen as M, which is simiar as the agebraic Bethe ansatz. Without osing generaity, we put L 2 = M. In order to construct the eigenvaues of the fused transfer matrices, we introduce three z(u) functions Here z m (u) is defined as z 1 (u) = z 1 (u)+x 1 (u), z 2 (u) = z 2 (u), z 3 (u) = z 3 (u). (5.3) z m (u) = sinh(2u) sinh(2u + 3η) sinh(2u+(m 1)η)sinh(2u+mη) K(m) (u)b 0 (u) Q(m 1) (u+η)q (m) (u η), Q (m 1) (u)q (m) (u) m = 1,2,3, (5.4) 4 It is remarked that for eiptica integrabe modes asymptotic behaviors such as (4.8)-(4.11) wi be repaced by the associated quasi-periodicities of the fused transfer matrices (for an exampe, see [31] (or [32]) for the XYZ cosed chain (or the XYZ open chain)). 14

15 with the notations Q (0) (u) = b 0 (u), Q (3) (u) = 1 and x 1 (u) is defined as x 1 (u) = sinh(2u)sinh(2u+3η)a 0 (u)b 0 (u) f 1(u)Q (2) ( u η), (5.5) Q (1) (u) where K (m) (u) are the decompositions of the quantum determinant and f 1 (u) is a function which wi be determined ater. The nested functiona T Q ansatz is expressed as Λ(u) = 3 z i1 (u) i 1 =1 = sinh(2u+3η) sinh(2u+η) K(1) (u)a 0 (u) Q(1) (u η) Q (1) (u) + sinh(2u)sinh(2u+3η) sinh(2u+η)sinh(2u+2η) K(2) (u)b 0 (u) Q(1) (u+η)q (2) (u η) Q (1) (u)q (2) (u) + sinh(2u) sinh(2u+2η) K(3) (u)b 0 (u) Q(2) (u+η) Q (2) (u) + sinh(2u)sinh(2u+3η)a 0 (u)b 0 (u) f 1(u)Q (2) ( u η), (5.6) Q (1) (u) [ ] Λ 2 (u) = ρ 2 (2u η) z i1 (u) z i2 (u η) x 1 (u)z 2 (u η) Λ 3 (u) = 1 i 1 <i 2 3 = ρ 2 (2u η)b 0 (u η) { sinh(2u 2η)sinh(2u+3η) K (1) (u)k (2) (u η)a 0 (u) Q(2) (u 2η) sinh(2u) sinh(2u η) Q (2) (u η) + sinh(2u 2η)sinh(2u+3η) K (1) (u)k (3) (u η)a 0 (u) Q(1) (u η)q (2) (u) sinh(2u + η) sinh(2u) Q (1) (u)q (2) (u η) + sinh(2u 2η)sinh(2u+3η) sinh(2u+η)sinh(2u+2η) K(2) (u)k (3) (u η)b 0 (u) Q(1) (u+η) Q (1) (u) + sinh(2u 2η) sinh(2u) Q(2) (u)k (3) (u η) Q (2) (u η) sinh(2u)sinh(2u+3η)a 0 (u)b 0 (u)f 1 (u) Q(2) ( u η) Q (1) (u) }, (5.7) 3 ρ 2 (2u kη)z 1 (u)z 2 (u η)z 3 (u 2η), (5.8) k=1 15

16 where the non-negative integer L 1 is L 1 = N +M +6. (5.9) Because the eigenvaues Λ m (u) are the trigonometric poynomias, the residues of right hand sides of Eqs.(5.6)-(5.8) shoud be zero, which gives the constraints of the Bethe roots {λ (r) } thus the Bethe ansatz equations. The Bethe ansatz equations obtained from the reguarity of Λ(u) shoud be the same as that obtained from the reguarity of Λ 2 (u). The function Q (r) (u) has two zero points, λ (r) and λ (r) obtained from these two points aso shoud be the same, which requires rη. The Bethe ansatz equations K (r) (u)k (r) ( u rη) = K (r+1) (u)k (r+1) ( u rη), r = 1,2. (5.10) We note that Λ 3 (u) is a trigonometric poynomia automaticay. The fact that Λ 3 (u) shoud be the quantum determinant requires K (1) (u)k (2) (u η)k (3) q {K (u)} q {K + (u)} (u 2η) = 3 (5.11) k=1sinh(2u+kη)sinh(2u (k +1)η). The consistency of Bethe ansatz equations aso require that the function f 1 (u) has the crossing symmetry f 1 (u) = f 1 ( u η). (5.12) Furthermore, the eigenvaues Λ m (u) shoud satisfy the functiona reations (4.2). This gives other constraints of the function f 1 (u). Considering a the above requirements, we parameterize the function f 1 (u) as f 1 (u) = ( 1 4 )3 hsinh(2u)sinh 2 (2u+η)sinh(2u+2η)sinh(2u η)sinh(2u+3η), (5.13) where h is a constant which is determined by the asymptotic behaviors of the Λ m (u). Now, we are ready to give the Bethe ansatz equations as foowing 1+ sinh(2λ(1) ) sinh(2λ (1) +2η) K (2) (λ (1) )b 0 (λ (1) ) K (1) (λ (1) )a 0 (λ (1) ) Q (1) (λ (1) +η)q (2) (λ (1) η) Q (1) (λ (1) η)q (2) (λ (1) ) = hsinh 2 (2λ (1) )sinh 3 (2λ (1) +η)sinh(2λ (1) +2η)sinh(2λ (1) +3η) sinh(2λ (1) η) b 0 (λ (1) )Q (2) (λ (1) η) 4 3 K (1) (λ (1) )Q (1) (λ (1) η), = 1,...,M +N +6, (5.14) sinh(2λ (2) k +3η) sinh(2λ (2) k +η) K (2) (λ (2) k ) K (3) (λ (2) Q(1) (λ k ) Q (1) (λ (2) (2) k +η)q (2) (λ (2) k η) k )Q(2) (λ (2) 16 k +η) = 1, k = 1,...,M.(5.15)

17 It is easy to check that the Bethe ansatz equations (5.14)-(5.15) guarantee the reguarities of the ansatz Λ(u) given by (5.6) and the ansatz Λ 2 (u) given by (5.7). Moreover, the ansatz (5.6)-(5.8) indeed satisfy the function reations (4.2). The eft tasks are to determine the vaue of h in the ansatz (5.6)-(5.8), which can be done by anayzing the asymptotic behaviors, and to check the consistency of vaues at the specia points (4.5)-(4.6). Because there are three kinds of boundary refection matrices and they give the different behaviors, et us consider them one by one. For the refection matrices K ± given by (2.16) and (2.19), the decomposition K (i) (u) can be chosen as K (1) (u) = e 2η (e u sinh(ζ +u) ce 2u sinh2u) (e u η 2 sinh(ζ u+ η 2 )+c e 2(u η 2 ) sinh2(u η )), (5.16) 2 K (2) (u) = e 2η (e u+η sinh(ζ u η)+ce 2(u+η) sinh2(u+η)) (e u 3η 2 sinh(ζ +u+ 3η 2 ) c e 2(u+3η 2 ) sinh2(u+ 3η 2 )), (5.17) K (3) (u) = K (2) (u). (5.18) From the asymptotic behaviors (4.8)-(4.11) of the corresponding trigonometric poynomias, we arrive at the vaue of h h = c e (M+N)η+15η +c c 1c 2 c e (M+N)η 13η (c 1 c 2 +c 1c 2 e 2η ). (5.19) Then one can check that the vaues of ansatz (5.6)-(5.7) at the specia points (4.5)-(4.6) are the same as those of the corresponding fused transfer matrices, and we finish our construction. Taking thehomogeneous imit{θ j = 0,j = 1,...,N}, weconcude thatthet Qreation Λ(u) given by (5.6) is the eigenvaue of the transfer matrices t(u) of the trigonometric SU(3) open spin chain with the most genera off-diagona integrabe boundary conditions. The energy of the Hamitonian (2.24) reads E = sinhη nλ(u) u=0,{θj }=0, (5.20) u where the Bethe roots shoud satisfy the Bethe ansatz equations (5.14)-(5.15). Above resuts can be reduced to the diagona boundaries ones obtained by the agebraic Bethe ansatz [50, 51, 52]. 17

18 Numerica soutions of the BAEs (5.14)-(5.15) for sma size 5 with a random choice of η impy that the BAEs indeed give the compete soutions of the mode. Here we present the resut for the N = 2 case: the numerica soutions of the BAEs for the N = 2 case are shown in Tabe 1, whie the cacuated Λ(u) curves for the case of N = 2 are shown in Figure Λ(u) n=1 n=4 n=5 n=6 n=2 n=3 n=7 n=8 n= u Figure 1: (coor onine) Λ(u) vs. u for the case of N = 2 and {θ j = 0}. The curves cacuated from T Q reation (5.6) and the nested BAEs (5.14)-(5.15) are exacty the same as those obtained from the exact diagonaization of the transfer matrix t(u). 6 Diagona boundary case When the parameters c,c 1,c 2,c,c 1,c 2 in the refection matrices K ± (u) given by (2.16) and (2.19) vanish, the corresponding K-matrices become diagona ones. Let us denote them by K ± (u), namey, K (u) = e u sinh(ζ u) e u sinh(ζ u) 0 e u sinh(ζ +u), (6.1) K + (u) = M K ( u 3η/2) ζ ζ. (6.2) 5 It is sti an interesting open probem to investigate the root pattern of the BAEs (5.14)-(5.15) associated with the inhomogeneous T Q reation. Nevertheess, a standard method to study the thermodynamic imit was deveoped in [34] by considering a sequence of discrete η vaues at which the inhomogeneous terms in the BAEs vanish and in the thermodynamic imit these discrete η vaues become dense. 18

19 19 Tabe 1: Soutions of BAEs (5.14)-(5.15) for the case of N = 2,η = 0.3,{θ j = 0},ζ = 0.1,c = 1.0,c 1 = 0.5,ζ = 0.1,c = 0.5,c 1 = 0.7. The symbo n indicates the number of the eigenvaues, and E n is the corresponding eigenenergy. The energy E n cacuated from (5.20) is the same as that from the exact diagonaization of the Hamitonian (2.24). λ (1) 1 λ (1) 2 λ (1) 3 λ (1) 4 λ (1) 5 λ (1) 6 λ (1) i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i λ (1) 8 λ (1) 9 λ (1) 10 λ (2) 1 λ (2) 2 E n n i i i i i i i i i i i i i i i i i i i i i

20 Moreover, the corresponding decomposition {K (i) (u)} in (5.16)-(5.18), if denoted by K (i) (u), are given by K (1) (u) = e 3 2 η sinh(ζ +u)sinh(ζ u+ η ), (6.3) 2 K (2) (u) = e 3 2 η sinh(ζ u η)sinh(ζ +u+ 3η 2 ), (6.4) K (3) (u) = K (2) (u). (6.5) Then the corresponding T Q reations (5.6)-(5.7) become the usua homogeneous ones and now are given by Λ(u) = sinh(2u+3η) sinh(2u+η) K (1) (u)a 0 (u) Q (1) (u η) Q (1) (u) + sinh2usinh(2u+3η) sinh(2u+η)sinh(2u+2η) K (2) (u)b 0 (u) Q (1) (u+η) Q (2) (u η) Q (1) (u) Q (2) (u) + sinh2u sinh(2u+2η) K (3) (u)b 0 (u) Q (2) (u+η), (6.6) Q (2) (u) Λ 2 (u) = ρ 2 (2u η)b 0 (u η) { sinh(2u 2η)sinh(2u+3η) sinh2usinh(2u η) + sinh(2u 2η)sinh(2u+3η) sinh(2u+η)sinh2u K (1) (u) K (2) (u η)a 0 (u) Q (2) (u 2η) Q (2) (u η) K (1) (u) K (3) (u η)a 0 (u) Q (1) (u η) Q (2) (u) Q (1) (u) Q (2) (u η) + sinh(2u 2η)sinh(2u+3η) sinh(2u+η)sinh(2u+2η) K (2) (u) K (3) (u η)b 0 (u) Q (1) (u+η) Q (1) (u) Here the corresponding Q-functions are M k Q (k) (u) = =1 }. (6.7) sinh(u λ (k) )sinh(u+λ (k) +kη), 0 M 2 N, 0 M 1 M 2. (6.8) The resuting homogeneous reation (6.6) recovers that obtained by the agebraic Bethe ansatz method [52], whie the reference state is chosen as vac >= , (6.9) and the associated creation operators are the off-diagona matrix eements of the 3-rd row of the doube-row monodromy matrix T(u) given by (2.22). 20

21 7 Concusions In this paper, we study the exact soution of the anisotropic quantum spin chain with generic open boundary conditions and associated with SU q (3) agebra. After giving the off-diagona refection matrixes, by using the fusion technique, we obtain some cosed operator identities among the transfer matrices, the degenerate points and the corresponding asymptotic behaviors. Based on them, we construct the nested inhomogeneous T Q reations and the Bethe ansatz equations. These resuts can be generaized to the higher rank case. Moreover, when the boundary parameters take specia vaues corresponding to the diagona refection matrices, our resuts recover those previousy obtained by the conventiona Bethe ansatzs. Acknowedgments We woud ike to thank Y. Wang for his vauabe discussions and continuous encouragement. The financia supports from the Nationa Natura Science Foundation of China (Grant Nos , , , and ), BCMIIS and the Strategic Priority Research Program of the Chinese Academy of Sciences are gratefuy acknowedged. Appendix: Asymptotic behaviors for the other two cases refecting matrices When the refection matrices K ± are given by (2.17) and (2.19), the asymptotic behaviors of fused transfer matrices t m (u) are t 1 (u) u + = 1 4 N+1e(2N+4)u+3η [ c c 1c 2 c e2η e 2ηQ(2) + (c 1 c 2 +c 1c 2 e 4η )e Nη e ηq(2)] +, (A.1) t 1 (u) u = 1 [ c c 1c 2 e 2η e 2ηQ(2) 4 N+1e (2N+4)u 3η c +(c 1 c 2 +c 1 c 2e 4η )e Nη e ηq(2)] +, t 2 (u) u + = 1 c c 1c 2 4 2N+3e4(N+3)u+5η c [ c e 2η e 2ηQ(2) +(c 1 c 2 +c 1 c 2e 4η )e Nη e ηq(2)] +, (A.2) (A.3) 21

22 t 2 (u) u = 1 [ c c 1c 2 4 2N+3e 4(N+3)u η c c e2η e 2ηQ(2) +(c 1 c 2 +c 1 c 2e 4η )e Nη e ηq(2)] +, (A.4) where the operator Q (2) is Q (2) = N =1 E 22, E 22 = (A.5) The asymptotic behaviors of fused transfer matrices t m (u) associated with the refection matrices K ± given by (2.18) and (2.19) are t 1 (u) u + = 1 4 N+1e(2N+4)u+5η [ c c 1c 2 c eη e 2ηQ(3) +(c 1 c 2 +c 1 c 2e 2η )e Nη e ηq(3)] +, t 1 (u) u = 1 [ e η e 2ηQ(3) 4 N+1e (2N+4)u η c (A.6) +(c 1 c 2 +c 1 c 2e 2η )e Nη e ηq(3)] +, t 2 (u) u + = 1 c c 1c 2 4 2N+3e4(N+3)u+8η c [ c c 1c 2 c e η e 2ηQ(3) +(c 1 c 2 +c 1c 2 e 2η )e Nη e ηq(3)] +, t 2 (u) u = 1 [ c c 1c 2 4 2N+3e 4(N+3)u+2η c c eη e 2ηQ(3) +(c 1 c 2 +c 1c 2 e 2η )e Nη e ηq(3)] +, (A.7) (A.8) (A.9) where the operator Q (3) is Q (3) = N =1 E 33, E 33 = (A.10) Some remarks are in order. Simiary as the operator Q (1) given by (3.32), the above operator Q (2) (resp. Q (3) ) generates a remaining U(1)-symmetry for the corresponding mode respectivey. Moreover the eigenvaues of the operator are 0, 1,,N. Hence the whoe Hibert space can be decomposed into the subspaces abeed by its eigenvaue, on which the transfer matrices are invariant. We can cacuate the asymptotic behaviors of the corresponding transfer matrices on each subspace. 22

23 For the refection matrix K ± given by (2.17) and (2.19), the asymptotic behaviors of Λ m (u) can be obtained by acting the operator (A.1)-(A.4) on the subspace on which the eigenvaue of the operator Q (2) is M. After some cacuations, we arrive at [ Λ 1 (u) u + = e(2n+4)u+3η c c ] 1c 2 4 N+1 c e2mη+2η +(c 1 c 2 +c 1c 2 e 4η )e (N M)η Λ 1 (u) u = 1 4 N+1e (2N+4)u 3η [ c e 2Mη+2η +(c 1 c 2 +c 1c 2 e 4η )e (M N)η] +, Λ 2 (u) u + = 1 c c 1c 2 4 2N+3e4(N+3)u+5η c +(c 1 c 2 +c 1c 2 e 4η )e (M N)η] +, [ c e 2Mη+2η Λ 2 (u) u = 1 [ c c 1c 2 4 2N+3e 4(N+3)u η c c e2mη+2η +(c 1 c 2 +c 1 c 2e 4η )e (N M)η] +. +,(A.11) (A.12) (A.13) (A.14) For the refection matrix K ± given by (2.18) and (2.19), the asymptotic behaviors of Λ m (u) can be obtained by acting the operator (A.6)-(A.9) on the subspace on which the eigenvaue of the operator Q (3) is M. The finia resuts read [ Λ 1 (u) u + = e(2n+4)u+5η c c ] 1c 2 4 N+1 c e2mη+η +(c 1 c 2 +c 1c 2 e 2η )e (N M)η Λ 1 (u) u = 1 4 N+1e (2N+4)u η [ c c 1c 2 c e 2Mη+η +(c 1 c 2 +c 1 c 2e 2η )e (M N)η] +, Λ 2 (u) u + = 1 c c 1c 2 4 2N+3e4(N+3)u+8η c [ c e 2Mη+η +(c 1 c 2 +c 1c 2 e 2η )e (M N)η] +, Λ 2 (u) u = 1 [ c c 1c 2 4 2N+3e 4(N+3)u+2η c c e2mη+η +(c 1 c 2 +c 1c 2 e 2η )e (N M)η] +. +,(A.15) (A.16) (A.17) (A.18) For the refection matrix K ± given by (2.17) and (2.19), the eigenvaues of the transfer matrices t(u) and t 2 (u) are described by the ansatz (5.6) and (5.7) with the decomposition K (1) (u) = e 2η (e u sinh(ζ +u) csinh2u) 23

24 (e u η 2 sinh(ζ u+ η 2 )+c sinh2(u η 2 )), (A.19) K (2) (u) = e 2η (e u+η sinh(ζ u η)+csinh2(u+η)) K (3) (u) = K (2) (u), (e u 3η 2 sinh(ζ +u+ 3η 2 ) c sinh2(u+ 3η 2 )), (A.20) (A.21) and the parameter h takes the vaue of h = c e (M+N)η 12η +c c 1c 2 c e(m+n)η+16η (c 1 c 2 +c 1 c 2e 4η ). (A.22) For the refection matrix K ± given by (2.18) and (2.19), the eigenvaues of the transfer matrices t(u) and t 2 (u) are aso characterized by the ansatz (5.6) and (5.7) with the parametrization K (1) (u) = e 2η (e u sinh(ζ u) ce 2u sinh2u) (e u+η 2 sinh(ζ +u η 2 )+c e 2(u η 2 ) sinh2(u η 2 )), (A.23) K (2) (u) = e 2η (e u η sinh(ζ +u+η)+ce 2(u+η)) sinh2(u+η)) K (3) (u) = K (2) (u), (e u+3η 2 sinh(ζ u 3η 2 ) c e 2(u+3η 2 ) sinh2(u+ 3η 2 )), (A.24) (A.25) and the constant h is given by h = c c 1c 2 c e (M+N)η 11η +c c 1c 2 c e(m+n)η+17η (c 1 c 2 +c 1 c 2e 2η )e 2η. (A.26) References [1] R.J. Baxter, Exacty Soved Modes in Statistica Mechanics, Academic Press, [2] V.E. Korepin, N. M. Bogoiubov and A. G. Izergin, Quantum Inverse Scattering Method and Correation Function, Cambridge University Press, [3] R.J. Baxter, Phys. Rev. Lett. 26 (1971) 832; R.J. Baxter, Phys. Rev. Lett. 26 (1971) 834; R.J. Baxter, Ann. Phys. 70 (1967) 323. [4] L.A. Takhtadzhan and L. D. Faddeev, Rush. Math. Surveys 34 (1979)

25 [5] E.K. Skyanin and L.D. Faddeev, Sov. Phys. Dok. 23 (1978) 902. [6] F.C. Acaraz, M.N. Barber, M.T. Batcheor, R.J. Baxter and G.R.W. Quispe, J. Phys. A 20 (1987) [7] E.K. Skyanin, J. Phys. A 21 (1988) [8] J. Cao, W.-L. Yang, K. Shi and Y. Wang, Phys. Rev. Lett. 111 (2013) [9] J. Cao, H.-Q. Lin, K.-J. Shi and Y. Wang, Nuc. Phys. B 663 (2003) 487. [10] J. de Gier and P. Pyatov, J. Stat. Mech. (2004) P03002; A. Nichos, V. Rittenberg and J. de Gier, J. Stat. Mech. (2005) P03003; J. de Gier, A. Nichos, P. Pyatov and V. Rittenberg, Nuc. Phys. B 729 (2005) 387. [11] J. de Gier and F.H.L. Esser, Phys. Rev. Lett. 95 (2005) ; J. de Gier and F.H.L. Esser, J. Stat. Mech. (2006) P [12] Z. Bajnok, J. Stat. Mech. (2006) P [13] P. Baseihac, Nuc. Phys. B 754 (2006) 309. [14] P. Baseihac and K. Koizumi, J. Stat. Mech. (2007) P [15] P. Baseihac and S. Beiard, Lett. Math. Phys. 93 (2010) 213; P. Baseihac and S. Beiard, Nuc. Phys. B 873 (2013) 550. [16] S. Beiard and N. Crampé, SIGMA 9 (2013) 072. [17] S. Beiard, Nuc. Phys. B 892 (2015) 1. [18] S. Beiard and R.A. Pimenta, Nuc. Phys. B 894 (2015) 527. [19] J. Avan, S. Beiard, N. Grosjean and R.A. Pimenta, Nuc. Phys. B 899 (2015) 229. [20] E.K. Skyanin, Lecture Notes in Physics 226 (1985) 196; E.K. Skyanin, J. Sov. Math. 31 (1985) 3417; E.K. Skyanin, Prog. Theor. Phys. Supp. 118 (1995) 35. [21] H. Frahm, A. See and T. Wirth, Nuc. Phys. B 802 (2008)

26 [22] G. Niccoi, Nuc. Phys. B 870 (2013) 397; G. Niccoi, J. Phys. A 46 (2013) [23] S. Fadea, N. Kitanine and G. Niccoi, J. Stat. Mech. (2014) P [24] N. Kitanine, J.-M. Maiet and G. Niccoi, J. Stat. Mech. (2014) P [25] R.I. Nepomechie, J. Phys. A 34 (2001) 9993; R.I. Nepomechie, Nuc. Phys. B 622 (2002) 615; R.I. Nepomechie, J. Stat. Phys. 111 (2003) 1363; R.I. Nepomechie, J. Phys. A 37 (2004) 433. [26] W.-L. Yang and Y.-Z. Zhang, J. High Energy Phys. 12 (2004) 019; W.-L. Yang and Y.-Z. Zhang, J. High Energy Phys. 01 (2005) 021. [27] A. Doikou and P.P. Martins, J. Stat. Mech. (2006) P06004; A. Doikou, J. Stat. Mech. (2006) P [28] W.-L. Yang, R.I. Nepomechie and Y.-Z. Zhang, Phys. Lett. B 633 (2006) 664. [29] R.I. Nepomechie, J. Phys. A 46 (2013) [30] J. Cao, W.-L. Yang, K. Shi and Y. Wang, Nuc. Phys. B 875 (2013) 152. [31] J. Cao, S. Cui, W.-L. Yang, K. Shi and Y. Wang, Nuc. Phys. B 886 (2014) 185. [32] J. Cao, W.-L. Yang, K. Shi and Y. Wang, Nuc. Phys. B 877 (2013) 152. [33] J. Cao, W.-L. Yang, K. Shi and Y. Wang, J. High Energy Phys. 04 (2014) 143. [34] Y.-Y. Li, J. Cao, W.-L. Yang, K. Shi and Y. Wang, Nuc. Phys. B 884 (2014) 17. [35] X. Zhang, J. Cao, W.-L. Yang, K. Shi and Y. Wang, J. Stat. Mech. (2014) P [36] K. Hao, J. Cao, G.-L. Li, W.-L. Yang, K. Shi and Y. Wang, J. High Energy Phys. 06 (2014) 128. [37] Y. Wang, W.-L. Yang, J. Cao and K. Shi, Off-Diagona Bethe Ansatz for Exacty Sovabe Modes, Springer Press, [38] J.H.H. Perk and C.L. Schutz, Phys. Lett. A 84 (1981)

27 [39] J.H.H. Perk and C.L. Schutz, Famiies of commuting transfer matrices in q-state vertex modes, in Non-inear integrabe systems - cassica theory and quantum theory, eds. M. Jimbo and T. Miwa, Word Scientific, 1983, pp [40] C.L. Schutz, Physica A 122 (1983) 71. [41] J.H. H. Perk and H. Au-Yang, Yang-Baxter Equation, in Encycopedia of Mathematica Physics, eds. J.-P. Françoise, G.L. Naber and T.S. Tsun, Academic Press, [42] V.V. Bazhanov, Phys. Lett. B 159 (1985) 321. [43] M. Jimbo, Commun. Math. Phys. 102 (1986) 537. [44] A. Lima-Santos, Nuc. Phys. B 558 (1999) 637. [45] K. Shi, G. -L. Li, H. Fan and B. Y. Hou, High Energy Phys. and Nuc. Phys. 24 (2000) 11. [46] R. Maara and A. Lima-Santos, J. Stat. Mech. (2006) P [47] M. Karowski, Nuc. Phys. B 153 (1979) 244; P.P. Kuish, N.Yu. Reshetikhin and E.K. Skyanin, Lett. Math. Phys. 5 (1981) 393; P.P. Kuish and E.K. Skyanin, Lecture Notes in Physics 151 (1982) 61; A.N. Kiriov and N.Yu. Reshetikhin, J. Sov. Math. 35 (1986) 2627; A.N. Kiriov and N.Yu. Reshetikhin, J. Phys. A 20 (1987) [48] L. Mezincescu and R.I. Nepomechie, Nuc. Phys. B 372 (1992) 597. [49] Y.-K. Zhou, Nuc. Phys. B 458 (1996) 504. [50] H.J. de Vega and A. Gonzáez-Ruiz, Nuc. Phys. B 417 (1994) 553. [51] H.J. de Vega and A. Gonzáez-Ruiz, Mod. Phys. Lett. A 09 (1994) [52] G.-L. Li, R. H. Yue and B. Y. Hou, Nuc. Phys. B 586 (2000)