Quantum trace formulae for the integrals of the hyperbolic Ruijsenaars-Schneider model

We conjecture the quantum analogues of the classical trace formulae for the integrals of motion of the quantum hyperbolic Ruijsenaars-Schneider model. This is done by departing from the classical construction where the corresponding model is obtained from the Heisenberg double by the Poisson reduction procedure. We also discuss some algebraic structures associated to the Lax matrix in the classical and quantum theory which arise upon introduction of the spectral parameter.


Introduction
The Ruijsenaars-Schneider (RS) models [1,2] continue to provide an outstanding theoretical laboratory for the study of various aspects of Liouville integrability, both at the classical and quantum level, see, for instance, [5][6][7][8][9][10]. Also, new interesting applications of these type of models were recently found in conformal field theories [11].
In this work we study some aspects related to the quantum integrability of the RS model with the hyperbolic potential. Recall that the definition of quantum integrability relies on the existence of a quantisation map which maps a complete involutive family of classical integrals of motion into a set of commuting operators on a Hilbert space. In general, there are different ways to choose a functional basis for this involutive family which is mirrored by the ring structure of the corresponding commuting operators. In particular, a classical integrable structure, most conveniently encoded into a Lax pair (L, M ), produces a set of canonical integrals which are simply the eigenvalues of the Lax matrix. Their commutativity relies on the existence of the classical r-matrix [12]. Provided this matrix exists one can build up different classical involutive families represented, for instance, by elementary symmetric functions of the eigenvalues of L or, alternatively, by traces TrL k for k ∈ Z. Concerning the particular class of the RS hyperbolic models, the JHEP05(2019)069 quantisation of a family of elementary symmetric functions associated to a properly chosen L is well known and given by the Macdonald operators [2,13]. In this paper we conjecture the quantum analogues of TrL k built up in terms of the same L-operator that is used to generate Macdonald operators through the determinant type formulae [14,15]. In fact there appear two commuting families I ± k that are given by the quantum trace formulae as quantisation of the classical integrals TrL k . In particular, R andR are two quantum dynamical R-matrices that depend rationally on the variables Q i = e q i , where q i , i = 1, . . . , N are coordinates, and satisfy a system of equations of Yang-Baxter type. Also, R is a parametric solution of the standard quantum Yang-Baxter equation. 1 Departing from I ± k and introducing q = e − , we then find that these integrals are related to the Macdonald operators S k through the q-deformed analogues of the determinant formulae that in the classical case relate the coefficients of characteristic polynomial of L with invariants constructed out of TrL k . The commutativity of I ± k and their relation to Macdonald operators has been checked by explicit computation for sufficiently large values of N .
We arrive to this expression for I ± k through the following chain of arguments. It is known that the Calogero-Moser-Sutherland models and their RS generalisations can be obtained at the classical level through the hamiltonian or Poisson reduction applied to a system exhibiting free motion on one of the suitably chosen initial finite-or infinitedimensional phase spaces [16]- [23], [3][4][5][6]. For instance, the RS model with the rational potential is obtained by the hamiltonian reduction of the cotangent bundle T * G = G G , where G ia Lie group and G is its Lie algebra. In [19] the corresponding reduction was developed for the Lie group G = GL(N, C) by employing a special parametrisation for the Lie algebra-valued element = T QT −1 ∈ G , where Q is a diagonal matrix and T is an element of the Frobenius group F ⊂ G. An analogous parametrisation is used for the group element g = U P −1 T −1 ∈ G, where U is Frobenius and P is diagonal. If one writes Q i = q i and P i = exp p i , then (p i , q i ) is a system of canonical variables with the Poisson bracket {p i , q j } = δ ij . In the new variables the Poisson structure of the cotangent bundle is then described in terms of the triangular dynamical matrix r satisfying the classical Yang-Baxter equation (CYBE) and of another matrixr. The cotangent bundle is easily quantised, in particular, the algebra of quantum T -generators is T 1 T 2 = T 2 T 1 R 12 and its consistency is guaranteed by the fact that the matrix R, being a quantisation of r, is triangular, R 12 R 21 = 1, and obeys the quantum Yang-Baxter equation. The quantum L-operator is then introduced as L = T −1 gT and it is an invariant under the action of F . In [19] the same formula for I k as given above 2 was derived by eliminating from the commuting operators Trg k = TrT L k T −1 an element T .
To build up the hyperbolic RS model, one can start from the Heisenberg double associated to a Lie group G. As a manifold, the Heisenberg double is G × G and it has a well-defined Poisson structure being a deformation of the one on T * G [24]. However,

JHEP05(2019)069
an attempt to repeat the same steps of the reduction procedure meets an obstacle: since the action of G on the Heisenberg double is Poisson, rather than hamiltonian, the Poisson bracket of two Frobenius invariants, {L 1 , L 2 }, is not closed, i.e. it is not expressed via L's alone. Moreover, for the same reason, the Poisson bracket {p i , p j } does not vanish on the Heisenberg double. On the other hand, a part of the non-abelian moment map generates second class constraints and to find the Poisson structure on the reduced manifold one has to resort to the Dirac bracket construction. 3 In this paper we work out the Dirac brackets for Frobenius invariants and show in detail how the cancellation of the non-invariant terms happens on the constraint surface. This leads to the canonical set of brackets for the degrees of freedom (p i , q i ) on the reduced manifold, the physical phase space of the RS model. However, continuing along the same path as in the rational case [19] does not seem to yield {T, L} and {T, T } brackets. The variable T is not invariant with respect to the stability subgroup of the moment map and computation of such brackets requires fixing a gauge, which makes the whole approach rather obscure. Moreover, the very simple and elegant bracket {L 1 , L 2 } emerging on the reduced phase space looks the same as in the rational case, with one exception: now the r-matrix r 12 entering this bracket is not skewsymmetric, i.e. r 12 = r 21 . We then find a quantisation of r 12 : a simple quantum R-matrix R + satisfying R +12 R −21 = 1, where R −12 is another solution of the quantum Yang-Baxter equation. In the absence of the triangular property for R +12 , assuming, for instance, the same algebra for T 's as in the rational case -that is T 1 T 2 = T 2 T 1 R +12 -would be inconsistent. Thus, at this point we simply conjecture that the integrals of the hyperbolic model have absolutely the same form as in the rational case, with the exception that the rational R-matrices are replaced by their hyperbolic analogues, which we explicitly find. That this conjecture yields integrals of motion can then be verified by tedious but direct computation and indeed holds true. Working out explicit expressions for these integrals for small numbers N of particles we find the determinant formulae relating these integrals to the standard basis of Macdonald operators. The rest of the paper is devoted to the model whose formulation includes the spectral parameter. Neither for the rational nor for the hyperbolic case the spectral parameter is actually needed to demonstrate their Liouville integrability, but its introduction leads to interesting algebraic structures and clarifies the origin of the shifted Yang-Baxter equation [22] and its scale-violating solutions.
The paper is organised as follows. In the next section we show how to obtain the hyperbolic RS model by the Poisson reduction of the Heisenberg double. This includes the derivation of the Poisson algebra of the Lax matrix via the Dirac bracket construction. We also introduce the spectral parameter and build up the theory based on spectral parameter-dependent (baxterised) r-matrices. We also describe a freedom in the definition of r-matrices that does not change the Poisson algebra of L's. In section 3 we consider the corresponding quantum theory. Finding the hyperbolic quantum R-matrices R ± and R, we conjecture our main formula for the quantum integrals I ± k and explain how it is related to the basis of the Macdonald operators. The rest of the section is devoted to the JHEP05(2019)069 quantum baxterised R-matrices and the quantum L-operator algebra. We show that in spite of the fact that the constant R-matrices satisfy the usual system of quantum Yang-Baxter equations, their baxterised counterparts instead obey its modification that involves rescalings of the spectral parameter with the quantum deformation parameter q = e − . Some technical details are relegated to two appendices. All considerations are done in the context of holomorphic integrable systems. 2 The classical model from reduction

Moment map and Lax matrix
We start with recalling the construction of the classical Heisenberg double associated to the group G = GL(N, C). Let the entries of matrices A, B ∈ G generate the coordinate ring of the algebra of functions on the Heisenberg double. The Heisenberg double is a Poisson manifold with the following Poisson brackets Here and elsewhere in the paper we use the standard notation where the indices 1 and 2 denote the different matrix spaces. The matrix quantities r ± are the following r-matrices In the following we also need the split Casimir 3) whose action on the tensor product C N ⊗ C N is a permutation. In these formulae E ij stand for the standard matrix units. The r-matrices (2.2) satisfy the classical Yang-Baxter equation (CYBE) and have the following properties: r + − r − = C and r ±21 = −r ∓12 . The variables (A, B) can be interpreted as a pair of monodromies of a flat connection on a punctured torus around its two fundamental cycles [30]. The monodromies are not gauge invariants as they undergo an adjoint action of the group of residual gauge transformations which coincides with G A → hAh −1 , B → hBh −1 . (2.4)

JHEP05(2019)069
If G is a Poisson-Lie group with the Sklyanin bracket and it generates the following infinitesimal transformations of (A, B) To perform the reduction, we fix the moment map to the following value where n is the Lie algebra element n = e ⊗ e t − 1 , where e is an N -dimensional vector with all its entries equal to unity, e t = (1, . . . , 1), and γ is a formal parameter which will be eventually interpreted as the coupling constant. Fixing this value of the moment map is posteriorly motivated by the fact that the dynamical model arising on the reduced space will have a close connection to the RS model we are after. We are thus led to find all A, B that solve the following matrix equation where on the right-hand side we worked out the explicit form of the exponential exp(γn). In the following we adopt the concise notation To solve (2.10), we introduce a convenient representation for A and B: Here Q and P are two diagonal matrices and T, U ∈ G are two Frobenius matrices, i.e. they satisfy the Frobenius condition T e = e , U e = e (2.14) and, therefore, belong to the Frobenius subgroup F of G.

JHEP05(2019)069
Introducing W = T −1 U ∈ F , equation (2.10) takes the form where we used the fact that U ∈ F . Furthermore, we write This equation can be elementary solved for W −1 and we get where we introduced c j = (e t T ) j . The condition W −1 ∈ F gives a set of equations to determine the coefficients c j : Here V is a Cauchy matrix with entries We apply the inverse of V , to obtain the following formula for the coefficients c j where we substituted β from (2.11). Finally, inverting W −1 we find W itself .

(2.18)
It is obvious, that eq. (2.10) becomes equivalent to the following two constraints

JHEP05(2019)069
where T, U ∈ F , and the quantities W (Q ), c(Q ) are given by (2.18) and (2.17), respectively. Any solution of e t T = c t can be constructed as T = hT 0 , where T 0 is a particular solution of this equation and h is a Frobenius group element which satisfies the additional constraint e t h = e t . In fact, the subgroup of ⊂ F ⊂ G determined by the conditions constitutes the stability group 4 of the moment map determined by the element n. Note that dim C F = N 2 − N and dim C = (N − 1) 2 . Now we can define a family of G-invariant dynamical systems 5 taking the combination L = W (Q )P −1 as their Lax matrix. Explicitly, After specifying the proper reality conditions, this L becomes nothing else but the Lax matrix of the RS family with the hyperbolic potential. Note that on the constrained surface the A, B-variables take the following form The reduced phase space can be singled out by fixing the gauge to, for instance, h = 1.

Poisson structure on the reduced phase space
Now we turn to the analysis of the Poisson structure of the reduced phase space. We find from (2.1) the following formula Next, we need to determine the bracket between Q j and P i . We have Here the first bracket on the right-hand side vanishes because all Q j commute with A. 6 To compute the second bracket, we consider the variation of B = U P −1 T −1 U −1 δB T P = U −1 δU − P −1 δP . 4 We do not include in the one dimensional dilatation subgroup C * {h ∈ G : h = c1, c = 0}, because its action on the phase space is not faithful. 5 The systems whose hamiltonians are invariant under the action of G. 6 The spectral invariants of A are central in the Poisson subalgebra of A, the latter is described by the Semenov-Tian-Shansky bracket [24] given by the first line in (2.1).

JHEP05(2019)069
Note that this formula does not include the variation δT . This is because T is solely determined by A, so is its variation. The condition δU e = 0 allows one to find We thus have and similarly one can check the bracket {Q i , Q j } = 0. These formulae suggests to employ the exponential parametrisation for both P and Q , that is, to set An -invariant extension of the Lax matrix away from the reduced phase space is naturally given by the following Frobenius invariant where T is an element of the Frobenius group entering the factorisation (2.12). The Poisson bracket of Q j with components of L is computed in a straightforward manner which is perfectly compatible with the form (2.21) of the Lax matrix on the reduced space.
In matrix form the previous formula reads as As to the brackets between the entries of L, this time they cannot be represented in terms of L alone but also involve T . Ultimately, such a structure is a consequence of the fact that the action of the Poisson-Lie group G on the phase space is Poisson rather than hamiltonian, so that there is an obstruction for the Poisson bracket of two Frobenius invariants to also be such an invariant. In addition, computing the Dirac brackets of L one cannot neglect a non-trivial contribution from the second class constraints and, therefore, the analysis of the Poisson structure for L requires, as an intermediate step, to understand the nature of the constraints (2.10) imposed in the process of reduction. The same argument holds for the Poisson brackets between any of the Frobenius invariants W = T −1 U and P , showing as a particular case that P i 's have a non-vanishing Poisson algebra on the Heisenberg double. 7 We save the details of the corresponding analysis for appendix B and JHEP05(2019)069 present here the final result for the Poisson bracket between the entries of the Lax matrix on the reduced phase space Clearly, the bracket (2.26) has the same form as the corresponding bracket for the rational RS model [19] albeit with new dynamical r-matrices for which we got the following explicit where we introduced the notation Q ij = Q i − Q j . This structure can be obtained as well after the computation of the Dirac brackets of W and P on the reduced phase space using the decomposition L = W P −1 . Remarkably the imposition of Dirac constraints makes the Poisson subalgebra {P i } abelian, allowing the interpretation of components p i = log P i as particle momenta. Concerning the properties of the matrices (2.27) and the Lax matrix, we note the following: first, r is expressed via r andr as Second, the matrix r is degenerate, det r = 0, and it obeys the characteristic equation r 2 = −r. Moreover, in contrast to the rational case [19], r is not symmetric, rather it has the property It is elementary to verify that the quantities are in involution with respect to (2.27). This property of I k is, of course, inherited from the same property for TrB k on the original phase space (2.1). We refer to (2.36) as the classical trace formula.

Introduction of a spectral parameter
Here we introduce a Lax matrix depending on a spectral parameter and discuss the associated algebraic structures and an alternative way to exhibit commuting integrals.
To start with, we point out one important identity satisfied by the Lax matrix (2.21). According to the moment map equation (2.15), we have The inverse on the right-hand side of the last expression can be computed with the help of the well-known Sherman-Morrison formula and we get where we used the fact that W is a Frobenius matrix, so that W e = e. Here the vector c has components (2.17) and satisfies the relation e t T = c t . Multiplying both sides of (2.19) with P −1 we obtain the following identity for the Lax matrix (2.21). Evidently, we can consider as another Lax matrix since the evolution equation of the latter is of the Lax forṁ

JHEP05(2019)069
where M is defined by the hamiltonian flow of L. Note that one can add to M any function of L without changing the evolution equation for L , which defines a class of equivalent M 's. Now, it turns out that due to the special dependence of L on the momentum, M and M fall in the same equivalence class. To demonstrate this point, it is enough to consider the simplest hamiltonian H = Tr L for which the matrix M is given by It follows from (2.25) that for the flow generated by this hamiltonian Therefore, Hence, M is in the same equivalence class as M and, therefore, we can take the dynamical matrix M to be the same for both L and L . The above observation motivates to introduce a Lax matrix depending on a spectral parameter just as a linear combination of L and L . Namely, we can define where λ ∈ C is the spectral parameter. The matrix L(λ) has a pole at zero and the original matrix L is obtained from L(λ) in the limit λ → ∞, in particular, The evolution equation for L(λ) must, therefore, be of the forṁ where M is the expression (2.42). The next task is to compute the Poisson brackets between the components of (2.43). We aim at finding a structure similar to (2.26), namely,

JHEP05(2019)069
where r(λ, µ), r (λ, µ) andr(λ) are some spectral-parameter-dependent r-matrices. We show how to derive these r-matrices in appendix B. Our considerations are essentially based on the identity (2.39). To state the corresponding result, we need the matrix The minimal solution 9 for the spectral-dependent r-matrices realising the Poisson algebra (2.46) is then found to be (2.48) The matrices r and r are skew-symmetric in the sense that Further, one can establish implications of the Jacobi identity satisfied by (2.46) for these r-matrices. Introducing the dilatation operator acting on the spectral parameter we find that the r-matrix r(λ, µ) does not satisfy the standard CYBE but rather the following modification thereof Following [23], we refer to (2.50) as the shifted classical Yang-Baxter equation. This equation can be rewritten in the form of the standard Yang-Baxter equation for the matrix differential operator There are also two more equations involving the matrixr The explanation of its minimal character will be given later.

JHEP05(2019)069
and [r 12 (λ),r 13 One can check that relations (2.50), (2.52) and (2.53) guarantee the fulfilment of the Jacobi identity for the brackets (2.25) and (2.26). Note that r is scale-invariant: (D λ + D µ )r (λ, µ) = 0, implying that it depends on the ratio λ/µ. This property does not hold, however, for r andr. The solution we found for the spectral-dependent dynamical r-matrices is minimal in the sense that there is a freedom to modify these r-matrices without changing the Poisson bracket (2.46). First of all, there is a trivial freedom of shifting r and r as where f is an arbitrary function of the ratio of the spectral parameters. This redefinition affects neither the bracket (2.26) nor equations (2.50), (2.52), (2.53).
Second, one can redefiner and r as where s(λ) is an arbitrary matrix function of the spectral parameter. Owing to the structure of the bracket (2.46) this redefinition of the r-matrices produces no effect on the latter, as r remains unchanged, while the matrix s decouples from the right-hand side of the LL bracket (see (2.46)). For generic s(λ), redefinition (2.55) affects, 10 however, equations (2.50), (2.52), (2.53). In particular, there exists a choice of s(λ) which turns the shifted Yang-Baxter equations forr and r into the conventional ones, where the derivative terms on the right hand side of (2.50), (2.52) and (2.53) are absent. One can take, for instance, (2.56) With the last choice the matrixr(λ) becomes while for r(λ, µ) one finds

JHEP05(2019)069
and the modified r-matrix is (2.58) The modified r-matrix still solves the CYBE and obeys the same relation (2.32).
There is no symmetry operating on r-matrices that would allow one to remove the scale-non-invariant terms from these matrices. Clearly, the r-matrices satisfying the shifted version of the Yang-Baxter equations have a simpler structure than their cousins subjected to the standard Yang-Baxter equations. This fact plays an important role when it comes to quantisation of the corresponding model and the associated algebraic structures. We also point out that the r-matrices we found here through considerations in appendix B also follow from the elliptic r-matrices of [23] upon their hyperbolic degeneration, albeit modulo the shift symmetries (2.54) and (2.55).
From (2.46) one then finds which, upon taking the limit λ → ∞, yields the Lax equation (2.45) with M given by (2.42). The conserved quantities are, therefore, , which generates I k (λ) in the power series expansion over the parameter ζ, defines the classical spectral curve (2.59) 3 Quantum model

Quantum Heisenberg double
At the classical level we obtained the hyperbolic RS model by means of the Poisson reduction of the Heisenberg double. It is therefore natural to start with the quantum analogue of the Heisenberg double. The Poisson algebra (2.1) can be straightforwardly quantised in the standard spirit of deformation theory. We thus introduce an associative unital algebra generated by the entries of matrices A, B modulo the relations [33] and they can be regarded as the quantisation of the Poisson relations (2.1). The quantum R-matrices here are defined as follows: first, we consider the following well-known solution

JHEP05(2019)069
of the quantum Yang-Baxter equation Using this R one can construct two more solutions R ± of the quantum Yang-Baxter equation, namely, These solutions are, therefore, related as and they also satisfy where C is the split Casimir. In the limit → 0 the matrices R ± expand as where r ± are the classical r-matrices (2.2). Further, we point out thatR ± = CR ± satisfy the Hecke condition (3.7) The first, or alternatively, the last line in (3.1) is a set of defining relations for the corresponding subalgebra that describes a quantisation of the Semenov-Tian-Shansky bracket, the latter has a set of Casimir functions generated by C k = TrA k . In the quantum case an analogue TrA k can be defined by means of the quantum trace formula where D is a diagonal matrix D = diag(q, q 2 , . . . , q n ). The elements C k are central in the subalgebra generated by A. Indeed, by successively using the permutation relations for A, one gets We then multiply both sides of this relation by D 1 and take the trace in the first matrix space i.e. Tr q A k is central in the subalgebra generated by A. Analogously, the I k = Tr q B k are central in the algebra generated by B and, in particular, the I k form a commutative family.
In principle, we can start with (3.1) and try to develop a proper parametrisation of the (A, B) generators suitable for reduction. It is an interesting path that should lead to understanding how to implement the Dirac constraints at the quantum level. We will find, however, a short cut to the algebra of the quantum L-operator.

Quantum R-matrices and the L-operator
An alternative route to the quantum R-matrices and to the corresponding L-operator algebra is based on the observation that in the classical theory, the Poisson brackets between the entries of the Lax matrix have the same structure (2.26) for both rational and hyperbolic cases. As a consequence, the equations satisfied by the classical rational and hyperbolic r-matrices are also the same. This should also be applied to the equations obeyed by the corresponding quantum R-matrices. We thus assume that the matrices R andR for the hyperbolic RS model satisfy the system of equations and R 12R13R23 =R 23R13 P 3 R 12 P −1 3 , (3.10) and have the standard semi-classical limit where they match the classical r-matrices (2.27).
Here and in the following (Q i , P i ) satisfy the quantum algebra being the standard quantisation of the Poisson algebra on the reduced phase space (2.23), (2.30). In fact, it is not difficult to guess a proper solution for these R-matrices based on the analogy with the rational case. For R we can take where r is given on the first line of (2.27). In the following we adopt the notation R + = R.
Since the classical r-matrix satisfies the property r 2 = −r, the exponential in (3.13) can be easily evaluated and we find (3.14) A direct check shows that (3.14) is a solution of (3.9). In comparison to the rational model, a new feature is that there exists yet another solution R − of the Yang-Baxter equation, namely, i.e. precisely in the same way as their non-dynamical counterparts, cf. (3.4). Furthermore, the matrices R ± satisfy equation They are also of Hecke type and the matricesR ± = CR ± have the following property Concerning the generalisation of equation (3.10) to the hyperbolic case, we can imagine two different versions -one involving R + and another R − , that is, 19) It appears that there exists a unique matrixR which satisfies both these equations. It is given byR and its inverse isR The matrix (3.20) also obeys (3.11), we find

JHEP05(2019)069
They also have another important property, usually referred to as the zero weight condition [32], Finally, the quantum L-operator is literally the same as its classical counterpart (2.21), of course with the natural replacement of p i by the corresponding derivative where ω = e −γ and j is the operator j = e − ∂ ∂q j . 11 On smooth functions f (Q 1 , . . . , Q N ) it acts as It is a straightforward exercise to check that this L-operator satisfies the algebraic relations with the R-matrices given by (3.14), (3.15), (3.20) and (3.24). The consistency of these relations follow from (3.16) and (3.26). One can alternatively derive equations (3.29) by direct quantisation of (2.28)-(2.29), where the classical matrix is chosen to be r 12 or, equivalently, −r 21 whose consistency follows from the same R-matrices relations. The algebraic relation (3.30) is also known as the quantum Frobenius group condition [19]. Concerning commuting integrals, the Heisenberg double has a natural commutative family I k = Tr q B k . It is not clear, however, how these integrals can be expressed via L, because we are lacking an analogue of the quantum factorisation formula B = T LT −1 , where T and L would be subjected to well-defined algebraic relations. Instead, what we could do is to conjecture the same formula as was obtained for quantum integrals in the rational case [19], where now the R-matrices are those of the hyperbolic model. Interestingly, the existence of two R-matrices, R ± , should give rise to two families of commuting integrals I ± k . Borrowing the corresponding expression from the rational case [19], we conjecture the following quantum trace formulae as quantisation of the classical integrals (2.36). In (3.32) the number k on the right-hand side gives a number of L 1 's and t 2 stands for the transposition in the second matrix space.
In particular, 11 In fact, j = P −1 j , we use j to signify that we talk about a particular representation for L.

JHEP05(2019)069
is a one-dimensional projector and from (3.14), (3.15) and (3.20) we get Commutativity of I ± k is then verified by tedious but direct computation which we do not reproduce here, rather our goal is to present a formula which relates I ± k with the commuting family given by Macdonald operators.
We denote by {S k } a commutative family of Macdonald operators, where The Macdonald operators have the following generating function where ζ is a formal parameter, L is the Lax operator (3.28). Under the sign : : of normal ordering the operators p j and q j are considered as commuting and upon algebraic evaluation of the determinant all j are brought to the right. In the classical theory the normal ordering is omitted and the corresponding generating function yields classical integrals of motion that are nothing else but the spectral invariants of the Lax matrix. We found an explicit formula that relates the families {I ± k } and S k . To present it, we need the notion of a q-number [k] q associated to an integer k so that [k] 1 = k, which corresponds to the limit → 0. Then S k is expressed via I + m or I − m as These formulae can be inverted to express each integral I ± k as the determinant of a k × k matrix depending on S j , namely, (3.39)

Spectral parameter and quantum L-operator
The quantum L-operator depending on the spectral parameter is naturally introduced as a normal ordered version of its classical counterpart where b j are the same as in (3.28). This L-operator satisfies the following quadratic relation In (3.41) the quantum R-matrices are (3.43) Here R + and R − are the solutions (3.14) and (3.15) of the quantum Yang-Baxter equation, R is (3.20) and we have introduced the matrix X ≡ X 12 , and We also present the formula for the inverse ofR(λ)

JHEP05(2019)069
where R ± are the same as given by (3.24). We note also the relation Finally, in addition to (3.42) there is one more relation between R(λ, µ) and R(λ, µ), namely, An interesting observation is that the combination solves the usual quantum Yang-Baxter equation with the spectral parameter. However, the full R-matrix in (3.43) differs from R YB by the terms that violate scale invariance. As a result, this matrix obeys the shifted version of the quantum Yang-Baxter equation, namely, In addition, there are two more equations -the one involving both R andR, and the other involvingR only, It is immediately recognisable that equations ( which is different from the canonical classical r-matrices (2.48) by allowed symmetry shifts. Thus, (3.43) should be regarded as a quantisation of the classical r-matrices satisfying the shifted Yang-Baxter equation. In this respect it is interesting to point out that the corresponding quantisation of the r-matrices solving the usual CYBE remains unknown. Finally, the algebra (3.41) should be completed by the following additional relations encoding the commutation properties of L with Q where Ω 12 = 1 − (1 − q)C 12 .

JHEP05(2019)069
Now we derive a couple of important consequences of the algebraic relation (3.41). Namely, we establish the quantum Lax representation, similar to the rational case, and also prove the commutativity of the operators TrL(λ) for different values of the spectral parameter.
Following considerations of the dynamics in the classical theory, we take H= lim λ→∞ TrL(λ) as the hamiltonian. From (3.41) we get where (3.48) was used. A straightforward computation reveals that the traces on the left and the right-hand side of the last expression are equal and that, for instance, where in the last expression we commuted the entries of L ij to the left so that it formally coincides with its classical counterpart (2.42). In the limit λ → ∞, A priori the trace of the commutator on the right-hand side might not be equal to zero, because it involves matrices with operator-valued entries. An involved calculation that uses representation (3.40) shows that it nevertheless vanishes, 12 identically for λ and µ.

JHEP05(2019)069
Fortunately, there is a simple and transparent way to show the commutativity of traces of the Lax operator, which directly relies on the algebraic relations (3.56), thus bypassing the construction of the quantum Lax pair. Indeed, let us multiply both sides of (3.41) with P −1 2R 12 (λ)P 2 R −1 12 (λ, µ) and take the trace with respect to both spaces. We get From (3.50) we have so that the right-hand side of the above equation can be transformed as Further progress is based on the fact that the matricesR 12 (λ) andR −1 12 (λ) are diagonal in the second space. We represent it in factorised form (3.20) and (3.44). Therefore, Although this expression involves the shift operator, it commutes with any function of coordinates q j , because when pushed through (3.64), this function will undergo the shifts of q j in opposite directions which compensate each other. Similarly, Consider first the left-hand side of (3.62) Using the cyclic property of the trace in the second space, this expression is equivalent to

JHEP05(2019)069
Taking into account that L = W P −1 and the commutativity of P −1 j G j (λ)P j with any function of coordinates, we arrive at Now we look at the right-hand side of (3.62): using the cyclic property of the trace, the matrix R 12 (λ, µ) can be moved to the left where it cancels with its inverse. This manipulation is allowed because L 1 (λ) and L 2 (µ) produce together a factor P −1 1 P −1 2 with which R 12 (λ, µ) commutes due to the zero weight condition (3.27). Also, the individual entries of R 12 (λ, µ) are freely moved through P −1 1R 21 (µ)P 1 , because of the diagonal structure of the latter matrix in the first matrix space, analogous to the similar property of (3.64). Then, to eliminateR 21 (µ), one employs the same procedure as was used for the left-hand side of (3.62) and the final result is TrL(λ)TrL(µ). This proves the commutativity of traces of the Lax matrix for different values of the spectral parameter.
We finally remark that writing the analogue of (3.36) with spectral parameter dependent Lax operator [14,15] : the quantities S k (λ) are commuting integrals and they are related to Macdonald operators (3.35) by a simple coupling-and spectral parameter-dependent rescaling

Conclusions
We have discussed the hyperbolic RS model in the context of Poisson reduction of the Heisenberg double [22]: we derive its Poisson structure and show that only on the reduced phase space does the Poisson algebra of the Lax matrix close and take a form very similar to the Lax matrix of the rational RS model [19]. We find a quantisation of the L-operator algebra governed by new R matrices R ± , along with a quantisation of the classical integrals in the form of quantum trace formulae I k (see (3.32)). We show how these quantum integrals are related to the well-known Macdonald operators through determinant formulae. Along the way we present a second Lax matrix that we can use to introduce a spectral parameter in the model. At the classical level this yields r-matrices that satisfy the shifted Yang-Baxter equation due to scale-violating terms. We show that this L-operator algebra admits a quantisation as well, with new R matrices satisfying the shifted quantum Yang-Baxter equation.
A particularly interesting observation is that one cannot obtain the quantum L-operator algebra from the quantum Heisenberg double in the same way as was done for the quantum cotangent bundle. It would be interesting to pursue the question whether and how one can impose the Dirac constraints after quantisation in order to reconstruct JHEP05(2019)069 the quantum L-operator algebra. A first step in that direction could be finding an analytic proof that the Dirac bracket for L on the reduced phase space is closed for general N . Another interesting question is to find the relation between our quantum trace formulae and the commuting traces obtained by the fusion procedure [34,35] for the equations (3.29). In addition, it would be interesting to extend our results to the RS models with spin, in particular, to those discussed in [8], as well as to find an analogue of the formulae (3.32) for the model with elliptic potential or for other series of Lie algebras. Constructing the quantum spin versions of these models could further aid the understanding of the RS type models that appear in the study of conformal blocks as in [11].

A.1 Lax matrix and its Poisson structure
Consider the following matrix function on the Heisenberg double where T is the Frobenius solution of the factorisation problem (2.12). On the reduced space L turns into the Lax matrix of the hyperbolic RS model. For this reason we continue to call (A.1) the Lax matrix and below we compute the Poisson brackets between the entries of L considered as functions on the Heisenberg double. This will constitute the first step towards evaluation of the corresponding Dirac bracket. The standard manipulations give where we defined the following quantities By using (2.1), we get

JHEP05(2019)069
Here we introduced the dressed r-matriceš which have proved themselves to be a useful tool for the present calculation. The dressed r-matrices have essentially the same properties as their undressed counterparts, most importantly,ř because C 12 is an invariant element. Thus, for (A.2) we get Now we proceed with evaluation of T. Taking onto account that T satisfies (2.14), in components we have (δ ia δ kb ζ aj,bl + δ ij δ kb ζ ja,bl + δ ia δ kl ζ aj,lb + δ ij δ kl ζ ja,lb ) .
Here Q ij = Q i − Q j and we introduced the concise notation Using (2.1) and the fact that A = T Q T −1 , we find that With the help of (A.4) we find in components where C ij,kl = δ jk δ il are the entries of C 12 . Substitution of this tensor into (A.6) yields the following expression In the first line the summation can be extended to all values of a and b, because the expression which is summed vanishes for a = j and independently for b = l. For the same JHEP05(2019)069 reason, we have extended the summation over a in the second line, where we also substitute the explicit value for C ij,kl = δ jk δ il . In this way we find T ij,kl = ab − δ ia δ kbř−aj,bl + δ ij δ kbř−ja,bl + δ ia δ klř−aj,lb − δ ij δ klř−ja,lb This further yields the following expression Here the second line can be written in the concise form as the matrix element r Q ij,kl of the following matrix Therefore, Hence, where we introduced three r-matrices, a, b and c with entries Needless to say, the bracket thus obtained is compatible with the Frobenius condition (2.14), which means that for any values of the free indices. Now we turn our attention to B 12 , which in components reads as where we introduced the notation

JHEP05(2019)069
With the help of (2.1) we get and by using (A.4) obtain for components the following expression η aj,kl = Q ja (L ksř−aj,sl −ř −aj,ks L sl ) + L ks C aj,sl Q j .
With this expression at hand, we get Here the summation in the first line can be extended to include the term with a = j because the latter vanishes. The second line can be conveniently written as a matrix element of some r-matrix. Namely, In matrix form where a 12 is the same matrix as in (A.8) and we introduced We also need Sinceř −21 = −ř +12 , we have where the coefficients are (A.14)

JHEP05(2019)069
First, we note that these coefficients satisfy 15) which guarantees that spectral invariants of L are in involution on the Heisenberg double. Second, in (A.14) the apparent dependence on the variable T occurs in the single combination a 21 + b 12 − c 12 . To make further progress, consider a 21 = C 12 a 12 C 12 , as C 12 acts as the permutation. We have, written in components, (a 21 ) ij,kl = C im,kn (a 12 ) mr,ns C rj,sl = δ mk δ in δ mr ař −ma,ns δ js δ rl = δ kl ař −ka,ij = −δ kl ař +ij,ka . Therefore, The dependence on T disappears and we find a simple answer The only T -dependence is in the coefficient c 12 . This coefficient cannot be simplified or cancelled, so we leave it in the present form. Substituting in (A.14) the matrices (A.7), (A.11) and (A.16) and, performing necessary simplifications, we obtain our final result for the coefficients of the bracket (A.13) (A. 17) In fact, the identity matrix 1 ⊗ 1 appearing in k − and s + can be omitted as it cancels out in the expression (A.13). As was already mentioned, the only T -dependence left over is in the term c 12 , namely, This gives or in other words,

Equation (A.19) leads to the following relations between the coefficiencients
Notice that the fact that the right-hand side of the first two expressions is an invariant tensor. Relations (A.20) guarantee that the bracket (A.13) is skew-symmetric. Following similar steps, we can derive the Poisson brackets involving other Frobenius invariants on the Heisenberg double, namely W ij and P i coordinates. Introducing the notations which, for Frobenius elements g, h satisfies c hg 21 + c gh we can write where matrices r 12 andr 12 are defined in (2.27). The c hg -like terms in the brackets (A. 21) are not Frobenius invariants, despite the arguments of the brackets are so, as it happens for (A.13). These terms disappear after imposing Dirac constraints in the reduced phase space, as it will explicitly shown for the LL-bracket in A.2.

A.2 Dirac bracket
Here we outline the construction of the Dirac bracket between the entries of the Lax matrix (A.1). We argue that the contribution to the Dirac bracket coming from the second JHEP05(2019)069 class constraints has the same matrix structure as (A.13) and that this contribution precisely cancels all the terms c 12 in (A. 17), so that the resulting coefficients describing the Dirac bracket on the constraint surface are given by expressions (2.27) in the main text. We start with the Poisson algebra of the non-abelian moment map This is the Semenov-Tian-Shansky type bracket; it has N Casimir functions Tr(M k ) with k = 1, . . . , N . On the constraint surface S the moment map is fixed to the following value see (2.10). Substituting this expression into the right-hand side of (A.22) yields the following answer where Θ is the Heaviside step function For any X ∈ Mat(N, C) introduce the following quantities From these quantities we construct the projectors π (i) that have the following action on X

JHEP05(2019)069
In particular, π (0) projects on the Lie algebra of and π (3) -on the one-dimensional dilatation subalgebra C * . The completeness condition is From (A.24) it is readily seen that Analogously, we find Thus, projections π (0) (M) and π (3) (M) constitute (N − 1) 2 + 1 = N 2 − 2N + 2 constraints of the first class. Projections π (1) and π (2) yield a non-degenerate matrix of Poisson brackets and, therefore, represent 2(N − 1) constraints of the second class. This matrix should be inverted and used to define the corresponding Dirac bracket. Even simpler, the matrix (A.24) has rank 2(N − 1) and we can use any non-degenerate submatrix of this rank to define the corresponding Dirac bracket. Now we derive the Poisson relations between the moment map M and the Lax matrix given by (A.1). First, we compute Deriving this formula, we have used (2.7), as well as the fact that T ∈ F . Next, we obtain It is clear that the diagonal entries from this expression of L commute with all the constraints: {M ij , L kk } = 0, even without restricting to the constrained surface. On the constrained surface where M is given by (A.23), we have Here From the explicit expression (A.30) and the fact that T is an element of the Frobenius group, we further deduce that

JHEP05(2019)069
In other words, L commutes on the constraint surface with all constraints of the first class, independently on the value of T . With the help of (A.31) we obtain Taking into account that Now we come to the Dirac bracket construction. By picking a non-degenerate submatrix Ψ of the matrix M ij,kl , we invert it and define the corresponding Dirac bracket Here I = (ij) is a generalised index which we use to label matrix elements of M ij,kl that comprise the non-degenerate matrix Ψ IJ . To give an example, for N = 3 we can take as Ψ the following matrix In particular det Ψ = β 4 (ω + β) 2 (ω + 2β) 2 . Inverting Ψ, we find that that is, the contribution of the second class constraints has precisely the same structure as (A.13). Moreover, the corresponding coefficients are where c 12 is given by (A.18). Thus, in (A.32) all the terms c 12 cancel out. We have also performed a similar computation for N = 4, 5, 6, 7 with the same result. An analytic derivation for arbitrary N is still missing, although our findings leave little doubt that it holds true. In summary, on the reduced phase space the Dirac bracket between the components of the Lax matrix has the form (A.13) with the following coefficients (A.34) The coefficients have the following properties k ± 12 + k ± 21 = ±(C 12 − 1 ⊗ 1) , s − 12 = −s + 21 , (A. 35) which guarantee, in particular, skew-symmetry of (A.13). In addition, they satisfy the relation (A.15). In the main text we present the formula (A.13) in the r-matrix form (2.26) with the following identifications k + = r , s + = −r , k − = −r .
This ansatz depends on five undermined parameters: α, β, δ, δ 12 and δ 21 , which should eventually be expressed via λ and µ. We then plug this ansatz together with the expression (2.43) for the spectral-dependent Lax matrix into (B.11) and, by using the reduction formulae (B.7), (B.8) and (B.9), bring the resulting expression to the following