Abstract We present universal construction for the Calogero-Moser system with two types spins interaction of trigonometric potential based on the root system of semi-simple Lie algebra. In this formalism, we successfully build up the correct Lax pair as well as the $R$-matrix for this generalized Calogero-Moser models. Moreover using the property of root system, we make a concise explanation that in the quantized model, the $R$-matrix takes the same form as the classical one, which is the main new result of this paper. Keywords：root system;spin Calogero-Moser system;Lax pair;$R$-matrix;quantized model
The Calogero-Moser (CM) models are a collection of exactly integrable or solvable one-dimensional dynamical systems characterized by long-range pairwise interaction with the potential either given by rational, trigonometric, hyperbolic or elliptic function.[1] It turns out that both classical and quantum versions of these systems are connected with various branches of mathematical physics ranging from random matrix theory to two-dimensional gravity and play an important role in several areas of physical applications such as quantum Hall effect, solid state physics, and particle physics. Some new aspects have been demonstrated extensively in Refs.[2-3] that there exists profound interrelations between CM integrable models and Seiberg-Witten curve and differential in the framework of four-dimensional $\mathcal{N}=2$ supersymmetric Yang-Mills theory. In addition, a more delicate picture of relationships between integrable models and gauge theories in higher dimensions including five and six can be found in Refs.[4-5].
The big prominent advance was made by Bordner[6] who proposed a beautiful and universal formulation of Lax representation or generally called the root type Lax pair for different kinds of models of Calogero-Moser systems, which was closely related to the basis of the root systems of semi-simple Lie algebras.[7] This technical approach for CM integrable systems exhibits a novel feature that it provides us with a simple and effective method of deriving sufficiently many number of conserved quantities for all models in a unified description and allowing to prove in some cases the algebraic integrability of the relevant systems both at classical and quantum level. More importantly, the significance of this construction lies in its applications cover not only the usual Calogero-Moser systems but also the “twisted models" adopted by D'Hoker and Phong and their generalizations, namely “extended twisted models".[8-9] Furthermore, as pointed out by Takasaki,[10] a multi-spin generalization of this construction is also possible which is nothing but the genus-one case of Levin and Olshanetsky' framework.
Spin generalization is a direct and natural extension of the Calogero-Moser many-body system coupled to additional internal degrees of freedom with properties of the integrability related to the Dunkl operators and the presence of an algebra of symmetries which gives rise to a representation of the Yangian Y(gls).[11] Historically, the simplest sl$(n)$ spin generalization was first investigated by Krichever through introducing creation and annihilation analogue operators[12] and they discovered that the spin generalization just like the spinless case via the pole dynamics of the matrix KP hierarchy.[13] These models are slightly different from the ordinary one discussed by Bordner due to such spin generalization is determined by the semi-simple Lie algebras themselves rather than their root systems. Applying the abstract language of symplectic geometry, the classical spin variables can be formulated in the dual space of Lie algebra $\mathfrak{g}$ or a coadjoint orbit of the action of Lie group on the dual space of its Lie algebra. Meanwhile the Poisson brackets of spin variables are defined in terms of Kostant-Kirillov Poisson structure on the space of functions on $\mathfrak{g}^{*}$, which take the similar form as the Lie brackets of the Lie algebra basis.[14] This point of view is mainly inspired by the work of Kirillov in the study of coadjoint representation theory and later it also suggests us an insight and powerful tool for the studies of equivariant cohomology, localization principle,[15] and geometric quantization[16] and so on.
Along this way, Feher and Pusztai[17] introduced a new integrable model with two types of the interacting spin variables from the Hamiltonian reductions of the free geodesic motion on a non-compact simple Lie group $G$, which is analogous of the classical Calogero-Sutherland system. Their ideas are based mainly on the $G_{+}\times G_{+}$ symmetry of the extension of the symplectic form and Hamiltonian of the geodesic system defined on the cotangent bundle $T^{\ast}G$ and if one chooses generic values of the momentum map, we will get the desired new spin model. Recently, an intuitive geometrical picture was provided by Kharche, Levin, Olshanetsky, and Zotov[18] who gave an elegant Hitchin type description of this new integrable system. Roughly speaking, they constructed the Lax operators as the Higgs fields defined over a singular rational curve as well as the classical $R$-matrix depending on the spectral parameter. By means of the basis of meromorphic differentials on the curve, they obtained hierarchy of independent integrals of motion and proved that the number of these integrals is equal to the half of dimension of the fixed point set or usually called phase space.
In this paper we first review the basic knowledge of semi-simple Lie algebras and the associated root systems,some definitions and useful notations are given. Then the basic ingredients of the CM system with two types interacting spins via a trigonometric potential were introduced such as the Hamiltonian of the system and the Poisson brackets of the dynamic variables. This model is a straightforward generalization of the ordinary spin CM model and we present the Lax pair with spectral parameter expressed in the form of the root type matrix of semi-simple Lie algebra. Accordingly, we give a proof of the consistency between the Lax equation and the canonical equations of motion generated by Hamiltonian in detail. In order to illustrate the integrability of the CM system we construct the $R$-matrix and check that the $R$-matrix is equivalent to the Poisson brackets of Lax matrix. After these analysis, we come to the main discovery of this paper that we discuss the quantized model of the two type spins CM system and show that the $R$-matrix of this situation is the same as in the classical case. The final section of this paper is devoted to summary and comments.
2 Simple Lie Algebra and Root System
Let us briefly sketch some aspects of the concepts of semi-simple and simply-laced Lie algebra $\mathfrak{g}$ with rank $n$ which are very useful in the description of Calogero-Moser model. In fact we only need the data of its roots. Following the notations in Ref. [6], we denote the set of all roots by $\Delta$ and view the elements of $\Delta$ as real vectors in $n$-dimensional space which are normalised, without loss of generality, to 2, namely
$$ \Delta=\left\{\alpha,\beta,\gamma,\ldots \right\},\quad \alpha∈{R}^{n},\\ \alpha^{2}=\alpha\cdot\alpha=2, \quad \forall\alpha∈\Delta, $$ the dynamical variables are canonical coordinates $\left\{q_{i}\right\}$ associated with their canonical conjugate momenta $\left\{p_{i}\right\}$ satisfying the standard Poisson brackets:
$$ q_{1},\ldots,q_{n}\,, \quad p_{1},\ldots,p_{n}\,,\quad \left\{q_{i},p_{j}\right\}=\delta_{ij} \,,\\ \left\{q_{i},q_{j}\right\}=\left\{p_{i}\,,p_{j}\right\}=0\,, $$ it is convenient to denote them by $n$-dimensional vectors $q$ and $p$, that is
$$\displaylines{\quad q=(q_{1},\ldots,q_{n})\in \mathbb{R}^{n}\,,\quad p=(p_{1},\ldots,p_{n})\in \mathbb{R}^{n}\,,} $$ naturally the scalar products of $q$ and $p$ as well as the roots $\alpha\cdot q$, $p\cdot\beta$, etc. can be defined in the usual way.
In the subsequent part we intend to review the basic knowledge of simple Lie algebra and the associated root system.Let $\mathfrak{g}$ be a simple Lie algebra of rank $n$, $\mathfrak{b}$ a Cartan subalgebra,and $\Delta$ the associated root system, it is well known that the Cartan subalgebra, gives rise to a root space direct decomposition of the Lie algebra $\mathfrak{g}$,or in other words $\mathfrak{g}=\mathfrak{b}\oplus\bigoplus_{\alpha∈\Delta}\mathfrak{g}_{\alpha}$.We can choose a basis $\left\{E_{\alpha},H_{\mu}|\alpha∈\Delta,\mu=1,\ldots,n\right\}$ of $\mathfrak{g}$ with the following properties:
(i) $H_{\mu}~(\mu=1,\ldots,n)$ constitute an orthonormal basis of Cartan subalgebra $\mathfrak{b}$ equipped with the Killing form $B:\mathfrak{b}\times\mathfrak{b}\rightarrow\mathbb{C}$ induced from $\mathfrak{g}$ which satisfy $B(H_{\mu},H_{\nu})=\delta_{\mu\nu}$. Obviously by means of the Killing form there exists an isomorphism between $\mathfrak{b}^{*}=\mathrm{Hom}(\mathfrak{b},\mathbb{C})$ and $\mathfrak{b}$, hence for each $\alpha∈\mathfrak{b}^{*}$ one is able to determine an element $H_{\alpha}$. Now making use of the basis $H_{\mu}$ of $\mathfrak{b}$, this correspondence can be represented more explicitly in the form
$$ \alpha\mapsto H_{\alpha}=\sum_{\mu=1}^{n}\alpha(H_{\mu})H_{\mu}. $$ (ii) For each $\alpha∈\Delta$, the root subspace $\mathfrak{g}_{\alpha}$ is one-dimensional vector space together with the basis $E_{\alpha}$ such that $[E_{\alpha},E_{-\alpha}]=H_{\alpha}$ and for simplicity, we suppose that they fulfil the normalization condition $B(E_{\alpha},E_{-\alpha})=1$.
(iii) The Lie brackets of the basis elements take the form
$$ \left[H_{\mu},H_{\nu}\right]=0,\quad \left[E_{\alpha},E_{\beta}\right]=N_{\alpha,\beta} E_{\alpha+\beta} \quad (\alpha+\beta\neq 0), \\ \left[H_{\mu},E_{\alpha}\right]=\alpha(H_{\mu})E_{\alpha},\quad [E_{\alpha},E_{-\alpha}]=\sum_{\mu=1}^{n}\alpha(H_{\mu}) H_{\mu},$$ here the structure constants $N_{\alpha,\beta}$ are anti-symmetric with respect to the indices and vanish unless $\alpha+\beta$ is also a root included in $\Delta$. More importantly, it is necessary to stress the fact that these structure constants satisfy the general relations given by
$$ N_{\alpha+\beta,-\alpha}=N_{-\beta.\alpha+\beta} =N_{-\alpha,-\beta}=-N_{\alpha,\beta}. $$ Under these preparations, we are able to formulate the classical spin variables in the context of simple Lie algebra and root system, or more precisely, those spin variables are described as coordinates of the dual space $\mathfrak{g}^{*}=\mathrm{Hom}(\mathfrak{g},\mathbb{C})$. Let us introduce two sets of spin variables $F_{\alpha},F_{\mu}$ and $S_{\alpha},S_{\mu}$ respectively dual to the above basis $E_{\alpha}$ and $H_{\mu}$, which can be served as the coefficients of $E_{\alpha}$ and $H_{\mu}$ in the linear combination
$$ \sum_{\alpha∈\Delta}F_{-\alpha}E_{\alpha} +\sum_{\mu=1}^{n}F_{\mu}H_{\mu}, \sum_{\alpha∈\Delta}S_{-\alpha}E_{\alpha} +\sum_{\mu=1}^{n}S_{\mu}H_{\mu}, $$ furthermore we assume these spin variables satisfy the following generalized Poisson brackets
$$ \left\{F_{\alpha}\,,F_{-\alpha}\right\}=\sum_{\mu=1}^{n}\alpha(\mu)F_{\mu}\,,\quad \left\{S_{\alpha}\,,S_{-\alpha}\right\}=\sum_{\mu=1}^{n}\alpha(\mu)S_{\mu}\,,\\ \left\{F_{\alpha}\,,F_{\beta}\right\}=\frac{1}{4}(N_{\alpha,\beta}F_{\alpha+\beta}-N_{\alpha,-\beta}F_{\alpha-\beta}-N_{-\alpha,\beta}F_{\beta-\alpha}\\ \quad+N_{-\alpha,-\beta}F_{-\alpha-\beta})\,,\quad \alpha\pm\beta\neq 0\,,\\ \left\{S_{\alpha}\,,S_{\beta}\right\}=\frac{1}{4}(N_{\alpha,\beta}S_{\alpha+\beta}-N_{\alpha,-\beta}S_{\alpha-\beta}-N_{-\alpha,\beta}S_{\beta-\alpha}\\ \quad+N_{-\alpha,-\beta}S_{-\alpha-\beta})\,, \quad \alpha\pm\beta\neq 0\,,\\ \left\{F_{\mu}\,,F_{\alpha}\right\}=\frac{1}{2}\alpha(\mu)(F_{\alpha}-F_{-\alpha})\,,\quad \\\left\{S_{\mu}\,,S_{\alpha}\right\}=\frac{1}{2}\alpha(\mu)(S_{\alpha}-S_{-\alpha})\,,\\\left\{F_{\alpha}\,,S_{\beta}\right\}=0\,,\quad \left\{F_{\mu}\,,F_{\nu}\right\}=0\,,\quad \left\{S_{\mu}\,,S_{\nu}\right\}=0\,,\\ \left\{S_{\mu}\,,F_{\alpha}\right\}=0\,,\quad \left\{F_{\alpha}\,,S_{\mu}\right\}=0\,,$$ here $\alpha(\mu)$ is the short for $\alpha(H_{\mu})$. Comparing Eq.(8) with Eq.(5),
these expressions reveal their resemblances to the Lie brackets of the Lie algebra basis and using Eq.(6), it is easy for us to check that the Jacobi identities of the generalized Poisson bracket holds true.
3 Hamiltonian Dynamic System
As illustrated in classical and quantum mechanics, since quantization of a dynamic system requires full information about the canonical Poisson bracket relations between its variables, the Hamiltonian structure of nonlinear equations will be of utmost importance in future discussions. Thus let us restrict ourselves to the classical integrable system and our generalization of the Hamiltonian of the Calogero-Moser with two types of interacting spins is as follows[18-19]
$$ H=\frac{1}{2}p\cdot p-\frac{1}{4}\sum_{\alpha∈ \Delta} \Bigl(\frac{F_{-\alpha}F_{\alpha}-F_{\alpha}F_{\alpha}+S_{-\alpha}S_{\alpha} -S_{\alpha}S_{\alpha}}{\text{sin}^{2}(\alpha\cdot q)} +\frac{2\text{cos}(\alpha\cdot q)}{\text{sin}^{2}(\alpha\cdot q)}(F_{-\alpha}S_{\alpha}-F_{\alpha}S_{\alpha})\Bigr),$$ here $p$ and $q$ represent the positions and conjugate momenta of $n$ identical particles while $F_{\alpha},S_{\alpha}$ are the spin variables as explained above. For the purpose of obtaining the explicit equations of motion of this Hamiltonian which are governed by the generalized Poisson brackets Eq.(8), for $\beta\neq\pm\alpha$, a simple and direct calculation shows
$$ \left\{F_{\alpha},F_{-\beta}F_{\beta}\right\}=\frac{1}{4}F_{-\beta} (N_{\alpha,\beta}F_{\beta+\alpha}-N_{\alpha,-\beta}F_{\alpha-\beta} -N_{-\alpha,\beta}F_{\beta-\alpha} +N_{-\alpha,-\beta}F_{-\alpha-\beta}) +\frac{1}{4}(N_{\alpha,-\beta}F_{\alpha-\beta}-N_{\alpha,\beta}F_{\alpha+\beta} -N_{-\alpha,-\beta}F_{-\alpha-\beta} +N_{-\alpha,\beta}F_{\beta-\alpha})F_{\beta},$$ along with
$$ \sum_{β∈\Delta}\left\{F_{\mu},F_{-\beta}F_{\beta} -F_{\beta}F_{\beta}\right\} =\sum_{β∈\Delta}(-\beta(\mu))(F_{-\beta}-F_{\beta})F_{\beta}+\sum_{β∈\Delta}\beta(\mu)F_{-\beta} (F_{\beta}-F_{-\beta}) -\sum_{β∈\Delta}\beta(\mu)(F_{\beta}-F_{-\beta})F_{\beta} -\sum_{β∈\Delta}\beta(\mu)F_{\beta}(F_{\beta}-F_{-\beta}) =\sum_{β∈\Delta}\beta(\mu)(F_{\beta}-F_{-\beta})F_{-\beta}-\sum_{β∈\Delta}\beta(\mu)F_{\beta} (F_{-\beta}-F_{\beta}) +\sum_{β∈\Delta}\beta(\mu)(F_{-\beta}-F_{\beta})F_{-\beta} -\sum_{β∈\Delta}\beta(\mu) F_{\beta}(F_{\beta}-F_{-\beta})=0, \left\{F_{\alpha},F_{-\alpha}F_{\alpha}+F_{\alpha}F_{-\alpha}-F_{-\alpha}F_{-\alpha}\right\}= \sum_{\mu=1}^{n}\alpha(\mu)(F_{\mu}F_{\alpha}+F_{\alpha}F_{\mu}-F_{\mu}F_{-\alpha}-F_{-\alpha}F_{\mu}), $$ for convenience, we will define functions $u(\alpha\cdot q), v(\alpha\cdot q)$ as
$$ u(\alpha\cdot q)=\frac{1}{\text{sin}(\alpha\cdot q)},\quad \quad v(\alpha\cdot q)= \frac{\text{cos}(\alpha\cdot q)}{\text{sin}(\alpha\cdot q)}, $$ moreover, it should be remarked here that they obey the following crucial identities
$$ u'(\alpha\cdot q)=u'(-\alpha\cdot q),\quad \quad v'(\alpha\cdot q)=v'(-\alpha\cdot q), \\ u(\alpha\cdot q)u'(\beta\cdot q)-u(\beta\cdot q)u'(\alpha\cdot q) =u((\alpha+\beta)\cdot q)\Bigl(\frac{1}{\text{sin}^{2}(\alpha\cdot q)}-\frac{1}{\text{sin}^{2}(\beta\cdot q)}\Bigr), \\ v(\alpha\cdot q)v'(\beta\cdot q)-v(\beta\cdot q)v'(\alpha\cdot q) =v((\alpha+\beta)\cdot q)=\Bigl(\frac{1}{\text{sin}^{2}(\alpha\cdot q)} -\frac{1}{\text{sin}^{2}(\beta\cdot q)}\Bigr). $$ Let us now return to the original Hamiltonian Eq.(9) and adopt the convention of coordinates $q_{\mu}=q\cdot H_{\mu}$ as well as momenta $p_{\mu}=p\cdot H_{\mu}$, especially notice that $\left\{F_{\alpha}, F_{-\beta}\right\}=-\left\{F_{\alpha},F_{\beta}\right\}$ and $\left\{F_{\alpha}, F_{-\beta}F_{\beta}\right\}=-\left\{F_{\alpha}, F_{\beta}F_{\beta}\right\}$. Based on these preliminaries, we can derive the equations of motion of the generalized dynamic system straightforwardly
$$ \dot{q}_{\mu}=p_{\mu},\dot{p}_{\mu}=\frac{1}{2}\sum_{\alpha∈ \Delta}\alpha(\mu)\Bigl(\frac{\text{cos}(\alpha\cdot q)}{\text{sin}^{3}(\alpha\cdot q)}(F_{-\alpha}F_{\alpha}-F_{\alpha}F_{\alpha} +S_{-\alpha}S_{\alpha} -S_{\alpha}S_{\alpha})+\frac{\text{sin}^{2}(\alpha\cdot q)-2}{\text{sin}^{3}(\alpha\cdot q)} (F_{-\alpha}S_{\alpha}-F_{\alpha}S_{\alpha})\Bigr),\dot{F}_{\alpha}=\left\{F_{\alpha},H\right\}=-\frac{1}{4}\sum_{\beta\neq\pm\alpha}\frac{F_{-\beta}}{\text{sin}^{2} (\beta\cdot q)}(N_{\alpha,\beta}F_{\alpha+\beta} -N_{\alpha,-\beta} F_{\alpha-\beta}-N_{-\alpha,\beta}F_{\beta-\alpha}+N_{-\alpha,-\beta}F_{-\alpha-\beta})-\frac{1}{4}\sum_{\beta\neq\pm\alpha}\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)} \times(N_{\alpha,-\beta}F_{\alpha-\beta}-N_{\alpha,\beta}F_{\alpha+\beta}-N_{-\alpha,-\beta} F_{-\alpha-\beta}+N_{-\alpha,\beta}F_{\beta-\alpha})S_{\beta}-\frac{1}{2\text{sin}^{2}(\alpha\cdot q)}\sum_{\mu=1}^{n}\alpha(\mu)(F_{\mu}F_{\alpha}-F_{-\alpha}F_{\mu})-\frac{\text{cos}(\alpha\cdot q)}{2\text{sin}^{2}(\alpha\cdot q)} \times\sum_{\mu=1}^{n}\alpha(\mu)(F_{\mu}S_{\alpha}-F_{\mu}S_{-\alpha}), $$ and the similar expression for $\dot{S}_{\alpha}$. Also we have $\dot{F}_{\mu}=\dot{G}_{\mu}=0$, which implies the diagonal elements of the spin variables are conserved quantities.
4 Lax Pair and Lax Equation
For further analysis, it is useful to search for an appropriate Lax pair $(L(\lambda),M)$ which reflects the classical system is completely integrable in the Liouville sense. Generally speaking, the Lax pair may depend on an auxiliary parameter $\lambda$ usually referred to as the spectral parameter of the Lax equation. Indeed, due to the motion integrations $\dot{F}_{\mu}=\dot{G}_{\mu}=0$, or on the reduced phase space $F_{\mu}=f_{\mu},G_{\mu}=g_{\mu}$, here $f_{\mu},g_{\mu},\mu=1,2,...,n$ are some constants and if these constants satisfy the constraints
$$ \sum_{\mu=1}^{n}\alpha(\mu)f_{\mu}=0,\quad \quad \sum_{\mu=1}^{n}\alpha(\mu)g_{\mu}=0, $$ the previous Hamiltonian admits the following elegant Lax matrix representation
$$ L(\lambda)=\sum_{i=1}^{n}p_{\mu}H_{\mu}+\frac{1}{2}\sum_{\alpha∈\Delta}u(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})E_{\alpha} +\frac{1}{2}\sum_{\alpha∈\Delta}(v(\alpha\cdot q)+v(\lambda))(S_{-\alpha}-S_{\alpha})E_{\alpha},\\ M=-\frac{1}{2}\sum_{\alpha∈\Delta}u'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})E_{\alpha}-\frac{1}{2}\sum_{\alpha∈\Delta}v'(\alpha\cdot q)(S_{-\alpha}-S_{\alpha})E_{\alpha}, $$ here $\lambda$ is the spectral parameter and the functions $u(\alpha\cdot q),v(\alpha\cdot q)$ are defined as above.
${{Proof}}$ To easy of exposition, we define $\tilde{v}(\alpha\cdot q)=v(\alpha\cdot q)+v(\lambda)$.Let us take the time derivative of the matrix $L(\lambda)$ and the result can be divided into three parts $\dot{L}=A+B+C+D$, here symbols $A,B,C,D$ correspond to the terms of $t$-derivative of momentum $p_{\mu}$, $t$-derivative of functions $u(\alpha\cdot q),v(\alpha\cdot q)$ and $t$-derivative of spin variables $F_{\alpha},G_{\alpha}$ in $L(\lambda)$ respectively. Specifically, with the help of $\dot{F}_{\alpha}=-\dot{F}_{-\alpha}$ we find that
$$ A=\sum_{i=1}^{n}\dot{p}_{\mu}H_{\mu}=\frac{1}{2}\sum_{\alpha∈ \Delta}\Bigl[\frac{\text{cos}(\alpha\cdot q)}{\text{sin}^{3}(\alpha\cdot q)}(F_{-\alpha}F_{\alpha}-F_{\alpha}F_{\alpha} +S_{-\alpha}S_{\alpha}-S_{\alpha}S_{\alpha})H_{\alpha}-\frac{\text{sin}^{2}(\alpha\cdot q)-2}{2\text{sin}^{3}(\alpha\cdot q)} (F_{-\alpha}S_{\alpha}-F_{\alpha}S_{\alpha})H_{\alpha}\Bigr],B=\frac{1}{2}\sum_{\alpha∈\Delta}\dot{u}(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})E_{\alpha}+\dot{v}(\alpha\cdot q)(S_{-\alpha}-S_{\alpha})E_{\alpha} =\frac{1}{2}\sum_{\alpha∈\Delta}\alpha\cdot pu'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})E_{\alpha}+\alpha\cdot pv'(\alpha\cdot q)(S_{-\alpha}-S_{\alpha})E_{\alpha},C=\frac{1}{2}\sum_{\alpha∈\Delta}u(\alpha\cdot q)(\dot{F}_{-\alpha}E_{\alpha}-\dot{F}_{\alpha}E_{\alpha})=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}u(\alpha\cdot q)\frac{1}{\text{sin}^{2}(\beta\cdot q)}F_{-\beta}(N_{-\alpha,\beta}F_{\beta-\alpha}-N_{-\alpha,-\beta}F_{-\alpha-\beta}-N_{\alpha,\beta}F_{\alpha+\beta}+N_{\alpha,-\beta}F_{\alpha-\beta})E_{\alpha} -\frac{1}{4}\sum_{\alpha\neq\pm\beta} u(\alpha\cdot q)\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)}(N_{\alpha,-\beta}F_{\alpha-\beta}-N_{\alpha,\beta}F_{\alpha+\beta}-N_{-\alpha,-\beta}F_{-\alpha-\beta}+N_{-\alpha,\beta}F_{\beta-\alpha})S_{\beta}E_{\alpha},D=\frac{1}{2}\sum_{\alpha∈\Delta}\tilde{v}(\alpha\cdot q)(\dot{S}_{-\alpha}E_{\alpha}-\dot{S}_{\alpha}E_{\alpha}) =-\frac{1}{4}\sum_{\alpha\neq\pm\beta} \tilde{v}(\alpha\cdot q)\frac{1}{\text{sin}^{2}(\beta\cdot q)}\Bigl[S_{-\beta}(N_{-\alpha,\beta}S_{\beta-\alpha}-N_{-\alpha,-\beta}S_{-\alpha-\beta}-N_{\alpha,\beta}S_{\alpha+\beta} +N_{\alpha,-\beta}S_{\alpha-\beta}) E_{\alpha}-\frac{1}{4}\sum_{\alpha\neq\pm\beta}\tilde{v}(\alpha\cdot q)\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)}(N_{\alpha,-\beta}S_{\alpha-\beta}-N_{\alpha,\beta}S_{\alpha+\beta}-N_{-\alpha,-\beta}S_{-\alpha-\beta}+N_{-\alpha,\beta}S_{\beta-\alpha})F_{\beta}E_{\alpha}\Bigr]. $$ on the other hand, the commutation relation of the Lax pair $(L(\lambda), M)$ is expressible in the form of $[L(\lambda), M]=R+T+X+Y$, here $R$ is the commutator of $H_{\mu}$ and $M$, that is
$$ R=\Bigl[\sum_{i=1}^{n}p_{\mu}H_{\mu},M\Bigr]=-\frac{1}{2}\sum_{i=1}^{n}\sum_{\alpha\in\Delta}[u'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+v'(\alpha\cdot q)(S_{-\alpha}-S_{\alpha})]p_{\mu}\left[H_{\mu},E_{\alpha}\right]\\=-\frac{1}{2}\sum_{i=1}^{n}\sum_{\alpha\in\Delta}[u'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+v'(\alpha\cdot q)(S_{-\alpha}-S_{\alpha})]p_{\mu}\alpha(\mu)E_{\alpha}\\=-\frac{1}{2}\sum_{i=1}^{n}\sum_{\alpha\in\Delta}\alpha\cdot p[u'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+v'(\alpha\cdot q)(S_{-\alpha}-S_{\alpha})]E_{\alpha}\,,$$ while $T=T_{1}+T_{2}+T_{3}+T_{4}$ stands for the terms of Cartan part ($\alpha+\beta=0$) in the commutator $[L(\lambda), M]$. At first, due to the symmetry with respect to $\alpha$ and $-\alpha$ we conclude
$$ T_{1}=-\frac{1}{4}\sum_{\alpha\in\Delta}u(\alpha\cdot q)u'(-\alpha\cdot q)(F_{-\alpha}-F_{\alpha})(F_{\alpha}-F_{-\alpha})\left[E_{\alpha},E_{-\alpha}\right] =\frac{1}{4}\sum_{\alpha\in \Delta}\frac{\text{cos}(\alpha\cdot q)}{\text{sin}^{3}(\alpha\cdot q)} (F_{-\alpha}F_{\alpha} -F_{-\alpha}F_{-\alpha}-F_{\alpha}F_{\alpha}+F_{\alpha}F_{-\alpha})H_{\alpha} =\frac{1}{2}\sum_{\alpha\in \Delta}\frac{\text{cos}(\alpha\cdot q)}{\text{sin}^{3}(\alpha\cdot q)} (F_{-\alpha}F_{\alpha}-F_{\alpha}F_{\alpha})H_{\alpha}\,,$$ here we have already taken the advantage of $H_{\alpha}+H_{-\alpha}=0$. In addition, since $v'(\alpha\cdot q)$ is an even function, applying the symmetry between $\alpha$ and $-\alpha$ in a similar way we assert
$$ \sum_{\alpha∈\Delta}v'(-\alpha\cdot q)(S_{-\alpha}-S_{\alpha})(S_{\alpha}-S_{-\alpha})H_{\alpha} =\frac{1}{2}\sum_{\alpha∈\Delta}[v'(-\alpha\cdot q)(S_{-\alpha}-S_{\alpha})(S_{\alpha}-S_{-\alpha})H_{\alpha} +v'(\alpha\cdot q)(S_{\alpha}-S_{-\alpha})(S_{-\alpha}-S_{\alpha})H_{-\alpha}] = \frac{1}{2}\sum_{\alpha∈\Delta}v'(\alpha\cdot q)(S_{-\alpha}-S_{\alpha})(S_{\alpha}-S_{-\alpha})(H_{\alpha}+H_{-\alpha})=0\,, $$ which leads to
$$ T_{2}=-\frac{1}{4}\sum_{\alpha∈\Delta}(v(\alpha\cdot q)+v(\lambda))v'(-\alpha\cdot q)(S_{-\alpha}-S_{\alpha})(S_{\alpha}-S_{-\alpha})\left[E_{\alpha},E_{-\alpha}\right]=\frac{1}{2}\sum_{\alpha∈ \Delta}\frac{\text{cos}(\alpha\cdot q)}{\text{sin}^{3}(\alpha\cdot q)}(S_{-\alpha}S_{\alpha}-S_{\alpha}S_{\alpha})H_{\alpha}, $$ analogously we have
$$ T_{3}=-\frac{1}{4}\sum_{\alpha\in\Delta}u(\alpha\cdot q)v'(-\alpha\cdot q)(F_{-\alpha}-F_{\alpha})(S_{\alpha}-S_{-\alpha})\left[E_{\alpha},E_{-\alpha}\right] =-\frac{1}{4}\sum_{\alpha\in\Delta}u(\alpha\cdot q)v'(-\alpha\cdot q)(F_{-\alpha}S_{\alpha}+F_{\alpha} S_{-\alpha} -F_{\alpha}S_{\alpha}-F_{-\alpha}S_{-\alpha})H_{\alpha} =-\frac{1}{4}\sum_{\alpha\in\Delta}(u(\alpha\cdot q)v'(-\alpha\cdot q)-u(-\alpha\cdot q)v'(\alpha\cdot q))(F_{-\alpha}S_{\alpha}-F_{\alpha}S_{\alpha})H_{\alpha} =\frac{1}{2}\sum_{\alpha\in\Delta}\frac{1}{\text{sin}^{3}(\alpha\cdot q)}(F_{-\alpha}S_{\alpha}-F_{\alpha}S_{\alpha})H_{\alpha}\,,$$ as well as
$$ T_{3}+T_{4}=-\frac{1}{2}\sum_{\alpha∈\Delta}\frac{\text{sin}^{2}(\alpha\cdot q)-2}{\text{sin}^{3}(\alpha\cdot q)}(F_{-\alpha}S_{\alpha}-F_{\alpha}S_{\alpha})H_{\alpha}.$$ Now let us turn the attention to $X$ as given in the commutator which are the terms of off-Cartan part ($\alpha+\beta\neq 0$) of $[L(\lambda), M]$, firstly for pure spin variables $F_{\alpha}$ we simply have
$$ X_{1}=-\frac{1}{4}\sum_{\alpha,\beta\in\Delta,\alpha+\beta\neq 0}u(\alpha\cdot q)u'(\beta\cdot q)(F_{-\alpha}-F_{\alpha})(F_{-\beta}-F_{\beta})\left[E_{\alpha},E_{\beta}\right]\\=-\frac{1}{4}\sum_{\alpha,\beta\in\Delta,\alpha\neq\pm\beta}u(\alpha\cdot q)u'(\beta\cdot q)N_{\alpha,\beta}(F_{-\alpha}F_{-\beta}+F_{\alpha}F_{\beta}-F_{-\alpha}F_{\beta}-F_{\alpha}F_{-\beta})E_{\alpha+\beta}\,, $$ making use of the symmetry between $\alpha$ and $\beta$ one obtains
$$ \sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}u(\alpha\cdot q)u'(\beta\cdot q)N _{\alpha,\beta}F_{-\alpha}F_{-\beta}E_{\alpha+\beta} =\frac{1}{2}\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}[u(\alpha\cdot q)u'(\beta\cdot q)N_{\alpha,\beta}F_{-\alpha}F_{-\beta} +u(\beta\cdot q)u'(\alpha\cdot q)N _{\beta,\alpha}F_{-\beta}F_{-\alpha}],E_{\alpha+\beta}=\frac{1}{2}\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}[u(\alpha\cdot q)u'(\beta\cdot q)-u(\beta\cdot q)u'(\alpha\cdot q)] N _{\alpha,\beta}F_{-\alpha}F_{-\beta}E_{\alpha+\beta}=\frac{1}{2}\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}u((\alpha+\beta)\cdot q)\Bigl(\frac{1}{\text{sin}^{2}(\alpha\cdot q)}-\frac{1}{\text{sin}^{2}(\beta\cdot q)}\Bigr) N _{\alpha,\beta}F_{-\alpha}F_{-\beta}E_{\alpha+\beta}=\frac{1}{2}\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}u(\alpha\cdot q)\frac{1}{\text{sin}^{2}(\beta\cdot q)}(N _{\beta,\alpha-\beta}F_{-\beta}F_{\beta-\alpha}E_{\alpha}-N _{\alpha-\beta,\beta}F_{\beta-\alpha}F_{-\beta}E_{\alpha})=\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}u(\alpha\cdot q)\frac{1}{\text{sin}^{2}(\beta\cdot q)}N _{-\alpha,\beta}F_{-\beta}F_{\beta-\alpha}E_{\alpha}, $$ a similar deduction for the other terms in Eq.(25) we arrive at
$$ X_{1}=-\frac{1}{4}\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}u(\alpha\cdot q)\frac{1}{\text{sin}^{2}(\beta\cdot q)}(N_{-\alpha,\beta}F_{-\beta}F_{\beta-\alpha}E_{\alpha}-N_{\alpha,\beta}F_{-\beta}F_{\alpha+\beta}E_{\alpha}-N_{-\alpha,-\beta}F_{-\beta}F_{-\alpha-\beta}E_{\alpha}+N_{\alpha,-\beta}F_{-\beta}F_{\alpha-\beta}E_{\alpha}), $$ now the analogous expression for the single spin variable $S_{\alpha}$ turns out to be
$$ X_{2}=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}\tilde{v}(\alpha\cdot q)v'(\beta\cdot q)(S_{-\alpha}-S_{\alpha})(S_{-\beta}-S_{\beta})\left[E_{\alpha},E_{\beta}\right]=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}\Bigl[v(\alpha\cdot q)\frac{1}{\text{sin}^{2}(\beta\cdot q)}(N_{-\alpha,\beta}S_{-\beta}S_{\beta-\alpha}E_{\alpha}-N_{\alpha,\beta}S_{-\beta}S_{\alpha+\beta}E_{\alpha}-N_{-\alpha,-\beta}S_{-\beta}S_{-\alpha-\beta}E_{\alpha}+N_{\alpha,-\beta}S_{-\beta}S_{\alpha-\beta}E_{\alpha})-\frac{1}{4}v(\lambda)\frac{1}{\text{sin}^{2}(\beta\cdot q)}(N_{-\alpha,\beta}S_{-\beta}S_{\beta-\alpha}E_{\alpha}-N_{\alpha,\beta}S_{-\beta}S_{\alpha+\beta}E_{\alpha}-N_{-\alpha,-\beta}S_{-\beta}S_{-\alpha-\beta}E_{\alpha}+N_{\alpha,-\beta}S_{-\beta}S_{\alpha-\beta}E_{\alpha})\Bigr]=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}\tilde{v}(\alpha\cdot q)\frac{1}{\text{sin}^{2}(\beta\cdot q)}(N_{-\alpha,\beta}S_{-\beta}S_{\beta-\alpha}E_{\alpha}-N_{\alpha,\beta}S_{-\beta}S_{\alpha+\beta}E_{\alpha}-N_{-\alpha,-\beta}S_{-\beta}S_{-\alpha-\beta}E_{\alpha}+N_{\alpha,-\beta}S_{-\beta}S_{\alpha-\beta}E_{\alpha}), $$ besides, the terms of mixed spin variables $F_{\alpha}$ and $G_{\beta}$ are given by
$$ X_{3}=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}u(\alpha\cdot q)\tilde{v}'(\beta\cdot q)(F_{-\alpha}-F_{\alpha})(S_{-\beta}-S_{\beta})\left[E_{\alpha},E_{\beta}\right]=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}u(\alpha\cdot q)v'(\beta\cdot q)N _{\alpha,\beta}(F_{-\alpha}S_{-\beta}+F_{\alpha}S_{\beta}-F_{-\alpha}S_{\beta}-F_{\alpha}S_{-\beta})E_{\alpha+\beta}X_{4}=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}\tilde{v}(\alpha\cdot q)u'(\beta\cdot q)(F_{-\alpha}-F_{\alpha})(S_{-\beta}-S_{\beta})\left[E_{\alpha},E_{\beta}\right]=-\frac{1}{4}\sum_{\alpha\neq\pm\beta}\tilde{v}(\alpha\cdot q)u'(\beta\cdot q)N _{\alpha,\beta}(S_{-\alpha}F_{-\beta}+S_{\alpha}F_{\beta}-S_{-\alpha}F_{\beta}-S_{\alpha}F_{-\beta})E_{\alpha+\beta},(31)$$ thus we should exploit the functions $u(\alpha\cdot q)v'(\beta\cdot q),\tilde{v}(\alpha\cdot q)u'(\beta\cdot q)$ mentioned above and pay attention that $v(\beta\cdot q)=u(\beta\cdot q)\text{cos}(\beta\cdot q)$ which immediately leads to
$$ u((\alpha+\beta)\cdot q)v'(-\beta\cdot q)-\tilde{v}(-\beta\cdot q)u'((\alpha+\beta)\cdot q)=(u((\alpha+\beta)\cdot q)u'(-\beta\cdot q)-u(-\beta\cdot q)u'((\alpha+\beta)\cdot q))\text{cos}(\beta\cdot q)+u((\alpha+\beta)\cdot q)u(-\beta\cdot q)\text{sin}(\beta\cdot q)-v(\lambda)u'((\alpha+\beta)\cdot q)=u(\alpha\cdot q)\Bigl(\frac{1}{\text{sin}^{2}((\alpha+\beta)\cdot q)}-\frac{1}{\text{sin}^{2}(\beta\cdot q)}\Bigr)\text{cos}(\beta\cdot q)-\frac{1}{\text{sin}((\alpha+\beta)\cdot q)}-v(\lambda)u'((\alpha+\beta)\cdot q)=-u(\alpha\cdot q)\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)}+\frac{\text{cos}(\beta\cdot q)-\text{sin}(\alpha\cdot q)\text{sin}((\alpha+\beta)\cdot q)}{\text{sin}(\alpha\cdot q)\text{sin}^{2}((\alpha+\beta)\cdot q)}-v(\lambda)u'((\alpha+\beta)\cdot q)=-u(\alpha\cdot q)\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)}+\frac{\text{cos}(\alpha\cdot q)\text{cos}((\alpha+\beta)\cdot q)}{\text{sin}(\alpha\cdot q)\text{sin}^{2}((\alpha+\beta)\cdot q)}+v(\lambda)\frac{\text{cos}((\alpha+\beta)\cdot q)}{\text{sin}^{2}((\alpha+\beta)\cdot q)}=-u(\alpha\cdot q)\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)} +\tilde{v}(\alpha\cdot q)\frac{\text{cos}((\alpha+\beta)\cdot q)}{\text{sin}^{2}((\alpha+\beta)\cdot q)}, $$ therefore from Eq.(32) one is capable of evaluating the addition of $X_{3}, X_{4}$ as
$$ X_{3}+X_{4} =\frac{1}{4}\sum_{\alpha\neq\pm\beta} u(\alpha\cdot q)\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)}(N_{\alpha,-\beta}F_{\alpha-\beta} -N_{\alpha,\beta} F_{\alpha+\beta}-N_{-\alpha,-\beta}F_{-\alpha-\beta} +N_{-\alpha,\beta}F_{\beta-\alpha})S_{\beta}E_{\alpha} +\sum_{\alpha\neq\pm\beta}\tilde{v}(\alpha\cdot q)\frac{\text{cos}(\beta\cdot q)}{\text{sin}^{2}(\beta\cdot q)} (N_{\alpha,-\beta}S_{\alpha-\beta}-N_{\alpha,\beta} S_{\alpha+\beta}-N_{-\alpha,-\beta}S_{-\alpha-\beta} +N_{-\alpha,\beta}S_{\beta-\alpha})F_{\beta}E_{\alpha}\,,$$ through the above calculations, we verify that the consistency of Lax pair is essentially the same as for the canonical equations of the Hamiltonian system. That is complete the proof.
5 R-matrix
We have seen that the Lax pair method is an excellent and conceptual tool for the investigation of soliton theory, or more generally the issue of integrability in the classical mechanic system. Another essential and more particularly heuristic ingredient of the integrable system is inspired from the point of view of the $R$-matrix[20] from which, like the previously known Lax pair, it means that the system possesses sufficiently many number of conserved quantities in soliton theory. The significance of this method lies in it works not only for the classical case but also is a central object in the modern approach to the studies of quantum integrable systems deriving in particular a basic equation satisfied by it, namely the Yang-Baxter equation and perhaps even their supersymmetric generalizations of the associated system. Now considering the Calogero-Moser with two type spins model discussed above, we should emphasize that it is only the reduced system or in the equivalent form of $\dot{F}_{\mu}=\dot{G}_{\mu}=0$ that the original system is integrable. Actually, according to the condition Eq. (15), we establish the following $R$-matrix, which can be formulated naturally via the linear combinations of tensor product of root matrix
$$ a(\lambda,\eta)=-\frac{1}{2}(\text{cot}(\lambda+\eta)+\text{cot}(\lambda-\eta)), \\ b(\alpha\cdot q,\lambda,\eta)=-\frac{1}{2}(\text{cot}(\lambda-\eta)+\text{cot}(\alpha\cdot q)), \\ c(\alpha\cdot q,\lambda,\eta)=-\frac{1}{2}(\text{cot}(\lambda+\eta)+\text{cot}(\alpha\cdot q)), $$ obviously we have
$$ b(\alpha\cdot q,\lambda,\eta)=-b(-\alpha\cdot q,\eta,\lambda), \\ c(\alpha\cdot q,\eta,\lambda)=c(\alpha\cdot q,\lambda,\eta). $$ With the help of this $R$-matrix, it is possible to generate a hierarchy of integrable equations which clearly indicates a prominent role played in the case of classical integrable systems. For the sake of simplicity, we reformulate the matrix $L(\lambda)=\sum_{\mu}L_{\mu}H_{\mu}+\sum_{\alpha∈\Delta}L_{\alpha}(\lambda)E_{\alpha}$ as obviously we have
$$L_{\mu}=p_{\mu},L_{\alpha}(\lambda)=\frac{1}{2}u(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+\frac{1}{2}\tilde{v}(\alpha\cdot q,\lambda)(S_{-\alpha}-S_{\alpha}),(39) $$ and it is easily to deduce the commutation relations ($\alpha\neq\pm\beta$) between dynamic functions $L_{\mu},L_{\alpha}$ from the Poisson bracket (8)
obviously we have
$$ \left\{L_{\mu},L_{\nu}\right\}=0,\quad \quad \left\{L_{\mu},L_{\alpha}(\lambda)\right\}=\frac{1}{2}\alpha(\mu)[u'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha}) +\tilde{v}'(\alpha\cdot q,\lambda)(S_{-\alpha}-S_{\alpha})], \left\{L_{\alpha}(\lambda),L_{-\alpha}(\eta)\right\}=-\frac{1}{2}\sum_{\mu=1}^{n}\alpha(\mu) [u(\alpha\cdot q)u(-\alpha\cdot q)F_{\mu} +\tilde{v}(\alpha\cdot q,\lambda)\tilde{v}(-\alpha\cdot q,\eta)S_{\mu}],\left\{L_{\alpha}(\lambda),L_{\beta}(\eta)\right\}=\frac{1}{4}u(\alpha\cdot q)u(\beta\cdot q)(N_{\alpha,\beta} F_{\alpha+\beta}-N_{-\alpha,\beta}F_{\beta-\alpha}-N_{\alpha,-\beta}F_{\alpha-\beta}+N_{-\alpha,-\beta}F_{-\alpha-\beta}) +\frac{1}{4}\tilde{v}(\alpha\cdot q,\lambda)\tilde{v}(\beta\cdot q,\eta)(N_{\alpha,\beta}S_{\alpha+\beta}-N_{-\alpha,\beta}S_{\beta-\alpha}-N_{\alpha,-\beta}S_{\alpha-\beta}+N_{-\alpha,-\beta}S_{-\alpha-\beta}),$$ noting $L_{1}=L\otimes1$, $L_{2}=1\otimes L$ together with the antisymmetric property of $N_{\alpha,\beta}$ we have obviously we have
$$\left\{L_{1}(\lambda),L_{2}(\eta)\right\}=\sum_{\mu,\alpha∈\Delta}\left\{L_{\mu},L_{\alpha}(\eta)\right\}H_{\mu}\otimes E_{\alpha}+\sum_{\alpha,β∈\Delta}\left\{L_{\alpha}(\lambda),L_{\beta}(\eta)\right\}E_{\alpha}\otimes E_{\beta}+\sum_{\mu,\alpha∈\Delta}\left\{L_{\alpha}(\lambda),L_{\mu}\right\}E_{\alpha}\otimes H_{\mu}=\frac{1}{2}\sum_{\alpha∈\Delta}[u'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+\tilde{v}'(\alpha\cdot q,\eta)(S_{-\alpha}-S_{\alpha})]H_{\alpha}\otimes E_{\alpha}-\frac{1}{2}\sum_{\alpha∈\Delta}[u'(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+\tilde{v}'(\alpha\cdot q,\lambda)\times(S_{-\alpha}-S_{\alpha})]E_{\alpha}\otimes H_{\alpha}-\frac{1}{4}\sum_{\alpha\neq\pm\beta}u(\alpha\cdot q)u(\beta\cdot q)[N_{\alpha,\beta}(F_{-\alpha-\beta}-F_{\alpha+\beta})-N_{-\alpha,\beta}(F_{\alpha-\beta}-F_{\beta-\alpha})]E_{\alpha}\otimes E_{\beta}-\frac{1}{4}\sum_{\alpha\neq\pm\beta}\tilde{v}(\alpha\cdot q,\lambda)\tilde{v}(\beta\cdot q,\eta)[N_{\alpha,\beta}(S_{-\alpha-\beta}-S_{\alpha+\beta})-N_{-\alpha,\beta}(S_{\alpha-\beta}-S_{\beta-\alpha})]E_{\alpha}\otimes E_{\beta}-\frac{1}{2}\sum_{\mu,\alpha∈\Delta}\alpha(\mu)[u(\alpha\cdot q)u(-\alpha\cdot q)F_{\mu}+\tilde{v}(\alpha\cdot q,\lambda)\tilde{v}(-\alpha\cdot q,\eta)S_{\mu}]E_{\alpha}\otimes E_{-\alpha}, $$ noting that on the reduced phase space (15), the last term in Eq.(41) vanishes.As an aside, we find obviously we have
$$ \left[r_{12}(\lambda,\eta),L_{1}(\lambda)\right]=\sum_{\mu,\alpha∈\Delta}\left[a(\lambda,\mu)H_{\mu}\otimes H_{\mu},L_{\alpha}(\lambda)E_{\alpha}\otimes I\right]+\sum_{\mu,\alpha∈\Delta}[b(\alpha\cdot q,\lambda,\eta)E_{\alpha}\otimes E_{-\alpha}+c(\alpha\cdot q,\lambda,\eta)E_{\alpha}\otimes E_{\alpha},L_{\mu}H_{\mu}\otimes I+L_{\alpha}(\lambda)E_{\alpha}\otimes I],\left[r_{21}(\eta,\lambda),L_{2}(\eta)\right]\hspace {25.5mm}=\sum_{\mu,\alpha∈\Delta}[a(\eta,\lambda)H_{\mu}\otimes H_{\mu},L_{\alpha}(\eta)I\otimes E_{\alpha}]+\sum_{\mu,\alpha∈\Delta}[b(\alpha\cdot q,\eta,\lambda)E_{-\alpha}\otimes E_{\alpha}+c(\alpha\cdot q,\eta,\lambda)E_{\alpha}\otimes E_{\alpha},L_{\mu}I\otimes H_{\mu}+L_{\alpha}(\eta)I\otimes E_{\alpha}],$$ one can then write obviously we have
$$ \sum_{\mu,\alpha∈\Delta}(\left[a(\lambda,\eta)H_{\mu}\otimes H_{\mu},L_{\alpha}(\lambda)E_{\alpha} \otimes I\right]-\left[a(\eta,\lambda)H_{\mu}\otimes H_{\mu},L_{\alpha}(\eta)I\otimes E_{\alpha}\right]) =\sum_{\mu,\alpha∈\Delta}a(\lambda,\eta)L_{\alpha}(\lambda)\left[H_{\mu},E_{\alpha}\right] \otimes H_{\mu} +\sum_{\mu,\alpha∈\Delta}a(\eta,\lambda)L_{\alpha}(\eta)H_{\mu}\otimes [H_{\mu},E_{\alpha}] =\sum_{\mu,\alpha∈\Delta}a(\lambda,\eta)L_{\alpha}(\lambda)\alpha(\mu)E_{\alpha}\otimes H_{\mu}+\sum_{\mu,\alpha∈\Delta}a(\eta,\lambda)L_{\alpha}(\eta)\alpha(\mu)H_{\mu}\otimes E_{\alpha}=\sum_{\alpha∈\Delta}a(\lambda,\eta)L_{\alpha}(\lambda)E_{\alpha}\otimes H_{\alpha} +\sum_{\alpha∈\Delta}a(\eta,\lambda)L_{\alpha}(\eta)H_{\alpha}\otimes E_{\alpha}=\frac{1}{2}\sum_{\alpha∈\Delta}a(\lambda,\eta)[u(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+\tilde{v}(\alpha\cdot q,\lambda)(S_{-\alpha}-S_{\alpha})]E_{\alpha}\otimes H_{\alpha}-\frac{1}{2}\sum_{\alpha∈\Delta}a(\eta,\lambda)[u(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+\tilde{v}(\alpha\cdot q,\eta)(S_{-\alpha}-S_{\alpha})]H_{\alpha}\otimes E_{\alpha}, $$ now in terms of $b(\alpha\cdot q,\lambda,\eta)$ the commutator assumes the form
$$\sum_{\mu,\alpha,β∈\Delta}\left[b(-\alpha\cdot q,\eta,\lambda)E_{\alpha}\otimes E_{-\alpha},L_{\mu}I\otimes H_{\mu} +L_{\beta}(\eta)I\otimes E_{\beta} \right]=\sum_{\mu,\alpha∈\Delta}b(-\alpha\cdot q,\eta,\lambda)L_{\mu}\alpha(\mu)E_{\alpha}\otimes E_{-\alpha}-\sum_{\alpha∈\Delta}b(-\alpha\cdot q,\eta,\lambda)\times L_{\alpha}(\eta)E_{\alpha}\otimes H_{\alpha}+\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}b(-\alpha\cdot q,\eta,\lambda)N_{-\alpha,\beta}L_{\beta}(\eta)E_{\alpha}\otimes E_{-\alpha+\beta}=-\sum_{\mu,\alpha∈\Delta}b(\alpha\cdot q,\lambda,\eta)L_{\mu}\alpha(\mu)E_{\alpha}\otimes E_{-\alpha}-\frac{1}{2}\sum_{\alpha∈\Delta}b(-\alpha\cdot q,\eta,\lambda)[u(\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+\tilde{v}(\alpha\cdot q,\eta)(S_{-\alpha}-S_{\alpha})]E_{\alpha}\otimes H_{\alpha}-\frac{1}{2}\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}b(\alpha\cdot q,\lambda,\eta)N_{\alpha,\beta}[u((\alpha+\beta)\cdot q)(F_{-\alpha-\beta}-F_{\alpha+\beta})+\tilde{v}((\alpha+\beta)\cdot q,\eta)(S_{-\alpha-\beta}-S_{\alpha+\beta})]E_{\alpha}\otimes E_{\beta}, $$ and for $c(\alpha\cdot q,\lambda,\eta)$ term we also have
$$\sum_{\mu,\alpha,β∈\Delta}\left[c(\alpha\cdot q,\eta,\lambda)E_{\alpha}\otimes E_{\alpha},L_{\mu}I\otimes H_{\mu} +L_{\beta}(\eta)I\otimes E_{\beta} \right]=-\sum_{\mu,\alpha∈\Delta}c(\alpha\cdot q,\eta,\lambda)L_{\mu}\alpha(\mu)E_{\alpha}\otimes E_{\alpha}-\frac{1}{2}\sum_{\alpha∈\Delta}c(\alpha\cdot q,\eta,\lambda)\times[u(-\alpha\cdot q)(F_{-\alpha}-F_{\alpha})+\tilde{v}(-\alpha\cdot q,\eta)(S_{-\alpha}-S_{\alpha})]E_{\alpha}\otimes H_{\alpha}+\frac{1}{2}\sum_{\alpha,β∈\Delta,\alpha\neq\pm\beta}c(\alpha\cdot q,\eta,\lambda)N_{-\alpha,\beta}[u((\beta-\alpha)\cdot q)(F_{\alpha-\beta}-F_{\beta-\alpha})+\tilde{v}((\beta-\alpha)\cdot q,\eta)(S_{\alpha-\beta}-S_{\beta-\alpha})]E_{\alpha}\otimes E_{\beta}, $$ now in order to prove $\{L_{1}(\lambda),L_{2}(\eta)\}=[R_{12}(\lambda,\eta),L_{1}(\lambda)]$ $-\left[R_{21}(\eta,\lambda),L_{2}(\eta)\right]$, a comparison of the previous equations Eqs. (41)-(47), it remains to guarantee the identities below
$$ u'(\alpha\cdot q)=-a(\lambda,\eta)u(\alpha\cdot q)-b(-\alpha\cdot q,\eta,\lambda)u(\alpha\cdot q)-c(\alpha\cdot q,\eta,\lambda)u(-\alpha\cdot q), \\ \tilde{v}'(\alpha\cdot q,\lambda)=-a(\lambda,\eta)\tilde{v}(\alpha\cdot q,\lambda)-b(-\alpha\cdot q,\eta,\lambda)\tilde{v}(\alpha\cdot q,\eta)-c(\alpha\cdot q,\eta,\lambda)\tilde{v}(-\alpha\cdot q,\eta), \\ \frac{1}{2}u(\alpha\cdot q)u(\beta\cdot q)=b(-\beta\cdot q,\lambda,\eta)u((\alpha+\beta)\cdot q)-b(\alpha\cdot q,\lambda,\eta)u((\alpha+\beta)\cdot q), \\ \frac{1}{2}\tilde{v}(\alpha\cdot q)\tilde{v}(\beta\cdot q)=b(-\beta\cdot q,\lambda,\eta)\tilde{v}((\alpha+\beta)\cdot q,\lambda)-b(\alpha\cdot q,\lambda,\eta)\tilde{v}((\alpha+\beta)\cdot q,\eta), \\ \frac{1}{2}u(\alpha\cdot q)u(\beta\cdot q)=-c(\beta\cdot q,\lambda,\eta) u((\alpha-\beta)\cdot q) -c(\alpha\cdot q,\lambda,\eta)u((\beta-\alpha)\cdot q), \\ \frac{1}{2}\tilde{v}(\alpha\cdot q)\tilde{v}(\beta\cdot q)=-c(\beta\cdot q,\lambda,\eta)\tilde{v}((\alpha-\beta)\cdot q,\lambda)-c(\alpha\cdot q,\lambda,\eta)\tilde{v}((\beta-\alpha)\cdot q,\eta), $$ the proofs of these equations are elementary and direct, so we do not want to explain them here any more.
6 Quantization Case
In this section we mainly illustrate the new result of this paper that although the classical and quantum treatments differ in many aspects, the procedure above indeed is applicable to quantized counterpart of Eq.(9), which plays a vital role in the formulation of the quantum inverse scattering problem and to explore the quantum integrability of the corresponding dynamic system, we give the following more detailed description. It is well known that in order to quantize the classical system one should replace the Poisson bracket $\{ , \}$ with the Dirac commutator bracket $-i[\quad]$. In this situation the dynamic functions $L_{\mu},L_{\alpha}$ all being operators and because of the non-commutativity of the operators, the commutation relation between $R$-matrix and Lax matrix $L_{1,2}$ is a slightly different from classical case. To be specific, the commutators in quantized case are expressible in the form of
$$ \sum_{\mu,\alpha,β∈\Delta}\left[b(\alpha\cdot q,\lambda,\eta)E_{\alpha} \otimes E_{-\alpha},L_{\mu}H_{\mu}\otimes 1\right] =\sum_{\mu,\alpha∈\Delta}b(\alpha\cdot q,\lambda,\eta)L_{\mu}\left[E_{\alpha},H_{\mu}\right]\otimes E_{-\alpha} +i\sum_{\mu,\alpha∈\Delta}b'(\alpha\cdot q,\lambda,\eta)\alpha(\mu)H_{\mu}E_{\alpha}\otimes E_{-\alpha} =-\sum_{\mu,\alpha∈\Delta}b(\alpha\cdot q,\lambda,\eta)L_{\mu}\alpha(\mu)E_{\alpha}\otimes E_{-\alpha}+i\sum_{\alpha∈\Delta}b'(\alpha\cdot q,\lambda,\eta)H_{\alpha}E_{\alpha}\otimes E_{-\alpha}, $$ $$\sum_{\mu,\alpha,β∈\Delta}\left[b(-\alpha\cdot q,\eta,\lambda)E_{\alpha} \otimes E_{-\alpha},L_{\mu}I\otimes H_{\mu} \right]=\sum_{\mu,\alpha∈\Delta}b(-\alpha\cdot q,\eta,\lambda)L_{\mu}\alpha(\mu) E_{\alpha}\otimes \left[E_{-\alpha},H_{\mu} \right]-i\sum_{\mu,\alpha∈\Delta}b'(-\alpha\cdot q,\eta,\lambda)\alpha(\mu) E_{\alpha} \otimes H_{\mu}E_{-\alpha} =\sum_{\mu,\alpha∈\Delta}b(-\alpha\cdot q,\eta,\lambda)L_{\mu} \alpha(H_{\mu})E_{\alpha}\otimes E_{-\alpha} +i\sum_{\alpha∈\Delta}b'(-\alpha\cdot q,\eta,\lambda)E_{\alpha} \otimes H_{-\alpha}E_{-\alpha}. $$ To proceed we firstly recall that $H_{\mu},E_{\alpha}$ constitute a complete basis of the simple Lie algebra $\mathfrak{g}$ and from $\left[H_{\mu},E_{\alpha}\right]=\alpha(\mu)E_{\alpha}$ we have
$$ \mathrm{Tr}(H_{\alpha}E_{\alpha})=\sum_{\mu=1}^{n} \alpha(\mu)\mathrm{Tr}(H_{\mu}E_{\alpha})=\sum_{\mu=1}^{n}\mathrm{Tr}(H_{\mu}\left[H_{\mu},E_{\alpha}\right])=0, $$ which means $H_{\alpha}E_{\alpha}$ belonging to $\mathfrak{g}$, hence for every $\alpha$ one can rewrite $H_{\alpha}E_{\alpha}$ explicitly as
$$ H_{\alpha}E_{\alpha}=\sum_{\mu=1}^{n}C_{\mu} H_{\mu}+\sum_{β∈\Delta}C_{\beta}E_{\beta}, $$ here $C_{\mu},C_{\beta}$ are some coefficients and there must $\exists \nu$ satisfies $\alpha(\nu)\neq 0$ which immediately leads to
$$ \alpha(h_{\nu})H_{\alpha}E_{\alpha}=\left[H_{\nu},H_{\alpha}E_{\alpha}\right]=\sum_{β∈\Delta}C_{\beta}\beta(\nu)E_{\beta}, $$ combining Eqs. (52) with (53) we find $C_{\mu}=0$ for all $\mu$.
In a similar spirit, for $\forall β∈\Delta,\beta\neq \alpha$, there exist $\nu$ satisfies $(\alpha-\beta)(\nu)\neq 0$, once again applying $[H_{\nu},H_{\alpha}E_{\alpha}]=\sum_{β∈\Delta}C_{\beta}[H_{\nu},E_{\beta}]$ we infer $(\alpha-\beta)(\nu)C_{\beta}=0$, which implies $C_{\beta}=0$. It is then seen that for every $\alpha$ one obtains $H_{\alpha}E_{\alpha}=C_{\alpha}E_{\alpha}$ for some constant $C_{\alpha}$ and imposition of this idea on other cases we confirm
$$ E_{\alpha}H_{\alpha}=\tilde{C}_{\alpha}E_{\alpha},\quad H_{-\alpha}E_{-\alpha}=C_{-\alpha}E_{-\alpha}, \\ E_{-\alpha}H_{-\alpha}=\tilde{C}_{-\alpha}E_{-\alpha}, $$ here $\tilde{C}_{\alpha}$ are the corresponding coefficients. Thus if expanding the left hand side of commutation relation $[H_{\alpha}E_{\alpha},E_{-\alpha}]=C_{\alpha}H_{\alpha}$ we have
$$ H^{2}_{\alpha}-\sum_{\mu=1}^{n}\alpha^{2}(\mu) E_{-\alpha}E_{\alpha}=C_{\alpha}H_{\alpha}, $$ as well as
$$ H^{2}_{-\alpha}-\sum_{\mu=1}^{n}\alpha^{2}(\mu)E_{\alpha} E_{-\alpha}=-C_{-\alpha}H_{\alpha}, \\ H^{2}_{\alpha}-\sum_{\mu=1}^{n}\alpha^{2}(\mu)E_{\alpha}E_{-\alpha}=\tilde{C}_{\alpha}H_{\alpha}, \\ H^{2}_{-\alpha}-\sum_{\mu=1}^{n}\alpha^{2} (\mu)E_{-\alpha}E_{\alpha}=-\tilde{C}_{-\alpha}H_{\alpha}, $$ putting all these results together one deduces that
$$ C_{\alpha}+C_{-\alpha}=\sum_{\mu=1}^{n}\alpha^{2}(\mu),\quad \quad C_{\alpha}-\tilde{C}_{\alpha}=\sum_{\mu=1}^{n}\alpha^{2}(\mu)\,,C_{-\alpha}-\tilde{C}_{-\alpha}=\sum_{\mu=1}^{n}\alpha^{2}(\mu),\quad -\tilde{C}_{-\alpha}-\tilde{C}_{\alpha} =\sum_{\mu=1}^{n}\alpha^{2}(\mu)\,,$$ with the solution of the constants
$$ C_{\alpha}=C_{-\alpha}=\frac{1}{2}\sum_{\mu=1}^{n}\alpha^{2}(\mu),\quad \\ \tilde{C}_{-\alpha}=\tilde{C}_{\alpha}=-\frac{1}{2}\sum_{\mu=1}^{n}\alpha^{2}(\mu),\quad \forall\alpha∈ \Delta, $$ now from the previous detailed analysis one claims that
$$ H_{\alpha}E_{\alpha}\otimes E_{-\alpha}=C_{\alpha}E_{\alpha}\otimes E_{-\alpha}=E_{\alpha}\otimes C_{-\alpha} E_{-\alpha} \\ =E_{\alpha}\otimes H_{-\alpha}E_{-\alpha}, $$ and assemble all these relations together with Eqs. (49), (50) we have
$$ \sum_{\mu,\alpha,β∈\Delta}\left[b(\alpha\cdot q,\lambda,\eta)E_{\alpha}\otimes E_{-\alpha},L_{\mu}H_{\mu}\otimes 1\right]\\ =\sum_{\mu,\alpha,β∈\Delta}\left[b(-\alpha\cdot q,\eta,\lambda)E_{\alpha}\otimes E_{-\alpha},L_{\mu}I\otimes H_{\mu} \right], $$ also we get
$$ \sum_{\mu,\alpha,β∈\Delta}\left[c(\alpha\cdot q,\lambda,\eta)E_{\alpha}\otimes E_{\alpha},L_{\mu}H_{\mu}\otimes 1 \right]\\ =\sum_{\mu,\alpha,β∈\Delta}\left[c(\alpha\cdot q,\eta,\lambda)E_{\alpha}\otimes E_{\alpha},L_{\mu}I\otimes H_{\mu} \right], $$ now by using of these identities, the other terms in Dirac commutator $[L_{1}(\lambda),L_{2}(\eta)]$ and $[R_{12}(\lambda,\eta),L_{1}(\lambda)],$ $[R_{21}(\eta,\lambda),L_{2}(\eta)]$ are the same as in classical Poisson bracket $\{L_{1}(\lambda),L_{2}(\eta)\}$ and $[R_{12}(\lambda,\eta),L_{1}(\lambda)],[R_{21} (\eta,\lambda),$ $L_{2}(\eta)]$, therefore it is necessary to comment on that in the quantized version of the two types spin Calogero-Moser Hamiltonian system, the $R$-matrix coincides with the original form as in classical case.
7 Summary
The constructions of universal Lax pair and $R$-matrix for the Calogero-Moser system with two type interacting spins formulated in terms of root system of any semi-simple Lie algebras are illustrated in this paper. The primary difference compared to the ordinary spin Calogero-Moser model is that we generalize the Poisson brackets of spin dynamic variables, which are deeply related to the structure constants of root system of the semi-simple Lie algebra. Under this postulation, a detailed exposition of consistency of the Lax pair with spectral parameter is given in the case of the trigonometric type of interaction potential in a direct way, of course, one can apply the similar calculations to the hyperbolic potential straightforwardly. Thus based on the notion of Lax pair, we explore the universal R-matrix in classical case and argue that it is equivalent to the Poisson brackets between original Lax matrix with different spectral parameters. For the quantized model as opposed to the classical case, taking into account of the completion of the root system basis, we demonstrate that the $R$-matrix behaves the same manner as in the classical situation, which is the major conclusion of this paper. We should mention here that the essential ingredient of the construction of Lax pair for the two type interacting spins heavily relying on the elegant property of trigonometric function, nevertheless, this idea is not available to the rational and elliptic interactions. Thus it is plausible to pose new expectations of the further work that how to give an analogous and precise description of two type spins Calogero-Moser system with rational or elliptic potential. Along this line, it appears that this viewpoint of universal construction may be employed to the other CM model such as the generalization of supersymmetric spin Calogero-Moser systems with two type interacting spins or deformed quantum Calogero-Moser systems. Further investigation in this direction is a challenging problem.
S.Khoroshkin, M.Matushko, E.Sklyanin , On Spin Calogero-Moser System at Infinity},arXiv:math-ph/1608. 00599. [Cited within: 1]
I.Krichever, O.Babelon, E.Billey, M.Talon , Spin generalization of the Calogero-Moser system and the Matrix KP equation, arXiv:hep-th/9411160. [Cited within: 1]
V. Pashkov and A.Zabrodin, Spin Generalization of the Calogero-Moser Hierarchy and the Matrix KP Hierarchy,arXiv:math-ph/1712. 06943. [Cited within: 1]
R. J.Szabo, Equivariant Cohomology and Localization of Path Integrals, Springer New York (2000). URL [Cited within: 1] ABSTRACT Publisher’s description: This book, addressing both researchers and graduate students, reviews equivariant localization techniques for the evaluation of Feynman path integrals. The author gives the relevant mathematical background in some detail, showing at the same time how localization ideas are related to classical integrability. The text explores the symmetries inherent in localizable models for assessing the applicability of localization formulae. Various applications from physics and mathematics are presented.
N. M.J.Woodhouse. , Geometric quantization, Oxford University Press New York ( 1997). [Cited within: 1]
L.Feher and B.G Pusztai, Lett. Math. Phys. 79(2007) 263. [Cited within: 1]
S.Kharchev, A.Levin, M.Olshanetsky, A.Zotov , Quasi-Compact Higgs Bundles and Calogero-Sutherland Systems with Two Types Spins, arXiv: 1712. 08851. [Cited within: 2]
,Department of Physics, Zhejiang University, Hangzhou 310027, China
