Symmetry, Integrability and Geometry: Methods and Applications SIGMA 7 (2011), 012, 22 pages Spin Chains with Non-Diagonal Boundaries and Trigonometric SOS Model with Reflecting End? Ghali FILALI † and Nikolai KITANINE ‡ † Universite´ de Cergy-Pontoise, LPTM UMR 8089 du CNRS, 2 av. Adolphe Chauvin, 95302 Cergy-Pontoise, France E-mail: ghali.filali@u-cergy.fr ‡ Universite´ de Bourgogne, Institut de Mathe´matiques de Bourgogne UMR 5584 du CNRS, 9 av. Alain Savary - B.P. 47 870, 21078 Dijon, France E-mail: Nikolai.Kitanine@u-bourgogne.fr Received October 28, 2010, in final form January 11, 2011; Published online January 27, 2011 doi:10.3842/SIGMA.2011.012 Abstract. In this paper we consider two a priori very different problems: construction of the eigenstates of the spin chains with non parallel boundary magnetic fields and compu- tation of the partition function for the trigonometric solid-on-solid (SOS) model with one reflecting end and domain wall boundary conditions. We show that these two problems are related through a gauge transformation (so-called vertex-face transformation) and can be solved using the same dynamical reflection algebras. Key words: algebraic Bethe ansatz; spin chains; dynamical reflection algebra; SOS models 2010 Mathematics Subject Classification: 82B20; 82B23 1 Introduction In 1988 E. Sklyanin proposed a way [25] to generalize the algebraic Bethe ansatz [7] to the open integrable systems. In particular, this approach permitted to construct the eigenstates for the quantum XXZ spin chain with external boundary magnetic fields parallel to the z axis (diagonal boundary terms). General algebraic framework of this method (quantum inverse scattering method) worked also for the non-parallel boundary magnetic fields, i.e. the double row monodromy matrix, commuting transfer matrices and the trace identities were obtained in the most general case. However, the conservation of the z component of the total spin turned out to be essential for the eigenstates construction. It is important to mention that the algebraic Bethe ansatz establishes a clear relation between the quantum spin chains and two-dimensional models in statistical mechanics. The periodic XXZ spin chain was solved using the transfer matrix of the six-vertex model, while the open chain with diagonal boundary terms was solved using the transfer matrix of the six-vertex model with reflecting ends. This equivalence with two-dimensional models in statistical mechanics also turned out to be essential for the computation of the correlation functions. It was shown by V. Korepin [20] that the partition function of the six-vertex model with domain wall boundary conditions (DWBC) is the key element for the study of the correlation functions. The determinant representation for this partition function found by A. Izergin [14] is crucial for the computation of the correlation functions starting from the algebraic Bethe ansatz [15, 18, 19]. ?This paper is a contribution to the Proceedings of the International Workshop “Recent Advances in Quantum Integrable Systems”. The full collection is available at http://www.emis.de/journals/SIGMA/RAQIS2010.html 2 G. Filali and N. Kitanine For the open case the corresponding partition function was computed in [27]. This represen- tation was used first to compute the scalar products and norms of the Bethe vectors [28] and then to study the correlation functions of the open spin chains with external boundary magnetic fields parallel to the z axis [16, 17]. As we mentioned above, the construction of the eigenstates for the spin chains with arbitrary boundary magnetic fields is more complicated. Recent advances in the study of these systems permitted to apply the usual Bethe ansatz technique to the non-diagonal case [21, 3, 30]. This case has many interesting applications, for example it can be used to study out-of-equilibrium systems such as asymmetric simple exclusion processes (ASEP). The Bethe ansatz solution permitted to deal with problems which were out of reach for the usual techniques, such as the relaxation dynamics [5, 6]. The algebraic Bethe ansatz technique developed in [3, 30] is based on a gauge transformation (so-called vertex-face transformation) which diagonalizes the boundary matrices and relates the spin chain to a trigonometric solid-on-solid (SOS) model with reflecting ends. Once again, it can be demonstrated that these two models can be described using the same algebraic structures. The dynamical Yang–Baxter algebra proved to be the main tool to solve the SOS models [12, 8, 10, 9]. In particular, it was used to compute the partition function of the SOS model with DWBC [24, 23]. Unfortunately, the final results obtained in these papers cannot be written as a single determinant as in the six vertex case and thus, this result cannot be easily used for the computation of the correlation functions. In this paper we consider one possible connection between the spin chains with non-diagonal boundary terms and trigonometric SOS models in a systematic way. In the first part of our paper we review the construction of the eigenstates and we give a more simple picture than [3, 30] using two types of SOS reflection algebras. Finally, we consider the SOS model with one reflecting end and we use this algebraic framework to compute the partition functions for different types of domain wall boundary conditions generalizing the results of our paper [11]. The main advantage of our result is that all these partition functions can be written as a single determinant. The main aim of this paper is to relate the trigonometric SOS model with reflecting end and spin chains with non-parallel boundary magnetic fields in a simple and comprehensible way and to show the algebraic structures which can be used in the future to compute the correlation functions and form factors. It is important to note that the Bethe ansatz cannot be applied for arbitrary boundary fields. It works only if the parameters satisfy certain boundary constrains. It seems that it is absolutely necessary to have at least one relation between the boundary parameters in order to construct the Bethe ansatz. Here we consider a special case when two boundary matrices can be diagonalized simultaneously by the vertex-face transformation. It requires two constraints [30] instead of the only one obtained initially [21]. We would like to underline that this case is more complicated than the general one, as in this situation the Bethe ansatz solution is not complete (it is a degenerated situation where completeness conjecture [22] does not work). The paper is organized as follows. In Section 2 we review the quantum inverse scattering method for the spin chains with non-diagonal boundary terms. Following [25] we construct the “vertex type” double row monodromy matrix as well as a family of commuting transfer matrices. We show how the Hamiltonian can be obtained from these transfer matrices. In Sec- tion 3 we introduce the vertex-face transformation which diagonalizes the boundary matrices and transform the vertex type monodromy matrices into “SOS type” objects, satisfying dynam- ical Yang–Baxter algebra relations. We construct two double row monodromy matrices which satisfy the dynamical reflection relations or dual dynamical reflection relations. In Section 4 we use these algebras to construct eigenstates of the spin chain with non diagonal boundary terms. Finally, in the last section we use these two reflection algebras to compute the partition functions of the SOS model with reflecting end and domain wall boundary conditions. Some Spin Chains with Non-Diagonal Boundaries 3 details on the algebraic Bethe ansatz, Bethe states and relations between the two dynamical reflection algebras are given in the Appendices. 2 XXZ spin chain with non diagonal boundary We consider the open XXZ spin chain with the most general non-diagonal boundary terms H = N−1∑ i=1 ( σxi σ x i+1 + σ y i σ y i+1 + cosh ησ z i σ z i+1 ) + sinh η sinh ζ sinh δ [ − cosh ζ cosh δσz1 + sinh τσ x 1 − i cosh τσ y 1 ] + sinh η sinh ζ sinh δ [ − cosh ζ cosh δσzN + sinh τσ x N − i cosh τσ y N ] . (2.1) The Hamiltonian acts in the Hilbert space H = ⊗Nm=1Hm, Hm ∼ C 2. Two boundary magnetic fields are described here by six parameters δ, δ, ζ, ζ, τ and τ . This Hamiltonian can be obtained using the Quantum Inverse Scattering Method (QISM) [7] or more precisely its boundary version [25] starting from the Yang–Baxter equation and usual six-vertex trigonometric R-matrix. The key element of the QISM is the quantum R-matrix R : C −→ End(V ⊗ V ), V ∼ C2 satisfying the Yang–Baxter equation. R12(λ1 − λ2)R13(λ1 − λ3)R23(λ2 − λ3) = R23(λ2 − λ3)R13(λ1 − λ3)R12(λ1 − λ2). The solution corresponding to the XXZ spin chain is the six-vertex R-matrix R(λ) =     sinh(λ+ η) 0 0 0 0 sinhλ sinh η 0 0 sinh η sinhλ 0 0 0 0 sinh(λ+ η)     . This R-matrix is symmetric, namely PR12P = R12, where P is the permutation operator on V ⊗ V and satisfies the following important relations: • initial condition R(0) = sinh ηP, • unitarity R12(λ)R21(−λ) = − sinh(λ− η) sinh(λ+ η) Id, (2.2) • Z2 symmetry σy1σ y 2R12(λ)σ y 1σ y 2 = R12(λ), • crossing symmetry −σy1R t1 12(−λ− η)σ y 1 = R21(λ). (2.3) 4 G. Filali and N. Kitanine Here R21 = P12R12P12, t1 denotes the matrix transposition on the first space of the tensor product, and σx,y,z are the usual Pauli matrices. The bulk monodromy matrix T0(λ) ∈ End(V0 ⊗ H) of the inhomogeneous system is con- structed as: T0(λ) = R01(λ− ξ1) · · ·R0N (λ− ξN ). (2.4) In these last expressions, R0m denotes the R-matrix in End(V0 ⊗ Hm), and ξ1, ξ2, . . . , ξN are arbitrary complex parameters (inhomogeneity parameters) attached to the different sites of the chain of length N . It is easy to show that it satisfies the Yang–Baxter algebra relations: R12(λ1 − λ2)T1(λ1)T2(λ2) = T2(λ2)T1(λ1)R12(λ1 − λ2). This monodromy matrix can be used to construct the eigenstates of the periodic XXZ chain [7]. For the open chain a more complicated construction is needed to reconstruct the Hamiltonian and, eventually, to solve the model [25]. Taking account boundaries require the boundary K± matrix which satisfy the reflection equation [4, 25]: R12(λ12)K−(λ1)R21(λ12)K−(λ2) = K−(λ2)R12(λ12)K−(λ1)R21(λ12) (2.5) and its dual: R12(λ21)K t1 + (λ1)R21(−λ12 − 2η)K t2 + (λ2) = K t2 + (λ2)R12(−λ12 − 2η)K t1 + (λ1)R21(λ21), here we used the following short notations λ12 = λ1 − λ2, λ12 = λ1 + λ2. We consider the most general solution of these equations [13]: K−(λ) ≡ K−(λ; δ, ζ, τ) =      cosh(δ + ζ)e−λ − cosh(δ − ζ)eλ 2 sinh(δ + λ) sinh(λ+ ζ) e−τ sinh(2λ) 2 sinh(δ + λ) sinh(λ+ ζ) −eτ sinh(2λ) 2 sinh(δ + λ) sinh(λ+ ζ) cosh(δ + ζ)eλ − cosh(δ − ζ)e−λ 2 sinh(δ + λ) sinh(λ+ ζ)      and K+(λ) = K−(−λ − η; δ, ζ, τ). The complete inhomogeneous spin chain algebra is actually described by the double row monodromy matrix: (U−)0(λ) = T0(λ)(K−)0(λ)T̂0(λ) = ( A−(λ) B−(λ) C−(λ) D−(λ) ) [0] , or, equivalently, by the dual one: (U t0+ )0(λ) = T t0 0 (λ)(K t0 + )0(λ)T̂ t0 0 (λ) = ( A+(λ) C+(λ) B+(λ) D+(λ) ) [0] , where, due to the unitarity (2.2) and the crossing symmetry (2.3): T̂0(λ) ≡ RN0(λ+ ξN ) · · ·R10(λ+ ξ1) = γ(λ)σ y 0T t0 0 (−λ− η)σ y 0 = γ̂(λ)T −1 0 (−λ), with normalization coefficients: γ(λ) = (−1)N , γ̂(λ) = γ(λ) N∏ i=1 sinh(λ+ ξi − η) sinh(λ+ ξi + η). Spin Chains with Non-Diagonal Boundaries 5 The double row monodromy matrix satisfies the reflection algebra relation R12(λ12)(U−)1(λ1)R21(λ12)(U−)2(λ2) = (U−)2(λ2)R12(λ12)(U−)1(λ1)R21(λ12). The corresponding relation for the dual double monodromy matrix is R12(λ21) ( U t1+ ) 1(λ1)R21(−λ12 − 2η) ( U t2+ ) 2(λ2) = ( U t2+ ) 2(λ2)R12(−λ12 − 2η) ( U t1+ ) 1(λ1)R21(λ21). The transfer matrix TXXZ(λ) ∈ EndH can be constructed from both double row monodromy matrices TXXZ(λ) = tr0{(K+)0(λ) (U−)0(λ)} = tr0 {( Kt0− ) 0(λ) ( U t0+ ) 0(λ) } . (2.6) Here the trace is taken over the auxiliary space V0. It was shown by Sklyanin [25] that the transfer matrices commute for any value of the spectral parameter, [TXXZ(λ),TXXZ(µ)] = 0. In the homogeneous limit (ξm = 0 for m = 1, . . . , N), the Hamiltonian of the open XXZ spin chain with most general boundary fields (2.1) can be obtained as the following derivative of the transfer matrix (2.6): H = c1 d dλ TXXZ(λ) λ=0 + const, (2.7) where c1 = −8 sinh(δ) sinh(δ − η) sinh(ζ) sinh(ζ − η). This procedure permits to construct a commuting family of the conserved charges. However to construct the eigenstates in the framework of the algebraic Bethe ansatz one needs a reference state |0〉 which is an eigenstate for the operators A± andD± and annihilated by the operators C±, C±(λ)|0〉 = 0, A±(λ)|0〉 = a±|0〉, D±(λ)|0〉 = d±|0〉. Evidently, for the spin chains with diagonal boundary terms a ferromagnetic state with all the spins up is a state with such properties. But spin chain where the total spin component Sz is not conserved (like the periodic XYZ spin chain which is related to the elliptic 8-vertex model) do not posses a similar reference state. Indeed, to diagonalize the Hamiltonian one needs a gauge transformation which maps the underlying vertex model into a SOS model as it was pointed out by Faddeev and Takhtadjan [26]. In our case, even if the open XXZ bulk Hamiltonian still conserve this U(1) symmetry, the boundary terms break it. Thus our next goal is to construct a gauge transformation (vertex-face) which leads to a SOS model with a reflecting end where the boundary terms are diagonal. 3 Vertex-face transformation The idea to use a gauge transformation to diagonalize the XXZ Hamiltonian with non-diagonal boundary terms was first applied in [3], then in a more algebraic form the vertex-face transfor- mation for this case was studied in [30]. Here we propose a simpler form for this transformation and a more explicit construction of the eigenstates. 6 G. Filali and N. Kitanine 3.1 Dynamical Yang–Baxter algebra We introduce the following local gauge transformation S(λ; θ) ≡ S(λ; θ, ω) = eλ/2 ( e−(λ+θ+ω) e−(λ−θ+ω) 1 1 ) . This matrix satisfies two important properties. • It diagonalizes the boundary matrix K−(λ), K−(λ) = S −1(λ; δ − ζ, τ)K−(λ)S(−λ; δ − ζ, τ), where the diagonal matrix K−(λ) is K−(λ) ≡ K−(λ; δ, ζ) = S −1(λ; δ − ζ, τ)K−(λ)S(−λ; δ − ζ, τ) =     sinh(δ − λ) sinh(δ + λ) 0 0 sinh(ζ − λ) sinh(ζ + λ)     . (3.1) • It is a vertex-face transformation, namely: R12(λ12)S1(λ1; θ)S2(λ2; θ − ησ z 1) = S2(λ2; θ)S1(λ1; θ − ησ z 2)R12(λ12; θ), (3.2) where we introduced the dynamical R-matrix R(λ; θ) =        sinh(λ+ η) 0 0 0 0 sinhλ sinh(θ − η) sinh θ sinh η sinh(θ − λ) sinh θ 0 0 sinh η sinh(θ + λ) sinh θ sinhλ sinh(θ + η) sinh θ 0 0 0 0 sinh(λ+ η)        . (3.3) This transformation can be understood as mapping the vertex configuration into a face con- figuration in the dual lattice [2] and the R-matrix (3.3) contains the statistical weights of the trigonometric SOS model. It is a very well studied object generating the Felder’s dynamical Yang–Baxter algebra [8]. It solves the dynamical Yang–Baxter equation (DYBE), R12(λ12; θ − ησ z 3)R13(λ13; θ)R23(λ23; θ − ησ z 1) = R23(λ23; θ)R13(λ13; θ − ησ z 2)R12(λ12; θ), and satisfies the following properties which will be essential to construct the double row mono- dromy matrices: • Ice rule [σz1 + σ z 2 ,R12(λ; θ)] = 0. This symmetry is responsible of the six vertex texture of the statistical weight: Rµναβ = 0 unless α+ β = µ+ ν. It is easy to see that this relation induces a similar relation for the transposed R-matrix: [σz1 − σ z 2 ,R t1 12(λ; θ)] = 0, which translate into: (Rt)µναβ = 0 unless α− β = µ− ν. Spin Chains with Non-Diagonal Boundaries 7 • Unitarity: R12(λ; θ)R21(−λ; θ) = − sinh(λ− η) sinh(λ+ η) Id. • Crossing symmetries. The crossing relations for the dynamical R-matrix are not as simple as for the vertex type R-matrix, here we write it in the following compact form1: −σy1R t1 12(−λ− η; θ + ησ z 1)σ y 1 sinh(θ − ησz2) sinh θ = R21(λ; θ), (3.4) −σy1R t1 21(−λ− η; θ − ησ z 1)σ y 1 sinh(θ + ησz2) sinh θ = R12(λ; θ). (3.5) We should stress that these two relations are not equivalent at all since the dynamical R-matrix are not symmetric, but there is a weaker condition • Parity: R21(λ; θ) = σ x,y 1 σ x,y 2 R12(λ; θ)σ x,y 1 σ x,y 2 = R12(λ;−θ). (3.6) Now we can easily introduce the dynamical monodromy matrix as an ordered product of the dynamical R-matrices T0(λ; θ) = R01 ( λ− ξ1; θ − η N∑ i=2 σzi ) · · ·R0N (λ− ξN ; θ). (3.7) T0(λ; θ) satisfies the zero weight condition: [T0(λ; θ), σ z 0 + S z] = 0. It is easy to show that this monodromy matrix satisfies the dynamical Yang–Baxter relation R12(λ12; θ − ηSz)T1(λ1; θ)T2(λ2; θ − ησz1) = T2(λ2; θ)T1(λ1; θ − ησ z 2)R12(λ12; θ), where the z component of the total spin Sz = N∑ i=1 σzi . This dynamical monodromy matrix is related to the usual one (2.4) by the vertex-face trans- formation S, S−({ξ}; θ)S0(λ; θ − ηSz)T0(λ; θ) = T0(λ)S0(λ; θ)S−({ξ}; θ − ησz0), where S−({ξ}; θ) = SN (ξN ; θ) · · ·S1 ( ξ1; θ − η N∑ i=2 σzi ) . We need to introduce a second vertex-face transformation to proceed with the construction of the double row dynamical monodromy matrix (it will be necessary to construct a dynamical analog of the U+ double row monodromy matrix in the next section). As for the vertex case, dual algebras in dynamical context are closely related to the antipode of the corresponding dynamical quantum group, which, due to the complicated form of the crossing symmetries, is not represented by the simple matrix transposition or inversion. 1We assume in this paper the following normal ordering: the σz0 in the argument of the R01 matrix (which does not commute with it) is always on the right of all other operators involved in the definition of R. 8 G. Filali and N. Kitanine The matrix R(λ; θ) can be also obtained by the second vertex-face transformation defined by the same matrix S R12(λ12)S2(λ2; θ)S1(λ1; θ + ησ z 2) = S1(λ1; θ)S2(λ2; θ + ησ z 1)R12(λ12; θ), which leads to the following dynamical Yang–Baxter equation R12(λ12; θ)R13(λ13; θ + ησ z 2)R23(λ23; θ) = R23(λ23; θ + ησ z 1)R13(λ13; θ)R12(λ12; θ + ησ z 3). (3.8) Now we can introduce the “crossed” L-operator, which will be used to construct the dynamical analog of the double row monodromy matrix U+ Lt112(λ; θ) = R t1 12(λ; θ + ησ z 1) sinh(θ − ησz2) sinh θ , and its inverse (in a convenient for us normalization) L̂t121(λ; θ) = R t1 21(λ; θ − ησ z 1) sinh(θ + ησz2) sinh θ . Due to the dynamical Yang–Baxter equation (3.8) and the crossing symmetry (3.4), this new L-operator also satisfies a dynamical Yang–Baxter type equation: R12(λ21; θ + ησ z 3)L t1 13(λ13; θ)L t2 23(λ23; θ − ησ z 1) = L t2 23(λ23; θ)L t1 13(λ13; θ − ησ z 1)R12(λ21; θ). It is easy to derive the properties of this L-operator, such as unitarity L̂t121(−λ− η; θ)L t1 12(λ− η; θ) = − sinh(λ− η) sinh(λ+ η) Id, ice rule [σz1 − σ z 2 ,L t1 12(λ; θ)] = [σ z 1 − σ z 2 , L̂ t1 21(λ; θ)] = 0, and parity σx,y1 σ x,y 2 L t1 12(λ; θ)σ x,y 1 σ x,y 2 = L̂ t1 21(λ; θ) = L t1 12(λ;−θ), from the corresponding symmetries of the dynamical R-matrix. Starting from this L-operator we can construct a dual monodromy matrix Vt00 (λ; θ) = L t0 0N ( λ− ξN ; θ + η N−1∑ i=1 σzi ) · · · Lt001(λ− ξ1; θ). (3.9) The Yang–Baxter relation for this matrix can be written in the following form R12(λ21; θ + ηSz)V t1 1 (λ1; θ)V t2 2 (λ2; θ − ησ z 1) = V t2 2 (λ2; θ)V t1 1 (λ1; θ − ησ z 2)R12(λ21; θ). Once again, it is possible to construct a gauge transformation which relates this matrix to the transposed vertex monodromy matrix S+({ξ}; θ)S˜(λ+ η; θ + ηSz)0V t0 0 (λ; θ) = T t0 0 (λ)S˜(λ+ η; θ)0S+({ξ}; θ − ησ z 0), where S+({ξ}; θ) = S1(ξ1; θ) · · ·SN ( ξN ; θ + η N−1∑ i=1 σzi ) , and S˜0(λ; θ) = σ y 0S0(λ; θ)σ y 0 . Spin Chains with Non-Diagonal Boundaries 9 3.2 Dynamical reflection algebra Now we can proceed with a construction of the double row monodromy matrix. The key ingre- dient is once again the reflection equation. It easy to check that for the diagonal K matrix (3.1) satisfy the following dynamical reflection relation for θ = δ − ζ R12(λ12; θ)(K−)1(λ1)R21(λ12; θ)(K−)2(λ2) = (K−)2(λ2)R12(λ12; θ)(K−)1(λ1)R21(λ12; θ). (3.10) This is essentially the reflection equation introduced in [4], with the dynamical R-matrix. This equation can be checked directly but also it can be derived from the vertex reflection equation using the gauge transformation S: Proposition 3.1. The vertex type reflection equation (2.5) for K− is equivalent to the SOS type reflection equation for K− if we the SOS parameter is chosen as θ = δ − ζ, ω = τ . Using (3.2) the left hand side of the reflection equation (2.5) can be rewritten as R12(λ12)K−(λ1)R21(λ12)K−(λ2) = S2(λ2; θ)S1(λ1; θ − ησ z 2)R12(λ12; θ)K1(λ1) ×R21(λ12; θ)K2(λ2)S −1 1 (−λ1; θ − ησ z 2)S −1 2 (−λ2; θ). To obtain this formula we used the fact that K− is diagonal. Similar calculation for the right hand side give us: K−(λ2)R12(λ12)K−(λ1)R21(λ12) = S2(λ2; θ)S1(λ1; θ − ησ z 2) K2(λ2)R12(λ12; θ) ×K1(λ1; )R21(λ12; θ)S −1 1 (−λ1; θ − ησ z 2)S −1 2 (−λ2; θ). Which finally leads to the equation (3.10). In a similar way one can prove that the dynamical double row monodromy matrix defined as (U−)0(λ; θ) = T0(λ; θ)(K−)0(λ)T̂0(λ; θ) = ( A−(λ; θ) B−(λ; θ) C−(λ; θ) D−(λ; θ) ) [0] , (3.11) where T̂0(λ; θ) ≡ RN0(λ+ ξN ; θ) · · ·R10 ( λ+ ξ1; θ − η N∑ i=2 σzi ) = γ(λ)σy0T t0 0 (−λ− η; θ + ησ z 0)σ y 0 sinh(θ − η Sz) sinh(θ) = γ̂(λ)T −10 (−λ; θ), satisfies the dynamical reflection algebra relation, R12(λ12; θ − η Sz)(U−)1(λ1; θ)R21(λ12θ − ηSz)(U−)2(λ2; θ) = (U−)2(λ2; θ)R12(λ12; θ − η Sz)(U−)1(λ1; θ)R21(λ12; θ − η Sz). (3.12) This double row monodromy matrix has a clear interpretation in statistical mechanics, it is an essential tool to study the trigonometric SOS model with a reflecting end. We will return to this correspondence in the last section of this paper. It can be easily shown that the vertex-face transformation S gives a relation between this dynamical double row monodromy matrix and the usual one S−({ξ}; θ)S0(λ; θ − ηSz)(U−)0(λ; θ) = (U−)0(λ)S−({ξ}; θ)S0(−λ; θ − ηSz). (3.13) 10 G. Filali and N. Kitanine 3.3 Dual reflection algebra In the vertex case we constructed two double row monodromy matrices: one around the K− matrix and second one around the K+ matrix. The same procedure is possible in the dynamical case. The only difference is that because of a more complicated crossing relations (3.4) and (3.5) instead of the transposed monodromy matrix we will use here the dual monodromy matrix Vt0(λ) (3.9). We start once again with a diagonal solution K+(λ) of the dual reflection equation with the dynamical R-matrix, R12(λ21; θ) ( Kt1+ ) 1(λ1)R21(−λ12 − 2η; θ) ( Kt2+ ) 2(λ2) = ( Kt2+ ) 2(λ2)R12(−λ12 − 2η; θ) ( Kt1+ ) 1(λ1)R21(λ21; θ). It can be expressed in terms of the K− matrix, K+(λ) ≡ K+(λ; δ, ζ) = K−(−λ− η; δ, ζ). The crucial point here is that both the vertex boundary matrix K+ and the dynamical bounda- ry K+ are related each other by mean of the second Vertex-face transformation: Kt+(λ) = S˜ −1(λ+ η; δ − ζ, τ)Kt+(λ)S˜(−λ− η; δ − ζ, τ). Note that the dynamical parameter θ should not a priori coincide with the parameter θ from the previous section as it is chosen here θ = δ − ζ to diagonalize the boundary matrix K+. As in the previous section we can construct the following double row monodromy matrix (U+) t0 0 (λ; θ) = V t 0(λ; θ) ( Kt0+ ) 0(λ)V̂ t0 0 (λ; θ) = ( A+(λ; θ) C+(λ; θ) B+(λ; θ) D+(λ; θ) ) [0] , where V̂t00 (λ; θ) ≡ L̂ t0 10(λ+ ξ1; θ) · · · L̂ t0 N0 ( λ+ ξN ; θ + η N−1∑ i=1 σzi ) = γ˜(λ) ( Vt00 )−1 (−λ− 2η; θ). The normalization factor in the last equation can be easily computed using the unitarity of the operator L, γ˜(λ) = γ(λ) N∏ i=1 sinh(λ+ ξi) sinh(λ+ ξi + 2η). This double row monodromy matrix satisfies the following dynamical reflection equation R12(λ21; θ + ηSz)(U t1 + )1(λ1; θ)R21(−λ12 − 2η; θ + η S z))(U t2+ )2(λ2; θ) = (U t2+ )2(λ2; θ)R12(−λ12 − 2η; θ + η S z)(U t1+ )1(λ1; θ)R21(λ21; θ + η S z). (3.14) Its relation with the corresponding vertex double row monodromy matrix can be easily written in the following form S+({ξ}; θ)S˜0(λ+ η; θ + η Sz) ( U t+ ) 0(λ; θ) = ( U t+ ) 0(λ)S+({ξ}; θ)S˜0(−λ− η; θ + ηS z). In the next section we show how to use these matrices to construct eigenstates of the open XXZ spin chain with non-diagonal boundary terms. Spin Chains with Non-Diagonal Boundaries 11 4 Algebraic Bethe ansatz Our goal is to construct eigenstates of the spin chain Hamiltonian (2.1) which is equivalent to the construction of eigenstates of the transfer matrix TXXZ(λ) (2.6). We have seen that usual algebraic Bethe ansatz procedure cannot be applied to this case as there is no simple reference state. Thus our strategy is to consider dynamical monodromy matrix instead of vertex one using the vertex-face correspondence. The advantage of this approach is the fact that this procedure allows us (under certain constraints) to diagonalize the boundary matrices, and hence, it becomes possible to use the algebraic Bethe ansatz techniques. Using (3.13) the transfer matrix (2.6) can rewritten in terms of the dynamical algebra gene- rators TXXZ(λ)S−({ξ}, θ) = S−({ξ}, θ) × Tr0 { K+(λ)S0(λ; θ − ηSz)(U−)0(λ; θ)S −1 0 (−λ; θ − ηS z) } . (4.1) This trace can be rewritten in terms of the operator entries of the double row monodromy matrix Tr0{K+(λ)S0(λ; θ − ηSz)(T−)0(λ; θ)S −1 0 (−λ; θ − ηS z)} = {S−10 (−λ; θ − ηS z)K+(λ)S0(λ; θ − ηSz)} + +A−(λ; θ) + {S−10 (−λ; θ − ηS z)K+(λ)S0(λ; θ − ηSz)} − −D−(λ; θ) + {S−10 (−λ; θ − ηS z − 2η)K+(λ)S0(λ; θ − ηSz)} − +B−(λ; θ) + {S−10 (−λ; θ − ηS z + 2η)K+(λ)S0(λ; θ − ηSz)} + −C−(λ; θ). The next crucial step is to restrict our analysis to a subspace with a fixed z component of the total spin (in the face picture, it has nothing to do with the XXZ spin). This restriction follows from the results of [3, 21], it was shown that the number of Bethe roots for all the eigenstates is the same. It means that in the framework of the algebraic Bethe ansatz the number of creation operators and hence the total spin should be fixed. Unfortunately it also means that in this approach the algebraic Bethe ansatz cannot lead to the complete description of the eigenstates. Thus we consider the action of the trace in (4.1) on the states |ψ〉 with a fixed z component of the total spin Sz|ψ〉 = s|ψ〉. To apply the algebraic Bethe ansatz we need to diagonalize both boundary matrices simul- taneously, it means that the non diagonal terms should be zero, { S−10 (−λ; δ − ζ − ηs, τ)K+(λ)S0(λ; δ − ζ − η(s− 2), τ) }− + = 0, { S−10 (−λ; δ − ζ − ηs, τ)K+(λ)S0(λ; δ − ζ − η(s+ 2), τ) }+ − = 0. These lead to two constraints on the boundary parameters: cosh(δ − ζ) = cosh(δ − ζ − ηs+ τ − τ − η), (4.2) cosh(δ − ζ) = cosh(δ − ζ − ηs− τ + τ + η), (4.3) which can be solved by imposing: τ = τ + η + ipin, δ − ζ = δ − ζ − ηs+ 2ipim. (4.4) With these constraints we obtain the following result: TXXZ(λ)S−({ξ}; θ)|ψ〉 = S−({ξ}; θ)Tr0(K˜+(λ; δ, ζ)T−(λ; θ))|ψ〉, 12 G. Filali and N. Kitanine where K˜+(λ; δ, ζ) = sinh(δ − ζ − ησz) sinh(δ − ζ) K−(−λ− η; δ, ζ). Therefore the eigenstates of the vertex transfer matrix TXXZ can be obtained from the eigen- states of the SOS transfer matrix TSOS1(λ; θ) = Tr0 ( K˜+(λ; δ, ζ)U−(λ; θ) ) . (4.5) Using the reflection algebra relation (3.12) which contains the commutation relations for the operators A−(λ; θ), B−(λ; θ), C−(λ; θ), and D−(λ; θ), we can actually construct two sets of eigenstates of TSOS1(λ; θ). The advantage of the SOS picture is the existence of the reference state |0〉 with all the spins up, which is an eigenstate of the operators A−(λ; θ) and D−(λ; θ) and it is annihilated by the operator C−(λ; θ), C−(λ; θ)|0〉 = 0. Then the eigenstates are constructed using the algebraic Bethe ansatz technique: Proposition 4.1. Let θ = δ − ζ, and δ, ζ, δ and ζ satisfy the boundary constraints (4.4) with total spin s being even if N is even and odd if N is odd, |s| < N . Then the state |ψ1−({λ}M )〉 = B−(λ1; θ) · · · B−(λM ; θ)|0〉, M = N − s 2 , (4.6) is an eigenstate of TSOS1(µ; θ) with eigenvalue Λ1(µ; {λ}; δ, ζ, δ, ζ) = sinh(ζ − µ) sinh(δ + µ) sinh(δ − µ) sinh(2µ+ 2η) sinh(ζ − µ− η) sinh(δ − µ− η) sinh(δ + µ) sinh(2µ+ η) × M∏ i=1 sinh(µ+ λi) sinh(µ− λi − η) sinh(µ+ λi + η) sinh(µ− λi) N∏ i=1 sinh(µ+ ξi + η) sinh(µ− ξi + η) + sinh(ζ + µ+ η) sinh(δ + µ+ η) sinh(ζ − µ− η) sinh 2µ sinh(ζ − µ− η) sinh(δ + µ) sinh(ζ + µ) sinh(2µ+ η) × M∏ i=1 sinh(µ+ λi + 2η) sinh(µ− λi + η) sinh(µ+ λi + η) sinh(µ− λi) N∏ i=1 sinh(µ+ ξi) sinh(µ− ξi), provided the λi, i = 1, . . . ,M satisfy the Bethe equation: y1(λi, {λ}; δ, ζ, δ, ζ) = y1(−λi − η, {λ}; δ, ζ, δ, ζ), i = 1, . . . ,M, (4.7) with y1(λi, {λ}; δ, ζ, δ, ζ) = sinh(ζ + λi) sinh(δ − λi) sinh(ζ − λi) sinh(δ + λi) × M∏ k 6=i sinh(λi + λk) sinh(λi − λk − η)× N∏ j=1 sinh(λi + ξj + η) sinh(λi − ξj + η). The proof of this proposition is a standard algebraic Bethe ansatz computation, some details on the commutation relations leading to this result are given in the Appendix A. It is possible to construct a second set of Bethe states starting from the second reference state |0〉 with all the spins down. It is also an eigenstate of the operators A−(λ; θ) and D−(λ; θ), while B−(λ; θ)|0〉 = 0. Spin Chains with Non-Diagonal Boundaries 13 Proposition 4.2. |ψ2−({λ}M )〉 = C−(λ1; θ) · · · C−(λM ; θ)|0〉, M = N + s 2 , (4.8) is an eigenstate of TSOS1(µ; θ) with eigenvalue: Λ2(µ; {λ}; δ, ζ, δ, ζ) = Λ1(µ; {λ}; ζ, δ, ζ, δ), (4.9) provided the λi, i = 1, . . . ,M satisfy the Bethe equations: y2(λi, {λ}; δ, ζ, δ, ζ) = y2(−λi − η, {λ}; δ, ζ, δ, ζ), i = 1, . . . ,M, (4.10) with y2(λi, {λ}; δ, ζ, δ, ζ) = y1(λi, {λ}; ζ, δ, ζ, δ). (4.11) These two set of Bethe states are actually related each other. More details on this relations are given the Appendix B. Now it is easy to construct the eigenstates of the XXZ transfer matrix; Proposition 4.3. Let θ = δ − ζ, and δ, ζ, δ and ζ satisfy the boundary constraints (4.4) with total spin s being even if N is even and odd if N is odd, |s| < N . Then the states S−({ξ}; θ)|ψ 1 −({λ}MB )〉, with MB = N−s2 are eigenstates of TXXZ(λ) provided λ1, . . . , λMB satisfy the Bethe equa- tions (4.7). A second set of Bethe states can be constructed as S−({ξ}; θ)|ψ 2 −({λ}MC )〉, with MC = N+s2 , if λ1, . . . , λMC satisfy (4.10). The corresponding eigenvalues of the Hamiltonian is then given by (2.7) E = M∑ j=1 ε(λj) + c1N coth η, ε(λj) = c1 sinh η sinh(λ+ η) sinhλ . Thus we obtained a simple algebraic expression for the eigenstates of the spin chain with non- diagonal boundaries in terms of the SOS operators and vertex-face transformation S. This method however has two clear limitations: first, it works in its present form only if the con- straints (4.4) are satisfied, and second, it permits to construct only states with fixed number of operators B−(λ; θ) (or C−(λ; θ)) and this set of states is evidently incomplete (it is clear from the SOS picture). Thus the problem of the complete description of the spectrum remains open. The construction of the same eigenstates can be equivalently performed starting from the U+(λ; θ) reflection algebra. Following the same lines as before we obtain: TXXZ(λ)S+({ξ}; θ) = S+({ξ}; θ) × Tr0 { S˜0(λ+ η; θ + ηSz)(U t0 + )0(λ; θ)S˜ −1 0 (−λ− η; θ + ηS z)Kt−(λ) } . Again, it is possible to write the trace in terms of the operator entries of the double row monod- romy matrix the U+(λ, θ). It is easy to check that if we consider the action of this trace on the states with z component of total spin s and if the constraints (4.4) are satisfied the non-diagonal terms in this expressions vanish TXXZ(λ)S+({ξ}; θ)|ψ〉 = S+({ξ}; θ)Tr0 ( U t0+ (λ; θ)K˜ t0 − (λ; δ, ζ) ) |ψ〉, 14 G. Filali and N. Kitanine where K˜−(λ; δ, ζ) = sinh(δ − ζ − ησz) sinh(δ − ζ) K+(−λ− η; δ, ζ). It means that we replace once again the vertex transfer matrix with non-diagonal boundary matrices by a SOS transfer matrix with diagonal boundaries. TSOS2(λ, θ) = Tr0 ( U t0+ (λ; θ)K˜ t0 − (λ; δ, ζ) ) . Proposition 4.4. Let θ = δ − ζ, and δ, ζ, δ and ζ satisfy the boundary constraints (4.4) with total spin s being even if N is even and odd if N is odd, |s| < N . |ψ1+({λ}M )〉 = B+(λ1; θ) · · · B+(λM ; θ)|0〉, M = N − s 2 , is an eigenstate of TSOS2(µ; θ) with eigenvalue Λ2(µ; δ, ζ, δ, ζ) provided the λi, i = 1, . . . ,M satisfy the Bethe equations y2(λi, {λ}; δ, ζ, δ, ζ) = y2(−λi − η, {λ}; δ, ζ, δ, ζ), i = 1, . . . ,M. (4.12) Proposition 4.5. |ψ2+({λ})〉 = C+(λ1; θ) · · · C+(λM ; θ)|0〉, M = N + s 2 , is an eigenstate of TSOS2(µ; θ) with eigenvalue Λ1(µ; δ, ζ, δ, ζ) provided the λi, i = 1, . . . ,M satisfy the Bethe equations y1(λi, {λ}; δ, ζ, δ, ζ) = y1(−λi − η, {λ}; δ, ζ, δ, ζ), i = 1, . . . ,M. (4.13) Finally we can once again construct the corresponding eigenstates for the XXZ transfer matrix using the vertex-face transformation S+. Proposition 4.6. Let θ = δ − ζ, and δ, ζ, δ and ζ satisfy the boundary constraints (4.4) with total spin s being even if N is even and odd if N is odd, |s| < N . Then the states S+({ξ}; θ)|ψ 1 +({λ}MB )〉, with MB = N−s2 are eigenstates of TXXZ(λ) provided λ1, . . . , λMB satisfy the Bethe equa- tions (4.12). A second set of Bethe states can be constructed as S+({ξ}; θ)|ψ 2 +({λ}MC )〉, with MC = N+s2 , if λ1, . . . , λMC satisfy (4.13). It is easy to check that Λ1(µ; δ, ζ, δ, ζ) = Λ2(µ; δ, ζ, δ, ζ), y1(λi, {λ}; δ, ζ, δ, ζ) = y2(λi, {λ}; δ, ζ, δ, ζ). It means that the eigenstates constructed by the first method (Proposition 4.3) and the states constructed using the dual algebra (Proposition 4.6) are exactly the same states. However it was already shown for diagonal boundary conditions [16] that it is important to be able to construct the same states by two different methods to proceed with the computation of the scalar products and correlation functions. Spin Chains with Non-Diagonal Boundaries 15 5 SOS model with reflecting end In the previous section we used the trigonometric dynamical reflection algebra to construct the eigenstates of the spin chains with non-diagonal boundary fields. However this algebra has another clear interpretation, it describes a trigonometric solid-on-solid (SOS) model with one reflecting end (in the same way as usual reflection algebra with diagonal matrix K describes a six-vertex model with reflecting end). The SOS model is a two dimensional statistical mechanics lattice model which can be defined in terms of a height function. Every square of the lattice is characterized by a height θ and its values for two adjacent squares differ by η. There are 6 possible face configurations θ − η θ − 2η θ θ − η θ + η θ + 2η θ θ + η θ − η θ θ θ + η θ + η θ θ θ − η θ + η θ θ θ + η θ − η θ θ θ − η and the corresponding statistical weights Rabcd are given by the dynamical R-matrix (3.3). We consider this model with a reflecting end, which means that each horizontal line makes a U-turn on the left (right) side of the lattice. It produces two following configurations charac- terized by the boundary matrix K−(λ) (Kt+(λ)): Figure 1. Boundary configuration with external height θ. It was seen in the vertex case that it is crucial to study the partition function of this model with domain wall boundary conditions. It is, in particular, necessary to study the correlation functions of the spin chains with non-diagonal boundary terms. These domain wall boundary conditions can be easily written for both cases: with left or right reflection end, here we give the example of the left reflection end (given by the matrix K−), see Fig. 2. Here the heights decrease from left to rights on the upper boundary, the heights grow from left to right on the lower boundary (we will call this situation case I). As left external height is fixed these two conditions determine completely the configuration on the right boundary (heights decreasing in the upward direction). Equivalently one can consider heights increasing on the upper boundary and decreasing on the lower one (case II). The same two situations (case I and case II) are possible for the reflection on the right boundary (with boundary weights given by the matrix Kt+). Thus we have four possible configurations and for all of them the partition functions can be written in terms of the double row monodromy matrix. Namely for the left reflecting end the 16 G. Filali and N. Kitanine θ-Nηθ-(N-1)η θ θ θ θ θ+Nη θ-η θ-(N-1)η θ-(N-2)η θ-(N-3)η θ+(N-3)η θ+(N-2)η θ+(N-1)η θ+(N-1)ηθ+η Figure 2. Domain wall boundary conditions. case I partition function can be written as ZB−N,2N ({λ}, {ξ}, δ, ζ) = 〈0¯| N∏ i=1 B−(λi; θ)|0〉, while for the case II we obtain ZC−N,2N ({λ}, {ξ}, δ, ζ) = 〈0| N∏ i=1 C−(λi; θ)|0¯〉, where we used our standard notations θ = δ − ζ. If the reflecting end is on the right, we obtain the following representations: ZB+N,2N ({λ}, {ξ}, δ, ζ) = 〈0¯| N∏ i=1 B+(λi; θ)|0〉, ZC+N,2N ({λ}, {ξ}, δ, ζ) = 〈0| N∏ i=1 C+(λi; θ)|0¯〉. All these partition functions can be computed following the same lines. The detailed derivation of the function ZB−N,2N is given in [11], it is shown that this function is uniquely defined by the following properties • For each parameter λi the normalized partition function Z˜N,2N ({λ}, {ξ}, δ, ζ) = exp ( (2N + 2) N∑ i=1 λi ) × sinh(δ + λi) sinh(ζ + λi)Z B− N,2N ({λ}, {ξ}, δ, ζ), is a polynomial of degree at most 2N + 2 in e2λi . Spin Chains with Non-Diagonal Boundaries 17 • For N = 1 the partition function is ZB−1,2 (λ, ξ, δ, ζ) = sinh η sinh(δ − ζ − η) sinh2(δ − ζ) ( sinh(δ − λ) sinh(δ + λ) sinh(λ− ξ) sinh(δ − ζ + λ+ ξ) + sinh(ζ − λ) sinh(ζ + λ) sinh(λ+ ξ) sinh(δ − ζ − λ+ ξ) ) . • ZB−N,2N ({λ}, {ξ}, δ, ζ) is symmetric in λi. • ZB−N,2N ({λ}, {ξ}, δ, ζ) is symmetric in ξi. • Crossing symmetry ZB−N,2N (−λi − η, {λ}, {ξ}, δ, ζ) = − sinh(2(λi + η)) sinh(λi + ζ) sinh(2λi) sinh(λi − ζ + η) × sinh(λi + δ) sinh(λi − δ + η) ZB−N,2N (λi, {λ}, {ξ}, δ, ζ). • Recursive relations ZB−N,2N ({λ}, {ξ}, δ, ζ) ∣ ∣ ∣ λ1=ξ1 = sinh η sinh(ζ − λ1) sinh(ζ + λ1) × N∏ i=1 sinh(λi + ξ1) sinh(δ − ζ + (N − 2i)η) sinh(δ − ζ + (N − 2i+ 1)η) × N∏ i=2 sinh(λ1 − ξi + η) sinh(λ1 + ξi + η) sinh(λi − ξ1 + η) × ZB−(N−1),2(N−1)({λ}2,...,N , {ξ}2,...,N , δ, ζ), ZB−N,2N ({λ}, {ξ}, δ, ζ) ∣ ∣ ∣ λN=−ξ1 = sinh η sinh(δ − λN ) sinh(δ + λN ) × N∏ i=1 sinh(λi − ξ1) sinh(δ − ζ + (N − 2i)η) sinh(δ − ζ + (N − 2i+ 1)η) × N∏ i=2 sinh(λN + ξi + η) sinh(λN − ξi + η) sinh(λi−1 + ξ1 + η) × ZB−(N−1),2(N−1)({λ}1,...,N−1, {ξ}2,...,N , δ, ζ). It means that if we find a function satisfying all these properties it is the partition function. The following proposition can be easily proved by the direct verification of these properties: Proposition 5.1. ZB−N,2N ({λ}, {ξ}, δ, ζ) = γ(λ) detMij N∏ i=1 ( sinh(δ − ζ + η(N − 2i) sinh(δ − ζ + η(N − i)) ) × N∏ i,j=1 sinh(λi + ξj) sinh(λi − ξj) sinh(λi + ξj + η) sinh(λi − ξj + η) ∏ 1≤i