arxiv:math-ph/ v3 22 Oct 2004

Transcription

1 General boundary conditions for the sl(n) and sl( N) open spin chains arxiv:math-ph/040601v3 Oct 004 D. Arnaudon a1, J. Avan b1, N. Crampé a1, A. Doikou a1, L. Frappat ac1, E. Ragoucy a1 a Laboratoire d Annecy-le-Vieux de Physique Théorique LAPTH, CNRS, UR 5108, Université de Savoie B.P. 110, F Annecy-le-Vieux Cedex, France b Laboratoire de Physique Théorique et odélisation Université de Cergy, 5 mail Gay-Lussac, Neuville-sur-Oise F Cergy-Pontoise Cedex c ember of Institut Universitaire de France Abstract Two types of boundary conditions ( soliton preserving and soliton non-preserving ) are investigated for the sl(n) and sl( N) open spin chains. The appropriate reflection equations are formulated and the corresponding solutions are classified. The symmetry and the Bethe Ansatz equations are derived for each case. The general treatment for non-diagonal reflection matrices associated to soliton preserving case is worked out. The connection between the soliton non-preserving boundary conditions and the twisted (super) Yangians is also discussed. SC: 81R50, 17B37 PACS: 0.0.Uw, Fd, Pq math-ph/ LAPTH-1050/04 June arnaudon@lapp.inp3.fr, avan@ptm.u-cergy.fr, crampe@lapp.inp3.fr, doikou@lapp.inp3.fr, frappat@lapp.inp3.fr, ragoucy@lapp.inp3.fr

2 Introduction The possibility of constructing and(at least partially) solving by algebraic and/or analytical methods, one-dimensional interacting quantum spin chains, is one of the maor achievements in the domain of quantum integrable systems. Its main tool is the quantum R-matrix, obeying a cubic Yang-Baxter equation, the coproduct properties of which allow the building of an L-site transfer matrix with identical exchange relations and the subsequent derivation of quantum commuting Hamiltonians [1]. The same structure is instrumental in formulating the quantum inverse scattering procedure, initiated by the Leningrad school []. A subsequent development was the definition of exactly solvable open spin chains with non-trivial boundary conditions. These are characterised by a second obect: the reflection matrix K, obeying a quadratic consistency equation with the R matrix, with the generic abstract form RKRK = KRKR [3 7]. Using again coproduct-like properties of this structure one constructs suitable transfer matrices yielding (local) commuting spin chain Hamiltonians by combining K and semitensor products of R [4]. any efforts have been devoted to this issue [8 3], based on the pioneering approach of Sklyanin and we here aim at treating a particular, but very significant, class of examples for this problem. To better characterise the type of spin chain which we will be considering here it is important to recall that both R and K matrices have an interpretation in terms of diffusion theory for particlelike obects identified in several explicit cases with exact eigenstates of some quantum integrable field theories such as sine Gordon [4], non-linear Schrödinger equation [4, 5] or principal chiral model [6]. R describes the basic -body scattering amplitudes and K describes a 1-body scattering or reflection on a boundary. The Yang-Baxter equation (YBE) and reflection equation (RE) then characterise consistent factorisability of any L-body amplitude in terms of 1- and -body scattering amplitudes, regardless of the order of occurrence of the 1- and -body events in the diffusion process. As a consequence, when one describes the scattering theory of a model with more than one type of particle involved, one is led to introduce several operators of R and K type. The case which we examine here corresponds to a situation where the states involved can be split into resp. particles and anti-particles with a suitable representation of CP transformation acting on the states. Within the context of integrable field theories it is ustified to denote them respectively solitons S and antisolitons A. Assuming that the -body diffusion conserves the soliton or antisoliton nature of the particles, but that the reflection may change it, one should therefore consider four types of R matrices (resp. RSS SS, RSA AS, RAS SA, RAA ) connected by CP operations; and four reflection matrices (resp. KS S, KS A, KA S, KA A ). It then becomes possible to define several non-equivalent constructions of commuting transfer matrices. As a consequence, one sees that a variety of spin chain models can be built using a Sklyanin-like procedure, depending on which transfer matrix is being constructed and which reflection matrices are used to build it. Locality arguments also come into play, leading to more complicated combinations of transfer matrices as we shall presently see. We shall here describe the construction, and present the resolution by analytical Bethe Ansatz methods [16,7], of open spin chains based on the simplest rational R matrix solutions of the Yang- Baxter equations for underlying sl(n) and sl( N) Lie (super) algebras. These solutions, with a 1

3 rational dependence on the spectral parameter, are instrumental in defining the Yangian [8]. We shall consider two types of associated reflection matrices K to build two distinct types of integrable spin chains: one which entails purely soliton- and antisoliton-preserving reflection amplitude by two matrices KS S and KA A (hereafter denoted soliton-preserving case or SP); the other which entails the two soliton-non-preserving reflection amplitudes KS A and KS A (hereafter denoted soliton-nonpreserving case or SNP). Closed spin chains based on sl( N) superalgebras were studied in e.g. [9] and, in the case of alternating fundamental-conugate representations of sl( N) in [30]. Open spin chains based on sl(1 ) have been studied in details in e.g. [1,]. The plan of our presentation is as follows: We first define the relevant algebraic obects, R matrix and K matrix, together with their compatibility (Yang Baxter and reflection) equations and their relevant properties. In particular we introduce the various reflection equations which arise in the SP and SNP cases. Let us emphasise that all notions introduced in the sl(n) case will be straightforwardly generalised to the sl( N) case, albeit with a graded tensor product. In a second part we define the commuting transfer matrices which can be built in both cases, and the local Hamiltonians which can be built from them. Locality requirement leads to considering a product of two transfer matrices, resp. soliton-antisoliton and antisoliton-soliton, in the SNP case. Once again this construction will be valid, with suitable modifications, for the sl( N) case. In a third part we discuss the symmetries of these transfer matrices induced by their respective YBE and RE structures, in particular focusing on the connection between the SNP case and twisted Yangians. We then start the discussion of the analytical Bethe Ansatz formulation for the sl(n) spin chains in the SNP case. The derivation of suitable new fusion formulae explicited in Appendix A and B makes it possible to get a set of Bethe equations. In Section 5 we consider the case of sl( N) super algebra as underlying algebra. Contrary to the previous case it is first needed to establish a classification for the reflection matrices based on the rational (super Yangian) quantum R-matrix solution, both for SP and SNP conditions. We then establish the Bethe equations for both SP and SNP cases. In the SP case in addition we consider spin chains built from general K matrix solutions, in the SNP case we restrict ourselves to diagonal K matrices. 1 Yang Baxter and reflection equations The R and K matrices obey sets of coupled consistency equations together with characteristic properties which we now describe. 1.1 The R matrix We will consider in a first stage the sl(n) invariant R matrices R 1 (λ) = λi+ip 1 (1.1)

4 where P is the permutation operator P 1 = N E i E i. (1.) i, E i are the elementary matrices with 1 in position (i,) and 0 elsewhere. We define a transposition t which is related to the usual transposition T by (A is any matrix): A t = V 1 A T V where (1.3) V = antidiag(1,1,...,1), or ( ) for which V = θ 0 = 1 V = antidiag 1,...,1, 1,..., 1 }{{}}{{}, for which V = θ 0 = 1. N/ N/ The second case is forbidden for N odd. This R matrix satisfies the following properties: (i) Yang Baxter equation [1, 31 33] (ii) Unitarity R 1 (λ 1 λ ) R 13 (λ 1 ) R 3 (λ ) = R 3 (λ ) R 13 (λ 1 ) R 1 (λ 1 λ ) (1.4) where R 1 (λ) = P 1 R 1 (λ)p 1 = R t 1t 1 (λ) = R 1 (λ). (iii) Crossing-unitarity where ρ = N and It obeys R 1 (λ) R 1 ( λ) = ζ(λ) (1.5) R t 1 1 (λ) R t 1 ( λ iρ) = ζ(λ+iρ) (1.6) ζ(λ) = (λ+i)( λ+i), ζ(λ) = (λ+iρ)( λ+iρ). (1.7) [A 1 A, R 1 (λ)] = 0 for any matrix A. (1.8) The R matrix can be interpreted physically as a scattering matrix [4, 33, 34] describing the interaction between two solitons that carry the fundamental representation of sl(n). To take into account the existence, in the general case, of anti-solitons carrying the conugate representation of sl(n), we shall introduce another scattering matrix, which describes the interaction between a soliton and an anti-soliton. This matrix is derived as follows R 1(λ) = R 1 (λ) := R t 1 1 ( λ iρ) (1.9) = R t 1 ( λ iρ) =: R 1 (λ) = R 1 (λ) (1.10) 3

5 In the case N = and for θ 0 = 1 (sp() case), R is proportional to R, so that there is no genuine notion of anti-soliton. This reflects the fact that the fundamental representation of sp() = sl() is self-conugate. This does not contradict the fact that for N = and for θ 0 = +1 (so() case), there exists a notion of soliton and anti-soliton. The equality between R 1(λ) and R 1 (λ) in (1.9) reflects the CP invariance of R, from which one also has R 1 = R 1, i.e. the scattering matrix of anti-solitons is equal to the scattering matrix of solitons. In (1.9), R1 (λ) = ( λ iρ)i + iq 1, in which Q 1 is proportional to a proector onto a one-dimensional space. It satisfies The R matrix (1.10) also obeys (i) A Yang Baxter equation (ii) Unitarity (iii) Crossing-unitarity Q = ρ Q and P Q = Q P = θ 0 Q. (1.11) R 1 (λ 1 λ ) R 13 (λ 1 ) R 3 (λ ) = R 3 (λ ) R 13 (λ 1 ) R 1 (λ 1 λ ) (1.1) R 1 (λ) R 1 ( λ) = ζ(λ) (1.13) R t 1 1 (λ) R t 1 ( λ iρ) = ζ(λ+iρ). (1.14) Remark: The crossing-unitarity relation written in the literature usually involves a matrix = V T V. In our case, turns out to be 1 for two reasons: (i) the factors q k of the quantum (trigonometric) case degenerate to 1 in the Yangian (rational) case; (ii) the signs usually involved in the super case are in this paper (section 5) taken into account in the definition of the super-transposition (5.3). 1. The K matrix The second basic ingredient to construct the open spin chain is the K matrix. We shall describe in what follows two different types of boundary conditions, called soliton preserving (SP) [8 11] and soliton non preserving (SNP) [1 14] Soliton preserving reflection matrices In the case of soliton preserving boundary conditions, the matrix K is a numerical solution of the reflection (boundary Yang Baxter) equation [3] R ab (λ a λ b ) K a (λ a ) R ba (λ a +λ b ) K b (λ b ) = K b (λ b ) R ab (λ a +λ b ) K a (λ a ) R ba (λ a λ b ), (1.15) 4

6 and it describes the reflection of a soliton on the boundary, coming back as a soliton. Another reflection equation is required for what follows, in particular for the fusion procedure described in the appendices R ab (λ a λ b ) Kā(λ a ) R ba (λ a +λ b ) K b (λ b ) = K b (λ b ) R ab (λ a +λ b ) Kā(λ a ) R ba (λ a λ b ). (1.16) Kā is a solution of the anti-soliton reflection equation obtained from (1.15) by CP conugation and actually identical to (1.15) due to the CP invariance of the R-matrix. It describes the reflection of an anti-soliton on the boundary, coming back as an anti-soliton. Equation (1.16) appears as a criterion for a consistent choice of a couple of solutions K a and Kā of (1.15), yielding the commutation of transfer matrices (hereafter to be defined). Graphically, (1.15) and (1.16) are represented as follows: a a a b b a b b These K matrices (solutions of the soliton preserving reflection equation (1.15)) have been classified for sl(n) Yangians in [35]. This classification can be recovered as a particular case of our proposition 5.1, where the sl( N) Yangians are studied. Yang Baxter and reflection equations will indeed take the same form albeit with a graded tensor product in the superalgebraic case (see section 5 for more details). Proposition 5. then provides the classification of pairs {K a (λ),kā(λ)} which obey (1.15) and the compatibility equation (1.16). 1.. Soliton non-preserving reflection matrices In the context of soliton non-preserving boundary conditions one has to consider [1 14] the case where a soliton reflects back as an anti-soliton. The corresponding reflection equation has the form R ab (λ a λ b ) K a (λ a ) R ba (λ a +λ b ) K b (λ b ) = K b (λ b ) R ab (λ a +λ b ) K a (λ a ) R ba (λ a λ b ). (1.17) Note that equation (1.17) is satisfied by the generators of the so-called twisted Yangian [36,37], which will be discussed in section 3.. Similarly to the previous case, one introduces Kā, describing an anti-soliton that reflects back as a soliton, satisfying (1.17) and the consistency condition R ab (λ a λ b ) Kā(λ a ) R ba (λ a +λ b ) K b (λ b ) = K b (λ b ) R ab (λ a +λ b ) Kā(λ a ) R ba (λ a λ b ). (1.18) Graphically, (1.17) and (1.18) are represented as follows: 5

7 a a a b b a b b The K matrices corresponding to the soliton non-preserving reflection equation are classified in proposition 5.3, in the case of sl( N) Yangians, where once again similar Yang Baxter and reflection equations occur. Proposition 5.4 then provides the classification of pairs { K a (λ), Kā(λ)} which obey (1.17) and the compatibility equation (1.18). The transfer matrix We are now in a position to build open spin chains with different boundary conditions from the obects K, K, R and R [4]. Our purpose is to determine the spectrum and the symmetries of the transfer matrix for the case where soliton non-preserving boundary conditions are implemented. We first recall the general settings for the soliton preserving case..1 Soliton preserving case Let us first define the transfer matrix for the well-known boundary conditions, i.e. the soliton preserving ones. The starting point is the construction of the monodromy matrices T a (λ) = T a (λ) K a (λ) ˆT a (λ), (.1) T a (λ) = Tā(λ) K ā (λ) ˆTā(λ) (.) and two transfer matrices (soliton soliton and anti-soliton anti-soliton) with t(λ) = Tr a K + a (λ) T a(λ), t(λ) = Tr a K + ā (λ) T a(λ), (.3) T a (λ)=r al (λ)...r a1 (λ), ˆTa (λ)=r 1a (λ)...r La (λ), Tā(λ)= R al (λ)... R a1 (λ), ˆTā(λ)= R 1a (λ)... R La (λ). (.4) The numerical matrices K a (λ), K ā (λ) are solutions of (1.15), (1.16) and K + a satisfies a reflection equation dual to (1.15), R ab ( λ a +λ b ) K + a (λ a) t R ba ( λ a λ b iρ) K + b (λ b) t = K + b (λ b) t R ab ( λ a λ b iρ) K + a (λ a) t R ba ( λ a +λ b ). (.5) 6

8 The solutions of (.5) take the form K + a (λ) = K t a( λ iρ), where K a (λ) is a solution of (1.15). Similarly K + ā (λ) satisfies a reflection equation dual to (1.15), the solutions of which being of the form K + ā (λ) = K t ā ( λ iρ). In addition, K+ ā satisfies also a compatibility condition dual to (1.16). Actually, the dual reflection equations happen to be the usual reflection equations after a redefinition λ c λ c iρ. From their explicit expression, one can deduce that the monodromy matrices T (λ) and T(λ) obey the following equations: R ab (λ a λ b ) T a (λ a ) R ba (λ a +λ b ) T b (λ b ) = T b (λ b ) R ab (λ a +λ b ) T a (λ a ) R ba (λ a λ b ) (.6) R ab (λ a λ b ) T a (λ a ) R ba (λ a +λ b ) T b (λ b ) = T b (λ b ) R ab (λ a +λ b ) T a (λ a ) R ba (λ a λ b ) (.7) R ab (λ a λ b ) T a (λ a ) R ba (λ a +λ b ) T b (λ b ) = T b (λ b ) R ab (λ a +λ b ) T a (λ a ) R ba (λ a λ b ), (.8) which ust correspond to the soliton preserving reflection equations (1.15) and (1.16). One can recognise, in the relation (.6), the exchange relation of reflection algebras based on R-matrix of Y(sl N ). The form of K ± determines the precise algebraic structure which is involved, see section 3. As usual in the framework of spin chain models, the commutativity of the transfer matrices (.3) [t(λ),t(µ)] = 0 (.9) [ t(λ), t(µ)] = 0 [t(λ), t(µ)] = 0 (.10) is ensured by the above exchange relations. The first commutator guarantees the integrability of the model, whose Hamiltonian is given by H = 1 d dλ t(λ) λ=0, (.11) the locality being ensured because R(0) = P. The commutators (.10) will be needed so that the fusion procedure be well-defined (see appendices). The transfer matrix for K + (u) = K (u) = 1 satisfies a crossing-like relation (see for details [14,7]) t(λ) = t( λ iρ). (.1) The eigenvalues of the transfer matrices as well as the corresponding Bethe Ansatz equations have been derived for diagonal K matrices in [8,16].. Soliton non-preserving case This case was studied in [14] for the sl(3) chain only. Here we generalise the results for any sl(n). One introduces two monodromy matrices T a (λ) = T a (λ) K a (λ) ˆTā(λ), T a (λ) = Tā(λ) K ā (λ) ˆT a (λ), (.13) 7

9 and two transfer matrices (anti-soliton soliton and soliton anti-soliton) defined by where now t(λ) = Tr a K+ a (λ) T a (λ), t(λ) = Tr a K+ā (λ) T a (λ), (.14) T a (λ)=r al (λ) R al 1 (λ)...r a (λ) R a1 (λ), ˆTā(λ)=R 1a (λ) R a (λ)...r L 1a (λ) R La (λ), Tā(λ)= R al (λ)r al 1 (λ)... R a (λ)r a1 (λ), ˆTa (λ)= R 1a (λ)r a (λ)... R L 1a (λ)r La (λ). (.15) Note that, in this case, the number of sites is L because we want to build an alternating spin chain, which is going to ensure that the Hamiltonian of the model is local. This construction is similar to the one introduced in [38], where however a different notion of R was used. The numerical matrices K a, K ā are solutions of (1.17), (1.18). The numerical matrices K + a and K + ā are solutions of the following reflection equations: The commutators R ab ( λ a +λ b ) K + a (λ a) t Rba ( λ a λ b iρ) K + b (λ b) t = K + b (λ b) t Rab ( λ a λ b iρ) K + a (λ a ) t R ba ( λ a +λ b ) (.16) R ab ( λ a +λ b ) K + ā (λ a) t Rba ( λ a λ b iρ) K + b (λ b ) t = K + b (λ b ) t Rab ( λ a λ b iρ) K + ā (λ a) t R ba ( λ a +λ b ) (.17) R ab ( λ a +λ b ) K + ā (λ a ) t R ba ( λ a λ b iρ) K + b (λ b) t = K + b (λ b) t R ab ( λ a λ b iρ) K + ā (λ a) t Rba ( λ a +λ b ) (.18) [t(λ),t(µ)] = 0, [ t(λ), t(µ)] = 0 and [t(λ), t(µ)] = 0 (.19) are ensured by the above exchange relations for K + and the relations for the monodromy matrices, namely R ab (λ a λ b ) T a (λ a ) R ba (λ a +λ b ) T b (λ b ) = T b (λ b ) R ab (λ a +λ b ) T a (λ a ) R ba (λ a λ b ), (.0) R ab (λ a λ b ) T a (λ a ) R ba (λ a +λ b ) T b (λ b ) = T b (λ b ) R ab (λ a +λ b ) T a (λ a ) R ba (λ a λ b ), (.1) R ab (λ a λ b ) T a (λ a ) R ba (λ a +λ b ) T b (λ b ) = T b (λ b ) R ab (λ a +λ b ) T a (λ a ) R ba (λ a λ b ). (.) The relation (.0) has to be compared with the exchange relation of twisted Yangians based on R-matrix of Y(sl N ). We will come back to this point in section 3. In the SNP case, one can show [14] that the transfer matrices for K + (u) = K (u) = 1 exhibit a crossing symmetry, namely t(λ) = t( λ iρ), t(λ) = t( λ iρ). (.3) Starting from (.15) and (.14), the Hamiltonian of the alternating open spin chain is derived as [14, 17] H = 1 d dλ t(λ) t(λ). (.4) λ=0 8

10 The locality is again ensured by R(0) = P, and integrability is guaranteed by (.19). We can explicitly write the open Hamiltonian in terms of the permutation operator and the R(0) and R(0) matrices. Let us first introduce some notations: R i = R i (0), R i = d dλ R i (λ) λ=0 and Ř i = P i d dλ R i(λ). (.5) λ=0 After some algebraic manipulations, in particular taking into account that Tr 0 P 0i R0i I, we obtain the following expression for the Hamiltonian (.4) (for a detailed proof see [14]) H + + L L 1 L 1 R 1 R 1 + L 1 R +1 + Ř + R +1 + R +1 + R 1 R 1 + R 1 + R 1 R+1 + R +1 + R 1 R 1 + Ř 1 +1 R 1 + R 1 R Tr 0 Ř 0 L R L 1 L P 0 L 1 R0 L 1 RL 1 L + R 1 Ř 1 R 1, (.6) which is indeed local including terms that describe interaction up to four first neighbours. It is easily shown, acting on (.15) by full transposition, that t(λ) t t 1...t L (λ) provided that ( K a )t K ā. Eigenvectors of H in (.6) are determined by sole evaluation of eigenvectors of t(λ). We shall therefore only need to consider diagonalisation of t(λ) in what follows. 3 Symmetry of the transfer matrix In the two (SP and SNP) boundary cases, the use of exchange relations for the monodromy matrices allows us to determine the symmetry of the transfer matrix. For simplicity, we fix K + (λ) (or K + (λ)) to be I, leaving K (λ) (or K (λ)) free. 3.1 Soliton preserving boundary conditions In this case, the general form for K (λ) [35] is conugated (through a constant matrix) to the following diagonal matrix: K (λ) = iξi+λe with E = diag(+1,...,+1, 1,..., 1). (3.1) }{{}}{{} m n In the particular case of diagonal solutions, one recovers the scaling limits of the solutions obtained in the quantum case in [39]. The monodromy matrix then T(λ) generates a B(N,n) reflection algebra as studied in [40]. Taking the trace in space a of the relation (.6), we obtain: (λ a λ b )[t(λ a), T(λ b )] = (iλ a N)[T(λ a ), T(λ b )] (3.) 9

11 From the asymptotic behaviour (λ ) of the R matrix we deduce T (λ) = λ L+1 ( ( R 0i (λ) = λ I+ i ) λ P 0i, (3.3) E+ i λ ( ξ + ) L ) B 0 +O( 1 λ ) where B 0 = P 0 E 0 +E 0 P 0. We canwrite B 0 = N α,β=1 E αβ b αβ, with E αβ theelementary matrices acting in space 0, and b αβ realizing the gl(m) gl(n) algebra (in the space ). Picking up the coefficient of λ L b in the relation (3.), one then concludes that [ ] t(λ a ), i.e. the transfer matrix commutes with the gl(m) gl(n) algebra. 3. Soliton non-preserving boundary conditions L b αβ (3.4) = 0, (3.5) We have already mentioned that the monodromy matrix T(λ) satisfies one of the defining relations (.0) for the twisted Yangians Y ± (N) [36,37]. The general form for K (λ) is given in proposition 5.3: it is constant and obeys ( K ) t = ǫ K with ǫ = ±1. We need to investigate the asymptotic behaviour (λ ) of the R and R matrices given by (1.1) and (1.10) respectively: ( R 0i (λ) = λ I+ i ) λ P 0i and R0i (λ) = λ (I+ i λ P ) 0i, (3.6) where we have introduced P 0i = ρi Q 0i. (3.7) Accordingly, the monodromy matrices (.15) take the following form and finally T 0 (λ) = ( λ ) L ( I+ i λ ˆT 0(λ) = ( λ ) L ( I+ i λ L i=1 (P 0,i + P 0,i 1 )+O( 1 λ ) ), L ( P0,i +P 0,i 1 )+O( 1 ) λ ) i=1 (3.8) where T 0 (λ) K ˆT 0(λ) = λ 4L ( K + i λ L i=1 S 0i +O( 1 λ ) ), (3.9) S 0,i = P 0,i K + K P0,i, S 0,i 1 = K P 0,i 1 + P 0,i 1 K. (3.10) 10

12 Similarly to the previous case, one can show from equation (.0) that (this time for K = 1 only): [ t(λ a ), L i=1 S 0i ] = 0 (3.11) In this particular case S realises the so(n) (or sp(n)) generators. When N is even, as already mentioned, there are two possibilities for the proector Q. ore specifically the choice θ 0 = 1 in (1.3) corresponds to the so(n) case, whereas θ 0 = 1 corresponds to the sp(n) case. Remark: The same construction (open chain with twisted boundary conditions) can be done starting from the so(n), sp(n) and osp( N) R-matrix [41]. However, since the fundamental representations of these algebras are self-conugated, solitons and anti-solitons define the same obect. Hence, the twisted boundary conditions should be equivalent to open chains with ordinary boundary conditions. Indeed, it has been shown in [41] that boundary reflection equations (defining the boundary algebra) and twisted reflection equations (defining the twisted Yangian) are identical. 4 Spectrum of the transfer matrix Our purpose is to determine the spectrum of the transfer matrix for the sl(n) case. 4.1 Treatment of non-diagonal reflection matrices (SP case) In the soliton preserving case, the classification of reflection matrices associated to the Yangian Y(gl(N)) has been computed in [35]. It can been recovered from the sl( N) case given in proposition 5.1. Using this classification, it is easy to show the following proposition [4]: Proposition 4.1 Let K(λ) be any diagonalizable reflection matrix. D(λ), the corresponding diagonal reflection matrix, can be written as D(λ) = U 1 K(λ)U (4.1) where U, the diagonalization matrix, is constant. Let t K (λ) = Tr a (T a (λ)k a (λ)ˆt a (λ)) and t D (λ) = Tr a (T a (λ)d a (λ)ˆt a (λ)) be the corresponding transfer matrices (we set K + = I). Then, t K (λ) and t D (λ) have the same eigenvalues, their eigenvectors (say v K and v D respectively) being related through v K = U 1 U...U L v D (4.) Proof: The fact that the diagonalization matrix is a constant (in λ) is a consequence of the classification (see [35] and proposition 5.1). Using the property (1.8), one can show that t K (λ) = U 1 U...U L t D (λ)(u 1 U...U L ) 1, (4.3) which is enough to end the proof. The general treatment (including the super case) for diagonal reflection matrices is done in section 11

13 From the above property, this treatment (for diagonal matrices) is enough to obtain the spectrum for all the transfer matrices associated to all the reflection matrices (provided they are diagonalisable). As an illustration of this proposition, we compute the eigenvalues associated to the non-diagonal reflection matrix [3, 10] λ+iξ 0 0 kλ. 0 cλ+iξ.. 0 K(λ) = with c = 4k +1 (4.4). 0.. cλ+iξ 0 kλ 0 0 λ+iξ It is easy to see that K(λ) is diagonalized by the constant matrix The corresponding diagonal matrix is given by k k 0 0 ξ ξ U = (4.5) c 1 c c c D(λ) = c diag( λ+iξ,λ+iξ,...,λ+iξ ), with ξ = ξ c, (4.6) in accordance with the classification of reflection matrices. The t K (λ) eigenvalues as well as the Bethe equations are identical to the ones of t D (λ). They can be deduced from the general treatment given in section 5.4.4, taking formally n 1 = n = N = 0 and specifying m 1 = 1. They can also be viewed as the scaling limit of quantum groups diagonal solutions [15]. The eigenvectors are related using the formula (4.), with the explicit form (4.5) for U. This general procedure can be applied for an arbitrary spin chain, provided the diagonalization matrix is independent from the spectral parameter and commutes, see eq. (1.8), with the R-matrix under consideration. When K + (λ) is not I, this technics can also be used if K + (λ) and K (λ) can be diagonalized in the same basis [43]. Note however that the classification does not ensure the full generality of such an assumption. 4. Pseudo-vacuum and dressing functions We present below the case when soliton non-preserving boundary conditions are implemented with the simplest choice K ± = I. Results for more general choices of K can be deduced from the superalgebra case treated in section

14 We first derive the pseudo-vacuum eigenvalue denoted as Λ 0 (λ), with the pseudo-vacuum being 1 L ω + = + i where + = 0. CN. (4.7) i=1 0 Note that the pseudo-vacuum is an exact eigenstate of the transfer matrix (.3), and Λ 0 (λ) is given by the following expression with N Λ 0 (λ) = (a(λ) b(λ)) L g 0 (λ)+(b(λ) b(λ)) L l=1 g l (λ)+(ā(λ)b(λ)) L g N 1 (λ). (4.8) a(λ) = λ+i, b(λ) = λ, ā(λ) = a( λ iρ), b(λ) = b( λ iρ) (4.9) and g l (λ) = λ+ i(ρ θ 0) λ+ iρ, 0 l < N 1 gn 1(λ) = 1, for N odd g l (λ) = g N l 1 ( λ iρ). (4.10) We remind that ρ = N. We make at this point the assumption that any eigenvalue of the transfer matrix can be written as N Λ(λ) = (a(λ) b(λ)) L g 0 (λ)a 0 (λ)+(b(λ) b(λ)) L l=1 g l (λ)a l (λ)+(ā(λ)b(λ)) L g N 1 (λ)a N 1 (λ) (4.11) where the so-called dressing functions A i (λ) need now to be determined. We immediately get from the crossing symmetry (.3) of the transfer matrix: A 0 (λ) = A N 1 ( λ iρ), A l (λ) = A N l 1 ( λ iρ). (4.1) oreover, we obtain from the fusion relation (A.10) the following identity, by a comparison of the forms (4.11) for the initial and fused auxiliary spaces: Gathering the above two equations (4.1), (4.13) we conclude A 0 (λ+iρ)a N 1 (λ) = 1. (4.13) A 0 (λ)a 0 ( λ) = 1. (4.14) Finally from equations (B.13) important relations between the dressing functions are entailed for both soliton preserving and soliton non-preserving boundary conditions. In particular, N 1 l=0 A l (λ+i(n 1) il) = 1. (4.15) 13

15 Taking into account the constraints (4.1), (4.14) and (4.15) one derives the dressing functions: A 0 (λ) = A l (λ) = (1) (l) (l+1) λ+λ (1) i λ+λ (1) + i λ+λ (l) + il +i λ+λ (l) + il λ+λ (l+1) λ+λ (l+1) λ λ (1) i, λ λ (1) + i λ λ (l) + il +i λ λ (l) + il + il i + il + i λ λ (l+1) λ λ (l+1) + il i + il + i, 1 l < N 1 (4.16) together with the property A l (λ) = A N 1 l ( λ iρ), and, for N = n+1: A n (λ) = (n) λ+λ (n) + in +i λ+λ (n) + in λ λ (n) + in +i λ λ (n) + in λ+λ (n) + in i λ+λ (n) + in + i λ λ (n) + in i λ λ (n) + in +, (4.17) i Note that the dressing does not depend on the value of θ 0. The numbers (l) in (4.17) are related as customary to the eigenvalues of diagonal generators S l of the underlying symmetry algebra (determined in the previous section), namely S 1 = 1 (0) (1), S l = (l 1) (l) with S l = 1 (E ll E l l), 1 l < N 1 (4.18) with (0) = L. Recallthatforthesl(N)casethecorresponding numbers (l), seee.g. [11], aregivenbythefollowing expressions E ll = (l 1) (l), l = 1,...,N (4.19) with (0) = L and (N) = 0. If we now impose (l) = (N l) and consider the differences E ll E l l, we end up with relations (4.18), in accordance with the folding of sl(n) leading to so(n) and sp(n) algebras. 4.3 Bethe Ansatz equations From analyticity requirements one obtains the Bethe Ansatz equations which read as: 14

16 4.3.1 sl(n + 1) algebra (1) e 1 (λ (1) i ) L = e (λ (1) i λ (1) ) e (λ (1) i +λ (1) ) () (l) 1 = e (λ (l) i λ (l) ) e (λ (l) i +λ (l) ) l =,...,n 1, τ=±1 e 1 (λ (1) i λ () ) e 1 (λ (1) i +λ () ), (l+τ) e 1 (λ (l) i λ (l+τ) ) e 1 (λ (l) (n) e 1(λ (n) i ) = e 1 (λ (n) i λ (n) ) e 1 (λ (n) i +λ (n) )e (λ (n) i λ (n) ) e (λ (n) i +λ (n) ) (n 1) where we have introduced e 1 (λ (n) i i +λ (l+τ) ) λ (n 1) ) e 1 (λ (n) i +λ (n 1) ), (4.0) e x (λ) = λ+ ix λ ix. It is interesting to note that equations (4.0) are exactly the Bethe Ansatz equations of the osp(1 N 1) case (see e.g. [14,19]) sl(n) algebra (1) e 1 (λ (1) i ) L = (l) 1 = (n) e θ0 (λ (n) i ) = e (λ (1) i λ (1) ) e (λ (1) i +λ (1) ) () e (λ (l) i λ (l) ) e (λ (l) i +λ (l) ) l =,...,n 1, (n 1) e (λ (n) i e 1 (λ(n) i λ (n) ) e (λ (n) i +λ (n) ) τ=±1 e 1 (λ (1) i λ () ) e 1 (λ (1) i +λ () ), (l+τ) e 1 (λ (l) i λ (l+τ) ) e 1 (λ (l) i +λ (l+τ) ) λ (n 1) ) e 1 (λ(n) i +λ (n 1) ). (4.1) The Bethe Ansatz equations are essentially the same as the ones obtained from the folding of the usual sl(n) Bethe equations (see e.g. [8,11]) for (l) = (N l). It can be realised from the study of the underlying symmetry of the model that this folding has algebraic origins, as mentioned previously. 15

17 5 sl( N) superalgebra In this section, we generalise the previous approach on (SP and SNP) boundary conditions to the Z -graded case based on the sl( N) superalgebra. 5.1 Notations The Z -gradation, depending on a sign θ 0 = ±, is defined to be ( 1) [i] = θ 0 for i an sl() index and ( 1) [i] = θ 0 an sl(n) index. The sl( N) invariant R matrix reads R 1 (λ) = λi+ip 1, (5.1) where P is from now on the super-permutation operator (i.e. X 1 PX 1 P) such that P = +N i, ( 1) [] E i E i (5.) The usual super-transposition T is defined for any matrix A = i Ai E i, by A T = i ( 1) [i][]+[] A i E i = i ( A T ) i Ei. (5.3) As for the sl(n) case, we will use a super-transposition t of the form A t = V 1 A T V. (5.4) The convention for θ 0 andthe expression of V are chosen accordingly to the selected Dynkin diagram. Let us recall that for a basic Lie superalgebra, unlike the Lie algebraic case, there exist in general many inequivalent simple root systems (i.e. that are not related by a usual Weyl transformation), and hence many inequivalent Dynkin diagrams. This situation occurs when a simple root system contains at least one isotropic fermionic root. For each basic Lie superalgebra, there is a particular Dynkin diagram which can be considered as canonical: it contains exactly one fermionic root. Such a Dynkin diagram is called distinguished. In the case of sl( N), it has the following form: }{{} N 1 }{{} 1 In the case of sl( n) superalgebras, there exists a symmetric Dynkin diagram with two isotropic fermionic simple roots in positions n and +n: }{{} n 1 }{{} 1 }{{} n 1 16

18 The generalization of the Weyl group for a basic Lie superalgebra gives a method for constructing all the inequivalent simple root systems and hence all the inequivalent Dynkin diagrams. For more details, see [44 47]. (i) Distinguished Dynkin diagram basis In this case, we consider that the sl() part occupies the upper part of the matrices and corresponds to bosonic degrees of freedom, whereas the sl(n) part occupies the lower part and corresponds to fermionic degrees. ore precisely, the gradation takes the form: { ( 1) [i] 1 for 1 i = (5.5) 1 for +1 i +N, and the matrix V reads ( V 0 V = 0 V N ). (5.6) In the above formula, V (resp. V N ) is the (resp. N N) matrix given in (1.3) for θ 0 = +1. (ii) Symmetric Dynkin diagram basis As in the sl(n) case, one has to take or N even. Note that for the odd odd case no symmetric Dynkin diagram exists and consequently no twisted super-yangian [48]. Here, we choose N to be even: N = n. The sl() part lies in the middle part of the matrices and corresponds to fermionic degrees of freedom, whereas the sl(n) part occupies the upper and lower part and is associated to bosonic degrees of freedom. Correlatively, θ 0 = 1 in this case. The gradation is given by { ( 1) [i] 1 for 1 i n and +n+1 i +N = (5.7) 1 for n+1 i +n, while ( V = antidiag 1,...,1 }{{} n+ ), 1,..., 1 }{{} n. (5.8) We will mostly use the distinguished Dynkin diagram basis in the soliton preserving case, and the symmetric one in the soliton non-preserving case. In both cases, the R-matrix obeys the properties stated in section 1.1, with R 1 (λ) = R t 1 1 ( λ iρ) and ρ = θ 0 ( N). The K-matrices will obey the defining relations stated in section 1., and the properties of the transfer matrix (section ) also hold; the tensor product is now Z -graded. 5. Classification of reflection matrices for Y( N) This section is devoted to the classification of reflection matrices for the super-yangian Y( N) based on sl( N), both for soliton preserving (prop. 5.1 and 5.) and for soliton non-preserving boundary conditions (prop. 5.3 and 5.4). 17

19 5..1 Soliton preserving reflection Proposition 5.1 Any bosonic invertible solution of the soliton preserving reflection equation (RE) R 1 (λ 1 λ )K 1 (λ 1 )R 1 (λ 1 +λ )K (λ ) = K (λ )R 1 (λ 1 +λ )K 1 (λ 1 )R 1 (λ 1 λ ) where R 1 (λ) = λi+ip 1 is the super-yangian R-matrix, takes the form K(λ) = U (iξi+λe)u 1 where either (i) E is diagonal and E = I (diagonalisable solutions) (ii) E is strictly triangular and E = 0 (non-diagonalisable solutions) and U is an element of the group GL() GL(N). The classification is done up to multiplication by a function of the spectral parameter. Proof: Firstly it is obvious that for any solution K(λ) to the RE, and for any function f(λ), the product f(λ)k(λ) is also a solution to the RE, so that the classification will be done up to multiplication by a function of λ. Expanding the reflection equation, one rewrites it as: [K (λ 1 ), K (λ )] = i(λ 1 +λ )(K (λ 1 )K 1 (λ ) K (λ )K 1 (λ 1 )) +i(λ 1 λ )(K 1 (λ 1 )K 1 (λ ) K (λ )K (λ 1 )) (5.9) One then considers the RE with λ 1 and λ exchanged, and sums these two. After multiplication by P 1, one gets (for λ 1 λ ): [K 1 (λ 1 ), K 1 (λ )] = [K (λ 1 ), K (λ )] (5.10) the only solution of which is [K(λ 1 ), K(λ )] = 0. In other words, the matrices K(λ) at different values of λ s are diagonalisable (or triangularisable) in the same basis and must satisfy (λ 1 +λ )(K (λ 1 )K 1 (λ ) K (λ )K 1 (λ 1 ))+(λ 1 λ )(K 1 (λ 1 )K 1 (λ ) K (λ )K (λ 1 )) = 0 (5.11) By setting λ = λ 1, we get K(λ)K( λ) = k(λ)i for some function k(λ). If one now considers the case of invertible matrices, and since we are looking for solutions up to a multiplicative function, we can take K(λ)K( λ) = I, a condition which is generally assumed for reflection matrices. We first consider the case where these matrices can be diagonalised: K(λ) = U D(λ)U 1, where U is a group element of GL(m) GL(n). Proecting the RE on the basis element E ii E, one gets ( ) ( ) (λ 1 +λ ) d (λ 1 )d i (λ ) d (λ )d i (λ 1 ) = (λ 1 λ ) d (λ )d (λ 1 ) d i (λ 1 )d i (λ ) (5.1) where D(λ) = diag(d 1 (λ),d (λ),...,d m+n (λ)). Since K(λ) is supposed invertible, all the d s are not zero, and we consider q i (λ) = d i(λ) (5.13) d (λ) 18

20 which obeys ( ) ( ) (x+y) q(y) q(x) +(x y) q(x)q(y) 1 = 0. (5.14) The solution to this equation is q(x) = x+iξ x iξ so that, considering q 1 (λ), we get where ξ is some complex parameter (including ξ = ), d (λ) = λ+iξ λ iξ d 1 (λ), (5.15) Requiring q i (λ) to obey the equation (5.14) shows that one must have d (λ) = ǫ λ+iξ with ǫ = ±1, (5.16) where we have used the invariance under multiplication by a function. This yields the form (i), with E = diag(ǫ 1,...,ǫ m+n ). We now turn to the case K(λ) = U T(λ)U 1 where T(λ) is triangular. The proection of the RE on E ii E shows that the diagonal part of T(λ) is still of the form (i). We distinguish two cases: E has two different eigenvalues (which are ±1), or E is proportional to I (and then the diagonal of T(λ) is also proportional to I). If E has two eigenvalues, we proect, in the first auxiliary space, on two diagonal elements E and E kk associated to these eigenvalues: ( ) ( ) (λ 1 +λ ) (iξ ±λ )T(λ 1 ) (iξ ±λ 1 )T(λ ) = (λ 1 λ ) T(λ 1 )T(λ ) (iξ ±λ 1 )(iξ ±λ ) (5.17) The difference and the sum of these equations read: From equation (5.18), one gets λ T(λ 1 ) λ 1 T(λ ) = λ λ 1, λ 1 +λ 0 (5.18) ( ) ) iξ(λ 1 +λ ) T(λ 1 ) T(λ ) = (λ 1 λ ) (T(λ 1 )T(λ )+ξ λ 1 λ (5.19) T(λ 1 ) I λ 1 = T(λ ) I λ = T 0, i.e. T(λ) = I+λT 0 (5.0) where T 0 is a triangular matrix. Plugging this solution in equation (5.19), we obtain (iξ 1)(λ 1 +λ )T 0 = λ 1 λ (T 0 I)+(ξ +1)I, λ 1,λ (5.1) whose only (constant) solution is of the form (i) with iξ = 1. We are thus left with the case where the diagonal of T(λ) is proportional to the identity matrix: T(λ) = I + S(λ) with S(λ) strictly triangular. Proecting once more on a diagonal element in the first auxiliary space, we obtain ( ) λ S(λ 1 ) λ 1 S(λ ) = (λ 1 λ )S(λ 1 )S(λ ) (5.) S(λ 1 ) λ 1 (I+S(λ 1 )) = S(λ ) λ (I+S(λ )) = σ (5.3) 19

21 where σ is strictly triangular. We therefore have T(λ) = I + λσ(i λσ) 1. With this form for T(λ), the RE rewrites (σ 1 σ )σ 1 σ = 0, whose solution (for strictly triangular matrices) is given by σ = 0. Using this property, we get the solution (ii). Note that the solutions given in this proposition are all of the form K(λ) = iξi+λe with E = I or E = 0 (taking E = UEU 1 ). Proposition 5. Given a solution K(λ) = iξi+λe to the soliton preserving RE R 1 (λ 1 λ ) K 1 (λ 1 ) R 1 (λ 1 +λ ) K (λ ) = K (λ ) R 1 (λ 1 +λ ) K 1 (λ 1 ) R 1 (λ 1 λ ), (5.4) and a solution K(λ) = iξ I+λE to the anti-soliton preserving RE identical to (5.4), the compatibility condition R 1 (λ 1 λ ) K 1(λ 1 ) R 1 (λ 1 +λ ) K (λ ) = K (λ ) R 1 (λ 1 +λ ) K 1(λ 1 ) R 1 (λ 1 λ ) (5.5) is solved by E = E t and ξ +ξ = θ 0 N Str E. Proof: Straightforwardly, equation (5.5) is equivalent to E t Q 1 E = E Q 1 E t (5.6) (ξ +ξ )[E,Q 1 ] = [E,Q 1 E Q 1 ] (5.7) The first equation yields E = E t. Using Q 1 E Q 1 = θ 0 N Q 1 StrE one gets the relation between ξ and ξ. 5.. Soliton non-preserving reflection Proposition 5.3 Any bosonic invertible solution of the soliton non-preserving RE R 1 (λ 1 λ ) K 1 (λ 1 ) R 1 (λ 1 +λ ) K (λ ) = K (λ ) R 1 (λ 1 +λ ) K 1 (λ 1 ) R 1 (λ 1 λ ) (5.8) where R 1 (λ) = λi+ip 1 is the super-yangian R-matrix, is a constant (up to a multiplication by a scalar function) matrix such that K t = ± K. Proof: Writing the R and R matrices in terms of I, P 1 and Q 1, and taking the part of (5.8) which is symmetric in the exchange of λ 1 and λ, yields the following equation K 1 (λ 1 )Q 1 K1 (λ )+ K 1 (λ )Q 1 K1 (λ 1 ) = K (λ 1 )Q 1 K (λ )+ K (λ )Q 1 K (λ 1 ) (5.9) In the same way, exchanging the role of spaces 1 and from the original equation, one gets K 1 (λ 1 )Q 1 K (λ )+ K (λ 1 )Q 1 K1 (λ ) = K (λ )Q 1 K1 (λ 1 )+ K 1 (λ )Q 1 K (λ 1 ) (5.30) Transposing both equations (5.9) and (5.30) in space 1 and eliminating P 1, one gets after some algebra K t (λ ) = f(λ 1,λ ) K(λ 1 ) λ 1,λ (5.31) from which the final result follows. 0

22 Proposition 5.4 Given a solution K 1 to the soliton non-preserving RE R 1 (λ 1 λ ) K 1 (λ 1 ) R 1 (λ 1 +λ ) K (λ ) = K (λ ) R 1 (λ 1 +λ ) K 1 (λ 1 ) R 1 (λ 1 λ ), (5.3) and a solution K 1 to the CP-conugate RE identical to (5.3), the compatibility condition R 1 (λ 1 λ ) K 1(λ 1 ) R 1 (λ 1 +λ ) K (λ ) = K (λ ) R 1 (λ 1 +λ ) K 1(λ 1 ) R 1 (λ 1 λ ) (5.33) ) 1. is solved by K 1 ( K1 Proof: Straightforward. 5.3 Pseudo-vacuum and its dressing We can determine explicitly the eigenvalue Λ 0 (λ) of the transfer matrix (defined as in section ) acting on the pseudo-vacuum ω + (which is always bosonic in our conventions). We take here K ± = I (resp. K± = I), whilst cases with non trivial K ± (resp. K± ) are studied in section (resp ). Λ 0 (λ) is given by the following expression +N Λ 0 (λ) = α(λ) L g 0 (λ)+β(λ) L ( 1) [l+1] g l (λ)+γ(λ) L ( 1) [+N] g +N 1 (λ) (5.34) l=1 where for (using the notation given in (4.9)): (i) Soliton preserving boundary conditions with L sites (distinguished Dynkin diagram) α(λ) = a (λ), β(λ) = γ(λ) = b (λ) (5.35) and g l (λ) = g l (λ) = i( N) λ(λ+ ) (λ+ il i(l+1) )(λ+ ), l = 0,..., 1 λ(λ+ i( N) ), l =,...,+N 1 (5.36) (λ+ i( l 1) )(λ+ i( l) ) (ii) Soliton non preserving boundary conditions with L sites (symmetric Dynkin diagram) ( ), ( ), ( α(λ) = a(λ) b(λ) β(λ) = b(λ) b(λ) γ(λ) = ā(λ)b(λ)) (5.37) and g l (λ) = λ+ i (ρ+1) λ+ iρ, +N 1 0 l < g+n 1(λ) = 1, if +N odd g l (λ) = g N+ l 1 ( λ iρ). (5.38) We remind that θ 0 = 1 in that case. 1

23 From the exact expression for the pseudo-vacuum eigenvalue, we introduce the following assumption for the structure of the general eigenvalues: +N Λ(λ) = α(λ) L g 0 (λ)a 0 (λ)+β(λ) L ( 1) [l+1] g l (λ)a l (λ) l=1 + γ(λ) L ( 1) [+N 1] g +N 1 (λ)a +N 1 (λ) (5.39) where the dressing functions A i (λ) need to be determined. The basic constraints that they have to satisfy are the fusion and crossing equations as well as analyticity requirements. 5.4 sl( N) with soliton preserving boundary conditions From the analyticity of Λ(λ), one gets A l ( il ) = A l 1( il ), l = 1,..., 1, A l ( il ) = A l 1( il ), l = N +1,..., 1 (5.40) Gathering the constraints (.1), (A.10) and (5.40), one can determine the dressing functions, i.e. A 0 (λ) = A l (λ) = A l (λ) = (1) (l) (l+1) (l) λ+λ (1) i λ+λ (1) + i λ+λ (l) + il +i (l+1) λ+λ (l) + il λ+λ (l+1) λ+λ (l+1) λ λ (1) i λ λ (1) + i λ λ (l) + il +i λ λ (l) + il + il i + il + i λ+λ (l) +i il i λ+λ (l) +i il λ+λ (l+1) λ+λ (l+1) λ λ (l+1) λ λ (l+1) +i il + i +i il i + il i + il + i λ λ (l) +i il i λ λ (l) +i il λ λ (l+1) λ λ (l+1) l = 1,..., 1 +i il + i +i il i l =,...,+N 1 (5.41)

24 5.4.1 Bethe Ansatz equations for the distinguished Dynkin diagram From analyticity requirements one obtains the Bethe Ansatz equations, (1) e 1 (λ (1) i ) L = e (λ (1) i λ (1) ) e (λ (1) i +λ (1) ) () (l) 1 = e (λ (l) i λ (l) ) e (λ (l) i +λ (l) ) 1 = τ=±1 e 1 (λ (1) i λ () ) e 1 (λ (1) i +λ () ), (l+τ) l =,..., 1,+1,...,+N ( 1) (+1) e 1 (λ () i e 1 (λ () i (+N ) 1 = e 1 (λ (+N 1) i (+N 1) λ ( 1) ) e 1 (λ () i +λ ( 1) ) λ (+1) ) e 1 (λ () i +λ (+1) ) e (λ (+N 1) i e 1 (λ (l) i λ (l+τ) ) e 1 (λ (l) λ (+N ) ) e 1 (λ (+N 1) i +λ (+N ) ) i +λ (l+τ) ) λ (+N 1) ) e (λ (+N 1) i +λ (+N 1) ) (5.4) We recover here for = and N = 1 the Bethe Ansatz equation of the supersymmetric t J model which corresponds to the sl( 1) case [49] Bethe Ansatz equations for arbitrary Dynkin diagrams We wrote above only the dressing functions that correspond to the distinguished Dynkin diagram. It is however possible to construct the g and dressing functions for all the inequivalent Dynkin diagrams of sl( N). The inequivalent Dynkin diagrams of the sl( N) superalgebras contain only bosonic root of same square length ( white dots ), usually normalized to, and isotropic fermionic roots ( grey dots ). A given diagram is completely characterized by the p-uple of integers 0 < n 1 <... < n p < + N labelling the positions of the grey dots of the diagram. Formally, we define n 0 = 0 and n p+1 = +N although there is actually no root at these positions. Such a diagram defined by the p-uple (n i ) i=1...p corresponds to the superalgebra sl( N) with = i odd i p+1 n i i even i<p+1 n i and N = i even i p+1 n i i odd i<p+1 n i. (5.43) The g functions have a form similar to (5.36), with a change of increasing or decreasing behaviour of the poles each time a grey (fermionic) root is met. Indeed ( ) λ λ+ i( N) g l (λ) = ( λ+ i δ )( l λ+ i (δ ), l = 0,...,+N 1 (5.44) l +1) 3

25 where δ 0 = 0 whilst the δ l for l = 1,...,+N 1 are obtained by iteration δ l 1 if l = n i for some i δ l = δ l 1 +1 if n i < l < n i+1 for some i δ l 1 1 if n i 1 < l < n i for some i The Bethe Ansatz equations read, for l = 1,...,+N 1 and i = 1,..., (l) +N 1 (k) e (1 α l,α l ) e αl,α k (λ (l) i λ (k) ) e αl,α k (λ (l) i +λ (k) 1 (λ (1) i ) L l = 1 ) = 1 l 1 k=1 (5.45) (5.46) where α l,α k is the scalar product of the simple roots associated to the chosen Dynkin diagram Bethe Ansatz equations for the symmetric Dynkin diagram We give the useful example of the symmetric Dynkin diagram for which N is even, with the indices ordered as in (5.7). The g functions are in this case g l (λ) = g l (λ) = g l (λ) = i( N) λ(λ+ ) (λ+ il i(l+1) )(λ+ ), l = 0,...,N/ 1 λ(λ+ i( N) ) (λ+ i(n l 1) )(λ+ i(n l) ), l = N/,...,+N/ 1 λ(λ+ i( N) ), l = +N/,...,+N 1 (5.47) (λ+ i(l ) )(λ+ i(l +1) ) and it is straightforward to get the A i s. The Bethe Ansatz equations take the form: (1) e 1 (λ (1) i ) L = e (λ (1) i λ (1) ) e (λ (1) i +λ (1) ) () (l) 1 = e (λ (l) i λ (l) ) e (λ (l) i +λ (l) ) 1 = τ=±1 e 1 (λ (1) i λ () ) e 1 (λ (1) i +λ () ), (l+τ) l =,...,+N, l N, N + (l+1) e 1 (λ (l) i λ (l+1) l = N, N + (+N ) 1 = e 1 (λ (+N 1) i (+N 1) e (λ (+N 1) i ) e 1 (λ (l) i +λ (l+1) ) (l 1) e 1 (λ (l) e 1 (λ (l) i λ (l+τ) i λ (l 1) ) e 1 (λ (l) λ (+N ) ) e 1 (λ (+N 1) i +λ (+N ) ) i +λ (l+τ) ) ) e 1 (λ (l) i +λ (l 1) ) λ (+N 1) ) e (λ (+N 1) i +λ (+N 1) ) (5.48) 4

26 The indices of the e s in the products are the entries of the Cartan matrix corresponding to the chosen Dynkin diagram, in accordance with the results obtained in the closed chain case [0, 9, 50] Non trivial soliton preserving boundary conditions We come back to the distinguished Dynkin diagram basis and implement non trivial soliton preserving boundary conditions. From the classification given in section 5.., we know that K (λ) is always conugated (by a constant matrix U) to a diagonal matrix of the form K(λ) = diag(α,...,α,β,...,β,β,...,β,α,...,α) (5.49) }{{}}{{}}{{}}{{} m 1 m n n 1 As in the section 4.1, it is easy to see that the spectrum and the symmetry of the model depend only on the diagonal (5.49), and not on U. Indeed, when considering two reflection matrices related by a constant conugation, the corresponding transfer matrices are also conugated. Thus, property 5.1 ensures that it is enough to consider diagonal K(λ) matrices to get the general case. Such a property, which relies on the form of the R-matrix, is a priori valid only in the rational sl(n) and sl( N) cases. For a diagonal solution (5.49) with m 1 +m =, n 1 +n = N, α(λ) = λ+iξ, β(λ) = λ+iξ, and the free boundary parameter ξ, one can compute the new form g l (λ) of the g-functions entering the expression of Λ 0 (λ), the new pseudo-vacuum eigenvalue. They take the form: g l (λ) = ( λ+iξ)g l (λ), g l (λ) = (λ+iξ +im 1 )g l (λ), l = 0,...,m 1 1 l = m 1,...,+n 1 g l (λ) = ( λ+iξ im +in )g l (λ), l = +n,...,+n 1 (5.50) where g l (λ) are given by (5.36). The dressing functions (5.41) keep the same form, but the Bethe Ansatz equations (5.4) are modified (by K (λ)), so that the value of the eigenvalues Λ(λ) are different from the ones obtained when K(λ) = I. The modifications induced on Bethe Ansatz equations (5.4) are the following: The factor e 1 ξ+m 1 (λ) appears in the LHS of the m th 1 Bethe equation. The factor e 1 ξ+m 1 m n (λ) appears in the LHS of the (+n ) th Bethe equation. 5.5 sl( N) with soliton non-preserving boundary conditions From equations of the type (B.13) for the supersymmetric case relations between the dressing functions are entailed for both soliton preserving and soliton non-preserving boundary conditions. In particular, for the case that corresponds to the symmetric Dynkin diagram one obtains (N = n, while can be even or odd), n 1 n 1 A l (λ il) A +n+l (λ+i i(n 1)+il) = l=0 l=0 1 l=0 A n+l (λ i(n 1)+il). (5.51) In fact the latter equation is only necessary for the soliton non-preserving boundary conditions. 5

27 5.5.1 Dressing functions As already mentioned we consider here the R matrix, that corresponds to the symmetric Dynkin diagram. Note that the sl(n ) case is isomorphic to sl( n) and entails the same dressing functions and BAE. From the constraints (4.1), (4.14), (5.51), we conclude that the dressing functions take the form: A 0 (λ) = A l (λ) = A l (λ) = (1) (l) (l+1) (l) λ+λ (1) i λ+λ (1) + i λ+λ (l) + il +i (l+1) λ+λ (l) + il λ+λ (l+1) λ+λ (l+1) λ λ (1) i, λ λ (1) + i λ λ (l) + il +i λ λ (l) + il + il i + il + i λ+λ (l) +in il i λ+λ (l) +in il λ+λ (l+1) λ+λ (l+1) +in il + i +in il i λ λ (l+1) λ λ (l+1) + il i + il + i λ λ (l) +in il i λ λ (l) +in il λ λ (l+1) λ λ (l+1), l = 1,...,n 1 (5.5) +in il + i +in il i, n l < n+ 1 and A l (λ) = A +n 1 l ( λ iρ), and for = m+1 A k (λ) = (k) λ+λ (k) +in ik i λ+λ (k) +in ik λ+λ(k) +in ik + i λ+λ (k) +in ik i λ λ (k) +in ik i λ λ (k) +in ik λ λ (k) +in ik + i λ λ (k) +in ik, k = m+n (5.53) i Recall that the sl(m + 1 n + 1) case is not examined because there is no symmetric Dynkin diagram and consequently the folding of the algebra can not be implemented. In any case it is known [48] that the twisted super-yangian does not exist for sl(m+1 n+1) Bethe Ansatz equations From the analyticity requirements one obtains the Bethe Ansatz equations which read as: 6