Abstract

We consider the classical integrable 1+1 trigonometric ${\rm gl}_{N}$ Landau-Lifshitz models constructed by means of quantum $R$-matrices satisfying also the associative Yang-Baxter equation. It is shown that 1+1 field analogue of the trigonometric Calogero-Moser-Sutherland model is gauge equivalent to the Landau-Lifshitz model which arises from the Antonov-Hasegawa-Zabrodin trigonometric non-standard $R$-matrix. The latter generalizes the Cherednik’s 7-vertex $R$-matrix in ${\rm GL}_{2}$ case to the case of ${\rm GL}_{N}$. Explicit change of variables between the 1+1 models is obtained.

1 Introduction

In this paper we consider two types of models arising as 1+1 field generalizations of classical integrable finite-dimensional systems. At the level of finite-dimensional mechanics the first one is the Calogero-Moser-Sutherland many-body system [10]. In the trigonometric case it is defined by the Hamiltonian

$$\begin{array}[]{c}\displaystyle{H^{\hbox{\tiny{CMS}}}=\frac{1}{2}\sum\limits_{i=1}^{N}\frac{p_{i}^{2}}{2}-\sum\limits_{i<j}\frac{c^{2}}{4\sinh^{2}(\frac{q_{i}-q_{j}}{2})}\,,}\end{array}$$ (1.1)

where $p_{i}$ and $q_{i}$, $i=1,...,N$ are canonically conjugated momenta and positions of particles

$$\begin{array}[]{c}\displaystyle{\{q_{i},p_{j}\}=\delta_{ij}\,,\qquad\{p_{i},p_{j}\}=\{q_{i},q_{j}\}=0}\end{array}$$ (1.2)

and $c\in\mathbb{C}$ is a coupling constant. The second type model at the level of classical mechanics is an integrable trigonometric top [14] of the Euler-Arnold type [4]. The latter means that the Hamiltonian is given as

$$\begin{array}[]{c}\displaystyle{H^{\hbox{\tiny{top}}}=\frac{1}{2}\,{\rm tr}(SJ(S))\,,\qquad S\in{\rm Mat}(N,\mathbb{C})\,,}\end{array}$$ (1.3)

where the matrix elements of $S$ are dynamical variables and $J(S)$ is some special linear functional. Together with the Poisson-Lie brackets on ${\rm gl}^{*}_{N}$ Lie coalgebra

$$\begin{array}[]{c}\displaystyle{\left\{S_{ij},S_{kl}\right\}=\frac{1}{N}\,\Big{(}S_{il}\delta_{kj}-S_{kj}\delta_{il}\Big{)}}\end{array}$$ (1.4)

the Hamiltonian (1.3) provides equations of motion in the Euler-Arnold form

$$\begin{array}[]{c}\displaystyle{\{S,H^{\hbox{\tiny{top}}}\}\equiv{\dot{S}}=[S,J(S)]\,.}\end{array}$$ (1.5)

Both models (1.1) and (1.3) are described by the Lax equations

$$\begin{array}[]{c}\displaystyle{{\dot{L}}(z)=[M(z),L(z)]\,,\quad L(z),M(z)\in{\rm Mat}(N,\mathbb{C})\,,}\end{array}$$ (1.6)

where $z\in\mathbb{C}$ is a spectral parameter. The dimension of the phase space of the integrable top depends on the choice of the coadjoint orbit, i.e. on some fixation of eigenvalues of $S$ which are the Casimir functions of (1.4). In the case of the coadjoint orbit of minimal dimension (when $N-1$ eigenvalues of $S$ coincide) the phase space has dimension $2N-2$, which is equal to the one for the Calogero-Moser-Sutherland model in the center of mass frame. In this particular case (and for some special $J(S)$) two models can be shown to be gauge equivalent, that is there exist a gauge transformation matrix $g(z)=g(z,q_{1},...,q_{N})\in{\rm Mat}(N,\mathbb{C})$, which maps one Lax matrix to another one:

$$\begin{array}[]{c}\displaystyle{L^{\hbox{\tiny{top}}}(z)=g(z)L^{\hbox{\tiny{CMS}}}(z)g^{-1}(z)\,.}\end{array}$$ (1.7)

This allows to compute explicit change of variables $S=S(p,q,\nu)$, which provides the Poisson canonical map between phase spaces of both models endowed with the Poisson brackets (1.2) and (1.4) respectively. The gauge transformation (1.7) can be considered as the classical analogue of the IRF-Vertex correspondence between dynamical and non-dynamical quantum $R$-matrices: [9]:

$$\begin{array}[]{c}\displaystyle{R^{\hbox{\tiny{Vertex}}}_{12}(\hbar,z_{1}-z_{2})=g_{2}(z_{2},q)\,g_{1}(z_{1},q-\hbar^{(2)})\,R^{\hbox{\tiny{IRF}}}_{12}(\hbar,z_{1}-z_{2}|\,q)g^{-1}_{2}(z_{2},q-\hbar^{(1)})g_{1}^{-1}(z_{1},q)\,,}\\ \\ \displaystyle{g_{1}(z_{1},q-\hbar^{(2)})=P_{2}^{-\hbar}g_{1}(z_{1},q)P_{2}^{\hbar}\,,\quad P_{2}^{\hbar}=\sum\limits_{k=1}^{N}1_{N}\otimes E_{kk}\exp(\hbar\partial_{q_{k}})\,,}\end{array}$$ (1.8)

where $g(z,q)=g(z,q_{1},...,q_{N})\in{\rm Mat}(N,\mathbb{C})$ is some special matrix (the intertwining matrix) providing the IRF-Vertex transformation. The matrix $g(z)$ from (1.8) is, in fact, exactly the one, which is used in (1.7). The non-dynamical (vertex type) $R$-matrix, by definition, satisfies the quantum Yang-Baxter equation:

$$\begin{array}[]{c}\displaystyle{R_{12}^{\hbar}R_{13}^{\hbar}R_{23}^{\hbar}=R_{23}^{\hbar}R_{13}^{\hbar}R_{12}^{\hbar}\,,\qquad R^{\hbar}_{ab}=R^{\hbox{\tiny{Vertex}}}_{ab}(\hbar,z_{a}-z_{b})\,.}\end{array}$$ (1.9)

The classical IRF-Vertex relation between many-body systems and top-like models was proposed in [17] (see also [2, 18]) in the framework of Hitchin approach to integrable systems. It was called the symplectic Hecke correspondence since it changes a certain characteristic class of underlying bundles. For the trigonometric $R$-matrices in ${\rm GL}_{N}$ case the relation (1.8) was described in [3], and the vertex-type $R$-matrix was evaluated in the special trigonometric limit starting from the elliptic $R$-matrix. The resultant $R$-matrix was called the non-standard trigonometric $R$-matrix. It generalizes the Cherednik’s 7-vertex $R$-matrix [11] in ${\rm GL}_{2}$ case. At the classical level the relation (1.8) was described in [14].

The purpose of the paper is to describe the IRF-Vertex type relation between the classical 1+1 field generalizations of the trigonometric models (1.1) and (1.3). In 1+1 case the Lax equation (1.6) turns into the Zakharov-Shabat equation written for the fields depending on the time variable $t$ and the space variable $x$:

$$\begin{array}[]{c}\displaystyle{\partial_{t}{U}(z)-k\partial_{x}{V}(z)+[{U}(z),{V}(z)]=0\,,\qquad{U}(z),{V}(z)\in{\rm Mat}(N,\mathbb{C})\,,}\end{array}$$ (1.10)

where $k\in\mathbb{C}$ is a constant parameter. The limit to the finite-dimensional mechanics (when all fields are independent of $x$) corresponds to $k\rightarrow 0$. Then the zero curvature equation (1.10) becomes the Lax equation (1.6). The field analogue for integrable finite-dimensional many-body systems was first proposed by A. Mikhailov for the Toda model [20]. The particles momenta and positions of particles become fields with the Poisson brackets

$$\begin{array}[]{c}\displaystyle{\{q_{i}(x),p_{j}(y)\}=\delta_{ij}\delta(x-y)\,,\qquad\{p_{i}(x),p_{j}(y)\}=\{q_{i}(x),q_{j}(y)\}=0\,.}\end{array}$$ (1.11)

Hereinafter we assume that all fields are periodic functions on a circle, i.e. $q_{i}(x+2\pi)=q_{i}(x)$ and similarly for all other fields. The field generalization for the Calogero-Moser-Sutherland model was introduced in [15] and [17]. In paper [1] the final explicit description was suggested for ${\rm sl}_{N}$ (i.e. $N$-body) case. Let us also mention the field analogue of the Ruijsenaars-Schneider model [24], which generalizes results of [1] to semi-discrete Zakharov-Shabat equations.

The filed generalization of the top-like models is given by the Landau-Lifshitz type models of 1-dimensional magnets [16]. In ${\rm sl}_{2}$ case equation of motion takes the form

$$\begin{array}[]{l}\displaystyle{\partial_{t}{S}=[{S},J({S})]+[{S},\partial_{x}^{2}{S}]\,,\quad S=S(t,x)\in{\rm Mat}(2,\mathbb{C})\,.}\end{array}$$ (1.12)

Integrability of the latter model was proved in [22]. For our purpose we need some higher rank ${\rm gl}_{N}$ extension of the Landau-Lifshitz equation. One possible higher rank generalization of (1.12) was proposed in [13]. In this paper we deal with another construction for ${\rm gl}_{N}$ case suggested in [6]. It is based on the associative Yang-Baxter equation [12]:

$$\begin{array}[]{c}\displaystyle{R^{\hbar}_{12}R^{\eta}_{23}=R^{\eta}_{13}R^{\hbar-\eta}_{12}+R^{\eta-\hbar}_{23}R^{\hbar}_{13}\,,\qquad R^{x}_{ab}=R^{x}_{ab}(z_{a}-z_{b})}\end{array}$$ (1.13)

written for the vertex type $R$-matrices. Namely, it was shown in [6] that the coefficients of expansion of $R$-matrix satisfying (1.13) can be used for construction of the $U$-$V$ pair with spectral parameter providing some higher rank generalizations of (1.12) through the Zakharov-Shabat equation (1.10). Classification of trigonometric ${\rm GL}_{N}$ solutions of (1.13) was suggested in [21, 23]. It includes (the properly normalized) non-standard $R$-matrix from [3]. We use this $R$-matrix and apply the construction from [6] to define the higher rank Landau-Lifshitz model. Then we show that there exists a gauge transformation $G\in{\rm Mat}(N,\mathbb{C})$, which relates it with the 1+1 Calogero-Moser-Sutherland model via the gauge transformation:

$$\begin{array}[]{l}\displaystyle{U^{\hbox{\tiny{LL}}}(z)=G(z)U^{\hbox{\tiny{2dCMS}}}(z)G^{-1}(z)+k\partial_{x}G(z)G^{-1}(z)\,,}\\ \\ \displaystyle{V^{\hbox{\tiny{LL}}}(z)=G(z)V^{\hbox{\tiny{2dCMS}}}(z)G^{-1}(z)+\partial_{t}G(z)G^{-1}(z)\,.}\end{array}$$ (1.14)

The existence of the relation (1.14) was argued in [17], and the explicit changes of variables were found in [5] for ${\rm sl}_{2}$ models (in the rational, trigonometric and elliptic cases) and in [7] for ${\rm sl}_{N}$ rational models. In this resect the aim of this paper is to extend the results of [7] to the trigonometric models.

The phenomenon of gauge equivalence is known for other models. For example, the equivalence exists between 1+1 Heisenberg magnet and the nonlinear Schrodinger equation [25]. Another example is the gauge transformation with the classical $r$-matrix structure in WZNW (and Toda) theories described in [8].

2 Antonov-Hasegawa-Zabrodin non-standard $R$-matrix

Any non-dynamical $R$-matrix in the fundamental representation of ${\rm GL}_{N}$ is represented in the form

$$\begin{array}[]{c}\displaystyle{R_{12}^{\hbar}(z_{1},z_{2})=\sum\limits_{i,j,k,l=1}^{N}R_{ij,kl}(\hbar,z_{1},z_{2})E_{ij}\otimes E_{kl}\in{\rm Mat}(N,\mathbb{C})^{\otimes 2}\,,}\end{array}$$ (2.1)

where $z_{1},z_{2}$ are spectral parameters, $\hbar$ is the Planck constant, the set $E_{ij}\in{\rm Mat}(N,\mathbb{C})$, $i,j=1,...,N$ is the basis of matrix units in ${\rm Mat}(N,\mathbb{C})$ and $R_{ij,kl}(\hbar,z_{1},z_{2})$ is a set of functions. By definition of a quantum $R$-matrix, it satisfies the quantum Yang-Baxter equation (1.9). In what follows we also assume

$$\begin{array}[]{c}\displaystyle{R_{12}^{\hbar}(z_{1},z_{2})=R_{12}^{\hbar}(z_{1}-z_{2})\,.}\end{array}$$ (2.2)

In [3] the following $R$-matrix (up to a choice of normalization coefficients) was computed by a special trigonometric limit:

$$\begin{array}[]{c}\displaystyle{R_{ij,kl}^{\eta}(z)=\delta_{ij}\delta_{kl}\delta_{ik}\frac{N}{2}(\operatorname{coth}(Nz/2)+\operatorname{coth}(N\eta/2))+}\\ \\ \displaystyle{+\delta_{ij}\delta_{kl}\varepsilon(i\neq k)\frac{Ne^{(i-k)\eta-\operatorname{sgn}(i-k)N\eta/2}}{2\sinh(N\eta/2)}+\delta_{il}\delta_{kj}\varepsilon(i\neq k)\frac{Ne^{(i-k)z-\operatorname{sgn}(i-k)Nz/2}}{2\sinh(Nz/2)}+}\\ \\ \displaystyle{+N\delta_{i+k,j+l}e^{(i-j)z+(j-k)\eta}(\varepsilon(i<j<k)-\varepsilon(k<j<i))+}\\ \\ \displaystyle{+Ne^{-N\Lambda}\delta_{i+k,j+l+N}\left(\delta_{iN}e^{-jz-l\eta}-\delta_{kN}e^{lz+j\eta}\right)\,,}\end{array}$$ (2.3)

where $\Lambda\in\mathbb{C}$ is a free constant. It is called the non-standard trigonometric $R$-matrix. Here the following notation is used:

$$\begin{array}[]{c}\displaystyle{\varepsilon(\hbox{A})=\left\{\begin{array}[]{l}1\,,\hbox{if A is true}\,,\\ 0\,,\hbox{if A is false}\,.\end{array}\right.}\end{array}$$ (2.4)

In the $N=2$ case it is the 7-th vertex $R$-matrix proposed by I. Cherednik [11]:

$$\begin{array}[]{c}R^{\hbar}(z)=\left(\begin{array}[]{cccc}\coth(z)+\coth(\hbar)&0&0&0\\ 0&\sinh^{-1}(\hbar)&\sinh^{-1}(z)&0\\ 0&\sinh^{-1}(z)&\sinh^{-1}(\hbar)&0\\ -4\,e^{-2\Lambda}\sinh(z+\hbar)&0&0&\coth(z)+\coth(\hbar)\end{array}\right)\,.\end{array}$$ (2.5)

The $R$-matrix (2.3) satisfies the skew-symmetry

$$\begin{array}[]{c}\displaystyle{R^{\hbar}_{12}(z)=-R_{21}^{-\hbar}(-z)=-P_{12}R_{12}^{-\hbar}(-z)P_{12}\,,\qquad P_{12}=\sum\limits_{i,j=1}^{N}E_{ij}\otimes E_{ji}\,,}\end{array}$$ (2.6)

where $P_{12}$ is the matrix permutation operator. Also, it satisfies the unitarity property

$$\begin{array}[]{c}\displaystyle{R^{\hbar}_{12}(z)R^{\hbar}_{21}(-z)=f^{\hbar}(z)\,\,1_{N}\otimes 1_{N}}\end{array}$$ (2.7)

with the normalization factor

$$\begin{array}[]{c}\displaystyle{f^{\hbar}(z)=\frac{N^{2}}{4}\left(\frac{1}{\sinh^{2}(N\hbar/2)}-\frac{1}{\sinh^{2}(Nz/2)}\right)}\end{array}$$ (2.8)

and the following set of local expansions:

$$\begin{array}[]{c}\displaystyle{R^{\hbar}_{12}(z)=\frac{1}{\hbar}\,1_{N}\otimes 1_{N}+r_{12}(z)+\hbar\,m_{12}(z)+O(\hbar^{2})\,,}\end{array}$$ (2.9)

$$\begin{array}[]{c}\displaystyle{R^{\hbar}_{12}(z)=\frac{1}{z}\,P_{12}+R^{\hbar,(0)}_{12}+O(z)\,,}\end{array}$$ (2.10)

$$\begin{array}[]{c}\displaystyle{R^{\hbar,(0)}_{12}=\frac{1}{\hbar}\,1_{N}\otimes 1_{N}+r^{(0)}_{12}+O(\hbar)\,,\qquad r_{12}(z)=\frac{1}{z}\,P_{12}+r^{(0)}_{12}+O(z)\,.}\end{array}$$ (2.11)

The skew-symmetry (2.6) yields

$$\begin{array}[]{c}r_{12}(z)=-r_{21}(-z)\,,\qquad m_{12}(z)=m_{21}(-z)\,,\\ \\ R^{\hbar,(0)}_{12}=-R^{-\hbar,(0)}_{21}\,,\qquad r_{12}^{(0)}=-r_{21}^{(0)}\,.\end{array}$$ (2.12)

The coefficient $r_{12}(z)$ from (2.9) is the classical $r$-matrix:

$$\begin{array}[]{c}\displaystyle{r_{ij,kl}(z)=\delta_{ij}\delta_{kl}\delta_{ik}\frac{N}{2}\operatorname{coth}(Nz/2)+}\\ \\ \displaystyle{+\delta_{ij}\delta_{kl}\varepsilon(i\neq k)\left((i-k)-\frac{N\operatorname{sgn}(i-k)}{2}\right)+\delta_{il}\delta_{kj}\varepsilon(i\neq k)\frac{Ne^{(i-k)z-\operatorname{sgn}(i-k)Nz/2}}{2\sinh(Nz/2)}+}\\ \\ \displaystyle{+Ne^{(i-j)z}\delta_{i+k,j+l}(\varepsilon(i<j<k)-\varepsilon(k<j<i))+Ne^{-N\Lambda}\delta_{i+k,j+l+N}\left(e^{-jz}\delta_{iN}-e^{lz}\delta_{kN}\right)\,.}\end{array}$$ (2.13)

The next coefficient $m_{12}(z)$ in the quasi-classical limit (2.9) is also used in what follows:

$$\begin{array}[]{c}\displaystyle{m_{ij,kl}(z)=\delta_{ij}\delta_{kl}\delta_{ik}\frac{N^{2}}{12}+\delta_{ij}\delta_{kl}\varepsilon(i\neq k)\left(\frac{(i-k)^{2}}{2}+\frac{N^{2}}{12}-\frac{N}{2}|i-k|\right)+}\\ \\ \displaystyle{+N(j-k)e^{(i-j)z}\delta_{i+k,j+l}\Big{(}\varepsilon(i<j<k)-\varepsilon(k<j<i)\Big{)}-}\\ \\ \displaystyle{-Ne^{-N\Lambda}\delta_{i+k,j+l+N}\left(le^{-jz}\delta_{iN}+je^{lz}\delta_{kN}\right)\,.}\end{array}$$ (2.14)

Some more explicit formulae are presented in the Appendix.

$R$-matrix satisfies the quantum Yang-Baxter equation (1.9). But it is even more important for our purpose that it satisfies the associative Yang-Baxter equation (1.13). This statement follows from the classification [21, 23] of trigonometric solutions of (1.13). See also [14] for a review.

3 Calogero-Moser-Sutherland model and trigonometric top

Here we briefly review the relation between the finite-dimensional Calogero-Moser-Sutherland model (1.1) and the trigonometric top. More details can be found in [14], where this relation was described at the level of relativistic models, i.e. between the trigonometric Ruijsenaars-Schneider model and the relativistic (the one governed by quadratic $r$-matrix structure) trigonometric top.

The Calogero-Moser-Sutherland model.

The Lax pair has the form:

$$\begin{array}[]{c}\displaystyle{L^{\hbox{\tiny{CSM}}}_{ij}(z)=}\\ \\ \displaystyle{=\delta_{ij}\left[-p_{i}+\frac{1}{N}\sum\limits_{k=1}^{N}p_{k}+\frac{c}{2}\coth\left(\frac{Nz}{2}\right)\right]+\left(1-\delta_{ij}\right)\frac{c}{2}\left[\coth\left(\frac{q_{j}-q_{i}}{2}\right)+\coth\left(\frac{Nz}{2}\right)\right]}\end{array}$$ (3.1)

and

$$\begin{array}[]{c}\displaystyle{M^{\hbox{\tiny{CSM}}}_{ij}(z)=-\delta_{ij}\sum\limits_{k:k\neq i}\frac{c}{4\sinh^{2}(\frac{q_{i}-q_{k}}{2})}+\left(1-\delta_{ij}\right)\frac{c}{4\sinh^{2}(\frac{q_{i}-q_{j}}{2})}\,.}\end{array}$$ (3.2)

The equations of motion generated by the Hamiltonian (1.1) are represented in the Lax form (1.6). Hereinafter we work in the center of mass frame, i.e.

$$\begin{array}[]{c}\displaystyle{\sum\limits_{k=1}^{N}q_{k}=0\,.}\end{array}$$ (3.3)

Trigonometric top.

It was proved in [19] that for any $R$-matrix satisfying the associative Yang-Baxter equation (1.13) and the properties (2.6)-(2.12) one can construct an integrable top-like model of type (1.3)-(1.5). Namely, define the Lax pair as follows:

$$\begin{array}[]{c}\displaystyle{L^{\hbox{\tiny{top}}}(z,S)={\rm tr}_{2}\left(r_{12}(z)S_{2}\right),}\end{array}$$ (3.4)

$$\begin{array}[]{c}\displaystyle{M^{\hbox{\tiny{top}}}(z,S)={\rm tr}_{2}\left(m_{12}(z)S_{2}\right),}\end{array}$$ (3.5)

where ${{\rm tr}}_{2}$ means the trace over the second tensor component and $S_{2}=1_{N}\otimes S$. That is, having $r$-matrix in the form

$$\begin{array}[]{c}\displaystyle{r_{12}(z)=\sum\limits_{i,j,k,l=1}^{N}r_{ij,kl}(z)E_{ij}\otimes E_{kl}\in{\rm Mat}(N,\mathbb{C})^{\otimes 2}}\end{array}$$ (3.6)

and using ${\rm tr}(E_{kl}S)=S_{lk}$, from (3.4) one gets

$$\begin{array}[]{c}\displaystyle{L^{\hbox{\tiny{top}}}(z)=\sum\limits_{i,j,k,l=1}^{N}r_{ij,kl}(z)S_{lk}E_{ij}\in{\rm Mat}(N,\mathbb{C})\,.}\end{array}$$ (3.7)

Then the Lax equation (1.6) with the Lax pair (3.4)-(3.5) provides the Euler-Arnold equations (1.5) with $J(S)$ defined as

$$\begin{array}[]{c}\displaystyle{J(S)={\rm tr}_{2}\left(m_{12}(0)S_{2}\right).}\end{array}$$ (3.8)

Explicit expressions for (3.4)-(3.5) and (3.8) coming from the $R$-matrix (2.3) can be found in the Appendix.

Gauge equivalence.

It was mentioned in [14] that in order to establish the gauge equivalence (1.7) one should fix the constant $\Lambda$ entering (2.3), (2.13)-(2.14) and the formulae from the Appendix as

$$\begin{array}[]{c}\displaystyle{\Lambda=\imath\pi\,.}\end{array}$$ (3.9)

Introduce the matrix

$$\begin{array}[]{c}\displaystyle{g(z,q)=\Xi(z,q)D^{-1}(q)\in{\rm Mat}(N,\mathbb{C})\,,}\end{array}$$ (3.10)

where $\Xi(z,q)$ is the Vandermonde type matrix (except the last row)

$$\begin{array}[]{c}\displaystyle{\Xi_{ij}(z,q)=e^{(i-1)\left(z+q_{j}\right)}+(-1)^{N}e^{-\left(z+q_{j}\right)}\delta_{iN}\,,\qquad\sum\limits_{k=1}^{N}q_{k}=0}\end{array}$$ (3.11)

and $D$ is the following diagonal matrix:

$$\begin{array}[]{c}\displaystyle{D_{ij}(q)=\delta_{ij}\prod\limits_{k:k\neq i}^{N}\Big{(}e^{q_{i}}-e^{q_{k}}\Big{)}\,.}\end{array}$$ (3.12)

The matrix $g(z,q)$ is the intertwining matrix entering the IRF-Vertex relation (1.8) with the non-standard trigonometric $R$-matrix (2.3) in the r.h.s of (1.8), see [3].

Then the gauge equivalence (1.7) holds true, where on the ”top” side we deal with the model corresponding to the minimal coadjoint orbit. The latter means that the matrix $S$ is of rank one.

Then the relation (1.7) is valid with the following change of variables:

$$\begin{array}[]{c}\displaystyle{S_{ij}(p,q,c)=\frac{(-1)^{j}\sigma_{j}(e^{q})}{N}\sum_{m=1}^{N}\frac{{\tilde{p}}_{m}\Big{(}e^{(i-1)q_{m}}+(-1)^{N}\delta_{iN}e^{-q_{m}}\Big{)}-cN\delta_{iN}(-1)^{N}e^{-q_{m}}}{\prod\limits^{N}_{l:\,l\neq m}\Big{(}e^{q_{m}}-e^{q_{l}}\Big{)}}\,,}\end{array}$$ (3.13)

where

$$\begin{array}[]{c}\displaystyle{{\tilde{p}}_{m}=-p_{m}+(i-1)c-\frac{c}{2}\sum\limits_{l:\,l\neq m}^{N}\coth\Big{(}\frac{q_{m}-q_{l}}{2}\Big{)}}\end{array}$$ (3.14)

and $\sigma_{j}(e^{q})=\sigma_{j}(e^{q_{1}},...,e^{q_{N}})$ are elementary symmetric functions of the variables $e^{q_{i}}$, that is

$$\begin{array}[]{c}\displaystyle{\prod_{k=1}^{N}\left(\zeta-e^{x_{k}}\right)=\sum_{k=0}^{N}(-1)^{k}\zeta^{k}\sigma_{k}(e^{x})\,.}\end{array}$$ (3.15)

It is easy to see from (3.13) that the matrix $S$ is indeed a matrix of rank one. The set of its eigenvalues is as follows:

$$\begin{array}[]{c}\displaystyle{{\rm Spec}(S)=(0,...,0,c)\,,}\end{array}$$ (3.16)

so that ${\rm tr}(S)=c$ and

$$\begin{array}[]{c}\displaystyle{S^{2}=cS\,.}\end{array}$$ (3.17)

One can show that the Poisson brackets $\{S_{ij}(p,q,c),S_{kl}(p,q,c)\}$ between elements of the matrix (3.13) computed by means of the canonical brackets (1.2) reproduce the linear Poisson-Lie brackets (1.4), that is the map between two models is Poisson (canonical).

Similar results are known for the rational and elliptic models including relativistic models, where on the many-body side one deals with the Ruijsenaars-Schneider model, and the relativistic top is described by quadratic Poisson brackets of Sklyanin type. See [17, 2, 18, 14, 24] for details.

4 1+1 Landau-Lifshitz equations from $R$-matrices

In the field case the dynamical variables are again arranged into matrix

$$\begin{array}[]{c}\displaystyle{S=\sum\limits_{i,j=1}^{N}E_{ij}S_{ij}\,,\quad S_{ij}=S_{ij}(t,x)\,.}\end{array}$$ (4.1)

In what follows we use short notations for partial derivatives:

$$\begin{array}[]{c}\displaystyle{\partial_{x}\varphi(t,x)=\varphi_{x}\,,\quad\partial_{t}\varphi(t,x)=\varphi_{t}\,.}\end{array}$$ (4.2)

Having solution of the associative Yang-Baxter equation (1.13), which satisfies also the properties (2.6)-(2.12), one can derive a wide set of identities for the coefficients of expansions (2.9)-(2.11). The first example (see [19]) is given by the relation

$$\begin{array}[]{c}\displaystyle{[m_{13}(z),r_{12}(z)]=[r_{12}(z),m_{23}(0)]-[\partial_{z}m_{12}(z),P_{23}]+[m_{12}(z),r_{23}^{(0)}]+[m_{13}(z),r_{23}^{(0)}]\,.}\end{array}$$ (4.3)

It leads to the construction of the Lax pair for the integrable top (3.4)-(3.8). The next identity is of the form:

$$\begin{array}[]{c}\displaystyle{r_{12}(z)r_{13}(z)=r_{23}^{(0)}r_{12}(z)-r_{13}(z)r_{23}^{(0)}-\partial_{z}r_{13}(z)P_{23}+m_{12}(z)+m_{23}(0)+m_{13}(z)\,.}\end{array}$$ (4.4)

Using (4.3)-(4.4) the following construction of $U$–$V$ pair was suggested in [6]. Consider the pair of matrices:

$$\begin{array}[]{c}\displaystyle{U^{\hbox{\tiny{LL}}}(z)=L(S,z)={\rm tr}_{2}\Big{(}r_{12}(z)S_{2}\Big{)}\in{\rm Mat}(N,\mathbb{C})}\end{array}$$ (4.5)

and

$$\begin{array}[]{c}\displaystyle{V^{\hbox{\tiny{LL}}}(z)=V_{1}(z)-cV_{2}(z)-\frac{Nc^{2}}{12}\,1_{N}\in{\rm Mat}(N,\mathbb{C})\,,}\end{array}$$ (4.6)

where

$$\begin{array}[]{c}\displaystyle{V_{1}(z)=-c\partial_{z}L(S,z)+L(E(S)S,z)\,,}\\ \\ \displaystyle{V_{2}(z)=L(T,z)\,.}\end{array}$$ (4.7)

Here we use the notation

$$\begin{array}[]{c}\displaystyle{E(S)={\rm tr}_{2}\left(r_{12}^{(0)}S_{2}\right)\in{\rm Mat}(N,\mathbb{C})\,,}\end{array}$$ (4.8)

see (A.12)-(A.14) in the Appendix. The matrix $T$ entering the definition of $V_{2}(z)$ is defined as

$$\begin{array}[]{c}\displaystyle{T=-\frac{k}{c^{2}}\,[S,S_{x}]\,,\quad S_{x}=\partial_{x}S\,.}\end{array}$$ (4.9)

It solves equation

$$\begin{array}[]{c}\displaystyle{-kS_{x}=[S,T]}\end{array}$$ (4.10)

if $S$ satisfies the condition (3.17). In what follows we assume that $S$ is a rank one matrix, so that

$$\begin{array}[]{c}\displaystyle{S^{2}={\rm tr}(S)S\,,\quad c={\rm tr}(S)}\end{array}$$ (4.11)

in (3.17). In this case a set of additional relations appear. In particular, $SE(S)=0$ (see [6]). Finally, the Zakharov-Shabat equation is equivalent to the following Landau-Lifshitz type equation:

$$\begin{array}[]{c}\displaystyle{\partial_{t}S=2c[S,J(S)]+\frac{k^{2}}{c}[S,S_{xx}]-2k[S,E(S_{x})]\,.}\end{array}$$ (4.12)

In the $N=2$ case $E(S)=0$ and the last term vanishes.

Equation (4.12) is Hamiltonian. The Poisson brackets are defined similarly to (1.4) but for the loop algebra:

$$\begin{array}[]{c}\displaystyle{\left\{S_{ij}(x),S_{kl}(y)\right\}=\frac{1}{N}\,\Big{(}S_{il}\delta_{kj}-S_{kj}\delta_{il}\Big{)}\delta(x-y)\,.}\end{array}$$ (4.13)

The Hamiltonian

$$\begin{array}[]{c}\displaystyle{{\mathcal{H}}^{\hbox{\tiny{LL}}}=\oint dxH^{\hbox{\tiny{LL}}}(x)}\end{array}$$ (4.14)

with the density

$$\begin{array}[]{c}\displaystyle{H^{\hbox{\tiny{LL}}}(x)=Nc\operatorname{tr}(SJ(S))-\frac{Nk^{2}}{2c}\operatorname{tr}\left(\partial_{x}S\partial_{x}S\right)+Nk\operatorname{tr}\left(\partial_{x}SE(S)\right)}\end{array}$$ (4.15)

provides the equation (4.12) through the equation $\partial_{t}S=\{S,{\mathcal{H}}^{\hbox{\tiny{LL}}}\}$ and the brackets (4.13).

5 1+1 Calogero-Moser-Sutherland field theory

Following [1] let us introduce the $U$–$V$ pair with spectral parameter for 1+1 field generalization of the Calogero-Moser-Sutherland model. Introduce the set of functions

$$\begin{array}[]{c}\displaystyle{\alpha_{i}^{2}=kq_{ix}-c\,,\quad i=1,...,N\,.}\end{array}$$ (5.1)

and a function

$$\begin{array}[]{c}\displaystyle{\kappa=-\frac{1}{Nc}\sum\limits_{l=1}^{N}p_{l}\left(c-kq_{lx}\right)=\frac{1}{Nc}\sum\limits_{l=1}^{N}p_{l}\alpha_{l}^{2}\,.}\end{array}$$ (5.2)

The $U$–$V$ pair is as follows:

$$\begin{array}[]{c}\displaystyle{U^{\hbox{\tiny{2dCMS}}}_{ij}=\delta_{ij}\left[-p_{i}+\frac{1}{N}\sum_{k=1}^{N}p_{k}-\frac{\alpha_{i}^{2}}{2}\coth\left(\frac{Nz}{2}\right)-\frac{k\alpha_{ix}}{\alpha_{i}}\right]-}\\ \\ \displaystyle{-\left(1-\delta_{ij}\right)\left[\coth\left(\frac{q_{j}-q_{i}}{2}\right)+\coth\left(\frac{Nz}{2}\right)\right]\frac{\alpha_{j}^{2}}{2}}\end{array}$$ (5.3)

and

$$\begin{array}[]{c}\displaystyle{V^{\hbox{\tiny{2dCMS}}}_{ij}=\delta_{ij}\left[-\frac{q_{it}}{2}\coth\left(\frac{Nz}{2}\right)-\frac{Nc\alpha_{i}^{2}}{4}\left(\frac{1}{\sinh^{2}(\frac{Nz}{2})}+\frac{1}{3}\right)-\frac{\alpha_{it}}{\alpha_{i}}+\widetilde{m}_{i}^{0}-\frac{1}{N}\sum_{l=1}^{N}\tilde{m}_{l}^{0}\right]-}\\ \\ \displaystyle{-\left(1-\delta_{ij}\right)\Bigg{[}\frac{Nc}{2}\coth\left(\frac{Nz}{2}\right)\left(\coth\left(\frac{q_{j}-q_{i}}{2}\right)+\coth\left(\frac{Nz}{2}\right)\right)+\frac{Nc}{2\sinh^{2}\left(\frac{q_{i}-q_{j}}{2}\right)}-}\\ \\ \displaystyle{-\widetilde{m}_{ij}\left(\coth\left(\frac{q_{j}-q_{i}}{2}\right)+\coth\left(\frac{Nz}{2}\right)\right)\Bigg{]}\frac{\alpha_{j}^{2}}{2}\,,}\end{array}$$ (5.4)

where the set of functions $\tilde{m}_{i}^{0}$ in the diagonal part of $V$-matrix is of the form:

$$\begin{array}[]{c}\displaystyle{\tilde{m}_{i}^{0}=p_{i}^{2}+\frac{k^{2}\alpha_{ixx}}{\alpha_{i}}+2\kappa p_{i}-\sum_{l:l\neq i}^{N}\Bigg{[}\frac{\left(2\alpha_{l}^{4}+\alpha_{i}^{2}\alpha_{l}^{2}\right)}{4}\left(\frac{1}{\text{sinh}^{2}\left(\frac{q_{i}-q_{l}}{2}\right)}+\frac{1}{3}\right)+}\end{array}$$ (5.5)

$$+2k\alpha_{l}\alpha_{lx}\left(\coth\left(\frac{q_{i}-q_{l}}{2}\right)-\frac{(q_{i}-q_{l})}{6}\right)\Bigg{]}\,.$$

The off-diagonal part of $V$-matrix contains the following expressions:

$$\begin{array}[]{c}\displaystyle{\widetilde{m}_{ij}=p_{i}+p_{j}+2\kappa+}\\ \\ \displaystyle{+\frac{k\alpha_{ix}}{\alpha_{i}}-\frac{k\alpha_{jx}}{\alpha_{j}}-\sum\limits_{n:n\neq i,j}^{N}\frac{\alpha_{n}^{2}}{2}\left(\coth\left(\frac{q_{i}-q_{n}}{2}\right)+\coth\left(\frac{q_{n}-q_{j}}{2}\right)-\coth\left(\frac{q_{i}-q_{j}}{2}\right)\right)\,.}\end{array}$$ (5.6)

In the above formulae the notations (4.2) are used. The statement is that the Zakharov-Shabat equation (1.10) provides the following set of equations of motion:

$$\begin{array}[]{c}\displaystyle{\partial_{t}q_{i}=2p_{i}\left(c-kq_{ix}\right)-\frac{2}{Nc}\sum_{l=1}^{N}p_{l}\left(c-kq_{lx}\right)\left(c-kq_{ix}\right)\,,\quad i=1,...,N}\end{array}$$ (5.7)

and

$$\begin{array}[]{c}\displaystyle{\partial_{t}p_{i}=-2kp_{i}p_{ix}+\frac{2k}{Nc}\left\{\sum_{l=1}^{N}p_{i}p_{l}\left(c-kq_{lx}\right)\right\}_{x}+k\left\{\frac{k^{3}q_{ixxx}}{2\left(c-kq_{ix}\right)}+\frac{k^{4}q_{ixx}^{2}}{4\left(c-kq_{ix}\right)^{2}}\right\}_{x}+}\\ \\ \displaystyle{+\sum\limits_{l:l\neq i}^{N}\Bigg{[}k^{3}q_{lxxx}\left(\coth\left(\frac{\left(q_{i}-q_{l}\right)}{2}\right)-\frac{\left(q_{i}-q_{l}\right)}{6}\right)-\frac{3k^{2}q_{lxx}\left(c-kq_{lx}\right)}{2\text{sinh}^{2}\left(\frac{q_{i}-q_{l}}{2}\right)}-}\\ \\ \displaystyle{-\frac{3k^{2}q_{lxx}\left(c-kq_{lx}\right)}{6}-\frac{\left(c-kq_{lx}\right)^{3}\coth\left(\frac{\left(q_{i}-q_{l}\right)}{2}\right)}{2~{}\text{sinh}^{2}\left(\frac{q_{i}-q_{l}}{2}\right)}\Bigg{]}\,,\quad i=1,...,N\,.}\end{array}$$ (5.8)

The Hamiltonian description is given by

$$\begin{array}[]{c}\displaystyle{{\mathcal{H}}^{\hbox{\tiny{2dCMS}}}=\oint dxH^{\hbox{\tiny{2dCMS}}}(x)}\end{array}$$ (5.9)

with the density

$$\begin{array}[]{c}\displaystyle{H^{\hbox{\tiny{2dCMS}}}(x)=\sum_{i=1}^{N}p_{i}^{2}\left(c-kq_{ix}\right)-\frac{1}{Nc}\left(\sum_{i=1}^{N}p_{i}\left(c-kq_{ix}\right)\right)^{2}-}\\ \\ \displaystyle{-\sum_{i=1}^{N}\frac{k^{4}q_{ixx}^{2}}{4\left(c-kq_{ix}\right)}+\frac{k^{3}}{4}\sum_{i\neq j}^{N}\Big{[}q_{ix}q_{jxx}-q_{jx}q_{ixx}\Big{]}\left(\coth\left(\frac{q_{i}-q_{j}}{2}\right)-\frac{\left(q_{i}-q_{j}\right)}{6}\right)-}\\ \\ \displaystyle{-\frac{1}{8}\sum_{i\neq j}^{N}\Bigg{[}\left(c-kq_{ix}\right)^{2}\left(c-kq_{jx}\right)+\left(c-kq_{ix}\right)\left(c-kq_{jx}\right)^{2}-}\\ \\ \displaystyle{-ck^{2}\left(q_{ix}-q_{jx}\right)^{2}\Bigg{]}\Bigg{(}\frac{1}{\text{sinh}^{2}\left(\frac{q_{i}-q_{j}}{2}\right)}+\frac{1}{3}\Bigg{)}\,.}\end{array}$$ (5.10)

Equations of motion (5.7)-(5.8) are reproduced as Hamiltonian equations $\partial_{t}f=\{f,{\mathcal{H}}^{\hbox{\tiny{2dCMS}}}\}$ with the Poisson brackets (1.11).

It is important to mention that the above formulae are valid in the center of mass frame, i.e. the condition (3.3) holds true in the field case as well. At the same time the sum of momenta is not equal to zero.

6 Gauge equivalence and change of variables

Introduce the function

$$\begin{array}[]{c}\displaystyle{b(x,t)=\prod_{i<j}^{N}(e^{q_{i}}-e^{q_{j}})^{1/N}\prod_{l=1}^{N}\left(kq_{l,x}-c\right)^{1/(2N)}}\end{array}$$ (6.1)

and the matrix

$$\begin{array}[]{c}\displaystyle{G(z,q)=b(x,t)g(z,q)=b(x,t)\Xi(z,{q})D^{-1}(q)\in{\rm Mat}(N,\mathbb{C})\,.}\end{array}$$ (6.2)

Then the gauge equivalence (1.14) holds true relating the $U$–$V$ pairs (5.3)-(5.4) and (4.6)-(4.7) with the following change of variables:

$$\begin{array}[]{c}\displaystyle{S_{ij}(x)=\frac{(-1)^{j}\sigma_{j}(e^{q})}{N}\sum_{m=1}^{N}\frac{P_{m}\Big{(}e^{(i-1)q_{m}}+(-1)^{N}\delta_{iN}e^{-q_{m}}\Big{)}+N\alpha_{m}^{2}(-1)^{N}\delta_{iN}e^{-q_{m}}}{\prod\limits_{l:\,l\neq m}^{N}\left(e^{q_{m}}-e^{q_{l}}\right)}\,,}\end{array}$$ (6.3)

where the notations (3.15) for the elementary symmetric functions are used and

$$\begin{array}[]{c}\displaystyle{P_{m}=-p_{m}-\frac{k\alpha_{mx}}{\alpha_{m}}-(i-1)\alpha_{m}^{2}+\frac{(N-2)}{2}kq_{m,x}+\frac{\alpha_{m}^{2}}{2}\sum_{l\,:l\neq m}^{N}\coth\left(\frac{q_{m}-q_{l}}{2}\right)}\end{array}$$ (6.4)

for $m=1,...,N$.

Similarly to the finite-dimensional case the set of eigenvalues of the matrix $S$ has the form

$$\begin{array}[]{c}\displaystyle{Spec(S)=(0,...,0,c)}\end{array}$$ (6.5)

and

$$\begin{array}[]{c}\displaystyle{{\rm tr}(S)=c,~{}~{}~{}S^{2}=cS\,.}\end{array}$$ (6.6)

The Poisson brackets $\{S_{ij}(x),S_{kl}(y)\}$ being computed for the expressions (6.3) by means of the canonical brackets (1.11), provide the linear Poisson brackets (4.13). That is the map between two models is the Poisson map.

The Hamiltonians of two models (4.14) and (5.9) coincide with the change of variables (6.3):

$$\begin{array}[]{c}\displaystyle{{\mathcal{H}}^{\hbox{\tiny{2dCMS}}}={\mathcal{H}}^{\hbox{\tiny{LL}}}\,.}\end{array}$$ (6.7)

In fact, the densities of the Hamiltonians (4.15) and (5.10) also coincide up to a constant:

$$\begin{array}[]{c}\displaystyle{H^{\hbox{\tiny{2dCMS}}}(x)=H^{\hbox{\tiny{LL}}}(x)-\frac{N^{2}c^{3}}{12}\,.}\end{array}$$ (6.8)

The above statements are verified by technically complicated but straightforward calculations, which are similar to the finite-dimensional case discussed in [14] and [2]. The coincidence of the Hamiltonians follows from the gauge equivalence since the Hamiltonians are generated by traces of powers of monodromy of the connection $\nabla_{x}=k\partial_{x}-U$. The canonicity of the map between models can be verified directly. For this purpose one can represent the matrix $S$ in the form $S_{ij}=(1/N)a_{i}b_{j}$ and show that $\{b_{i}(x),a_{j}(y)\}=\delta_{ij}\delta(x-y)$ (for $i,j=1,...N-1$) follows from the canonical brackets (1.11). In this way one easily come to (4.13). This calculation is also similar to the finite-dimensional case, see e.g. [2] in the rational case. Details of the proof will be given elsewhere.

7 Appendix

Here we collect some useful explicit expressions entering expansions of the non-standard $R$-matrix, Lax equations and equations of motion of the trigonometric top and the corresponding Landau-Lifshitz model. We emphasize that the below given formulae contain dependence on the arbitrary constant $\Lambda$ entering $R$-matrix (2.3). In order to use these formulae for relation to the Calogero-Moser-Sutherland model one should fix $\Lambda=\imath\pi$ as in (3.9).

Expansions of the non-standard $R$-matrix.

We begin with the explicit expressions for $m_{12}(0)$ and $r_{12}^{(0)}$:

$$\begin{array}[]{c}\displaystyle{m_{ij,kl}(0)=\delta_{ij}\delta_{kl}\delta_{ik}\frac{N^{2}}{12}+\delta_{ij}\delta_{kl}\varepsilon(i\neq k)\left(\frac{(i-k)^{2}}{2}+\frac{N^{2}}{12}-\frac{N}{2}|i-k|\right)+}\\ \\ \displaystyle{+N(j-k)\delta_{i+k,j+l}\Big{(}\varepsilon(i<j<k)-\varepsilon(k<j<i)\Big{)}-Ne^{-N\Lambda}\delta_{i+k,j+l+N}\left(l\delta_{iN}+j\delta_{kN}\right)}\end{array}$$ (A.1)

and

$$\begin{array}[]{c}\displaystyle{r_{ij,kl}^{(0)}=\left(\delta_{ij}\delta_{kl}\varepsilon(i\neq k)+\delta_{il}\delta_{kj}\varepsilon(i\neq k)\right)\left((i-k)-\frac{N\operatorname{sgn}(i-k)}{2}\right)+}\\ \\ \displaystyle{+N\delta_{i+k,j+l}\Big{(}\varepsilon(i<j<k)-\varepsilon(k<j<i)\Big{)}+Ne^{-N\Lambda}\delta_{i+k,j+l+N}\left(\delta_{iN}-\delta_{kN}\right)\,.}\end{array}$$ (A.2)

$L$-matrix.

The Lax matrix $L^{\hbox{\tiny{top}}}(z,S)$ following from (3.4) has the following form. Its diagonal part:

$$\begin{array}[]{c}\displaystyle{L_{ii}^{\hbox{\tiny{top}}}(z,S)=\sum_{k=1}^{i-1}\frac{(2i-2k-N)}{2}S_{kk}+\sum_{k=i+1}^{N}\frac{(2i-2k+N)}{2}S_{kk}+\frac{NS_{ii}}{2}\coth\left(\frac{Nz}{2}\right)\,.}\end{array}$$ (A.3)

Next, for $i<j$:

$$\begin{array}[]{c}\displaystyle{L_{ij}^{\hbox{\tiny{top}}}(z,S)=N\sum\limits_{k=j+1}^{N}e^{z(i-j)}S_{i-j+k,k}+\frac{NS_{ij}e^{z(i-j)+\frac{Nz}{2}}}{2\sinh(Nz/2)}}\end{array}$$ (A.4)

and for $i>j$ we have:

$$\begin{array}[]{c}\displaystyle{L_{ij}^{\hbox{\tiny{top}}}(z)=Ne^{-N\Lambda}\delta_{iN}\sum_{k=j+1}^{N}e^{-jz}S_{k-j,k}-}\\ \\ \displaystyle{-N\sum_{k=1}^{j-1}e^{z(i-j)}S_{i-j+k,k}-NS_{i-j,N}e^{z(i-j)-\Lambda N}+\frac{NS_{ij}e^{z(i-j)-\frac{Nz}{2}}}{2\sinh(Nz/2)}\,.}\end{array}$$ (A.5)

$M$-matrix.

The $M$-matrix defined through (3.5) and (2.14) has the following explicit form. The diagonal part:

$$\begin{array}[]{c}\displaystyle{M_{ii}^{\hbox{\tiny{top}}}(z,S)=\sum_{k=1}^{i-1}\frac{1}{12}S_{kk}\left(6i^{2}-12ik-6iN+6k^{2}+6kN+N^{2}\right)+}\\ \\ \displaystyle{+\sum_{k=i+1}^{N}\frac{1}{12}S_{kk}\left(6i^{2}-12ik+6iN+6k^{2}-6kN+N^{2}\right)+\frac{1}{12}N^{2}S_{ii}\,.}\end{array}$$ (A.6)

For $i<j$:

$$\begin{array}[]{c}\displaystyle{M_{ij}^{\hbox{\tiny{top}}}(z,S)=N\sum_{k=j+1}^{N}(j-k)e^{z(i-j)}S_{i-j+k,k}}\end{array}$$ (A.7)

and for $i>j$:

$$\begin{array}[]{c}\displaystyle{M_{ij}^{\hbox{\tiny{top}}}(z,S)=Ne^{-N\Lambda}\delta_{iN}\sum_{k=j+1}^{N}(j-k)e^{-jz}S_{k-j,k}-N\sum_{k=1}^{j-1}(j-k)e^{z(i-j)}S_{i-j+k,k}-}\\ \\ \displaystyle{-NjS_{i-j,N}e^{z(i-j)-\Lambda N}\,.}\end{array}$$ (A.8)

$J(S)$-matrix.

The matrix $J(S)$ is defined through (3.8) and (A.1). Its diagonal part is as follows:

$$\begin{array}[]{c}\displaystyle{J_{ii}(S)=\sum_{k=1}^{i-1}\frac{1}{12}S_{kk}\left(6i^{2}-12ik-6iN+6k^{2}+6kN+N^{2}\right)+}\\ \\ \displaystyle{+\sum_{k=i+1}^{N}\frac{1}{12}S_{kk}\left(6i^{2}-12ik+6iN+6k^{2}-6kN+N^{2}\right)+\frac{1}{12}N^{2}S_{ii}\,.}\end{array}$$ (A.9)

For $i<j$ :

$$\begin{array}[]{c}\displaystyle{J_{ij}(S)=N\sum_{k=j+1}^{N}(j-k)S_{i-j+k,k}}\end{array}$$ (A.10)

and for $i>j$ :

$$\begin{array}[]{c}\displaystyle{J_{ij}(S)=Ne^{-N\Lambda}\delta_{iN}\sum_{k=j+1}^{N}(j-k)S_{k-j,k}-N\sum_{k=1}^{j-1}(j-k)S_{i-j+k,k}-jNe^{-N\Lambda}S_{i-j,N}\,.}\end{array}$$ (A.11)

$E(S)$-matrix.

The matrix $E(S)$ is defined as in (4.8) with $r_{12}^{(0)}$ given in (A.2). Its diagonal part has the form:

$$\begin{array}[]{c}\displaystyle{E_{ii}(S)=\sum_{k=1}^{i-1}\frac{1}{2}S_{kk}(2i-2k-N)+\sum_{k=i+1}^{N}\frac{1}{2}S_{kk}(2i-2k+N)\,.}\end{array}$$ (A.12)

The off-diagonal part with $i<j$ we have:

$$\begin{array}[]{c}\displaystyle{E_{ij}(S)=N\sum_{k=j+1}^{N}S_{i-j+k,k}+\left(i-j+\frac{N}{2}\right)S_{ij}}\end{array}$$ (A.13)

and finally for $i>j$:

$$\begin{array}[]{c}\displaystyle{E_{ij}(S)=Ne^{-N\Lambda}\delta_{iN}\sum_{k=j+1}^{N}S_{k-j,k}-N\sum_{k=1}^{j-1}S_{i-j+k,k}-Ne^{-N\Lambda}S_{i-j,N}+\left(i-j-\frac{N}{2}\right)S_{ij}\,.}\end{array}$$ (A.14)