arxiv: v1 [math-ph] 22 Dec 2014

Transcription

1 arxiv: v1 [math-ph] 22 Dec 2014 Bethe states of the XXZ spin- 1 chain with arbitrary 2 boundary fields Xin Zhang a, Yuan-Yuan Li a, Junpeng Cao a,b, Wen-Li Yang c,d1, Kangjie Shi c and Yupeng Wang a,b2 a Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing , China b Collaborative Innovation Center of Quantum Matter, Beijing, China c Institute of Modern Physics, Northwest University, Xian , China d Beijing Center for Mathematics and Information Interdisciplinary Sciences, Beijing, , China Abstract Based on the inhomogeneous T Q relation constructed via the off-diagonal Bethe Ansatz (ODBA, the Bethe-type eigenstates of the XXZ spin- 1 2 chain with arbitrary boundary fields are constructed. It is found that by employing two sets of gauge transformations, proper generators and reference state for constructing Bethe vectors can be obtained respectively. Given an inhomogeneous T Q relation for the eigenvalue, it is proven that the resulting Bethe state is an eigenstate of the transfer matrix, provided that the parameters of the generators satisfy the associated Bethe Ansatz equations (BAEs. PACS: Pq, Vf, Pm Keywords: Spin chain; Bethe Ansatz; T Q relation; Scalar product 1 Corresponding author: wlyang@nwu.edu.cn 2 Corresponding author: yupeng@iphy.ac.cn

2 1 Introduction In this paper we focus on constructing the Bethe-type eigenstates (Bethe states of the quantum XXZ spin- 1 chain with arbitrary boundary fields, defined by the Hamiltonian 2 H = = N 1 N 1 { σ x j σj+1 x +σ y } j σy j+1 +coshησz jσj+1 z + h1 σ 1 + h N σ N { } σ x j σj+1 x +σy j σy j+1 +coshησz j σz j (coshα sinhβ σ1 z sinhα coshβ +coshθ σ1 x +isinhθ σ y 1 (coshα + sinhβ + σn z coshθ + σn x isinhθ + σ y N, (1.1 sinhα + coshβ + where σj α (α = x,y,z is the Pauli matrix on the site j along the α direction and α ±, β ±, θ ± are the boundary parameters associated with the boundary fields. The model has played a fundamental role in the study of quantum integrable system[1, 2] with boundaries. Moreover, it has many applications in the non-perturbative analysis of quantum systems appearing in string and super-symmetric Yang-Mills (SYM theories [3] (and references therein, lowdimensional condensed matter physics [4] and statistical physics [5, 6]. However, the Bethe Ansatz solution of the model for generic values of boundary fields has challenged for many years since Sklyanin s elegant work [7], and many efforts had been made [8, 9, 10, 11, 5, 12, 13, 14, 15, 16, 17] to approach this nontrivial problem. The off-diagonal Bethe Ansatz (ODBA provides an efficient method [18] for solving the eigenvalue problem of integrable models with generic integrable boundary conditions. Several long-standing models [18, 19, 20, 21, 22, 23, 24] including the XXZ spin- 1 chain 2 have since been solved via this method. The central point is to construct a proper T Q relation [1], which immediately leads to the Bethe Ansatz solution for the eigenvalues, with an extra off-diagonal (or inhomogeneous term for the eigenvalues of the transfer matrix based on their functional relations. An interesting issue in this framework is how to retrieve the Bethe states from the obtained spectrum. Indeed, significant progress has been achieved in this aspect recently. For example, based on the inhomogeneous T Q relation obtained in [19], the Bethe states of the open XXX spin chain was conjectured in [25] and then proven in [26]. Alternatively, a set of eigenstates of the inhomogeneous XXZ transfer matrix was 2

3 derived in [16, 27] via the separation of variables (SoV method [28]. However, how to get the homogeneous limit (if there is any of those SoV states is still an open problem. For the open XXZ chain, when the boundary fields are all along the z-direction (or the diagonal boundaries, the corresponding Bethe states were constructed by the algebraic Bethe Ansatz method [7, 29]. The unparallel boundary fields break the U(1-symmetry (i.e, the total spin is not conserved any more. This makes the problem of constructing Bethe vectors rather unusual because of the absence of an obvious reference state. So far, the Bethe states could only be obtained for some constrained boundary parameters. When the boundary parameters obey a constraint [8, 9], which is already in U(1-symmetry-broken case, the associated Bethe states were constructed[9] within the framework of the generalized algebraic Bethe Ansatz [1, 30]. Very recently, based on small sites analysis of the model with triangular boundaries, the corresponding Bethe states are conjectured [31] and proven in [32]. In this paper we study the Bethe states of the transfer matrix for the quantum XXZ spin- 1 chain with arbitrary boundary fields based on the inhomogeneous T Q relation of 2 the eigenvalues obtained by ODBA. The paper is organized as follows. Section 2 serves as an introduction of our notations and the ODBA solutions of the model. In section 3, after introducing the gauge transformations, we compute the associated commutation relations among the matrix elements of the two gauged double-row monodromy matrices, and their actions on the associated vacuum states. In section 4, two particular gauge transformations are chosen according to the boundary parameters of K-matrices respectively. Based on the chosen parameters of the resulting transformations, the Bethe-type eigenstates of the transfer matrix are constructed. In section 5, we summarize our results and give the concluding remarks. Some useful formulae and technical proofs are given in Appendices A-C respectively. 2 ODBA solution For the XXZ spin chain with generic boundaries, the associated R-matrix and the reflection matrix K (u [33, 34] read R(u = 1 sinh(u + η sinhu 0 0 sinhu sinh(u + η, (2.1 3

4 K (u = ( K 11 (u K 12 (u K 21(u K 22(u, K 11 (u = 2(sinh(α cosh(β cosh(u+cosh(α sinh(β sinh(u, K 22 (u = 2(sinh(α cosh(β cosh(u cosh(α sinh(β sinh(u, K 12 (u = eθ sinh(2u, K 21 (u = e θ sinh(2u, (2.2 and K + (u = K ( u η, (2.3 (α,β,θ ( α +, β +,θ + where η is the crossing parameter, and α, β, θ are the boundary parameters associated with boundary fields (see (1.1. The R-matrix is a solution of 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 K (u satisfy the following reflection equations (RE R 12 (u 1 u 2 K1 (u 1R 21 (u 1 +u 2 K2 (u 2 = K2 (u 2 R 12 (u 1 +u 2 K1 (u 1 R 21 (u 1 u 2, (2.5 and R 12 (u 2 u 1 K 1 + (u 1 R 21 ( u 1 u 2 2ηK 2 + (u 2 = K 2 + (u 2 R 12 ( u 1 u 2 2ηK 1 + (u 1 R 21 (u 2 u 1. (2.6 Weintroducethe row-to-row (orone-rowmonodromymatricest 0 (uand ˆT 0 (u,which are 2 2 matrices with elements being operators acting on the tensor space V N, T 0 (u = R 0N (u θ N R 0N 1 (u θ N 1 R 01 (u θ 1, (2.7 ˆT 0 (u = R 10 (u+θ 1 R 20 (u+θ 2 R N0 (u+θ N. (2.8 Here {θ j j = 1,,N} are the inhomogeneous parameters. For open spin chains, one needs to consider the double-row monodromy matrix U 0 (u U 0 (u = T 0 (uk0 (uˆt 0 (u. (2.9 The double-row transfer matrix t(u is thus given by t(u = tr 0 (K 0 + (uu 0(u. (2.10 4

5 The QYBE (2.4 and REs (2.5 and (2.6 lead to the fact that the transfer matrices with different spectral parameters commute with each other [7]: [t(u, t(v] = 0. Then t(u serves as the generating functional of the conserved quantities of the corresponding system, which ensures the integrability of the open spin chain. The Hamiltonian (1.1 is expressed in terms of the transfer matrix (2.10 with the K- matrices (2.2 and (2.3 by H = lnt(u u=0,θj =0 N coshη tanhη. (2.11 u It was proven in [20] that for generic {θ j } the transfer matrix given by (2.10 for arbitrary boundary parameters satisfies the following operator identities t(θ j t(θ j η = a(θ j d(θ j η id, (2.12 t( u η = t(u, t(u+iπ = t(u, (2.13 t(0 = 2 3 sinhα coshβ sinhα + coshβ + coshη N sinh(η θ l sinh(η +θ l sinh 2 id, (2.14 η l=1 t( iπ 2 = 23 coshα sinhβ coshα + sinhβ + coshη N sinh( iπ +θ 2 l +ηsinh( iπ +θ 2 l η sinh 2 id, (2.15 η l=1 lim t(u = cosh(θ θ + e ±[(2N+4u+(N+2η] u ± 2 2N+1 sinh 2N η where the functions a(u and d(u are given by id+..., (2.16 a(u = 2 2sinh(2u+2η sinh(2u+η sinh(u α cosh(u β sinh(u α + cosh(u β + Ā(u,(2.17 N sinh(u θ l +ηsinh(u+θ l +η d(u = a( u η, Ā(u = sinh 2. (2.18 η l=1 The above operator relations lead to that the corresponding eigenvalue of the transfer matrix, denoted by Λ(u, enjoys the following properties Λ(θ j Λ(θ j η = a(θ j d(θ j η, j = 1,...,N, (2.19 Λ( u η = Λ(u, Λ(u+iπ = Λ(u, (2.20 N Λ(0 = 2 3 sinh(η θ l sinh(η +θ l sinhα coshβ sinhα + coshβ + coshη sinh 2, η ( l=1

6 Λ( iπ 2 = 23 coshα sinhβ coshα + sinhβ + coshη N sinh( iπ +θ 2 l +η sinh( iπ +θ 2 l η sinh 2, (2.22 η l=1 lim Λ(u = cosh(θ θ + e ±[(2N+4u+(N+2η] u ± 2 2N+1 sinh 2N η +... (2.23 Λ(u, as an entire function of u, is a trigonometric polynomial of degree 2N +4. (2.24 Each solution of (2.19-(2.24 can be given in terms of the following inhomogeneous T Q relation [19, 35, 36] 3 Λ(u = a(u Q(u η +d(u Q(u+η Q(u Q(u 2c sinh(2usinh(2u+2η + Ā(uĀ( u η, (2.25 Q(u where c is a constant depending on the boundary parameters c = cosh(α +β +α + +β + +(1+Nη cosh(θ θ +, (2.26 and the Q-function is given by Q(u = N with the parameters {λ j } satisfying the associated BAEs sinh(u λ j sinh(u+λ j +η, (2.27 a(λ j Q(λ j η+d(λ j Q(λ j +η + 2c sinh2λ j sinh(2λ j +2ηĀ(λ jā( λ j η = 0, j = 1,...,N. (2.28 We shall show in Section 4 that for each solution of (2.19-(2.24, one can construct the corresponding Bethe-type eigenstate (see (4.10 below of the transfer matrix (2.10 with the eigenvalue given by (2.25. Therefore the relations (2.19-(2.24 (or the inhomogeneous T Q relation (2.25 indeed completely characterize the spectrum of the transfer matrix. Some remarks are in order. There exist various possible ways [19] to parameterize the solution of(2.19-(2.24, but they are all equivalent to each other because of the finite number of solutions. For generic boundary parameters, the minimal degree of the Q-polynomial is 3 The T Q relation (2.25 corresponds to the case of M = 0 in [19]. A generalization to other cases is straightforward. 6

7 N, while the degree of the Q-polynomial may be reduced to a small value in case of the inhomogeneous term (or the third term in (2.25 vanishing. In this case the T Q relation becomes a homogeneous one (the well-known Baxter s T Q relation. This happens in case of U(1 symmetry or in degenerate cases [9], for which the transfer matrix can be diagonalized in smaller blocks. 3 Gauge transformations and the associated operators A particular set of gauge transformation (the six-vertex version of the vertex-face correspondence, which have played a key role to construct the associated Bethe states, was proposed in [9]. Recently, such gauge transformation was adopted in constructing the SoV eigenstates for the open chains [27]. In this paper, we use two sets of such gauge transformation and the inhomogeneous T Q relation (2.25 to construct the Bethe states for the quantum XXZ spin- 1 chain with arbitrary boundary fields. 2 Following [9], let us introduce two column vectors as follows ( ( e [u+(α+mη] e [u+(α mη] X m (u α =, Y 1 m (u α = 1, (3.1 where α and m are two arbitrary complex parameters. For generic α and m, the two vectors are linearly independent. Thus one can introduce the following gauge matrices M m (u α = ( X m (u α, Y m (u α (, M 1 Y m (u = m (u α, (3.2 X m (u α M m (u α = ( X m+1 (u α, Y m 1 (u α ( Ỹm 1 (u α, M 1 m (u α =, (3.3 X m+1 (u α where M m (u α = ( Xm 1 (u α, Ŷ m+1 (u α X m (u α = eu+αη ( 1, e [u+(α+mη] 2sinhmη Y m (u α = eu+αη ( 1, e [u+(α mη] 2sinhmη ( Y, M 1 m (u α = m+1 (u α X m 1 (u α, (3.4, (3.5, (3.6 X m (u α = eη sinhmη sinh(m 1η X m(u α, Ỹ m (u α = eη sinhmη sinh(m+1η Y m(u α, (3.7 X m (u α = e η sinh(m+2η sinh(m+1η X m(u α, Ŷ m (u α = e η sinh(m 2η sinh(m 1η Y m(u α. (3.8 7

8 We remark that the vectors X m (u α and X m (u α only depend on α+m, while the vectors Y m (u α and Y m (u α only depend on α m, up to a scaling factor. These column and row vectors satisfy some intertwining relations [9], which are listed in Appendix A (see (A.1-(A.28 below. These relations allow us to introduce the following gauged operators and the associated K + -matrix ( A m (u α B m (u α U (m,α u = C m (u α D m (u α ( Y m (u αu(u = X m 2 ( u α Y m (u αu(uŷm( u α, (3.9 X m (u αu(u X m ( u α X m (u αu(uŷm+2( u α ( K + K + 11 (m,α u = (m,α u K+ 12 (m,α u K + 21 (m,α u K+ 22 (m,α u ( Y = m ( u αk + (ux m (u α Y m+2 ( u αk + (uy m (u α X m 2 ( u αk + (ux m (u α X m ( u αk +. (3.10 (uy m (u α With the help of the relations (A.29-(A.31, we can rewrite the transfer matrix (2.10 in terms of the above gauged operators and K-matrix, namely, t(u = tr { K + (uu (u } = K + 11(m,α ua m (u α+k + 21(m,α ub m (u α +K + 12(m,α uc m (u α+k + 22(m,α ud m (u α { } = tr U (m,α uk + (m,α u. (3.11 The QYBE (2.4, the RE (2.5 and the intertwining relations given in Appendix A allow us to derive the commutation relations among the matrix elements of U (m, α u. Here we present some relevant relations for our purpose: C m (u 1 αc m+2 (u 2 α = C m (u 2 αc m+2 (u 1 α, (3.12 [ Dm 2 (u 2 α, D m 2 (u 1 α ] = sinh(mη +u 1 +u 2 sinhmηsinh(u 1 +u 2 +η C m 2(u 1 αb m (u 2 α sinh(mη +u 1 +u 2 sinhmηsinh(u 1 +u 2 +η C m 2(u 2 αb m (u 1 α, (3.13 D m 2 (u 2 αc m 2 (u 1 α = sinh(u 1 u 2 +ηsinh(u 1 +u 2 sinh(u 1 +u 2 +ηsinh(u 1 u 2 C m 2(u 1 αd m (u 2 α sinh(mη u 1 +u 2 sinh(u 1 +u 2 sinhmηsinh(u 1 u 2 sinh(u 1 +u 2 +η C m 2(u 2 αd m (u 1 α 8

9 sinh(mη +u 1 +u 2 sinhmηsinh(u 1 +u 2 +η C m 2(u 2 αa m (u 1 α, (3.14 [ Dm (u 2 α, A m (u 1 α ] = sinh(m+1ηsinh(mη u 1 +u 2 sinh(u 1 +u 2 +2η sinh(m+2ηsinh(m 1ηsinh(u 1 u 2 sinh(u 1 +u 2 +η [C m (u 1 αb m+2 (u 2 α C m (u 2 αb m+2 (u 1 α]. (3.15 The proof of the above relations is relegated to Appendix B. Let us introduce the left local vacuum state of the n-th site in the lattice as follows: ω;m,α n = X m+n N 1 (θ n α, n = 1,,N, (3.16 where the row vector X m (u is given by (3.5. Further, we introduce the following global vacuum state α+m = 2 N e N l=1 θ l αnη N sinh(m lη l=1 N ω;m n. (3.17 Theexplicitexpression(3.5oftherowvectorX m (uimpliesthattheaboveleftvacuumstate depends only on α+m. Following the method in [9, 37, 38], after some tedious calculation, we obtain the actions of the gauged operators C m (u α, A m (u α and D m (u α on the state as follows: α+m C m (u α =K 21 (m N,α u sinh(m+2η sinh(m+2 Nη N sinh(u θ j +ηsinh(u+θ j sinh 2 α+m+2, (3.18 η α+m D m (u α =K 22(m N,α u N +K 21 (m N,α u N n=1 sinh(u θ j +ηsinh(u+θ j +η sinh 2 α+m η sinh(u θ j +η α+m+1 B m+1 (u α, (3.19 α+m A m (u α = sinh(2u (m 1η sinh(2u+ηsinh(1 mη { N K sinh(u θ j +ηsinh(u+θ j +η 22(m N,α u sinh 2 α+m η N } + K 21 (m N,α u sinh(u θ j +η α+m+1 B m+1 (u α +F(u. (3.20 9

10 Here we have introduced the gauged K -matrix ( K (l K 11,α u = (l,α u K 12 (l,α u K 21 (l,α u K 22 (l,α u ( Y = l (u αk (u X l 2( u α Y l (u αk (uŷl ( u α X l (u αk (u X,(3.21 l ( u α X l (u αk (uŷl +2( u α with l = m N, and the gauged operator B m (u α is given by B m (u α = Y m N+1 ( u αˆt(uŷm+1( u α. (3.22 The extra term F(u in (3.20 actually vanishes at the points { θ j j = 1,,N}, namely, This fact gives rise to the following important relations F( θ j = 0, j = 1,...,N. (3.23 α+m A m ( θ j α = sinh((m 1η +2θ j sinh(m 1ηsinh(2θ j η α+m D m( θ j α. (3.24 The associated right vacuum state, which only depends on α+m, is given by [9] α+m = N X m+n n+1 (θ n α, (3.25 n=1 and the associated gauged operators are ( Am (u α B U (m,α u = m (u α, C m (u α D m (u α ( Ỹm 2 (u αu (ux m ( u α Ỹ m (u αu(uy m ( u α = X m (u αu (ux m ( u α Xm+2 (u αu (uy m ( u α. (3.26 The matrix elements of the above gauged monodromy matrix acting on the right vacuum state (3.25 were given in [9]. Here we present some relevant ones C m (u α α+m = K21 1η (l,α usinh(m+n sinh(m 1η N A m (u α α+m = K 11(l,α u sinh(u θ j sinh(u+θ j +η sinh 2 α+m 2, (3.27 η +K 21(l,α u N sinh(u θ j +ηsinh(u+θ j +η sinh 2 α+m η N sinh(u+θ j +η B m 1 (u α α+m 1, (

11 with l = m+n. Here another gauged K -matrix is (c.f., (3.21 ( K K (l,α u = 11 (l,α u K12 (l,α u K21 (l,α u K 22 (l,α u, ( Ỹl 2 (u αk (ux l ( u α Ỹ l (u αk (uy l ( u α = X l (u αk (ux l ( u α Xl+2 (u αk (uy l ( u α, (3.29 and the gauged operator B m (u α is given by 4 Bethe states B m (u α = Ỹm 1(uT(uY m+n 1 (u. (3.30 Up to now, the parameters α and m in the definitions of the gauged operator U (m,α u in (3.9 and the associated K-matrix K + (m,α u in (3.10 (resp. U (m,α u in (3.26 and the associated K-matrix K (m,α u in (3.29 are arbitrary. The works in [25, 26] shed light on the two important facts to construct the Bethe-type eigenstates of the U(1-symmetrybroken integrable models: (1 The inhomogeneous T Q relation has played a central role in constructing the Bethe states because it enables one in this case to tell the wanted term from the unwanted ones within the framework of the algebraic Bethe Ansatz method; (2 It also suggests that in order to construct the right Bethe states 4 of the transfer matrix (2.10, one may choose the two parameters α and m according to the boundary parameters α +, β + and θ + to construct the creation operator (resp. according to the boundary parameters α, β and θ to seek the associated reference state. For this purpose, let us choose the gauge parameters in (3.10 as follows { def αη = α (l η = η θ + +i π mod(2iπ, 2 (4.1 mη def = m (l η = α + +β + i π mod(2iπ. 2 In this particular choice of the gauged parameters, the corresponding gauged K-matrix K + (m,α u given by (3.10 becomes diagonal K + (m (l,α (l u = Diag(K + 11 (m(l,α (l u,k + 22 (m(l,α (l u, (4.2 where the non-vanishing matrix elements read K + 11(m (l,α (l 2e u u= cosh(α + +β + sinh(u+α ++ηcosh(u+β + +ηcosh(α + +β + η, (4.3 4 The generalization to construct the left Bethe states is straightforward. 11

12 2e u K + 22(m (l,α (l u= cosh(α + +β + sinh(u α ++ηcosh(u β + +ηcosh(α + +β + +η. (4.4 In this case the the transfer matrix (2.10 becomes t(u = K + 11 (m(l,α (l ua m (l(u α (l +K + 22 (m(l,α (l ud m (l(u α (l. (4.5 Direct calculation shows that the following identity holds K + 22(m (l,α (l u+ sinh((m(l 1η 2u sinh(2u+ηsinh(m (l 1η K+ 11(m (l,α (l u = 2e usinh(2u+2η sinh(2u+η sinh(u α +cosh(u β +. (4.6 Then let us choose the gauge parameters in (3.29 such that the following relation is satisfied (m (r +α (r η = θ +α +β Nη +iπ mod(2iπ. (4.7 InthiscasethecorrespondinggaugedK-matrixK (m (r +N,α (r ugivenby(3.29becomes up-triangular with the matrix element K 11 (m(r +N,α (r u fixed, namely, K 21 (m(r +N,α (r u = 0, K 11 (m(r +N,α (r u = 2e u sinh(u α cosh(u β. (4.8 Although neither the parameter α (r nor m (r is fixed by the up-triangularity condition of K (m (r,α (r u, the sum of the two parameters is unique as shown in (4.7. This allows us to define an unique reference state Ω, where the vacuum state α (r +m (r is defined by (3.25. Ω = α (r +m (r, (4.9 Following the method developed in [26], we propose that the Bethe-type eigenstate of the transfer matrix (2.10 for the present model is given by λ 1,,λ N = C m (l(λ 1 α (l C m (l +2(λ 2 α (l C m (l +2(N 1(λ N α (l Ω, (4.10 where the two parameters α (l and m (l are given by (4.1. The Bethe state λ 1,,λ N becomes an eigenstate of the transfer matrix t(u with an eigenvalue Λ(u given by (2.25 provided that the N parameters {λ j j = 1,,N} satisfy the BAEs (2.28. The proof is relegated to Appendix C. 12

13 From the definitions (3.1-(3.8 of the gauge matrices, it is clear that both the reference state Ω and the creation operator C m (l(u α (l have well-defined homogeneous limits: {θ j 0}. This implies that the homogeneous limit of the Bethe state (4.10 exactly gives rise to the corresponding Bethe state of the homogeneous XXZ spin- 1 chain with arbitrary 2 boundary fields, where the associated T Q relation and BAEs are given by (2.25 and (2.28 with {θ j = 0}. It would be interesting to study the relation between our Bethe states and the eigenstates proposed in [27] for which the homogeneous limit is still unclear. 5 Conclusions Inthispaper, wederived thebethestatesofthetransfer matrixofthexxzspin- 1 2 chainwith arbitrary boundary fields specified by the most generic non-diagonal K-matrices given by (2.2 and(2.3. It should be emphasized that constructing the Bethe state of U(1-symmetrybroken models had challenged for many years because of the lacking of the inhomogeneous T Q relations such as (2.25. The idea of this paper to construct the Bethe state is to search for two gauge transformations such that one makes the resulting K + -matrix to be diagonal and the other makes the resulting K -matrix up-triangular. Then we find that the two parameters m (l and α (l of the first gauge transformation must obey the following equations { sinh(α+ +β + = sinh(θ + +(α (l 1η +m (l η, sinh(α + +β + = sinh(θ + +(α (l 1η m (l η, while the parameters of the second gauge transformation have to satisfy the relation (5.1 sinh(α +β +sinh(θ +(m (r +α (r η+nη = 0. (5.2 The equation (5.1 is to determine the creation operator C m (l(u α (l, while the equation (5.2 is to choose the associated reference state (such as (4.9. It is found that besides the solution given by (4.1 and (4.7 there exist the other three solutions of (5.1 and (5.2. Each of the three solutions gives rise to a set of Bethe states with eigenvalues parameterized by a T Q relation of the form(2.25 by replacing α ±, β ± with ±α ±, ±β ±. Nevertheless, different types of inhomogeneous T Q relations [19, 14] only give different parameterizations of the eigenvalues of the transfer matrix but not new solutions. We note that for the degenerate caseconsidered in[9], thepresent methodmaynot workbutthebethestatescanbeobtained via algebraic Bethe Ansatz. 13

14 Acknowledgments The financial supports from the National Natural Science Foundation of China (Grant Nos , , , , the National Program for Basic Research of MOST (973 project under grant No.2011CB921700, BCMIIS and the Strategic Priority Research Program of the Chinese Academy of Sciences are gratefully acknowledged. Two of the authors (W.-L. Yang and K. Shi would like to thank IoP/CAS for the hospitality. Appendix A: Intertwining relations We list some intertwining relations (or face-vertex correspondence relations in [9] which are useful to construct the reference state and the commutation relations among the gauged operators: 5 R 12 (u 1 u 2 Xm+2 1 (u 1Xm+1 2 (u 2 = sinh(u 1 u 2 +η Xm+2 2 (u 2Xm+1 1 (u 1, R 12 (u 1 u 2 Xm 1 (u 1Ym 1 2 (u 2 = sinh(u 1 u 2 sinh(m 1η Ym 2 sinhmη (u 2Xm+1 1 (u 1 + sinh(mη +u 1 u 2 Xm 2 sinh mη (u 2Ym 1 1 (u 1, R 12 (u 1 u 2 Ym(u 1 1 Xm+1(u 2 2 = sinh(u 1 u 2 sinh(m+1η X 2 sinhmη m(u 2 Ym 1(u sinh(mη u 1 +u 2 Y 2 sinh mη m(u 2 Xm+1(u 1 1, R 12 (u 1 u 2 Ym 2 1 (u 1Ym 1 2 (u 2 = sinh(u 1 u 2 +η Ym 2 2 (u 2Ym 1 1 (u 1, R 12 (u 1 u 2 X m 1 2 (u 2 X m 1 (u 1 = sinh(u 1 u 2 +η X m 2 (u 2 X m 1 1 (u 1, R 12 (u 1 u 2 X 2 m 1 (u 2Ŷ 1 m+2 (u 1 = sinh(u 1 u 2 sinh(m+1η sinhmη X 2 m 2 (u 2Ŷ 1 m+1 (u 1 (A.1 (A.2 (A.3 (A.4 (A.5 + sinh(mη u 1 +u 2 Ŷm+2 2 sinh mη (u 2 X m 1 1 (u 1, (A.6 R 12 (u 1 u 2 Ŷ m+1(u 2 2 X m 2(u 1 1 = sinh(u 1 u 2 sinh(m 1η Ŷ 2 sinhmη m+2(u 2 X m 1(u sinh(mη +u 1 u 2 X 2 sinh mη m 2(u 2 Ŷ m+1(u 1 1, (A.7 5 In fact these vectors depend also on α but as this parameter will not vary in the following relations, in this appendix we omit this argument for simplicity temporarily. 14

15 R 12 (u 1 u 2 Ŷ m+1 2 (u 2Ŷ m 1 (u 1 = sinh(u 1 u 2 +η Ŷm 2 (u 2Ŷ m+1 1 (u 1, (A.8 X 1 m 1 (u 1X 2 m 2 (u 2R 12 (u 1 u 2 = sinh(u 1 u 2 +η X 2 m 1 (u 2X 1 m 2 (u 1, (A.9 X 1 m 1(u 1 Y 2 m(u 2 R 12 (u 1 u 2 = sinh(u 1 u 2 sinh(m+1η Y 2 sinhmη m+1(u 2 X 1 m(u 1 + sinh(mη +u 1 u 2 X 2 sinh mη m 1(u 2 Y 1 m(u 1, (A.10 Y 1 m+1 (u 1X 2 m (u 2R 12 (u 1 u 2 = sinh(u 1 u 2 sinh(m 1η X 2 m 1 sinhmη (u 2Y 1 m (u 1 + sinh(mη u 1 +u 2 Y 2 m+1 sinh mη (u 2X 1 m (u 1, (A.11 Y 1 m+1 (u 1Y 1 m+2 (u 2R 12 (u 1 u 2 = sinh(u 1 u 2 +η Y 2 m+1 (u 2Y 1 m+2 (u 1, X 1 m+1(u 1 X 2 m(u 2 R 12 (u 1 u 2 = sinh(u 1 u 2 +η X 2 m+1(u 2 X 1 m(u 1, (A.12 (A.13 X m+1 1 (u 1Ỹ m 2 2 (u 2R 12 (u 1 u 2 = sinh(u 1 u 2 sinh(m+1η Ỹm 1 2 sinhmη (u 2 X m+2 1 (u 1 + sinh(mη +u 1 u 2 X m+1 2 sinh mη (u 2Ỹ m 2 1 (u 1, (A.14 Ỹ 1 m 1 (u 1 X 2 m+2 (u 2R 12 (u 1 u 2 = sinh(u 1 u 2 sinh(m 1η sinhmη X 2 m+1 (u 2Ỹ 1 m 2 (u 1 + sinh(mη u 1 +u 2 Ỹm 1 2 sinh mη (u 2 X m+2 1 (u 1, (A.15 Ỹm 1(u 1 1 Ỹ m(u 2 2 R 12 (u 1 u 2 = sinh(u 1 u 2 +η Ỹ 2 m 1(u 2 Ỹ m(u 1 1, X 2 m (u 2R 12 (u 1 u 2 Xm 1 (u 1 = sinh(u 1 u 2 sinh(m 1η X 2 m 1 sinhmη (u 2Xm+1 1 (u 1, X 2 m (u 2R 12 (u 1 u 2 Ym 1 (u 1 = sinh(u 1 u 2 +η X 2 m+1 (u 2Ym+1 1 (u 1 (A.16 (A.17 + sinh(mη u 1 +u 2 Y 2 m+1 sinh mη (u 2Xm+1 1 (u 1, (A.18 Y 2 m(u 2 R 12 (u 1 u 2 Xm(u 1 1 = sinh(u 1 u 2 +η Y 2 m 1(u 2 Xm 1(u sinh(mη +u 1 u 2 X 2 sinh mη m 1(u 2 Ym 1(u 1 1, (A.19 Y 2 m (u 2R 12 (u 1 u 2 Ym 1 (u 1 = sinh(u 1 u 2 sinh(m+1η Y 2 m+1 sinhmη (u 2Ym 1 1 (u 1, 15 (A.20

16 X m+1 1 (u 1R 12 (u 1 u 2 Xm+1 2 (u 2 = sinh(u 1 u 2 sinh(m+1η Xm 2 sinhmη (u 2 X m+2 1 (u 1, X m+1 1 (u 1R 12 (u 1 u 2 Ym 1 2 (u 2 = sinh(u 1 u 2 +η Ym 2 2 (u 2 X m 1 (u 1 + sinh(mη +u 1 u 2 Xm 2 sinh mη (u 2Ỹ m 2 1 (u 1, Ỹm 1(u 1 1 R 12 (u 1 u 2 Xm+1(u 2 2 = sinh(u 1 u 2 +η X 2 m+2(u 2 Ỹ m(u sinh(mη u 1 +u 2 Y 2 sinh mη m(u 2 X m+2(u 1 1, Ỹm 1 1 (u 1R 12 (u 1 u 2 Ym 1 2 (u 2 = sinh(u 1 u 2 sinh(m 1η Ym 2 sinhmη (u 2Ỹ m 2 1 (u 1, X 1 m 1 (u 1R 12 (u 1 u 2 X 2 m 1 (u 2 = sinh(u 1 u 2 sinh(m+1η sinhmη X 2 m 2 (u 2X 1 m (u 1, (A.21 (A.22 (A.23 (A.24 (A.25 X 1 m 1(u 1 R 12 (u 1 u 2 Ŷ m+1(u 2 2 = sinh(u 1 u 2 +η Ŷ 2 m(u 2 X 1 m 2(u 1 + sinh(mη +u 1 u 2 X 2 sinh mη m 2(u 2 Y 1 m(u 1, (A.26 Y 1 m+1 (u 1R 12 (u 1 u 2 X 2 m 1 (u 2 = sinh(u 1 u 2 +η X 2 m (u 2Y 1 m+2 (u 1 + sinh(mη u 1 +u 2 Ŷm+2 2 sinh mη (u 2X 1 m (u 1, (A.27 Y 1 m+1 (u 1R 12 (u 1 u 2 Ŷ m+1 2 (u 2= sinh(u 1 u 2 sinh(m 1η Ŷm+2 2 sinhmη (u 2Y 1 m (u 1, (A.28 where Xm 1(u, X2 m (u are embedding vector in the 1-st and 2-nd tensor space, respectively. Moreover, the vectors also enjoy the following orthonormal relations: Y m (ux m (u = 1, Y m (uy m (u = 0, X m (ux m (u = 0, X m (uy m (u = 1, ( 1 0 X m (uy m (u+y m (ux m (u =, (A Ỹ m 1 (ux m+1 (u = 1, Ỹ m 1 (uy m 1 (u = 0, X m+1 (ux m+1 (u = 0, Xm+1 (uy m 1 (u = 1, ( 1 0 X m+1 (uỹm 1(u+Y m 1 (u X m+1 (u =, (A Y m+1 (u X m 1 (u = 1, Y m+1 (uŷm+1(u = 0, X m 1 (u X m 1 (u = 0, X m 1 (uŷm+1(u = 1, 16

17 X m 1 (uy m+1 (u+ŷm+1(ux m 1 (u = ( Appendix B: Commutation relations Using QYBE (2.4 and the RE (2.5, one may derive that. (A.31 R 12 (u 1 u 2 U 1 (u 1 R 21 (u 1 +u 2 U 2 (u 2 = U 2 (u 2 R 21 (u 1 +u 2 U 1 (u 1 R 12 (u 1 u 2. (B.1 MultiplyingtheaboveequationwithX 1 m+1 (u 1X 2 m (u 2fromtheleftand X 1 m+1 ( u 1 X 2 m+2 ( u 2 from the right, and using the relations (A.5 and (A.9, we arrive at (3.12. Similarly, multiplying(b.1withx 1 m 1 (u 1X 2 m 2 (u 2(X 1 m 1 (u 1X 2 m 2 (u 2fromtheleftandŶ1 m+1 ( u 1Ŷ 2 m ( u 2 ( X 1 m 1 ( u 1Ŷ 2 m+2 ( u 2fromtherightandusing theintertwining relations(a.1-(a.28, one can obtain the relation (3.13 (or (3.14. Using the similar method and the relation (3.13, one can further check (3.15. Appendix C: Proof the Bethe state There are several ways [25, 26, 31, 32] to show that the state λ 1,,λ N constructed by (4.10 is the eigenstate of the transfer matrix (2.10. Here we adopt the method developed in [26] to demonstrate it. For arbitrary parameters α, m let us introduce the following left states 6 parameterized by the N inhomogeneous parameters {θ j }: α,m;θ p1 θ pn = α+m D m ( θ p1 α D m ( θ pn α, 1 q 1 < q 2 <... < q n N, n = 0,1,,N. (C.1 The commutation relations (3.13, (3.14 and (3.18 imply that α,m;θ p1,,θ pn C m (u α=g(u,{θ p1,,θ pn } α,m+2;θ p1,,θ pn, (C.2 where g(u,{θ p1,,θ pn }= g 0 (u m,α n sinh(u+θ pj +ηsinh(u θ pj sinh(u θ pj +ηsinh(u+θ pj, (C.3 6 Such states were used as a basis to construct the SoV eigenstates of the XXZ open chain [27]. Here we use two different gauge transformations respectively for the left and right pseudo vacuum states to reach the Bethe states. 17

18 and The above equations lead to that g 0 (u m,α = K 21 (m N;α u sinh(m+2η sinh(m+2 Nη N sinh(u θ j +ηsinh(u+θ j sinh 2. (C.4 η α,m;θ p1,,θ pn C m ( θ pj α = 0, j 1,,n. (C.5 Let Ψ be an eigenstate of the transfer matrix t(u with an eigenvalue Λ(u, namely, t(u Ψ = Λ(u Ψ. (C.6 Following the method used in [18, 26], we introduce the following scalar products F n (θ p1,,θ pn = α (l,m (l ;θ p1,,θ pn Ψ, n = 0,,N, (C.7 which altogether uniquely determine the eigenstate Ψ. With the particular choice of the parameters (4.1, one can derive the following recursive relations (see (C.8 below by considering the quantity of α (l,m (l ;θ p1,,θ pn t( θ pn+1 Ψ, Λ( θ pn+1 F n (θ p1,,θ pn = K + 11(m (l,α (l θ pn+1 α (l,m (l ;θ p1,,θ pn A m (l( θ pn+1 α (l Ψ +K + 22(m (l,α (l θ pn+1 F n+1 (θ p1,,θ pn,θ pn+1. The relations (3.14, (3.15, (3.24 and (C.5 enable us to further simplify the above equation { Λ( θ pn+1 F n (θ p1,,θ pn = F n+1 (θ p1,,θ pn,θ pn+1 K + 22 (m(l,α (l θ pn+1 } sinh((m(l 1η+2θ pn+1 sinh(m (l 1ηsinh(2θ pn+1 η K+ 11 (m(l,α (l θ pn+1 (4.6 sinh( 2θ = 2e θp pn+1 +2η n+1 sinh( 2θ pn+1 +η sinh( θ p n+1 α + cosh( θ pn+1 β + F n+1 (θ p1,,θ pn,θ pn+1. Then the inhomogeneous T Q relation (2.25 implies that {F n (θ p1,,θ pn } satisfy the following recursive relations F n+1 (θ p1,,θ pn,θ pn+1 = 2e θp n+1 sinh(θpn+1 +α cosh(θ pn+1 +β Q( θ p n+1 η Q( θ pn+1 Ā( θ p n+1 F n (θ p1,,θ pn. 18 (C.8

19 Iterating the above recursive relation, we readily obtain n { F n (θ p1,,θ pn = 2e θp j sinh(θpj +α cosh(θ pj +β Q( θ } p j η Ā( θ pj F 0, Q( θ pj n = 0,1,,N, (C.9 where F 0 = α (l,m (l Ψ is an overall scalar factor and the function Ā(u is given by (2.18. The relation (C.2 implies that n { } α (l,m (l Q( θpj η ;θ p1,,θ pn λ 1,,λ N = G 0 α (l,m ;θ p1,,θ pn Ω Q( θ pj n = 0,,N, (C.10 where m = m (l +2N and G 0 is an overall factor which does not depend on n, N G 0 = g 0 (λ j m (l +2(j 1,α (l. (C.11 Hence it is sufficient to compute the quantity { α (l,m ;θ p1,,θ pn Ω }. Due to the fact that the particular choice (4.7 of the parameters m (r, α (r makes the matrix element K 21 (m(r + N,α (r u vanishes (see (4.8, we can derive the following relations from (3.27 and (3.28 C m (r(u α (r Ω =0, A m (r(u α (r Ω = K11 (m(r (r +N,α uā(u Ω. (C.13 (C.12 The definitions (3.9 and (3.26 of the two gauged double-row monodromy matrices and the relations(a.29-(a.31allowustoexpress theoperators C m (u α (l and D m (u α (l interms of some linear combinations of A m (r(u α (r, B m (r(u α (r, C m (r(u α (r and D m (r(u α (r respectively, namely, C m ( u α (l =X m ( u α (l X m (r( u α (r A m (r( u α (r Y m (r(u α (r X m (u α (l +X m ( u α (l Y m 2( u α (r C (r m (r( u α (r Y m (r(u α (r X m (u α (l +X m ( u α (l X m +2( u α (r B (r m (r( u α (r X m (r(u α (r X m (u α (l +X m ( u α (l Y m (r( u α (r D m (r( u α (r X m (r(u α (r X m (u α (l, (C.14 D m ( u α (l =X m ( u α (l X m (r( u α (r A m (r( u α (r Y m (r(u α (r Ŷm +2(u α (l +X m ( u α (l Y m 2( u α (r C (r m (r( u α (r Y m (r(u α (r Ŷm +2(u α (l +X m ( u α (l X m +2( u α (r B (r m (r( u α (r X m (r(u α (r Ŷm +2(u α (l +X m ( u α (l Y m (r( u α (r D m (r( u α (r X m (r(u α (r Ŷm +2(u α (l. (C.15 19

20 The vanishing condition (C.5 implies that α (l,m ;θ p1,,θ pn C m ( θ pn+1 α (l Ω = 0, n = 0,1,,N 1. (C.16 Keeping the relations (C.12 and (C.13 in mind and using the above equations and the explicit expressions (3.1, (3.5-(3.8, after a tedious calculation, we can derive the following recursive relations α (l,m ;θ p1,,θ pn+1 Ω = K 11 (m(r +N,α (r θ pn+1 Ā( θ p n+1 α (l,m ;θ p1,,θ pn Ω (4.8 = 2e θp n+1 sinh(θpn+1 +α cosh(θ pn+1 +β Ā( θ p n+1 α (l,m ;θ p1,,θ pn Ω, n = 0,1,,N 1. (C.17 Iterating the above recursive relations, we have n { α (l,m ;θ p1,,θ pn Ω = 2e θ pj sinh(θpj +α cosh(θ pj +β Ā( θ p j } α (l +m Ω, n = 0,1,,N. Submitting the above relations into (C.10, finally we obtain the expressions of the scalar products of the Bethe state (4.10 with each state (C.1 α (l,m (l ;θ p1,,θ pn λ 1,,λ N n { = 2e θp j sinh(θpj +α cosh(θ pj +β Q( θ } p j η Ā( θ pj Q( θ pj G 0 α (l +m Ω, n = 0,1,,N, (C.18 where the scalar factor G 0, which does not depend on n, is given by (C.11. Comparing the above expression with (C.9, we conclude that the Bethe state λ 1,...,λ N given by (4.10 is an eigenstate of the transfer matrix t(u with an eigenvalue Λ(u given by (2.25, provided that the parameters {λ j } satisfy the BAEs (2.28. References [1] R. J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic Press,

21 [2] V. E. Korepin, N. M. Bogoliubov and A. G. Izergin, Quantum Inverse Scattering Method and Correlation Function, Cambridge University Press, Cambridge, [3] N. Beisert et al., Lett. Math. Phys. 99 (2012, 3. [4] T. Giamarchi, Quantum Physics in One Dimension, Oxford University Press, Oxford, [5] J. de Gier and F.H.L. Essler, Phys. Rev. Lett. 95 (2005, [6] J. Sirker, R.G. Pereira and I. Affleck, Phys. Rev. Lett. 103 (2009, [7] E.K. Sklyanin, J. Phys. A 21 (1988, [8] R. I. Nepomechie, J. Phys. A 34 (2001, 9993; R. I. Nepomechie, Nucl. Phys. B 622 (2002, 615; R. I. Nepomechie, J. Stat. Phys. 111 (2003, 1363; R. I. Nepomechie, J. Phys. A 37 (2004, 433. [9] J. Cao, H.-Q. Lin, K.-J. Shi and Y. Wang, Nucl. Phys. B 663 (2003, 487. [10] W.-L. Yang, Y.-Z. Zhang and M. Gould, Nucl. Phys. B 698 (2004, 503. [11] J. de Gier and P. Pyatov, J. Stat. Mech. (2004, P03002; A. Nichols, V. Rittenberg and J. de Gier, J. Stat. Mech. (2005, P03003; J. de Gier, A. Nichols, P. Pyatov and V. Rittenberg, Nucl. Phys. B 729 (2005, 387. [12] A. Doikou and P. P. Martins, J. Stat. Mech. (2006, P06004; A. Doikou, J. Stat. Mech. (2006, P [13] Z. Bajnok, J. Stat. Mech. (2006, P [14] W.-L. Yang, R.I. Nepomechie and Y.-Z. Zhang, Phys. Lett. B 633 (2006, 664. [15] W. Galleas, Nucl. Phys. B 790 (2008, 524. [16] G. Niccoli, J. Stat. Mech. (2012, P [17] S. Belliard, N. Crampé and E. Ragoucy, Lett. Math. Phys. 103 (2013, 493. [18] J. Cao, W.-L. Yang, K.-J. Shi and Y. Wang, Phys. Rev. Lett. 111 (2013,

22 [19] J. Cao, W.-L. Yang, K.-J. Shi and Y. Wang, Nucl. Phys. B 875 (2013, 152. [20] J. Cao, W.-L. Yang, K.-J. Shi and Y. Wang, Nucl. Phys. B 877 (2013, 152. [21] J. Cao, S. Cui, W.-L. Yang, K.-J. Shi and Y. Wang, Nucl. Phys. B 886 (2014, 185. [22] Y. -Y. Li, J. Cao, W. -L. Yang, K. Shi and Y. Wang, Nucl. Phys. B 884 (2014, 17. [23] J. Cao, W.-L. Yang, K.-J. Shi and Y. Wang, JHEP 04 (2014, 143. [24] K. Hao, J. Cao, G. -L. Li, W. -L. Yang, K. Shi and Y. Wang, JHEP 06 (2014, 128. [25] S. Belliard and N. Crampé, SIGMA 9 (2013, 072. [26] X. Zhang, Y.-Y. Li, J. Cao, W. -L. Yang, K. Shi and Y. Wang, arxiv: v2. [27] S. Faldella, N. Kitanine and G. Niccoli, J. Stat. Mech. (2014, P [28] E. K. Sklyanin, Lect. Notes Phys. 226 (1985, 196; J. Sov. Math. 31 (1985, 3417; Prog. Theor. Phys. Suppl. 118 (1995, 35. [29] E. K. Sklyanin and L. D. Faddeev, Sov. Phys. Dokl. 23 (1978, 902. [30] L. A. Takhtadzhan and L. D. Faddeev, Rush. Math. Surveys 34 (1979, 11. [31] S. Belliard, arxiv: [32] N. Crampé, arxiv: [33] H.J. de Vega and A. González-Ruiz, J. Phys. A 26 (1993, L519. [34] S. Ghoshal and A.B. Zamolodchikov, Int. J. Mod. Phys. A 9 (1994, [35] R. I. Nepomechie, J. Phys. A 46 (2013, [36] J. Cao, W.-L. Yang, K. Shi and Y. Wang, arxiv: [37] H. Fan, B.-Y. Hou, K.-J. Shi and Z.-X. Yang, Nucl. Phys. B 478 (1996, 723. [38] W.-L. Yang and R. Sasaki, J. Math. Phys. 45 ( ; Nucl. Phys. B 679 (