Classical N-reflection equation and Gaudin models

We introduce the notion of N-reflection equation which provides a generalization of the usual classical reflection equation describing integrable boundary conditions. The latter is recovered as a special example of the N=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=2$$\end{document} case. The basic theory is established and illustrated with several examples of solutions of the N-reflection equation associated with the rational and trigonometric r-matrices. A central result is the construction of a Poisson algebra associated with a non-skew-symmetric r-matrix whose form is specified by a solution of the N-reflection equation. Generating functions of quantities in involution can be identified within this Poisson algebra. As an application, we construct new classical Gaudin-type Hamiltonians, particular cases of which are Gaudin Hamiltonians of BCL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$BC_L$$\end{document}-type.


Introduction
Classical integrable systems have been formulated in terms of the classical r-matrix in [1][2][3] for which the central equation is called the classical Yang-Baxter equation  Then, the study of boundary conditions which preserve the integrability of the model leads to the classical reflection equation [4] r ab (λ, ν) k a (λ)k b (ν) + k a (λ) r ba (ν, −λ) k b (ν) − k b (ν) r ab (λ, −ν) k a (λ) −k a (λ)k b (ν) r ba (−ν, −λ) = 0. (2) This equation appears naturally in classical integrable systems based on the BC L root system and can be interpreted via a Z 2 action on the A 2L root system. This point of view on integrable boundary conditions, sometimes called "folding" for short, has been used extensively, for example, in [5][6][7]. The search for integrable systems associated with the action of more complicated finite groups, and in particular of the cyclic group Z N , has attracted a lot of attention recently (see, e.g., [8][9][10][11][12][13][14][15][16]).
In this paper, we introduce a new equation which generalizes the classical reflection equation (2) as well as algebraic structures related to Z N (cyclotomic) models. In Sect. 2, we define the N -reflection equation and prove that it allows one to construct new solutions of the classical Yang-Baxter equation from old ones. Then, in Sect. 3, we show on some examples that this equation has interesting solutions and, in particular, solutions where the action on the spectral parameters is a Mobiüs transformation. In Sect. 4, we show how the N -reflection equation allows one to define a certain Poisson subalgebra of a linear Poisson algebra defined by a classical r -matrix and use this to obtain new integrable Gaudin models. We also show that the Hamiltonian equations of motion generated by the elements in involution in the Poisson subalgebra can be written in Lax form and give an explicit formula for the second matrix of the Lax pair.

Classical N-reflection equation
In this paper, we consider a solution r ab (λ, μ) of the classical Yang-Baxter equation (1) acting on C n ⊗ C n ⊗ C n . We now introduce the following generalization of the classical reflection equation associated with r . Definition 2.1 Let τ and g ( j) , j = 0, . . . , N −1 be functions of the spectral parameter, τ j be defined recursively by τ j (ν) = τ (τ j−1 (ν)) with τ 0 (ν) = ν and The classical N -reflection equation for τ , g ( j) and the matrix k is defined by Without loss of generality, we can set g (0) (ν) = 1.
For the sake of comparison with the usual reflection equation, let us write explicitly the classical 2-reflection equation by multiplying by k b (ν) on the right-hand side: Then, we have the following results establishing the connection with the usual reflection equation.
Proposition 2.1 Suppose the r -matrix depends only on the difference of the spectral parameters and is skew-symmetric, taking g (1) (ν) = −1 and τ (ν) = −ν, the classical 2-reflection Eq. (5) becomes the usual classical reflection equation In the other well-known case where the r -matrix depends only on the quotient of the spectral parameters and is skew-symmetric, the choice g (1) (ν) = −1 and τ (ν) = 1/ν gives the multiplicative version of the classical reflection equation.
However, let us emphasize that the classical 2-reflection equation introduced here generalizes the usual ones because of the introduction of the functions g (1) and τ . We will see on the examples below that we can get new interesting solutions from this generalized form. One remarkable feature is the possibility to act with Möbius transformations on the spectral parameter using τ .
Proof For convenience, for a given r -matrix, let us introduce the notation Note that with our definition of r , the classical N -reflection Eq. (4) may be rewritten compactly as In turn, replacing λ by τ (λ) in (9), we see that (9) implies, Now, using (10) in the form (11) and its counterpart under the exchange b ↔ c, μ ↔ ν, we can write The proposition follows from the fact that CY B E(r ) abc (λ, μ, ν) = 0.
This proposition allows us to construct a new solution r of the classical reflection equation for each solution k of the N -reflection equation and each solution r of the classical Yang-Baxter equation. It generalizes the results of [17], stated for the usual reflection equation. In general, the r-matrices r obtained by this construction are not skew-symmetric even if the starting r -matrix r is. However, they still allow for the construction of interesting algebraic structures and integrable models as shown in Sect. 4.

Remark 1
In the course of the proof, we used a compact form of the N -reflection equation which is worth pointing out separately It is convenient to introduce the notion of N-unitary relation which generalizes the so-called unitary relation (k(λ)k(−λ) ∝ I n or k(λ)k(1/λ) ∝ I n ) often required for the usual reflection equation.

Definition 2.2
The matrix k is said to satisfy the N -unitary relation if where f (ν) is a scalar function.
. This equation is automatically satisfied if the k-matrix satisfies the N-unitary relation.

Remark 3
If τ N (ν) = ν and g ( j) (ν) = ω j with ω N = 1, then a solution k satisfying the N-unitary relation (14) and the following symmetry relation is a solution of the N -reflection equation. This is a case of a reduction by a Z N action on the r -matrix. Equation (15) underlies the results obtained in [10] and is a particular case of the N -reflection equation.

Some solutions of the N-reflection equation
In this section, we present some solutions of the N -reflection equation for various values of N and for the standard rational and trigonometric r -matrices. Even in the case N = 2, some of these solutions are new compared to the usual reflection equation and, in fact, cannot be accommodated by it. As explained previously, all these solutions allow us to obtain new solutions of the nonstandard classical Yang-Baxter equation.

Rational r-matrix
In this subsection, we provide some solutions of the N -reflection equation associated with the rational r-matrix where P is the permutation operator of C n ⊗ C n . In this case, the following matrices are solutions of the usual reflection Eq. (2) where θ is a free parameter and G is a n × n-matrix satisfying G 2 = I n . This result generalizes to the case of the N -reflection equation as follows: where G is an n × n-matrix satisfying G N = I n , which is a solution of the N -reflection Proof We start by proving the N -unitary relation: We have used well-known relations for the roots of unity, such as Next, multiplying the N -reflection equation by k (ν) on the right and using the property of the permutation operator, the N -reflection equation is equivalent to By multiplying on the left by P ab and rearranging the sum, one gets equivalently We remark that the expression inside the square bracket in the L.H.S. of the previous relation can be written as follows by using the explicit form (18) Then, we recognize a telescopic sum in the L.H.S. of (21) and the N -reflection equation reduces to By considering separately the terms in the space a and the ones in the space b and by using the N -unitary relation proved previously, we obtain that relation (23) holds, which proves that the N -reflection equation is satisfied.
The importance of the functions τ and g can be seen from the existence of the following solutions of the N -reflection equation.

Proposition 3.2 The identity matrix is a solution of the 2-reflection equation with
where a, b and c are free parameters such that τ = 0, ∞ identically. In particular, The identity matrix is a solution of the 3-reflection equation with where a, b, c and d are parameters constrained by a 2 + ad + bc + d 2 = 0 and such that τ = 0, ∞ identically. In particular, τ 3 (ν) = ν.

Proof
The N -reflection equation for the rational r -matrix and for a k-matrix equal to the identity matrix reduces to one functional equation One can show by direct computation that for choices (24) and (25), the last equation holds. This finishes the proof.
Again, this proposition allows for the construction of new solutions of the nonstandard classical Yang-Baxter equation where the poles of the r -matrix are governed by the parameters of the function τ . For b = c = 0 and in the case N = 2, we recover the map τ (ν) = −ν usually used for the reflection equation in the rational case. Despite some effort, we have not been able to produce a solution where τ is a Möbius transformation and k is not the identity matrix in the rational case. However, below we present examples where this is possible in the trigonometric case.

Trigonometric r-matrix
In this subsection, we study the N -reflection equation in the case of the 4 × 4 trigonometric solution of the classical Yang-Baxter equation (1) defined by

Proposition 3.3 The 2-reflection equation with
has two solutions: the identity matrix and the following diagonal matrix The 3-reflection equation with (30) and the constraint a 2 + bc + ad + d 2 = 0 has five solutions: the identity matrix and the following diagonal matrices Proof The proposition is proven by direct computation.
In the particular case a = 0 and b = c, the 2-reflection equation becomes the usual one in the case of the trigonometric r -matrix since τ (ν) = 1/ν.

Poisson subalgebra and Lax pair of the equations of motion
Let us consider the following Poisson algebra, associated with a skew-symmetric matrix r solution of the classical Yang-Baxter equation and define

Proposition 4.1 The elements B(λ) form a Poisson subalgebra of the form
where r is given by (7) and

k is a solution of the N -reflection equation.
Proof This is shown by direct computation starting from (33) and using the definition of B(λ) (34). If one uses the N -reflection equation, one can bring all the terms produced in the right-hand side into the simple form given in (35). The skew-symmetry of the Poisson bracket (35) is evident, and the Jacobi identity is a consequence of the classical Yang-Baxter equation (1) satisfied by r . The latter fact is another motivation for introducing the N -reflection equation in the first place.
The Poisson algebra (35) provides an example of a linear Poisson structure based on a non-skew-symmetry r -matrix, as discussed in [18] (see also [19]).

Remark 4 If the hypotheses of Remark 3 hold, then
is a Z N action on the Poisson algebra (33) and B a (λ) contains the generators of the fixed-point subalgebra of this action.
Remark 5 Relation (34) looks similar to the reduction group scheme for the Lax operator considered in [20][21][22][23][24][25] in which k and τ provide a representation of the reduction group. In our case, the approach is different in that we do not assume a priori that k and τ are associated with a group. Instead, we look for conditions on them ensuring that the Poisson algebra for B(λ) closes. This leads to our N-reflection equation (4).
It is an open and interesting problem to study the interplay between the N-reflection equation and the reduction group approach.
It is well known [18] that the Poisson algebra (35) allows for the construction of Poisson commuting elements since, for p, q = 1, 2, . . . , one gets with The matrix M(λ, ν) satisfies the important relation This relation is a straightforward consequence of the N -reflection equation written as in (9). It is a generalization of the relation introduced in the context of the usual classical reflection equation [4] (which is recovered as explained after (5)). We recall that in that standard context, relation (40) can be taken as a definition of integrable boundary conditions. In [26], the Hamiltonian interpretation of this relation (and its generalization to the dynamical case) was presented in the case of the quadratic Poisson bracket and the usual reflection equation.
It is an open problem to generalize the results of [26] to the present case of the N -reflection equation.

Gaudin models
In this section, we apply the previous results to the construction of new integrable models which are generalizations of the Gaudin models [27]. The general procedure goes as follows. Suppose we are given a local representation of (33) in the form, 2 for m, p = 1, 2, . . . , L, Given a solution k of the N -reflection equation, we obtain a representation of relation (35) by setting where z m are free, mutually distinct, parameters. The Gaudin-type Hamiltonians are then defined by, These Hamiltonians are in involution due to relation (37) and allow us to define an integrable model. Therefore, each solution of the N -reflection equation allows us to get integrable Gaudin-type models based on a non-skew r -matrixr . In that sense, our construction produces explicit classes of models falling into the general scheme of [28] and generalizes the method used in [17].

Explicit example
We now present an explicit example associated with the 2-reflection equation in the rational case for n = 2. Let us introduce L copies of the su(2) Poisson algebra with generators {s + j , s − j , s z j } satisfying where s j are some parameters. Then, satisfies (41) with the rational r-matrix (16). From the solution of the 2-reflection equation given in Propositions 2.2, and 3.2 provides the following r-matrix .
Note that the r-matrix (46) is related to the rational r-matrix (16) via the transformation where p(μ) = b+cμ 2 2c(a−cμ) . This is an example of the equivalence relations discussed for instance in [29]. The Gaudin-type Hamiltonians (43) associated with this r-matrix read Note that in the case b = c = 0, we find the Gaudin Hamiltonians based on the BC L root system [30], consistently with the remark at the end of Sect. 3.1.
Similarly, from the solution of the 3-reflection equation given in Proposition 3.2, we obtain the following r-matrix with Proposition 2.2 (we recall that in this case This r-matrix (49) can be also obtained from the following transformation of the rational r-matrix (16) where (cμ+d) . In this case, the Gaudintype Hamiltonians (43) read In the case b = c = 0 and a d = ω = e 2iπ 3 , this reduces to a Z 3 -cyclotomic classical Gaudin Hamiltonians whose quantum counterparts were introduced in [10] and further developed in [14].

Conclusions and outlook
We introduced the N -reflection equation as a generalization of the usual reflection equation [4]. An important motivation is the possibility to define a consistent Poisson algebra (35) when one applies the following map: on the Poisson algebra The map (53) appears as a generalization of the maps producing the eigenspaces of the involution where k(λ) k(σ (λ)) = I and σ (λ) = −λ or σ (λ) = 1/λ. These are known to be related to the usual classical reflection equation (6) (or its multiplicative form). But in our case, we do not necessarily require a priori that there is an underlying map of order N which would be the generalization of the involution (56). Several consequences of this construction are noteworthy. Firstly, it gives a systematic way to construct many examples of non-skew-symmetry r -matrices. In this respect, we note that once such an r -matrix has been obtained, it can be taken as the starting point to repeat the procedure. It is an interesting question whether this process terminates or not and under which conditions on g ( j) , τ and k. The relationship between all the solutions of the classical Yang-Baxter equation that can be obtained in this way and the entire set of such solutions is also an open problem which would require an understanding of the full classification of non-skew-symmetric solutions of that equation. Such a classification is not known so far. Secondly, one can define new Gaudin models in a systematic way. Despite the possibility to relate the r -matrix of the new models to the standard rational r -matrix via rescaling and reparametrization, the nature of these models is quite different from the standard Gaudin models. Indeed, rescaling and reparametrization of an r -matrix are known to have consequences on the model, one of which is that reparametrization affects the skew-symmetry of the r -matrix, hence the (non) ultralocality of the model (in field theoretical language). In our case, these operations depend on the Möbius transformations acting on the spectral parameters and it is the first time that such transformations are used in this context to our knowledge. Thirdly, our results give a natural context for (40) in cases where τ is not of order 2. Such examples have appeared before in the context of linearizable boundary conditions [31]. However, there is no hint of an underlying Hamiltonian description of (40) in that context. We have provided such a Hamiltonian formulation in what we could call a time-independent setting. We hope that our results can be incorporated into a larger theory which would contain (40) as a special case, along the same lines as the results obtained in [26] in the standard reflection equation case. A major step in achieving this would be to find the appropriate dynamical generalization of our N -reflection equation. This is a completely open problem. A related open question is the problem of quantization of the N -reflection equation in the spirit of the quantization of the classical reflection equation [32,33]. In that context, it is known that the semiclassical limit of the quantum structure gives the appropriate Poisson algebra on the matrices k.