Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary
We solve the XXZ Gaudin model with generic boundary using the modified algebraic Bethe ansatz. The diagonal and triangular cases have been recovered in this general framework. We show that the model for odd or even lengths has two different behaviors. The corresponding Bethe equations are computed f...
Gespeichert in:
Datum: | 2017 |
---|---|
1. Verfasser: | |
Format: | Artikel |
Sprache: | English |
Veröffentlicht: |
Інститут математики НАН України
2017
|
Schriftenreihe: | Symmetry, Integrability and Geometry: Methods and Applications |
Online Zugang: | http://dspace.nbuv.gov.ua/handle/123456789/149276 |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Zitieren: | Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary / N. Crampe // Symmetry, Integrability and Geometry: Methods and Applications. — 2017. — Т. 13. — Бібліогр.: 44 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-149276 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1492762019-02-20T01:24:36Z Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary Crampe, N. We solve the XXZ Gaudin model with generic boundary using the modified algebraic Bethe ansatz. The diagonal and triangular cases have been recovered in this general framework. We show that the model for odd or even lengths has two different behaviors. The corresponding Bethe equations are computed for all the cases. For the chain with even length, inhomogeneous Bethe equations are necessary. The higher spin Gaudin models with generic boundary is also treated. 2017 Article Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary / N. Crampe // Symmetry, Integrability and Geometry: Methods and Applications. — 2017. — Т. 13. — Бібліогр.: 44 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 81R12; 17B80; 37J35 DOI:10.3842/SIGMA.2017.094 http://dspace.nbuv.gov.ua/handle/123456789/149276 en Symmetry, Integrability and Geometry: Methods and Applications Інститут математики НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
description |
We solve the XXZ Gaudin model with generic boundary using the modified algebraic Bethe ansatz. The diagonal and triangular cases have been recovered in this general framework. We show that the model for odd or even lengths has two different behaviors. The corresponding Bethe equations are computed for all the cases. For the chain with even length, inhomogeneous Bethe equations are necessary. The higher spin Gaudin models with generic boundary is also treated. |
format |
Article |
author |
Crampe, N. |
spellingShingle |
Crampe, N. Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary Symmetry, Integrability and Geometry: Methods and Applications |
author_facet |
Crampe, N. |
author_sort |
Crampe, N. |
title |
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary |
title_short |
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary |
title_full |
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary |
title_fullStr |
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary |
title_full_unstemmed |
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary |
title_sort |
algebraic bethe ansatz for the xxz gaudin models with generic boundary |
publisher |
Інститут математики НАН України |
publishDate |
2017 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/149276 |
citation_txt |
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary / N. Crampe // Symmetry, Integrability and Geometry: Methods and Applications. — 2017. — Т. 13. — Бібліогр.: 44 назв. — англ. |
series |
Symmetry, Integrability and Geometry: Methods and Applications |
work_keys_str_mv |
AT crampen algebraicbetheansatzforthexxzgaudinmodelswithgenericboundary |
first_indexed |
2025-07-12T21:13:56Z |
last_indexed |
2025-07-12T21:13:56Z |
_version_ |
1837477233591058432 |
fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 13 (2017), 094, 13 pages
Algebraic Bethe Ansatz for the XXZ Gaudin Models
with Generic Boundary
Nicolas CRAMPE
Laboratoire Charles Coulomb (L2C), UMR 5221 CNRS-Université de Montpellier,
Montpellier, France
E-mail: nicolas.crampe@umontpellier.fr
Received November 01, 2017, in final form December 06, 2017; Published online December 13, 2017
https://doi.org/10.3842/SIGMA.2017.094
Abstract. We solve the XXZ Gaudin model with generic boundary using the modified
algebraic Bethe ansatz. The diagonal and triangular cases have been recovered in this
general framework. We show that the model for odd or even lengths has two different
behaviors. The corresponding Bethe equations are computed for all the cases. For the chain
with even length, inhomogeneous Bethe equations are necessary. The higher spin Gaudin
models with generic boundary is also treated.
Key words: integrability; algebraic Bethe ansatz; Gaudin models; Bethe equations
2010 Mathematics Subject Classification: 81R12; 17B80; 37J35
1 Introduction
The algebraic Bethe ansatz is a powerful method to solve analytically numerous integrable
models [41]. Then, the introduction of integrable boundaries in this framework was initiated
in [39]. However, the problem to solve models with generic integrable boundaries thanks to
the algebraic Bethe ansatz has been overcome only recently [1, 5, 6, 8, 15]. The problem lies
in the fact that there is no simple particular eigenvector, which is normally the starting point
of the method. Therefore, other approaches have been used to deal with this problem: the
functional Bethe ansatz [21, 23, 34, 35], the coordinate Bethe ansatz [16, 17], the separation of
variables [19, 22, 37], the q-Onsager approach [2, 4] and the matrix ansatz [18, 25, 28, 29, 30].
The algebraic Bethe ansatz has been also used in this context [7, 10, 33, 38, 42] but constraints
on the parameters defining the boundaries were necessary. Then, inhomogeneous T-Q relations
were introduced in [11, 12, 13, 36] where the authors obtained the eigenvalues and the Bethe
equations for the XXZ model with generic boundaries. Finally, the modified algebraic Bethe
ansatz allows us to compute the associated eigenvectors [1, 5, 6, 8, 44]. Let us mention also
that this method has been also exploited to study twisted XXX spin chain in [9] or the rational
Gaudin magnets in arbitrarily oriented magnetic fields in [20].
In this letter, we solve the Gaudin model (introduced in [24]) with generic boundaries using
the modified algebraic Bethe ansatz. The Gaudin model is one of the simplest quantum inte-
grable system and we hope that its resolution will shed more light on the method. Let us also
mention that the algebraic Bethe ansatz was already applied to solve open Gaudin model using
the vertex-IRF correspondence of the XXZ spin chain [26, 43]. However, as pointed in [14], the
model considered in [26, 43] is obtained as a limit from the XXZ spin chain and is different of
the one studied here which is constructed directly from a classical r-matrix.
This letter is organised as follows. In Section 2, we recall well-known results about the Gaudin
algebra and model to fix the notations used in the following. Then, we provide technical results
in the Sections 3.1 and 3.2 to implement the algebraic Bethe ansatz. In Section 3.3, we recover
the results for the diagonal and triangular boundaries [27, 31, 32]. Section 3.4 contains new
mailto:nicolas.crampe@umontpellier.fr
https://doi.org/10.3842/SIGMA.2017.094
2 N. Crampe
results for the generic boundaries when the length of the chain is odd. Finally, we deal with the
even chain in Section 3.5 where an additional computation is necessary. We explain also why
the computation is different depending on the parity of the length of the chain. In Section 4,
we generalize these results to solve the higher spin Gaudin models with generic boundary.
2 Gaudin models and algebras
In this section, we recall different well-known notions to construct Gaudin models and prove
their integrability.
For a matrix r(x, y) ∈ End(CN ⊗ CN ) depending on 2 parameters x, y, one defines the
fundamental equation, called classical Yang–Baxter equation, given by
[r13(x1, x3), r23(x2, x3)] = [r21(x2, x1), r13(x1, x3)] + [r23(x2, x3), r12(x1, x2)].
In the previous relation, we have used the usual notations: r12(x) = r(x)⊗ I, r23(x) = I⊗ r(x),
. . . , where I is the identity matrix. A solution of the classical Yang–Baxter equation is called
a r-matrix. In addition, if an r-matrix r(x, y) depends only on the quotient x/y and satisfies
the supplementary relations r12(x, y) = r12(x/y) = −r21(y/x) = −r12(y/x)t1t2 then it is called
a skew-symmetric r-matrix. For a given skew-symmetric r-matrix r(x/y), one defines also the
reflection equation
r12(x/y)k1(x)k2(y)− k1(x)k2(y)r21(x/y) = k2(y)r12(xy)k1(x)− k1(x)r21(xy)k2(y).
A solution k(x) of the reflection equation is called k-matrix. From a skew-symmetric r-matrix
and an associated k-matrix, we can construct a new r-matrix [40]
r12(x, y) = r12(x/y)− k1(x)r12(1/(xy))k1(x)−1, (2.1)
which is not skew-symmetric.
To each r-matrix r(x, y), on can associate a Lie algebra R. Indeed, if K (x) ∈ End
(
CN
)
⊗R,
then the following relation
[K1(x),K2(y)] = [r21(y, x),K1(x)] + [K2(y), r12(x, y)] (2.2)
defines a Lie commutator for R. Indeed, the antisymmetry of the product is obvious and the
Jacobi identity is satisfied due to the classical Yang–Baxter equation. For any invertible matrix
M(x) ∈ End
(
CN
)
, there is a Lie algebra isomorphism given by
ΦM : R → RM ,
K0(x) 7→M0(x)−1K0(x)M0(x),
where RM is the Lie algebra defined by (2.2) with the r-matrix
M1(x)−1M2(y)−1r12(x, y)M1(x)M2(y).
From the defining relations (2.2), one can show that the transfer matrix
t(x) =
1
4
tr0
(
K0(x)2
)
(2.3)
satisfies [t(x), t(y)] = 0. Then, when a representation for R is chosen, the coefficients of t(x)
provide an integrable hierarchy. Let us remark that, for a given matrix M(x), R and RM give
the same hierarchy since it is easy to see that tr
(
K0(x)2
)
= tr
(
ΦM (K0(x))2
)
.
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary 3
In this letter, we focus on the skew-symmetric r-matrix associated to the affine Kac–Moody
algebra ŝl2:
r(x) =
1
x− 1
−1
2(x+ 1) 0 0 0
0 1
2(x+ 1) −2 0
0 −2x 1
2(x+ 1) 0
0 0 0 −1
2(x+ 1)
. (2.4)
The most general k-matrix associated to this r-matrix is given by
k(x) =
(
β + γ/x −α(β+ρ)
2 (x− 1/x)
β−ρ
2α (x− 1/x) β + γx
)
, (2.5)
where α, β, γ and ρ are scalar parameters. As explained previously, the r-matrix and the k-
matrix allow us to construct a new r-matrix (see relation (2.1)). Without loss of generality, we
can use the Lie algebra isomorphism ΦM to transform this r-matrix. Therefore, we use in this
paper the following r-matrix
r̃12(x, y) = M1(x)−1M2(y)−1
(
r12(x/y)− k1(x)r12(1/(xy))k1(x)−1
)
M1(x)M2(y), (2.6)
where r(x) is given by (2.4), k(x) is given by (2.5) and
M(x) =
(
β+ρ
2ρ
α
x
(β−ρ)x
2αρ 1
)
.
Explicitly, this r-matrix reads as
r̃(x, y) =
(
−1
2ω(x, y)σz b(x)σz + f(y, x)σ−
c(x)σz − f(1/y, 1/x)σ+ 1
2ω(x, y)σz
)
,
where σz =
(
1 0
0 −1
)
, σ+ =
(
0 1
0 0
)
and σ− =
(
0 0
1 0
)
and
δ(x) =
((β − ρ)x2 + 2γx+ ρ+ β)
2ρ(x− 1/x)
, b(x) = − α
δ(x)
, c(x) =
ρ2 − β2
4αρ2δ(1/x)
, (2.7)
ω(x, y) =
x+ y
x− y
+
xy + 1
xy − 1
and f(x, y) = ω(x, y)
δ(x)
δ(y)
. (2.8)
Let us emphasize that the r-matrix r̃(x, y) has a simpler form than the one without the isomor-
phism (i.e., r12(x/y) − k1(x)r12(1/(xy))k1(x)−1). Without the isomorphism, the 16 entries of
the r-matrix should be non-zero whereas r̃(x, y) has 6 vanishing entries and each two by two
block submatrices are diagonal or triangular.
As explained previously, the r-matrix r̃(x, y) allows one to define a Lie algebra by using
relation (2.2) with
K̃0(x) =
(
A(x) B̃(x)
C̃(x) −A(x)
)
. (2.9)
We denote by R̃ this Lie algebra. Let us mention that the Lie algebra R̃ for special choice of
the parameters in the k-matrix has been identified as Onsager algebras [3]1.
1To compare this letter with [3], we must use the additional property r12(x)t1 = σy
1r12(x)σy
1 satisfied by the
r-matrix (2.4).
4 N. Crampe
Finally, let us introduce the following representation2
K̃0(x) =
L∑
j=1
r̃0j(x, vj). (2.10)
By using the classical Yang–Baxter equation, it is obvious to show that (2.10) satisfies the
commutation relations of R̃. The parameters vi are called inhomogeneous parameters. The
Hamiltonians are usually defined by
H̃j = Resx=vj t(x) =
1
4
Resx=vj tr0
(
K̃0(x)2
)
= −vj
L∑
p=1
p 6=j
r̃jp(vj , vp) + Lvj
(
3v2j
v2j − 1
− ν(vj)
)
,
where
ν(z) =
z2(ρ2 − β2)− 2βγz + ρ2 − β2
2ρ2(z2 − 1)δ(z)δ(1/z)
. (2.11)
These Hamiltonians are integrable due to the commutation of the transfer matrix and [H̃j , H̃k]
= 0. They are related to the Gaudin models with boundary (or also called BC Gaudin models).
Indeed, by conjugating the previous Hamilonians by M1(v1) · · ·ML(vL), one gets the Hamilto-
nians
Hj = −vj
L∑
p=1
p 6=j
(
rjp
(
vj
vp
)
+ kj(vj)rjp
(
1
vjvp
)
kj(vj)
−1
)
+ Lvj
(
3v2j
v2j − 1
− ν(vj)
)
,
which are the usual form for the XXZ Gaudin model with boundary.
3 Modified algebraic Bethe ansatz
In this section, we want to diagonalize the matrix 1
4 tr0
(
K̃0(x)2
)
when the 4 parameters in the
k-matrix are generic. The diagonalisation when the k-matrix is diagonal was obtained in [27]
and for the triangular case in [31, 32]. For the generic case, the situation is more complicated.
Here, we use the modified algebraic Bethe ansatz, introduced and used in different contexts
[1, 5, 6, 8, 15], to compute the eigenvalues and the eigenvectors for the generic case. The first
two subsections are technical. Then, Section 3.3 contains the known results about the diagonal
or triangular boundaries recovered here and Sections 3.4 and 3.5 contains the new results for
the chain of length odd or even.
3.1 Algebraic relations
In this subsection, by using the commutation relations obtained from (2.2), we provide necessary
propositions to perform the algebraic Bethe ansatz. Firstly, let us mention that the commutation
relations between A(x), B̃(x) and C̃(x) given by (2.2) with the r-matrix (2.6) and K̃ (x) given
by (2.9) are quite complicated in comparison to the usual ones. This is why we introduce the
two following shifted generators
B(x, n) = B̃(x)− (2n− 1)b(x) and C(x, n) = C̃(x)− (2n− 1)c(x),
2We keep the same notation K̃ (x) for the algebraic element or its representation.
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary 5
where b(x) and c(x) are defined by (2.7). For these generators and the generators A(x), the
commutation relations become
A(x)A(y) = A(y)A(x), (3.1)
B(x, n)B(y, n+ 1) = B(y, n)B(x, n+ 1), (3.2)
C(x, n)C(y, n− 1) = C(y, n)C(x, n− 1), (3.3)
A(x)B(y,m) = B(y,m)A(x) + ω(x, y)B(y,m)− f(x, y)B(x,m), (3.4)
A(x)C(y,m) = C(y,m)A(x)− ω(x, y)C(y,m)− f(1/x, 1/y)C(x,m), (3.5)
C(x, n)B(y, n) = B(y, n+ 1)C(x, n+ 1)
+ 2f(x, y)A(x)− 2f(1/y, 1/x)A(y)− 8nb(y)c(x), (3.6)
where ω(x, y) and f(x, y) are defined by (2.8).
We introduce also a shifted transfer matrix
t(x, j) =
1
2
A(x)2 +
1
4
B(x, j + 1)C(x, j + 1) +
1
4
C(x, j)B(x, j)− 1
2
b(x)c(x). (3.7)
For j = 0 in the previous relation, the transfer matrix (2.3) is recovered: t(x) = 1
4 tr0 K̃0(x)2 =
t(x, 0).
Secondly, we present a lemma permitting to commute the shifted transfer matrix with the
generator B(x, p). This lemma is standard in the context of the algebraic Bethe ansatz for the
Gaudin model except for the modification of the shift for the transfer matrix.
Lemma 3.1. If the algebraic relations (3.1)–(3.6) are satisfied, the following relation holds
t(x, p− 1)B(z, p)−B(z, p)t(x, p) = B(z, p)
(
ω(x, z)(A(x) + ν(x)) + 2b(x)c(x)
)
−B(x, p)
(
f(1/z, 1/x)(A(z) + ν(z)) + 4(p− 1)c(x)b(z)
)
, (3.8)
where ν(z) is defined by (2.11).
Proof. By straightforward computations using the commutation relations (3.1)–(3.6), one gets
t(x, p− 1)B(z, p)−B(z, p)t(x, p) = ω(x, z)B(z, p)A(x)− f(1/z, 1/x)B(x, p)A(z)
+
1
2
(
ω(x, z)2 + f(z, x)f(1/z, 1/x)
)
B(z, p)
− 1
2
(
ω(z, x)f(1/z, 1/x) + ω(x, z)f(x, z) + 4(2p− 1)c(x)b(z)
)
B(x, p).
Then, by using the following functional relations
ω(x, z)2 + f(z, x)f(1/z, 1/x) = 2ν(x)ω(x, z) + 4b(x)c(x),
ω(z, x)f(1/z, 1/x) + ω(x, z)f(x, z) = 2ν(z)f(1/z, 1/x)− 4c(x)b(z)
one gets relation (3.8). �
Let us define the following product
B(z) = B(z1, 1)B(z2, 2) · · ·B(zM ,M),
where z = {z1, z2, . . . , zM}. Let us remark that, due to relation (3.2), B(z) is invariant under
any permutation of the zi which justifies the definition of z as a set. One introduces also, for
k = 1, 2, . . . ,M
B(zk, x) = B(z1, 1) · · ·B(zk−1, k − 1)B(x, k)B(zk+1, k + 1) · · ·B(zM ,M).
Finally, we are in position to compute the action of the transfer matrix with the product B(z):
6 N. Crampe
Proposition 3.2. If the algebraic relations (3.1)–(3.6) are satisfied, the following relation holds
t(x)B(z) =
1
2
B(z)B(x,M + 1)C(x,M + 1)
+ B(z)
Λ(x) +
M∑
p=1
ω(x, zp)
A(x) + ν(x) +
∑
q 6=p
(x− 1/x)ω(zp, zq)
zp − 1/zp
−
M∑
p=1
B(zp, x)f(1/zp, 1/x)
A(zp) + ν(zp) +
∑
q 6=p
ω(zp, zq)
,
where
Λ(x) =
1
2
A(x)2 + xA′(x) + ν(x)A(x)− x2 + 1
x2 − 1
A(x)− 1
2
b(x)c(x).
Proof. One gets, using a telescopic sum, the following relation
t(x, 0)B(z)− B(z)t(x,M) =
M∑
p=1
B1 · · ·Bp−1 (t(x, p− 1)Bp −Bpt(x, p))Bp+1 · · ·BM ,
where Bp = B(zp, p). Then, using the Lemma 3.1, we express t(x, p − 1)Bp − Bpt(x, p) in the
previous sum. Then, we use repetitively relation (3.4) to put the operators A on the right.
Then, we rearrange the terms using the following relations
ω(x, zp)f(x, zq)− f(1/zp, 1/x)f(zp, zq) = ω(zq, zp)f(1/zq, 1/x)− 4c(x)b(zq),∑
q>p
ω(x, zp)ω(x, zq) =
∑
q 6=p
x− 1/x
zp − 1/zp
ω(zp, zq)ω(x, zp).
Finally, t(x,M) is simplified using relation (3.7) and (3.6) for y → x to prove the proposition. �
Let us emphasize again that this type of computation and result are the same for the periodic
or diagonal Gaudin model. The only difference lies on the shift in the operators.
3.2 Explicit representation
Up to now, we have only used the commutation relations of the algebra to prove the previous
results. Now, we choose the representation for K̃ (x) given by (2.10). Therefore, we get explicitly
A(x) = −1
2
L∑
j=1
ω(x, vj)σ
z
j ,
B(x, n) = (1− 2n)b(x) + b(x)
L∑
j=1
σzj +
L∑
j=1
f(vj , x)σ−j , (3.9)
C(x, n) = (1− 2n)c(x) + c(x)
L∑
j=1
σzj −
L∑
j=1
f(1/vj , 1/x)σ+j .
Let us introduce the following vectors Ω =
(
1
0
)⊗L
, called usually pseudo-vacuum. We deduce
that
A(x)Ω = a(x)Ω and C(x, n)Ω = (L+ 1− 2n)c(x)Ω (3.10)
with a(x) = −1
2
L∑
j=1
ω(x, vj).
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary 7
The Bethe vectors are defined as follows
V(z) = B(z)Ω. (3.11)
These Bethe vectors are of course invariant under any permutation of the zi since B(z) and C(z)
are invariant. From the previous results (Proposition 3.2 and relations (3.10)), we get easily the
action of the transfer matrix on these Bethe vectors
t(x)V(z) =
L− 2M − 1
2
c(x)B(z)B(x,M + 1)Ω
+
λ(x) +
M∑
p=1
ω(x, zp)
a(x) + ν(x) +
∑
q 6=p
(x− 1/x)ω(zp, zq)
zp − 1/zp
V(z)
−
M∑
p=1
f(1/zp, 1/x)
a(zp) + ν(zp) +
∑
q 6=p
ω(zp, zq)
V(zp, x), (3.12)
where
λ(x) =
1
2
a(x)2 + xa′(x) + ν(x)a(x)− x2 + 1
x2 − 1
a(x)− 1
2
b(x)c(x).
The second term in (3.12) is the wanted term (i.e., is proportional to the Bethe vector) and the
third term is the unwanted term. These terms are already present for the periodic or triangular
case. For the generic boundary, we must deal with the first term in the previous relation (3.12).
There are different possibilities:
• c(x) = 0, i.e., ρ = ±β. This case corresponds to triangular (or diagonal) boundaries and
the Bethe ansatz works as usual. We get different sectors depending on the number M of
excitations (see Section 3.3).
• The length of chain is odd and M = L−1
2 . It corresponds to a case where the number
of excitations is fixed. From relation (3.12), one gets only the half of the spectrum. The
other half is obtained starting from t(x)V(z) (see Section 3.4).
• The length of the chain is even. In this case, we must compute B(z)B(x,M + 1)Ω. It
corresponds to the modified algebraic Bethe ansatz (see Section 3.5).
At first sight, the difference of the behavior for the lengths odd or even seems strange.
However, we can show that, in the case of odd length, two subspaces of (C2)⊗L are stabilized
by the transfer matrix. More precisely, if we denote by VJ the eigenspace of the total spin
Sz = 1
2
L∑
j=1
σzj with eigenvalue J , one gets
t(x)V1/2 ⊕ V3/2 ⊕ · · · ⊕ VL/2 ⊂ V1/2 ⊕ V3/2 ⊕ · · · ⊕ VL/2,
t(x)V−1/2 ⊕ V−3/2 ⊕ · · · ⊕ V−L/2 ⊂ V−1/2 ⊕ V−3/2 ⊕ · · · ⊕ V−L/2. (3.13)
Then it is natural that the spectrum for odd chain splits into two sectors. For the chain of even
length, this feature disappears: there is no stable subspace.
3.3 Triangular and diagonal boundaries
In this subsection, we deal with the triangular boundary (i.e., ρ = ±β). In this case, we can
easily see that one gets, for M = 0, 1, . . . , L,
t(x)VL/2−M ⊕ VL/2−M+1 ⊕ · · · ⊕ VL/2 ⊂ VL/2−M ⊕ VL/2−M+1 ⊕ · · · ⊕ VL/2.
8 N. Crampe
For M = 0, 1, . . . , L, the Bethe vectors V(z) defined by (3.11) is just in the sector VL/2−M ⊕
· · · ⊕ VL/2 and becomes an eigenvector of the transfer matrix if the Bethe roots zi satisfy the
Bethe equations:
a(zp) + ν(zp) +
M∑
q=1
q 6=p
ω(zp, zq) = 0 for p = 1, 2, . . . ,M. (3.14)
Explicitly, they read, for p = 1, . . . ,M ,
−1
2
L∑
j=1
(
zp + vj
zp − vj
+
zpvj + 1
zpvj − 1
)
+
1
2
(
γz − β
γz + β
+
βz − γ
βz + γ
)
+
M∑
q=1
q 6=p
(
zp + zq
zp − zq
+
zpzq + 1
zpzq − 1
)
= 0. (3.15)
The corresponding eigenvalue of the transfer matrix t(x) = 1
4 tr0
(
K̃0(x)2
)
is
λ(x) +
M∑
p=1
ω(x, zp)
a(x) + ν(x) +
M∑
q=1
q 6=p
(x− 1/x)ω(zp, zq)
zp − 1/zp
. (3.16)
Starting from the triangular boundary, we recover the diagonal one for α = 0. We see that the
eigenvalue (3.16) and the Bethe equations (3.15) do not depend on α. Then, as shown previously
for the Gaudin model [31, 32] (or for other models [7, 38]), the spectrum for triangular boundaries
is similar to the one with diagonal boundaries. Evidently, the Bethe eigenvectors depends on α,
via the function b(x) present in B(x, j) (see (3.9)). In particular, the function b(x) vanishes for
α = 0 and the Bethe vectors V(z) are now in the sector VL/2−M which is consistent with
t(x)VL/2−M ⊂ VL/2−M .
The last equation means that the total spin is conserved by the transfer matrix in the case of
diagonal boundary.
3.4 Generic boundaries for an odd chain
If the length of the chain is odd (i.e., there is an integer ` such that L = 2`+ 1), we can chose
the number of excitations as M = L−1
2 = ` and the coefficient in front of B(z)B(x,M + 1)Ω
vanishes. From the explicit expression (3.9) of B(x, j), we see that
V(z) ∈ V1/2 ⊕ V3/2 ⊕ · · · ⊕ VL/2 for M = `,
which is one of the subspaces stabilized by the transfer matrix, as explained previously. Then,
the eigenvalues for this sector are given by (3.16) with the Bethe equations (3.14) for M = `.
The spectrum and the associated eigenvectors for the other sector (3.13) are obtained from
the Bethe vectors
V(z) = C(z1, 0)C(z2,−1) · · ·C(z`,−`+ 1)
(
0
1
)⊗L
.
It is easy to show that V(z) belongs to the second sector. A relation similar to (3.12) by
replacing V(z) by V(z) can be proved. Then, we show that, starting from V(z), one gets
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary 9
exactly the expression (3.16) for the eigenvalues and (3.14) for the Bethe equations. Therefore,
the spectrums in both sectors are the same.
Let us also mention that the separation into two sectors is not new in this context. Indeed,
the complete spectrum of the XXZ spin chain with certain constraints between the parameters
of the boundaries is also obtained from two different Bethe vectors [42].
3.5 Generic boundaries for an even chain
For a chain with an even length, we must deal with the term B(z)B(x,M+1)Ω. For a generic M ,
this vector has no particular property but for M = L, special feature appears. This is the crucial
point of the modified algebraic Bethe ansatz.
Firstly, let us state the following lemma:
Lemma 3.3. The entries of the matrix
M(z) =
−→
L∏
q=1
(
(z − vq)(z − 1/vq)
(z − zq)(z − 1/zq)
B(zq, q)
)
×
(
(z − 1/z)
(z − x)(z − 1/x)
B(x, L+ 1)
)(
δ(z)B(z, L+ 1)
)−1
have only simple poles at z = 0,∞, z1, z2, . . . , zL, 1/z1, 1/z2, . . . , 1/zL, x, 1/x.
Proof. The proof of the existence of the poles ofM(z) at z = z1, z2, . . . , zL, 1/z1, 1/z2, . . . , 1/zL
and x, 1/x is straightforward. One gets
B̂(z) = δ(z)B(z, L+ 1) = α(2L+ 1)− α
L∑
j=1
σzj +
L∑
j=1
2z(vj − 1/vj)
(z − vj)(z − 1/vj)
δ(vj)σ
−
j .
This matrix B̂(z) is a lower triangular matrix with the diagonal entries in the set {α(L + 1),
α(L + 3), . . . , α(3L + 1)}. In particular, this shows that it is invertible which justifies the
definition of M(z). Due to the Cayley–Hamilton theorem, B̂(z)−1 is a polynomial of B̂(z).
Then, by using the fact that (σ−)2 = 0, we prove that the poles at z = vj and z = 1/vj
of B̂(z)−1 remain simple. Then, M(z) has no pole at these points. In addition, its becomes
a non-singular diagonal matrix for z = 0 which proves the simplicity of the pole of M(z) at
z = 0. The existence of the pole at z =∞ is proven using the property M(z) = −1/z2M(1/z)
which concludes the proof. �
Similar technical lemma has been proven and used in [15]. We are now in position to provide
the following proposition:
Proposition 3.4. The following relation holds
B(z)B(x, L+ 1)Ω = λ̂(x)V(z) +
L∑
p=1
(zp − 1/zp) Resx=zp
(
λ̂(x)
)
(zp − x)(zp − 1/x)
V(zp, x)
with z = {z1, z2, . . . , zL} and
λ̂(x) =
α(L+ 1)
δ(x)
L∏
q=1
(
(x− vq)(x− 1/vq)
(x− zq)(x− 1/zq)
)
.
10 N. Crampe
Proof. The residues of M(z)Ω are obtained by straightforward computations:
Resz=zp(M(z)Ω) =
1
δ(zp)
L∏
q=1
q 6=p
(
(zp − vq)(zp − 1/vq)
(zp − zq)(zp − 1/zq)
)(
(zp − vp)(zp − 1/vp)
(zp − x)(zp − 1/x)
)
V(zp, x),
Resz=x(M(z)Ω) =
1
δ(x)
L∏
q=1
(
(x− vq)(x− 1/vq)
(x− zq)(x− 1/zq)
)
V(z),
Resz=0(M(z)Ω) = − 1
α(L+ 1)
B(z)B(x, L+ 1)Ω.
By using M(z) = −1/z2M(1/z), we deduce that Resz=1/z0(M(z)Ω) = Resz=z0(M(z)Ω) where
z0 = 0, z1, z2, . . . , zL, x. The proof is concluded by using the fact that the sum over all the
residues (with the point at infinity) of a rational function vanishes. �
Therefore, by using relation (3.12) for M = L and the result of Proposition 3.4, we prove
that V(z) (for z = {z1, z2, . . . , zL}) is an eigenvector of the transfer matrix t(x) if the Bethe
roots satisfy the following Bethe equations
a(zp) + ν(zp) +
L∑
q=1
q 6=p
ω(zp, zq) =
(L+ 1)c(zp)
4zp
Resx=zp
(
λ̂(x)
)
, for p = 1, 2, . . . , L. (3.17)
To get the previous relation, we have used the functional relation
− (L+ 1)(z − 1/z)c(x)
2f(1/z, 1/x)(z − x)(z − 1/x)
=
(L+ 1)c(z)
4z
.
Explicitly, they become
−1
2
L∑
j=1
(
zp + vj
zp − vj
+
zpvj + 1
zpvj − 1
)
+ ν(zp) +
L∑
q=1
q 6=p
(
zp + zq
zp − zq
+
zpzq + 1
zpzq − 1
)
=
(ρ2 − β2)(L+ 1)2
16ρ2
(zp − vp)(zp − 1/vp)
δ(zp)δ(1/zp)(z2p − 1)
L∏
q=1
q 6=p
(zp − vq)(zp − 1/vq)
(zp − zq)(zp − 1/zq)
,
where ν(z) is given by (2.11). The corresponding eigenvalue is
λ(x)− L+ 1
2
c(x)λ̂(x) +
M∑
p=1
ω(x, zp)
a(x) + ν(x) +
∑
q 6=p
(x− 1/x)ω(zp, zq)
zp − 1/zp
. (3.18)
Then the l.h.s. of the Bethe equations is similar to the ones for the diagonal or triangular
boundaries. The r.h.s. corresponds to the inhomogeneous part of the TQ relation discovered
in [11, 12] to solve the XXZ spin chain with twist or boundaries.
4 Higher spin Gaudin models
In this section, we generalize the previous construction and consider the higher spin Gaudin
models. Indeed, it is well-known that we can find representations of the algebra R using
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary 11
different representations of sl(2) at each sites. Indeed, the following matrix K̃ (x) satisfies the
defining relation (2.2):
K̃0(x) =
L∑
j=1
L
(sj)
0j (x, vj), (4.1)
where sj = 1/2, 1, 3/2, . . . is the spin of the representation at the sites j and
L (s)(x, y) =
(
−ω(x, y)Sz 2b(x)Sz + f(y, x)S−
2c(x)Sz − f(1/y, 1/x)S+ ω(x, y)Sz
)
. (4.2)
In (4.2), Sz, S+ and S− are the usual matrices representing the sl(2) generators for spin s. The
only necessary knowledge about these matrices for the following constructions is the fact that
there exists a highest-weight vector ws such that
Szws = sws and S+ws = 0 for s = 1/2, 1, 3/2, . . .
For s1, s2, . . . , sL = 1
2 , the representation (4.1) becomes (2.10). Let us emphasize that the
representation of sl(2) may be different at each sites.
Now, we can generalize the results of Section 3. All the results of Section 3.1 remains
unchanged since it is independent of the representation.
The pseudo-vacuum becomes3
Ω = ws1 ⊗ ws2 ⊗ · · · ⊗ wsL .
Then, relation (3.12) is formally unchanged except the definition of the function a(x) which
becomes
a(x) = −
L∑
j=1
ω(x, vj)sj (4.3)
and the L in the first line must be replaced by L = 2
L∑
p=1
sp.
Then, there is again three possibilities to pursue the computation:
• The triangular case corresponding to c(x) = 0. The Bethe equations and the eigenvalues
are given respectively by (3.14) and (3.16) with a(x) given by (4.3). The number M ranges
now from 0 to L.
• The odd case M = L−1
2 . The discussion of Section 3.4 is still valid.
• The even case. To deal with this case, we must generalize the Lemma 3.3. This lemma is
still valid if we define M(z) as follows
M(z) =
L∏
q=1
(z − vq)2sq(z − 1/vq)
2sq
−→
L∏
q=1
B(zq, q)
(z − zq)(z − 1/zq)
× (z − 1/z)
(z − x)(z − 1/x)
B(x, L+ 1)
(
δ(z)B(z, L+ 1)
)−1
.
Let us emphasize that the first product is up to L, the number of sites, and the second up
to L. The demonstration is very similar to the one given previously except we must use
3In this section, we keep the same notations for all the objects which generalize the objects of the previous
sections.
12 N. Crampe
(S−)2s+1 = 0 for a spin s representation. Then, the Bethe equations and the eigenvalues
are given respectively by (3.17) and (3.18) with L replace by L, a(z) given by (4.3) and
λ̂(x) =
α(L+ 1)
δ(x)
L∏
q=1
(x− vq)2sq(x− 1/vq)
2sq
L∏
q=1
1
(x− zq)(x− 1/zq)
.
Acknowledgements
I thank P. Baseilhac, S. Belliard and V. Caudrelier for their interest. This work has been done
during the stay of the author at the “Laboratoire de Mathématiques et Physique Théorique
CNRS/UMR 7350, Université de Tours”. I thank warmly the LMPT for hospitality.
References
[1] Avan J., Belliard S., Grosjean N., Pimenta R.A., Modified algebraic Bethe ansatz for XXZ chain on the
segment – III: Proof, Nuclear Phys. B 899 (2015), 229–246, arXiv:1506.02147.
[2] Baseilhac P., Belliard S., The half-infinite XXZ chain in Onsager’s approach, Nuclear Phys. B 873 (2013),
550–584, arXiv:1211.6304.
[3] Baseilhac P., Belliard S., Crampe N., FRT presentation of the Onsager algebras, arXiv:1709.08555.
[4] Baseilhac P., Koizumi K., Exact spectrum of the XXZ open spin chain from the q-Onsager algebra repre-
sentation theory, J. Stat. Mech. Theory Exp. 2007 (2007), P09006, 27 pages, hep-th/0703106.
[5] Belliard S., Modified algebraic Bethe ansatz for XXZ chain on the segment – I: Triangular cases, Nuclear
Phys. B 892 (2015), 1–20, arXiv:1408.4840.
[6] Belliard S., Crampe N., Heisenberg XXX model with general boundaries: eigenvectors from algebraic Bethe
ansatz, SIGMA 9 (2013), 072, 12 pages, arXiv:1309.6165.
[7] Belliard S., Crampe N., Ragoucy E., Algebraic Bethe ansatz for open XXX model with triangular boundary
matrices, Lett. Math. Phys. 103 (2013), 493–506, arXiv:1209.4269.
[8] Belliard S., Pimenta R.A., Modified algebraic Bethe ansatz for XXZ chain on the segment—II: General
cases, Nuclear Phys. B 894 (2015), 527–552, arXiv:1412.7511.
[9] Belliard S., Pimenta R.A., Slavnov and Gaudin–Korepin formulas for models without U(1) symmetry: the
twisted XXX chain, SIGMA 11 (2015), 099, 12 pages, arXiv:1506.06550.
[10] Cao J., Lin H.-Q., Shi K.-J., Wang Y., Exact solution of XXZ spin chain with unparallel boundary fields,
Nuclear Phys. B 663 (2003), 487–519, cond-mat/0212163.
[11] Cao J., Yang W.-L., Shi K.-J., Wang Y., Off-diagonal Bethe ansatz and exact solution a topological spin
ring, Phys. Rev. Lett. 111 (2013), 137201, 5 pages, arXiv:1305.7328.
[12] Cao J., Yang W.-L., Shi K.-J., Wang Y., Off-diagonal Bethe ansatz solution of the XXX spin chain with
arbitrary boundary conditions, Nuclear Phys. B 875 (2013), 152–165, arXiv:1306.1742.
[13] Cao J., Yang W.-L., Shi K.-J., Wang Y., Off-diagonal Bethe ansatz solutions of the anisotropic spin- 1
2
chains
with arbitrary boundary fields, Nuclear Phys. B 877 (2013), 152–175, arXiv:1307.2023.
[14] Cirilo António N., Manojlović N., Nagy Z., Trigonometric s`(2) Gaudin model with boundary terms, Rev.
Math. Phys. 25 (2013), 1343004, 14 pages, arXiv:1303.2481.
[15] Crampe N., Algebraic Bethe ansatz for the totally asymmetric simple exclusion process with boundaries,
J. Phys. A: Math. Theor. 48 (2015), 08FT01, 12 pages, arXiv:1411.7954.
[16] Crampe N., Ragoucy E., Generalized coordinate Bethe ansatz for non-diagonal boundaries, Nuclear Phys. B
858 (2012), 502–512, arXiv:1105.0338.
[17] Crampe N., Ragoucy E., Simon D., Eigenvectors of open XXZ and ASEP models for a class of non-diagonal
boundary conditions, J. Stat. Mech. Theory Exp. 2010 (2010), P11038, 20 pages, arXiv:1009.4119.
[18] Crampe N., Ragoucy E., Simon D., Matrix coordinate Bethe ansatz: applications to XXZ and ASEP
models, J. Phys. A: Math. Theor. 44 (2011), 405003, 17 pages, arXiv:1106.4712.
[19] Faldella S., Kitanine N., Niccoli G., The complete spectrum and scalar products for the open spin-1/2
XXZ quantum chains with non-diagonal boundary terms, J. Stat. Mech. Theory Exp. 2014 (2014), P01011,
35 pages, arXiv:1307.3960.
https://doi.org/10.1016/j.nuclphysb.2015.08.006
https://arxiv.org/abs/1506.02147
https://doi.org/10.1016/j.nuclphysb.2013.05.003
https://arxiv.org/abs/1211.6304
https://arxiv.org/abs/1709.08555
https://doi.org/10.1088/1742-5468/2007/09/P09006
https://arxiv.org/abs/hep-th/0703106
https://doi.org/10.1016/j.nuclphysb.2015.01.003
https://doi.org/10.1016/j.nuclphysb.2015.01.003
https://arxiv.org/abs/1408.4840
https://doi.org/10.3842/SIGMA.2013.072
https://arxiv.org/abs/1309.6165
https://doi.org/10.1007/s11005-012-0601-6
https://arxiv.org/abs/1209.4269
https://doi.org/10.1016/j.nuclphysb.2015.03.016
https://arxiv.org/abs/1412.7511
https://doi.org/10.3842/SIGMA.2015.099
https://arxiv.org/abs/1506.06550
https://doi.org/10.1016/S0550-3213(03)00372-9
https://arxiv.org/abs/cond-mat/0212163
https://doi.org/10.1103/PhysRevLett.111.137201
https://arxiv.org/abs/1305.7328
https://doi.org/10.1016/j.nuclphysb.2013.06.022
https://arxiv.org/abs/1306.1742
https://doi.org/10.1016/j.nuclphysb.2013.10.001
https://arxiv.org/abs/1307.2023
https://doi.org/10.1142/S0129055X13430046
https://doi.org/10.1142/S0129055X13430046
https://arxiv.org/abs/1303.2481
https://doi.org/10.1088/1751-8113/48/8/08FT01
https://arxiv.org/abs/1411.7954
https://doi.org/10.1016/j.nuclphysb.2012.01.020
https://arxiv.org/abs/1105.0338
https://doi.org/10.1088/1742-5468/2010/11/P11038
https://arxiv.org/abs/1009.4119
https://doi.org/10.1088/1751-8113/44/40/405003
https://arxiv.org/abs/1106.4712
https://doi.org/10.1088/1742-5468/2014/01/P01011
https://arxiv.org/abs/1307.3960
Algebraic Bethe Ansatz for the XXZ Gaudin Models with Generic Boundary 13
[20] Faribault A., Tschirhart H., Common framework and quadratic Bethe equations for rational Gaudin magnets
in arbitrarily oriented magnetic fields, SciPost Phys. 3 (2017), 009, 24 pages, arXiv:1704.01873.
[21] Frahm H., Grelik J.H., Seel A., Wirth T., Functional Bethe ansatz methods for the open XXX chain,
J. Phys. A: Math. Theor. 44 (2011), 015001, 19 pages, arXiv:1009.1081.
[22] Frahm H., Seel A., Wirth T., Separation of variables in the open XXX chain, Nuclear Phys. B 802 (2008),
351–367, arXiv:0803.1776.
[23] Galleas W., Functional relations from the Yang–Baxter algebra: eigenvalues of the XXZ model with non-
diagonal twisted and open boundary conditions, Nuclear Phys. B 790 (2008), 524–542, arXiv:0708.0009.
[24] Gaudin M., Diagonalisation d’une classe d’Hamiltoniens de spin, J. Physique 37 (1976), 1089–1098.
[25] Gorissen M., Lazarescu A., Mallick K., Vanderzande C., Exact current statistics of the ASEP with open
boundaries, Phys. Rev. Lett. 109 (2012), 170601, 5 pages, arXiv:1207.6879.
[26] Hao K., Yang W.-L., Fan H., Liu S.-Y., Wu K., Yang Z.-Y., Zhang Y.-Z., Determinant representations
for scalar products of the XXZ Gaudin model with general boundary terms, Nuclear Phys. B 862 (2012),
835–849, arXiv:1205.0597.
[27] Hikami K., Separation of variables in the BC-type Gaudin magnet, J. Phys. A: Math. Gen. 28 (1995),
4053–4061, solv-int/9506001.
[28] Lazarescu A., Matrix ansatz for the fluctuations of the current in the ASEP with open boundaries,
J. Phys. A: Math. Theor. 46 (2013), 145003, 21 pages, arXiv:1212.3366.
[29] Lazarescu A., Mallick K., An exact formula for the statistics of the current in the TASEP with open
boundaries, J. Phys. A: Math. Theor. 44 (2011), 315001, 16 pages, arXiv:1104.5089.
[30] Lazarescu A., Pasquier V., Bethe ansatz and Q-operator for the open ASEP, J. Phys. A: Math. Theor. 47
(2014), 295202, 48 pages, arXiv:1403.6963.
[31] Manojlović N., Salom I., Algebraic Bethe ansatz for the XXZ Heisenberg spin chain with triangular bound-
aries and the corresponding Gaudin model, Nuclear Phys. B 923 (2017), 73–106, arXiv:1705.02235.
[32] Manojlović N., Salom I., Algebraic Bethe ansatz for the trigonometric s`(2) Gaudin model with triangular
boundary, arXiv:1709.06419.
[33] Melo C.S., Ribeiro G.A.P., Martins M.J., Bethe ansatz for the XXX-S chain with non-diagonal open
boundaries, Nuclear Phys. B 711 (2005), 565–603, nlin.SI/0411038.
[34] Murgan R., Nepomechie R.I., Bethe ansatz derived from the functional relations of the open XXZ chain for
new special cases, J. Stat. Mech. Theory Exp. 2005 (2005), P05007, 12 pages, hep-th/0504124.
[35] Nepomechie R.I., Bethe ansatz solution of the open XXZ chain with nondiagonal boundary terms,
J. Phys. A: Math. Gen. 37 (2004), 433–440, hep-th/0304092.
[36] Nepomechie R.I., An inhomogeneous T -Q equation for the open XXX chain with general boundary terms:
completeness and arbitrary spin, J. Phys. A: Math. Theor. 46 (2013), 442002, 7 pages, arXiv:1307.5049.
[37] Niccoli G., Non-diagonal open spin-1/2 XXZ quantum chains by separation of variables: complete spec-
trum and matrix elements of some quasi-local operators, J. Stat. Mech. Theory Exp. 2012 (2012), P10025,
42 pages, arXiv:1206.0646.
[38] Pimenta R.A., Lima-Santos A., Algebraic Bethe ansatz for the six vertex model with upper triangular
K-matrices, J. Phys. A: Math. Theor. 46 (2013), 455002, 13 pages, arXiv:1308.4446.
[39] Sklyanin E.K., Boundary conditions for integrable quantum systems, J. Phys. A: Math. Gen. 21 (1988),
2375–2389.
[40] Skrypnyk T., Generalized Gaudin spin chains, nonskew symmetric r-matrices, and reflection equation alge-
bras, J. Math. Phys. 48 (2007), 113521, 17 pages.
[41] Takhtadzhan L.A., Faddeev L.D., The quantum method for the inverse problem and the Heisenberg XY Z
model, Russian Math. Surveys 34 (1979), no. 5, 11–68.
[42] Yang W.-L., Zhang Y.-Z., On the second reference state and complete eigenstates of the open XXZ chain,
J. High Energy Phys. 2007 (2007), no. 4, 044, 11 pages, hep-th/0703222.
[43] Yang W.-L., Zhang Y.-Z., Sasaki R., An−1 Gaudin model with open boundaries, Nuclear Phys. B 729
(2005), 594–610, hep-th/0507148.
[44] Zhang X., Li Y.-Y., Cao J., Yang W.-L., Shi K.-J., Wang Y., Retrieve the Bethe states of quantum integrable
models solved via the off-diagonal Bethe Ansatz, J. Stat. Mech. Theory Exp. 2015 (2015), P05014, 18 pages,
arXiv:1407.5294.
https://doi.org/10.21468/SciPostPhys.3.2.009
https://arxiv.org/abs/1704.01873
https://doi.org/10.1088/1751-8113/44/1/015001
https://arxiv.org/abs/1009.1081
https://doi.org/10.1016/j.nuclphysb.2008.04.008
https://arxiv.org/abs/0803.1776
https://doi.org/10.1016/j.nuclphysb.2007.09.011
https://arxiv.org/abs/0708.0009
https://doi.org/10.1051/jphys:0197600370100108700
https://doi.org/10.1103/PhysRevLett.109.170601
https://arxiv.org/abs/1207.6879
https://doi.org/10.1016/j.nuclphysb.2012.05.019
https://arxiv.org/abs/1205.0597
https://doi.org/10.1088/0305-4470/28/14/023
https://arxiv.org/abs/solv-int/9506001
https://doi.org/10.1088/1751-8113/46/14/145003
https://arxiv.org/abs/1212.3366
https://doi.org/10.1088/1751-8113/44/31/315001
https://arxiv.org/abs/1104.5089
https://doi.org/10.1088/1751-8113/47/29/295202
https://arxiv.org/abs/1403.6963
https://doi.org/10.1016/j.nuclphysb.2017.07.017
https://arxiv.org/abs/1705.02235
https://arxiv.org/abs/1709.06419
https://doi.org/10.1016/j.nuclphysb.2004.12.008
https://arxiv.org/abs/nlin.SI/0411038
https://doi.org/10.1088/1742-5468/2005/05/P05007
https://arxiv.org/abs/hep-th/0504124
https://doi.org/10.1088/0305-4470/37/2/012
https://arxiv.org/abs/hep-th/0304092
https://doi.org/10.1088/1751-8113/46/44/442002
https://arxiv.org/abs/1307.5049
https://doi.org/10.1088/1742-5468/2012/10/P10025
https://arxiv.org/abs/1206.0646
https://doi.org/10.1088/1751-8113/46/45/455002
https://arxiv.org/abs/1308.4446
https://doi.org/10.1088/0305-4470/21/10/015
https://doi.org/10.1063/1.2816256
https://doi.org/10.1070/RM1979v034n05ABEH003909
https://doi.org/10.1088/1126-6708/2007/04/044
https://arxiv.org/abs/hep-th/0703222
https://doi.org/10.1016/j.nuclphysb.2005.09.003
https://arxiv.org/abs/hep-th/0507148
https://doi.org/10.1088/1742-5468/2015/05/P05014
https://arxiv.org/abs/1407.5294
1 Introduction
2 Gaudin models and algebras
3 Modified algebraic Bethe ansatz
3.1 Algebraic relations
3.2 Explicit representation
3.3 Triangular and diagonal boundaries
3.4 Generic boundaries for an odd chain
3.5 Generic boundaries for an even chain
4 Higher spin Gaudin models
References
|