Exact solution for the Bariev model with boundary fields arxiv:cond-mat/ v1 [cond-mat.str-el] 4 Sep 2000

Transcription

1 Exact solution for the Bariev model with boundary fields arxiv:cond-mat/ v1 [cond-mat.str-el] 4 Sep 000 A. Foerster a, X.-W. Guan a, J. Links c, I. Roditi b, H.-Q Zhou c, Abstract a Instituto de Fisica da UFRGS, Av. Bento Goncalves, 9500, Porto Alegre, , Brasil b Centro Brasileiro de Pesquisas Fisicas, Rua Dr. Xavier Sigaud 150, , Rio de Janeiro- RJ, Brasil c Centre for Mathematical Physics, Department of Mathematics, The University of Queensland, 407, Australia The Bariev model with open boundary conditions is introduced and analysed in detail in the framework of the Quantum Inverse Scattering Method. Two classes of independent boundary reflecting K-matrices leading to four different types of boundary fields are obtained by solving the reflection equations. The models are exactly solved by means of the algebraic nested Bethe ansatz method and the four sets of Bethe ansatz equations as well as their corresponding energy expressions are derived. Key words: Integrable spin chains; Algebraic Bethe ansatz; Yang-Baxter algebra; Reflection equations; PACS: w; Fd; Jm 1. Introduction Much work has been devoted, during the last years, towards a better understanding of integrable models of strongly correlated electrons. Such an interest arises from the efforts to unveil a theory that can provide a consistent picture for observed phenomena of condensed matter physics such as high-t c superconductivity. Two prototypes of integrable models for this purpose are the one-dimensional (1D Hubbard model solved by Lieb and Wu [1] (see also [] and the t-j model [3] (which is integrable for the supersymmetric coupling [4], mainly because of their relevance in the description of electronic mechanisms. This is justified on physical grounds since the electron hopping is strongly disturbed by the on-site Coulomb interaction in the 1D Hubbard model and Preprint submitted to Elsevier Preprint 1 December 017

2 by the spin fluctuations through the antiferromagnetic coupling for the supersymmetric t-j model. More recently, motivated by the inclusion of additional interactions whether through internal impurities or external boundary fields, a renewed interest in the study of integrable strongly correlated electron systems has taken place. Some examples of the latter are the 1D Hubbard open chain [5], the supersymmetric t-j model with boundary fields [6,7] and other integrable electronic systems with boundary interactions associated with Lie superalgebra symmetries [8]. The study of integrable systems with open boundary conditions (OBC [5,9 11] has turned into an active domain since the introduction by Sklyanin [1] of a generalization of the Quantum Inverse Scattering Method (QISM [13,14] which provides a systematic approach to handle the boundary problem by introducing a set of equations, called reflection equations (RE. The solutions of these equations, referred to as the boundary K-matrices, in turn introduce boundary interactions into the Hamiltonian of the system in such a manner that integrability is preserved. Indeed, the presence of boundary fields drastically changes surface critical properties, and in some occasions results in the appearance of so-called boundary bound states [15 17]. Nevertheless, in spite of these results, several models still lack a comprehensive study regarding OBC, one of those being the 1D Bariev model [18,19]. The 1D Bariev (interacting XY chain is also a Hubbard-like integrable model of special interest, as it exhibits the existence of hole pairs of Cooper type which are relevant to theories of superconductors. The integrability of the 1D Bariev model within the framework of the QISM was first established by Zhou in [0] (see also [1]. To achieve this goal one has to build non-trivial higher conserved currents that allow the identification of a quantum R-matrix satisfying the Yang-Baxter equation [13,14]. This R-matrix does not possess the difference property for the spectral parameter, a feature that is also shared by the R-matrix associated with the Hubbard model []. Soon after this parametrization of the model had been obtained, the algebraic Bethe ansatz solution for the periodic chain was calculated independently by Martins and Ramos [3] and by Zhou [4]. On the other hand, as far as we know, studies concerning OBC for the Bariev model appear incomplete. This can be understood in light of the fact that the R-matrix is presented in a complicated form with the absence of both additivity for the spectral parameter and crossing unitarity. As such, we must resort to searching for solutions of a much more intricate RE. Though some known solutions to the RE for the 1D Bariev open chain provide a certain class of boundary K-matrices [5], a thorough investigation is still lacking and, moreover, the algebraic Bethe ansatz solution has not been obtained yet. In this paper, we shall study in detail the OBC for the 1D Bariev chain within the framework of the QISM. We will show that besides the K-matrices

3 obtained in [5] there exists a new class of solutions of the RE which leads to pure magnetic boundary fields in the Hamiltonian and provides more general nesting boundary K-matrices for a hidden asymmetric XXZ open chain. Here it is worth noticing that this class of boundary fields may have a feasible realization by applying boundary external voltages in experiments on quantum wires [6]. By means of the algebraic Bethe ansatz technique, we obtain the eigenvalues of the transfer matrix of the 1D Bariev model for four different kinds of possible boundary fields. For all cases we present the Bethe ansatz equations and the energy spectrum in explicit form, which provides a first step towards investigating the surface critical properties for the 1D Bariev model with these different boundary fields. The paper is organized as follows. In section we construct the Bariev model with four kinds of boundary fields by means of the QISM adapted to special boundaryconditions.thebasicquantities, e.g.,ther-matrix,thematricesk ± defining the boundaries, the monodromy matrices and the transfer matrices are defined. In particular, we solve the RE and obtain two independent classes of solutions for the boundary K-matrices using the variable-separation method [7]. The relation between the transfer matrices and the Bariev Hamiltonians with different boundary fields is established. Sections 3 and 4 are devoted to the solution of the models through Bethe ansatz methods. In particular, in section 3 we discuss the action of the transfer matrix on the pseudovacuum state. In section 4 we perform the nested algebraic Bethe ansatz and find the Bethe ansatz equations, the eigenvalues of the transfer matrices and the energy spectrum of the system. Section 5 is reserved for our conclusions and discussions.. The Bariev open chain with boundary fields Let us begin by introducing the Bariev model with boundary fields whose Hamiltonian reads H= N 1 j=1 {( σ + j σ j+1 +σ j σ+ j+1 exp ( βτ + j+1 τ j+1 + ( τ + j τ j+1 +τ j τ+ j+1 exp ( βσ + j σ j where the four different boundary terms are given by } +b.t. (.1 b.t. (I,I = exp(β c [sinhβσ z Nτ z N +coshβ (σ z N +τ z N] 3

4 exp(β [sinhβσ1 z c τz 1 +coshβ (σz 1 +τz 1 ], (. 1 b.t. (II,II = 1 exp(βc σ z N + 1 exp(βc 1τ z 1, (.3 b.t. (I,II = exp(β c [sinhβσ z N τz N +coshβ (σz N +τz N ] + 1 exp(βc 1τ z 1, (.4 b.t. (II,I = 1 exp(βc σ z N exp(β c 1 [sinhβσ z 1τ z 1 +coshβ (σ z 1 +τ z 1]. (.5 Above σ i and τ i are the two commuting species of Pauli matrices acting on site i, β is a coupling constant and c 1, c are the boundary parameters characterizing the strength of the boundary fields. Notice that in cases (.,.3 the boundary terms at the left and right ends are consistent with the bulk symmetry which is a combination of the inversion j N j + 1 and the exchange σ τ. The cases (.4,.5 are interchanged under this symmetry. In order to derive these models, we first consider the following R-matrix [0] R(u,v = ρ ρ 0 0 ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ( ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ 1 whose elements (Boltzmann weights are given explicitly in Appendix A. It satisfies the non-additive Yang-Baxter equation R 1 (u,vr 13 (u,wr 3 (v,w = R 3 (v,wr 13 (u,wr 1 (u,v (.7 4

5 and enjoys the following properties R 1 (u,vr 1 (v,u=1, (.8 R t 1 1 ( v,ur t 1 (u, v=1, (.9 R t 1 ( u,vr t 1 1 (v, u=1. (.10 Above, t i denotes matrix transposition in the ith space and R 1 = P.R 1.P, where P isthepermutationoperator. Thematrix R isanew object introduced in Ref. [5] in order to handle integrable boundary terms. We may take either (.9 or (.10 as its definition. We remark here that neither the crossing unitarity nor the difference property is valid for this R-matrix, which turns the calculations into a complex problem. The Yang-Baxter algebra can also be verified for the monodromy matrix T R 1 (u,v 1 T (u T (v = T (v 1 T (ur 1 (u,v, (.11 which is defined, as usual, by T(u = L 0N (u L 01 (u. (.1 Above the subscripts 0 and N,...,,1 denote the auxiliary and quantum spaces, respectively and the local Lax operator is given by P j + Q + j f 1 (up j + τj f 1 (uσ Q j + j f 1 (uf (uσj τj f L j (u = 1 (u P j + τ j + P j + Q j f (uσ j τ j + f (uσ Q j j f 1 (uσ j + Q + j f1(uσ j + τj Pj Q + j f (upj τj f 1 (uf (uσ j + τ j + f (uσ j + Q j f (u P j τ j + P j Q (.13 where P ± j = 1 (1±σz j + u exp(β(1 σz j ; P ± j = 1 (1±σz j + u (1 σz j ; Q ± j = 1 (1±τz j+ u exp(β(1 τz j; Q ± j = 1 (1±τz j+ u (1 τz j; f 1 (u= 1+u exp(β, f (u = 1+u 5

6 Next, following the QISM adapted to the case of special boundary conditions, we define the doubled monodromy matrix as U (u = T(uK (ut 1 ( u, (.14 such that the transfer matrix is given by τ(u = Tr 0 K + (uu (u (.15 AboveT 1 istheinverseofthemonodromyt andk ± arethematricesdefining the boundaries. The requirement that they obey the RE [5] R 1 (u,v 1 K (ur 1 (v, u K (v = K 1 (vr 1 (u, v K (ur 1 ( v, u, (.16 R t 1t 1 (v,u K t + (v R 1 ( v,u 1 K t 1 + (u R 1 ( u,v K t + (v = 1 K t 1 + (urt 1t 1 ( u, v (.17 together with the Yang-Baxter algebra assure that the transfer matrix commutes for different spectral parameters, proving the integrability of the model. Therefore the transfer matrix (.15 may be considered as the generating fuction of infinitely many integrals of motion for the system. We emphasize that there is no isomorphism between the matrices K + (u and K (u, due to the absence of the crossing unitarity for the R-matrix. Therefore, we have to solve the two RE separately in order to fix the boundaries. In Appendix A, this calculation is presented in detail using the variable-separation prescription. This approach, which was already proposed in [7] to treat the D ( model with OBC, is extended here for the case of an R-matrix without the difference property. It is in fact a systematic method which simplifies drastically the calculation of the K ± matrices. After a tedious algebra, we find two different classes of boundary K ± matrices consistent with the integrability of the Hamiltonian (.1, listed below K1 ± (u K K ± (u = ± (u K3 ± (u 0, ( K4 ± (u Solutions for K, corresponding to the left boundary Case I K1 (u=(u c 1 (h u c 1, 6

7 Case II K (u=(c 1 u(h u+c 1, K3 (u=(c 1 u(h u+c 1, (.19 K4 (u=(u+c 1 (h u+c 1, K1 (u=(1+c 1 hu, K (u=(1 c 1 hu, K3 (u=(1+c 1 hu, (.0 K4 (u=(1 c 1 hu. Solutions for K +, corresponding to the right boundary Case I Case II K1 + (u=(c u+1(h c u+1, K + (u=h (c u 1(c u+1, K3 + (u=(c u 1(c u+1, (.1 K4 + (u=(c u 1(c u h, K1 + (u=h(hu c, K + (u=h(hu c, K3 + (u=(u+hc, (. K4 + (u=(u+hc where h = exp(β. These two classes of boundary K ± -matrices provide four possible choices of BC, according to the combination of the boundary pairs (K + (u,k (u, which are being labeled by (I,I, (II, II, (I,II and (II,I and originate the boundary terms b.t. (I,I, b.t. (II,II, b.t. (I,II and b.t. (II,I respectively. Taking into account the fact that TrK + (0 = 0 and TrK + (0 = 0 in the cases (I,I and (I,II and TrK + (0 = 0 in the cases (II,II and (II,I, we can show that the Hamiltonian (.1 with the boundary terms b.t. (I,I and b.t. (I,II is related to the transfer matrix matrix (.15 in the following way τ(u = C 1 u +C (H +const. u 3 +, (.3 while τ(u = C 3 u+c 4 (H +const. u + (.4 for the Hamiltonian (.1 with the boundary terms b.t. (II,II and b.t. (II,I. The 7

8 difference in the leading term expansions (.3,.4 is a consequence of the fact that for case I the matrix K + (u is quadratic in the variable u while for case II it is linear. Above C i, i = 1, 4, are some scalar functions of the boundary parameters. In the next sections, we shall be focusing in the solution of the eigenvalue problem of the transfer matrix (.15 and, consequently, in the diagonalization of the present model. 3. Algebraic Bethe ansatz I: Action of the operators on the pseudovacuum In order to carry out the algebraic Bethe ansatz for the Bariev model with boundaries, we first need to find the eigenvalues and eigenvectors of the transfer matrix (.15 τ Φ n = λ Φ n (3.5 This will be achieved by extending the Bethe ansatz techniques developed in [8],[9] for the Hubbard model to the present case. Following this prescription, the eigenvectors of the transfer matrix are given by Φ n = Φ n.f 0. (3.6 where the components of F are coefficients of an arbitrary linear combination of vectors Φ n (see next section for details and 0 is the pseudovacuum state, choosen here as the standard ferromagnetic one 0 = N j=1 0 j, (3.7 where 0 i = (3.8 i i whichcorrespondstothedoublyoccupiedstate.thenitiseasytoworkoutthe action of the monodromy matrices T,T 1 on this state. We begin by writing 8

9 the periodic monodromy matrix T(u as B(u B 1 (u B (u F(u C T(u = 1 (u A 11 (u A 1 (u E 1 (u. (3.9 C (u A 1 (u A (u E (u C 3 (u C 4 (u C 5 (u D(u and T(u = T 1 ( u as T(u = B(u B1 (u B (u F(u C 1 (u Ā11(u Ā1(u Ē1(u. (3.30 C (u Ā1(u Ā(u Ē(u C 3 (u C 4 (u C5 (u D(u By computing the action of T(u and T(u on the psedudovacuum we find B(u 0 = B(u 0 = 0, D(u 0 = D(u 0 = [u] N 0, A aa (u 0 = Āaa(u 0 = [u exp(β] N 0, A 1 (u 0 = Ā1(u 0 = 0, A 1 (u 0 = Ā1(u 0 = 0, B a (u 0 0, Ba (u 0 0, E a (u 0 0, Ē a (u 0 0, F(u 0 0, F(u 0 0, C i (u 0 = C i (u 0 = 0, i = 1,,5,a,b = 1,. (3.31 Similarly, the doubled monodromy matrix U = T(uK (u T(u(.14 can be arranged as U (u = B(u B1 (u B (u F(u C 1 (u Ã11(u Ã1(u Ẽ1(u, (3.3 C (u Ã1(u Ã(u Ẽ(u C 3 (u C 4 (u C5 (u D(u 9

10 Then using eq.(3.31 and the Yang-Baxter algebra T (ur 1 (u, u 1 T (u = 1 T (ur 1 (u, u T (u, (3.33 we can calculate the action of the doubled monodromy U on the pseudovacuum,whichduetotheexistenceoftwoclassesofboundaryk matrices(.19, (.0 yields (see Appendix B B(u 0 = W 1 (ub(u B(u 0, (3.34  aa (u 0 = Wa+1 (ua aa(uāaa(u 0, (3.35 ˆD(u 0 = W4 (ud(u D(u 0, (3.36 B a (u 0 0, Ẽ a (u 0 0, à ab (u 0 = 0, F(u 0, C i (u 0 = 0, i = 1,,5, a b = 1,, (3.37 where for case I: W 1 (u = 1, W (u = W 3 (u = h u(uc 1 1 (h u 1(h u c 1, W 4 (u = u (uc 1 1(h +1(uc 1 h (u 1(u h (h u c 1 (u c 1, and for case II: W 1 (u = 1, W (u = h u (h u 1, W 3 (u = hu(hu+c 1 (h u 1(1+c 1 hu, W 4 (u = u 3 (h +1(u+hc 1 (u 1(u h (1+c 1 hu, Actually, in order to considerably simplify the calculations, the following transformations  aa (u=ãaa(u ρ (u, u ρ 1 (u, u B(u, ˆD(u= D(u ρ 4(u, u ρ 1 (u, u B(u 10

11 1+u h u Â11(u h (1+u h u  (u, (3.38 were used in the above equations (3.34 (3.36. In fact, these types of transformations have also been performed in other models with OBC [5,9 11]. They simplify not only these expressions for the action of the operators on the pseudovacuum, but also simplify the fundamental commutation relations, allowing for an easy recognition of the wanted and unwanted terms, essential to the Bethe ansatz machinery (see next section for details. Then, after this transformation and taking into account the two sets of boundary matrices K + (.1,(., the transfer matrix (.15 can be recast as τ(u = W 1 + (u B(u+ W a+1 + (uâaa(u+w 4 + (uˆd(u, (3.39 a=1 where for case I: W + 1 (u = u (u c (h +1(h u c (uc +1(uh c +1(u 1(h u 1, W + (u = u(uc 1(u c (h +1 (h uc +1(u h (uc +1, W + 3 (u = h u(uc 1(u c (h +1 (h uc +1(u h (uc +1, W + 4 (u = (uc 1(uc h (uc +1(h uc +1, and for case II: W + 1 (u = u(h +1(hc u+1 h(u 1(h u 1(hu c, W + (u = (h +1(u+hc h(u h (hu c W + 3 (u = u(h +1(h+uc (u h (hu c W + 4 (u = (u+hc h(hu c, Notice from eqs. (3.34 (3.37 that the Lax operators (3.3 act as creation (annihilation operators on the pseudovacuum. Before ending this section, we would like to recall that the boundary matrix pairs (K + (u,k (u lead to four kinds of boundary fields preserving the integrability of the model. Consequently, the factors W i ± (k i = 1,...,4 also form a pair (W i + (k,wi (k leading to four possible choices consistent with the boundary matrix pairs. The 11

12 perfect factorization provide us with the nesting boundary K-matices appearing in (3.35 and (3.39 respectively. In the following section, we shall prove that they indeed constitute the solutions of the RE for the hidden asymmetric XXZ spin open chain. 4. Algebraic Bethe ansatz II: The nesting procedure The next step in the solution of the eigenvalue problem of the transfer matrix consists in the computation of the necessary commutation relations between the diagonal fields and the creation fields. From the RE (.16 and the definition (3.3, after tedious calculations, we obtain the main fundamental commutation relations. We display below just those needed for future calculations where B(u B a (v = (u+v(uvh +1 (u v(uvh 1 B a (v B(u+u.t., (4.40 ˆD(u B a (v = (u+v(uvh +1 (u vh (uv 1 B a (vˆd(u+u.t., (4.41 Â ad1 (u B d (v = (u+v(1+uvh (u v(uvh 1 r 1(u, v c 1b a 1 c r 1 ( v, u d 1d b 1 b B c (vâc 1 b 1 (u+u.t., (4.4 B a (u B b (v = { Bc (v B d (u 1+h u + 1+v F(v ξ 1 uv ( I Ã(u}.r( v, u 1+h v 1+u F(u ξ 1 uv ( I Ã(v + (h u+v 1+h u 1+v [ F(v B(u (u h v(uv 1 (h v +u 1+h v ] 1+u (h u+v 1+h u F(u B(v. ξ, ( v ξ = (0,h 1 Ã 11 (u,1,0, Ã(u = Ã1(u, (4.44 Ã 1 (u Ã(u 1

13 c(u, v d(u, v 0 r 1 (u, v =, ( ā(u, v c(u, v a( v, u c( v, u 0 r( v, u =, ( c( v, u d( v, u and ā(u, v = ρ 7(u, vρ 1 (u, v ρ (u, v, ρ 1(u, v ρ (u, v d(u, v = ρ 10(u, vρ 1 (u, v ρ (u, v ρ 1(u, v ρ (u, v c(u, v = ρ 1(u, vρ 15 (u, v ρ 1(u, v ρ (u, v a( v, u = ρ 14( v, uρ 10 ( v, u ρ 13( v, u, ρ 1 ( v, uρ 14 ( v, u d( v, u = ρ 14( v, uρ 7 ( v, u ρ 1( v, u, ρ 1 ( v, uρ 14 ( v, u c( v, u = ρ 14( v, uρ 15 ( v, u ρ 1 ( v, uρ 13 ( v, u. (4.47 ρ 1 ( v, uρ 14 ( v, u It remains to consider the commutation relations for the creation field F. Since they are extremely involved, we have collected them in Appendix B. We remark here that the fundamental commutation relations are much more complicated in comparsion with those obtained in the case of periodic boundary conditions, turning the calculations here much more cumbersome and difficult to treat. In all relations, we have omitted the unwanted terms (those that can not generate an eigenvector of the transfer matrix to save space. It turns out that the auxiliary matrices r 1 (u, v and r( v, u are nothing but the rational R matrices of an isotropic asymmetric six-vertex model. Indeed, by performing the parameterization u = exp(ik 1,v = exp(ik and h = exp(β = exp( iη, one may present these matrices as 13

14 a(k r( v, u = r(k 1 k = 1 k c(k 1 k 0, ( c(k 1 k d(k 1 k r 1 (u, v = P.r(k 1 +k η, (4.49 r 1 (u, v = P.r(k 1 k.p (4.50 with the Boltzmann weights sinη a(k = sin ( exp k +η sinη d(k = sin ( exp k +η c(k = sin k sin ( k +η exp ( i 1 k, (4.51 (i 1 k, (4.5 ( i 1 k. (4.53 In view of the commutation relation (4.43, we may phenomenologically write [8] the n-particle state as Φ n (k 1,...,k n = Φ n (k 1,...,k n F 0. (4.54 with a recursive relation of n-particle vector Φ n (k 1,...,k n = B a1 (k 1 Φ n 1 (k,...,k n n {[ + ξ F(k1 ] Φ n (k,...,k j 1,k j+1,...,k n + j= B(k j G j 1(k (n 1,...,k n } n {[ ξ F(k1 ] Φ n (k,...,k j 1,k j+1,...,k n j= [ I Ã(k j ] H (n j 1(k 1,...,k n }. (4.55 Above F = F a 1,...,a n are the coefficients of arbitrary linear combination of the vectors reflecting the spin degrees of freedom with a i = 1, (i = 1,...,n and ξ plays the role of forbidding two spin up or two spin down particles at the same site. F(u creates a local particle pair with opposite spins. From the commutation relation (4.43, we can also conclude that Φ n (k 1,,k n satisfies the symmetry relation 14

15 Φ n (k 1,,k j,k j+1,,k n = Φ n (k 1,,k j+1,k j,,k n.r(k j k j+1 (4.56 This symmetry giving a restriction to the functions H j 1(k (n 1,,k n and G (n j 1(k 1,,k n is very useful to deduce the coefficients and simplify the unwanted terms in the eigenvalue of the transfer matrix. In fact, after performing three-particle scattering, the explicit form of these coefficients can be fixed. However checking three-particle scattering is indeed a extremely tough calculation and we had to leave the coefficients to be determined later. Explicitly, we display the two-particle state Φ (k 1,k = B a1 (k 1 B a (k 1+e ik 1 1+e i(k η + 1 e i(k 1+k F(k 1 ξ (I Ã(k + cos(k 1 k +η 1+e ik 1 1+e i(k η sin( k 1 k F(k +η(1 e i(k 1 B(k. ξ. ( k Implementing a rescaling k k + η and seting c 1 = e iξ,c = e iξ +, the eigenvalue of the transfer matrix (.15 acting on the state (4.54 can be given as n τ(k Φ n (k 1,...,k n = W+ 1 (kw 1 (k + W a+1 + (kw a+1 (keink a=1 n cos ( ( k+k i cos k ki i=1 sin ( ( k+k i sin k ki +W + 4 (kw 4 (kei4n(k+η n i=1 Λ (1 aa (k,{k i},{λa j } cos ( k+k i cos ( k+k i sin ( k+k i ( cos k ki ( cos k ki sin ( k ki i=1 sin ( k+k i +η sin ( k k i Φ n (k 1,...,k n (4.58 provided that +η W + 1 (kw 1 (k W + (kw (kexp(ink k = k i = Λ (1 (k = k i,{k i },{λ j }, (4.59 i = 1, n. HereΛ (1 ( k,{ki },{λ j } istheeigenvalueofthereducedtransfermatrixτ (1 (k, which turns out to be an inhomogenous isotropic six-vertex model with open 15

16 boundary conditions, namely, τ (1 (kf a 1 a n = Λ (1 (k,{k i },{λ j }F a 1 a n, (4.60 where τ (1 (k = Tr 0 K (1 + (kt (1 (kk (1 (kt (1 1 ( k. (4.61 The nesting monodromy matrices T (1 (k and T (1 1 ( k read T (1 (k=r 1 (k +k 1 a 1e 1 ac 1 r 1 (k +k a e a 1 c r 1 (k +k n anen a n 1 c n, (4.6 T (1 1 ( k=r 1 (k k n b n 1d n b ne n r 1 (k k b 1d b e r 1 (k k 1 ad 1 b 1 e 1, (4.63 and the nesting K (1 ± (k may be chosen as Case I K (1 + (k = eiη 0 0 e iη, (4.64 in such a way that the factors W + can be rewritten as W + 1 (k = W + (k = W+ 3 (k = sin( k + η ξ +sin( k η ξ +cosη cos( k + η +ξ +cos( k η +ξ +sinksin(k +η, (4.65 sin( k + η +ξ +sin( k + η ξ +cosη cos( k η +ξ +cos( k + η +ξ +sin(k +η, (4.66 W + 4 (k = sin(k + η +ξ +sin( k + 3η +ξ + cos( k + η +ξ +cos( k η +ξ +, (4.67 and K (1 (k = 1 0, ( with W 1 (k = 1, (4.69 W (k = W 3 (k = sin(k + η +ξ sin( k η ξ sink, (

17 W 4 (k = sin( k + η +ξ sin( k + 3η +ξ cosη sin( k η ξ sin( k + η ξ sin(k +ηsin(k +η, (4.71 Case II K (1 + (k = cos(k +η ξ +e ik 0, (4.7 0 cos( k +η +ξ + with W 1 + (k = cos( k +ξ +cosηe ik sin( k ξ +sinksin(k +η, (4.73 W + (k = W+ 3 (k = cosη sin( k ξ +sin(k +η, (4.74 W 4 + (k = cos(k +η ξ +e iη sin( k ξ, ( and K (1 (k = cos(k +ξ e ik 0 0 cos( k ξ, (4.76 with W1 (k = 1, (4.77 W (k = W 3 (k = 1 sinkcos( k +ξ, (4.78 W 4 (k = cos( k +η ξ cosηe i(k+η cos( k +ξ sin(k +ηsin(k +η. (4.79 We notice that the nesting boundary K (1 ± (k-matrices also constitute four classes of BC consistent with the boundary matrix pairs (K + (k,k (k. So far, the eigenvalue problem of the 1D Bariev model with four kinds of possible boundary fields reduces to the solution of the nested auxiliary transfer matrix (4.60, which can be associated with an inhomogeneous isotropic six-vertex model with open boundary conditions (4.64 (4.76. Furthermore, we can check that the r-matrix (4.48 and the nesting K (1 ± -matrices (4.64 (4.76, which arises from the factorization of the spin degree of freedom, indeed satisfy the Yang-Baxter algebra r 1 (k λ 1 T (1 (k T (1 (λ = T (1 (λ 1 T (1 (kr 1 (k λ, (

18 and the reflection equations where r 1 (k λ = K (1 r t 1t = 1 K (1 (kr 1 (k +λ K (1 (λ (λr 1 (k +λ 1 (λ k K (1 t 1 K (1 (kr 1 (k λ, ( K (1 t 1 + (k r 1 ( k λ K (1 t + (λ + (λ r 1 ( k λ 1 K (1 t 1 + (kr t 1t (λ k, (4.8 1 r 1 (kr 1 ( k=1, (4.83 r t 1 1 ( kr t 1 (k=1, (4.84 r t 1 ( kr t 1 1 (k=1. (4.85 Following all steps in solving XXZ open chain in [1], one can completely diagonalize the nesting transfer matrix (4.61 with the help of the main ingredients (4.80 -(4.8 describing the open BC for the hidden asymmetric XXZ open chain. Therefore, what follows is the exact solution of the XXZ open chain characterizing the spin degree of freedom, i.e., Λ (1 (k,{k i },{λ j } Φ (1 (λ 1,,λ M sin(k +η M = sin(k +η W+ 11(kW11(kA (1 (k + sink M sin(k +η W+ (kw(kd (1 (k j=1 j=1 sin( k+λ j sin( k λ j η sin( k+λ j +ηsin( k λ j sin( k+λ j +ηsin( k λ j +η sin( k+λ j +ηsin( k λ j Φ (1 (λ l,,λ M, (4.86 provided that W + 11(kW 11(kA (1 (k W + (kw (kd (1 (k = M l = 1, l j sin( λ j+λ l +ηsin( λ j λ l +η sin( λ j+λ l sin( λ j λ l, (4.87 η where A (1 (k Φ (1 = Φ (1, (

19 { D (1 (k Φ (1 n = Here for case I, we have i=1 sin( k+λ i sin( k λ i sin( k+λ i +ηsin( k λ i +η } Φ (1. (4.89 W 11 + (k = W 11 (k = 1, (4.90 W + (k = W (k = 1, (4.91 while for case II, W + 11 (k = cos(k ξ +e ik, W 11 (k = cos(k +ξ e ik, (4.9 W + (k = cos(k +η +ξ +e iη, W (k = cos(k +η ξ e iη, (4.93 respectively. Finally, making a shift on the spin rapidity λ i λ i η, the eigenvalue of the transfer matrix (.15 with two classes of boundary K-matrices can be given explicitly (up to a common factor { n cos( k+k icos( k k i { τ(u Φ n (k 1,,k n = i=1 sin( k+k isin( k k i W + 1 (kw1 (k sin(k +η + sin(k +η W+ (kw (kw 11(kW + 11(ke ink M sin( k+λ l η sin( k λ } l η l=1 sin( k+λ l+η sin( k λ l+η n cos( k+k icos( k k i { + i=1 sin( k+k i +ηsin( k k i W + +η 4 (kw4 (ke i4n(k+η + sink sin(k +η W+ 3 (kw 3 (kw+ (kw (keink M sin( k+λ l+3η sin( k λ }} l+3η l=1 sin( k+λ l+η sin( k λ Φ l+η n (k 1,,k n, (4.94 where all factorized coefficients W ± i (i = 1,,4, 11 and have been given explicitly in the above equations, which provide four kinds of possible exact solutions. Finally, the expression above is in fact the eigenvalue of the transfer matrix, if the following Bethe ansatz equations are satisfied for the charge and spin rapidities k i, λ j. 19

20 M e ink i ζ 1 (k i,ξ + ζ (k i,ξ = ψ 1 (λ j,ξ + ψ (λ j,ξ = M l = 1, l j l=1 n sin( k i+λ l +η sin( k i λ l +η sin( k i+λ l η sin( k i λ l η sin( λ j+k l +η sin( λ j k l +η l=1 sin( λ j+k l η sin( λ j k l η sin( λ j+λ l +ηsin( λ j λ l +η, (4.95 sin( λ j+λ l ηsin( λ j λ l η, (4.96 j = 1,...,M, i = 1,...,n. The expressions for ζ a (k i,ξ ± and ψ a (λ j,ξ ±, a = 1, depend on which type of boundaries we are considering. Below we display all cases Case (I,I ζ 1 (k i,ξ + = sin(k i+η +ξ + sin( k i η ξ +, ζ (k i,ξ = sin( k i +η +ξ sin( k i η ξ, (4.97 ψ 1 (λ j,ξ + = ψ (λ j,ξ = 1, (4.98 Case (II,II ζ 1 (k i,ξ + = cos(k i ξ + cos( k i +ξ +, ζ (k i,ξ = e iki, (4.99 ψ 1 (λ j,ξ + = cos(λ j η ξ + e iλ j cos( λ, j+η +ξ + (4.100 ψ (λ j,ξ = cos(λ j η +ξ e iλ j cos( λ, j+η ξ + (4.101 Case (I,II Case (II,I ζ 1 (k i,ξ + = sin(k i+η +ξ + sin( k i η ξ +, ζ (k i,ξ = e iki, (4.10 ψ 1 (λ j,ξ + = 1, ψ (λ j,ξ = cos(λ j η +ξ e iλ j cos( λ j+η ξ +, (4.103 ζ 1 (k i,ξ + = cos(k i ξ + cos( k i +ξ +, ζ (k i,ξ = sin( sin( k i η k i +η +ξ ξ, (

21 ψ 1 (λ j,ξ + = cos(λ j η ξ + e iλ j cos( λ, ψ j+η (λ j,ξ = 1. ( ξ + Above M is the number of electrons with spin down and n is the total number of the electrons. Considering in addition the relations (.3 and (.4 one can find the energy spectrum for all cases Case (I,I E = coshβ c 1 coshβ +4 c n cosk i, (4.106 i=1 Case (II,II n E = exp(β(c 1 +c +4 cosk i, (4.107 i=1 Case (I,II E = exp(βc 1 coshβ +4 c n cosk i, (4.108 i=1 Case (II,I E = coshβ +exp(βc +4 c 1 n cosk i. (4.109 i=1 5. Conclusion and discussion In summary, we have analysed in detail the open boundary conditions for the Bariev chain with special boundary conditions. Two classes of boundary K + - matrices have been obtained by solving the RE. It has been found that these two classes of solutions of the RE lead to four kinds of integrable boundary fields for the charge and spin degrees of freedom separately. Through the nesting procedure, we have diagonalized exactly the two level transfer matrices with four kinds of possible boundary fields. The eigenvalues of the transfer matrices, the energy spectrum and Bethe ansatz equations have been derived for all cases. The boundary fields indeed contribute nontrivially to the groundstate properties aswell asthe low-lying spectrum. The functions ζ a (k i,ξ ± and ψ a (λ j,ξ ± will contribute nontrivial phase-shift factors to the density of the charge rapidity and spin rapidity. We notice that the Bethe equations (4.95 and (4.96 would be reduced to the purely doubling of the ones for the 1D 1

22 Bariev model with periodic BC [3][4] if the boundary parameters ξ ±. These boundary parameters ξ ±, denoting the impurity strength too, might change the band filling, the boundary surface energy, the mesoscopic effects as well. The results obtained provide us with a basis to investigate the surface critical properties, and correlation functions [14,30,31] for the resulting model. Acknowledgements X.W.G., A.F. and I.R. thank CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico for financial support. J.L. thanks the Australian Research Council for support. H.Q.Z. acknowledges the support from the NNSF of China. Appendix A. The Boltzmanweightsand theboundary K ± -matrices Before practising our ansatz for fixing the boundary K ± (u-matrix, we display the Boltzmann weights of the R-matrix (.6, which appeared first in [0] 1+h u ρ = 1+h v ρ 1+uvh 1, ρ 3 = (u vh 1+uvh ρ 1, 1+h u ρ 4 = 1+h v 1+u 1+v ρ (1+uv(1+uvh 1, ρ 5 = h 1+h v 1+u (u v ρ (1+uv(1+uvh 1, 1+h v ρ 6 = 1+u [ (u v h (u v ] ρ (1+uv(1+uvh 1, ρ 7 = 1+ ρ (1+uv(1+uvh 1, 1+u 1+v ρ 8 = ρ 1, ρ 9 = (u v(u h v 1+uv (1+uv(1+uvh ρ 1, [ (u v ] ρ 10 = 1+ ρ (1+uv(1+uvh 1, ρ 11 = (u vh 1+uv ρ 1, ρ 1 = h 1+h u 1+v (u v 1+h u ρ (1+uv(1+uvh 1, ρ 13 = 1+v (u v ρ (1+uv(1+uvh 1, ρ 14 = (u v(h u v (1+uv(1+uvh ρ h(u v 1, ρ 15 = (1+uv(1+uvh ρ 1, As mentioned in section, due to the absence of the isomorphism between K ± -matrices, we have to solve the RE (.16 (.17 separately to fix them. Substituting K (u (.18 into the RE (.16, we find, at first glance, that

23 the RE (.16 involve two variables u and v which make the functional equations involving Ka (u, a = 1,...,4 much more complicated. However, if we observe the structures of the R-matrix and the RE (.16, we may pick up some simpler functional equations of the RE which allow us to separate Ka (u,ka (v into the following forms K1 (v K (v = ρ (u,vρ 3 (v, uk1 (u+ρ 3 (u,vρ (v, uk (u ρ (u, vρ 3 ( v, uk1 (u+ρ 3 (u, vρ ( v, uk (u, K1 (v K3 (v = ρ (u,vρ 3 (v, uk1 (u+ρ 3 (u,vρ (v, uk (u ρ (u, vρ 3 ( v, uk1 (u+ρ 3 (u, vρ ( v, uk (u, K3 (v K4 (v = ρ 11(u, vρ 8 ( v, uk3 (u+ρ 8 (u, vρ 11 ( v, uk4 (u, ρ 11 (u,vρ 8 (v, uk3 (u+ρ 8 (u,vρ (v, uk4 (u Substituting the Boltzmann weights of the R-matrix into above equations, and analying the consistant conditions of these equations, we can easily conclude an ansatz K1 (u=(c 1 +c u(c 3 +c 4 u(c 5 +c 6 u, K (u=(c 1 c u(c 3 +c 4 u(c 5 +c 6 u, K3 (u=(c 1 +c u(c 3 c 4 u(c 5 +c 6 u, K4 (u=(c 1 +c u(c 3 c 4 u(c 5 c 6 u, (A.1 (A. (A.3 (A.4 with minimal coefficients c i to be determined. Running the RE again with this ansatz (A.1 (A.4, it is easily found that only one coefficient is free. Thus the two classes of boundary K -matrices can be immediately chosen as the forms (.19 and (.0. Simlarly, substituting K + -matrix (.18 into the RE (.17 we have K1 + (v K + (v = ρ (v,u ρ 16 (v, uk1 + (u+ρ 3 (v,u ρ 3 (v, uk + (u ρ 3 ( u, v ρ (u, vk1 + (u+ρ ( u, v ρ 16 (u, vk + (u, K1 + (v K3 + (v = ρ (v,u ρ 16 (v, uk1 + (u+ρ 3 (v,u ρ (v, uk3 + (u ρ 3 ( u, v ρ 3 (u, vk1 + (u+ρ ( u, v ρ 16 (u, vk3 + (u, K3 + (v K4 + (v = ρ 8( u, v ρ 15 (u, vk3 + (u+ρ 11 ( u, v ρ 7 (u, vk4 + (u, ρ 11 (v,u ρ 9 (v, uk3 + (u+ρ 8 (v,u ρ 15 (v, uk4 + (u here we first have to determine ρ s from the relation (.9 and (.10, i.e. ρ (x,y= (1+h xy (1+xy 1+h y 1+h x, (x y (h x y(h y x ρ 3 (x,y= (1+h xy (1+xy 1+h y 1+h x (x y (h x y(h y xh, 3

24 ρ 7 (x,y= (1+h xy (1+xy 1+y 1+x, (x y (h x y(h y x ρ 9 (x,y= (1+h xy (1+xy h 1+y 1+x, (x y (h x y(h y x ρ 15 (x,y= (1+h xy(1+xy (x y (xh y, ρ 16 (x,y= (1+h xy (1+xy (x y (yh xh. With the help of the new Boltzmann weights, the consistency of above equations indeed can be able to lead to an ansatz K1 + (u=(c 1 +c u(c 3 +c 4 uh (c 5 +c 6 u, K + (u=(c u h c 1 (c 4 uh +c 3 (c 5 +c 6 u, K3 + (u=(c 1 +c u(c 4 u c 3 (c 5 +c 6 u, K4 + (u=(c 1 +c u(c 4 u c 3 (c 6 u c 5 h. (A.5 (A.6 (A.7 (A.8 Running the second RE (.17 again, the solutions (.1 and (. would be fixed. Appendix B. Useful commutation relations For factorizing the transfer matrix (.15 acting on the pseudo-vacuum state, we need the following commutation relations C 1 (u B 1 (u= ρ (u, u [ B(uB(u A11 (uā11(u ], ρ 1 (u, u (B.1 C 1 (ub 1 (u= ρ [ (u, u B(u B(u Ā 11 (ua 11 (u ], ρ 1 (u, u (B. C (u B (u= ρ (u, u [ B(uB(u A (uā(u ], ρ 1 (u, u (B.3 C (ub (u= ρ [ (u, u B(u B(u Ā (ua (u ], ρ 1 (u, u (B.4 C 3 (u F(λ= ρ ] 4(u, u [ B(uB(u D(u D(u ρ 1 (u, u ρ [ (u, u C4 (uē1(u+c 5 (uē(u ], ρ 1 (u, u (B.5 C 3 (u F(u= ρ 4(u, u ρ (u, u C 1 (ub 1 (u+ ρ 8(u, u [Ā11 (ua 11 (u D(u D(u ] ρ (u, u 4

25 ρ 1(u, u ρ (u, u C 4(uĒ1(u ρ 7(u, u ρ (u, u C 5(uĒ(u, C 3 (u F(u= ρ 4(u, u ρ (u, u C (ub (u+ ρ 8(u, u ρ (u, u (B.6 [Ā (ua (u D(u D(u ] ρ 10(u, u ρ (u, u C 4(uĒ1(u ρ 1(u, u ρ (u, u C 5(uĒ(u, which can be derived directly from (3.33. In addition, one also can show the relations C i (λ B a (λ = 0 for i j, i = 1,,3, and a = 1,. Using these commutation relations and after lengthy algebra, one may manage to present the factorized forms of (3.35 and (3.39. In order to check the eigenvalue problem of two-particle excitation, the following commutation equations are necessary (B.7 ρ 1 (u,vρ 14 (v, u F(u B(v = ρ 14 (u, vρ 4 ( v, u F(v B(u +ρ 1 ( v, u [ ρ 3 (u, v B (v B 1 (u+ρ 1 (u, v F (vã11(u ] +ρ 13 ( v, u [ ρ 3 (u, v B 1 (v B (u+ρ 13 (u, v F (vã(u ] +ρ 14 ( v, u [ ρ 1 (u, v B(v F(u+ρ (u, v B 1 (vẽ1(u +ρ (u, v B (vẽ(u+ρ 4 (u, v F(v D(u ], (B.8 ρ 9 (u,vρ 14 (v, u D(u F(v+ ρ 5 (u,v [ ρ 1 (v, uã11(u F(v+ρ 11 (v, uẽ1(uẽ(v ] +ρ 6 (u,v [ ρ 13 (v, uã(u F(v+ρ 11 (v, uẽ(uẽ1(v ] +ρ 4 (u,v [ ρ 4 (v, u B(u F(v+ρ 8 (v, u B 1 (uẽ1(v +ρ 8 (v, u B (uẽ(v+ρ 1 (v, u F(u D(v ] = ρ 4 (u, vρ 1 ( v, u B(v F(u+ρ 1 (u, vρ 1 ( v, u F(v D(u +ρ 8 (u, v [ ρ 1 ( v, u B 1 (vẽ1(u+ρ 1 ( v, u B (vẽ(u ], (B.9 ρ 3 (u,v [ ρ 3 (v, uã11(u F(v+ρ 13 (v, uẽ1(uẽ(v ] +ρ (u,v [ ρ (v, u B(u F(v+ρ 1 (v, u B 1 (uẽ1(v +ρ 10 (v, u B (uẽ(v+ρ 8 (v, u F(u D(v ] = ρ 11 ( v, u [ ρ 13 (u, v B (v B 1 (u+ρ 11 (u, v F(vÃ11(u ] +ρ 8 ( v, u [ ρ (u, v B(v F(u+ρ 1 (u, v B 1 (vẽ1(u +ρ 10 (u, v B (vẽ(u+ρ 8 (u, v F(v D(u ], (B.10 ρ 3 (u,v [ ρ 3 (v, uã(u F(v+ρ 1 (v, uẽ(uẽ1(v ] +ρ (u,v [ ρ (v, u B(u F(v+ρ 7 (v, u B 1 (uẽ1(v 5

26 +ρ 1 (v, u B (uẽ(v+ρ 8 (v, u F(u D(v ] = ρ 11 ( v, u [ ρ 13 (u, v B (v B 1 (u+ρ 11 (u, v F(vÃ(u ] +ρ 8 ( v, u [ ρ (u, v B(v F(u+ρ 7 (u, v B 1 (vẽ1(u +ρ 1 (u, v B (vẽ(u+ρ 8 (u, v F(v D(u ], (B.11 ρ 3 (u,v [ ρ 3 (v, uã1(u F(v+ρ 1 (v, uẽ1(uẽ1(v ] +ρ (u,vρ 15 (v, u B (uẽ1(v = ρ 11 ( v, u [ ρ 1 (u, v B 1 (v B (u+ρ 11 (u, v F(vÃ1(u ] +ρ 15 (u, vρ 8 ( v, u B (vẽ1(u, (B.1 ρ 3 (u,v [ ρ 3 (v, uã1(u F(v+ρ 13 (v, uẽ(uẽ(v ] +ρ (u,vρ 15 (v, u B 1 (uẽ(v = ρ 11 ( v, u [ ρ 1 (u, v B 1 (v B 1 (u+ρ 11 (u, v F(vÃ1(u ] +ρ 15 (u, vρ 8 ( v, u B 1 (vẽ(u. (B.13 References [1] E.H. Lieb and F.Y. Wu, Phys. Rev. Lett. 0 ( [] F.H.L. Essler and V.E. Korepin, Exactly Solvable Models of Strongly Correlated Electrons (Singapore: World Scientific (1994 [3] P.W. Anderson, Science 35 ( ; F.C. Zhang and T.M. Rice, Phys. Rev. B37 ( [4] F.H.L. Essler and V.E. Korepin, Phys. Rev. B46 ( ; A. Foerster and M. Karowski, Phys. Rev. B46 ( ; Nucl. Phys. B396 ( [5] H.Q. Zhou, Phys. Rev. B 54 ( ; H. Asakawa and M. Suzuki, J. Phys. A: Math. Gen. 9, (1996 5; X.-W. Guan, M.-S. Wang and S.-D. Yang, Nucl. Phys. B 485 ( ; J. Phys. A: Math. Gen. 30 ( ; M. Shiroishi and M. Wadati, J. Phys. Soc. Jpn. 66 ( ; M. Shiroishi and M. Wadati, J. Phys. Soc. Jpn, 66 ( [6] F.H.L. Essler, J. Phys. A 9 ( ; G. Bedurftig, F.H.L. Essler and H. Frahm, Phys. Rev. Lett. 78 ( ; Nucl. Phys. B 489 ( [7] K. Hibberd, I. Roditi, J. Links and A. Foerster, Mod. Phys. Lett. A 15 ( ; H.-Q. Zhou X.-Y. Ge, J. Links and Mark D. Gould, Nucl. Phys. B 546 ( ; H. Fan, M. Wadati and X.-M. Wang, Phys. Rev. B 61 (

27 [8] H.-Q. Zhou, X.-Y. Ge and Mark D. Gould, J. Phys. A: Math. Gen. 3 (1999 L137; A. Foerster, K.E. Hibberd, J. Links, and I. Roditi, Phys. Lett. A 71 ( ; M.J. Martins and X.-W. Guan, Nucl. Phys. B 56 ( [9] L. Mezincescu and R.I. Nepomechie, J. Phys. A: Math. Gen. 4 (1991 L17; Int. J. Mod. Phys. A7 ( ; Int. J. Mod. Phys. A 7 ( ; R.I. Nepomechie, J. Phys. A: Math. Gen. 33 (000 L1 [10] A. Foerster and M. Karowski, Nucl. Phys. B 408 ( ; H.J. de Vega and A. Gonzalez Ruiz, J. Phys. A: Math. Gen. 6 (1993 L519. J.R. Links and M.D. Gould, Int. J. Mod. Phys. B 10 ( ; A.J. Bracken, X.Y. Ge, Y.Z. Zhang and H.Q. Zhou, Nucl. Phys. B 516 ( ; [11] H. Fan, B.Y. Hou, K.J. Shi and Z.X. Yang, Nucl. Phys. B 478 ( ; G.-L. Li, R.-H. Yue and B.-Y. Hou, Nested Bethe ansatz for Perk-Schutz model with open boundary conditions, preprint, 000; A.Lima-Santos, Nucl. Phys. B 558 ( [1] E.K. Sklyanin, J. Phys. A: Math. Gen. 1 ( [13] L.D. Faddeev, in 1984 Les Houches, in Recent advances in field theory and statistical mechanics eds. J.B. Zuber and R. Stora (North-Holland, Amsterdam, 198; P.P. Kulish and E.K. Sklyanin, in Lecture Notes in Physics, Vol. 151 (Springer, Berlin, 198 p. 61 [14] V.E. Korepin, N.M. Bogoliubov, A.G. Izergin, Quantum Inverse Scattering Method and Correlation Functions (Cambridge University Press, Cambridge, 1993 [15] G. Bedürftig and H. Frahm, J. Phys. A: Math. Gen. 3 ( ; Physica E 4 ( ; H. Asakawa and M. Suzuki, Physica A 36 ( ; H. Asakawa, Physica A 56 ( [16] F. C. Alcaraz, M.N. Barber, M.T. Batcheler, R.J. Baxter and G.R.W. Quispel, J. Phys. A: Math. Gen. 0 ( [17] S. Skorik and A. Kapustin, J. Phys. A9 ( ; S. Skorik and H. Saleur, J. Phys. A8 ( [18] R.Z. Bariev, J. Phys. A: Math. Gen. 4 (1991 L919 [19] R.Z. Bariev, A. Klümper, A. Schadschneider and J. Zittartz, J. Phys. A: Math. Gen. 6 ( [0] H.-Q. Zhou, Phys. Lett. A 1 ( ; J. Phys. A: Math. Gen. 9 ( [1] M. Shiroishi and M. Wadati, J. Phys. A: Math. Gen. 30 (

28 [] B.S. Shastry, Phys. Rev. Lett. 56 ( ,453; J. Stat. Phys. 50 ( [3] M.J. Martins and P. B. Ramos, J. Phys. A: Math. Gen. 30 (1997 L465 [4] H.-Q. Zhou, J. Phys. A: Math. Gen. 30 (1997 L43 [5] H.-Q. Zhou, Phys. Rev. B 53 ( ; Phys. Lett. A30( [6] A. Yacoby, R. Schuster and M. Meiblum, Phys. Rev. B 53 ( ; V. T. Petrashov, V. N. Antonov, P. Delsing and T. Claeson, Phys. Rev. Lett. 74 ( [7] M.J. Martins and X.-W. Guan, Nucl. Phys. B 583 ( [8] M.J. Martins and P.B. Ramos, Nucl. Phys. B 5 ( [9] X.-W. Guan, J. Phys. A: Math. Gen. 33 ( [30] A. Klümper and R.Z. Bariev, Nucl. Phys. B 458( ; Y. Fujii and M. Wadati, Chaos Soliton and Fract. 11 ( [31] T. Kojima, Finite XXZ critical chain with double boundaries, nlin.si/