Background
The (rational) Calogero-Moser system is the dynamical system which describes the evolution of $n$ particles on the line $\mathbb{C}$ which repel each other with force proportional to the cube of their distance. If the particles have (distinct!) position $q_i$ and momentum $p_i$, then the Hamiltonian which describes this system is $$ H=\sum_i p_i^2+\sum_{i\neq k}\frac{1}{(x_i-x_k)^2} $$ There are many interesting properties of this system, but one of the first interesting properties is that it is `completely integrable'. This means that solving it explicitly amounts to solving a series of straight-forward integrals.
The integrability can most easily be shown by showing that the phase space for this system includes into a symplectic reduction of a certain matrix space, and then noticing that the above Hamiltonian is a restriction of a integrable Hamiltonian on the whole space. This is done by assigning to any ensemble of points $q_i$ and momenta $q_i$ a pair of $n\times n$ matrices $X$ and $Y$, where $X$ is the diagonal matrix with $q_i$ on the diagonal entries, while $Y$ is given by $$ Y_{ii}=p_i, \; Y_{ik}=(x_i-x_k)^{-1}, \; i\neq k $$ This matrix assignment defines a map from the configuration space $CM_n$ of the CM system to the space of pairs of matrices. The space of pairs of matrices $(X,Y)$ is naturally a symplectic space from the bilinear form $(X,Y)\cdot (X',Y')=Tr(XY')-Tr(X'Y)$, and the action of $GL_n$ by simultaneous conjugation naturally has a moment map. Therefore, we sympletically reduce the space of pairs of matrices at a specific coadjoint orbit (not the origin) and get a new symplectic space $\overline{CM}_n$.
Composing the above matrix assignment with symplectic reduction, we get a map $CM_n\rightarrow \overline{CM}_n$. This map turns out to be a symplectic inclusion which has dense image. We also notice that the functions $Tr(Y^i)$, as $i$ goes from $1$ to $n$, descend to a Poisson-commuting family of functions on $\overline{CM}_n$, and because $\overline{CM}_n$ is $2n$ dimensional, each of the functions $Tr(Y^i)$ gives an integrable flow on $\overline{CM}_n$. Finally, we notice that $Tr(Y^2)$ restricts to $H$ on $CM_n$.
The Massive Version of the CM System
Now, make the following change to the system. To every particle, assign a number $m_i$ (the mass), which can be in $\mathbb{C}$, but I am interested in the case where the $m_i$ are positive integers. Define a the massive CM Hamiltonian as $$H_m=\sum_i\frac{p_i^2}{m_i}+\sum_{i\neq k}\frac{m_im_k}{(x_i-x_k)^2} $$ The physical meaning of this equation is that particles still have force proportional to the inverse of the cube of their distance, but the force is proportional to the mass of that particle; also, particles resist acceleration proportional to their mass. If the force were to drop off proportional to the inverse square of their distance, and attract instead of repel, this would model how massive particles move under the influence of gravity.
Questions
- Is this system integrable?
- Can it be realized in a similar matrix form?
- Does it have any interesting or new behavior than the usual CM system?
An Idea
It is almost possible to realize this Hamiltonian in a simple modification of the previous approach. Let $M$ denote the diagonal matrix with the $m_i$s on the diagonal. Then $$Tr(MYMY)=\sum_im_i^2p_i^2+\sum_{i\neq k}\frac{m_im_k}{(x_i-x_k)^2}$$ The functions $Tr( (MY)^i)$ should again be a Poisson commutative family. Rescaling the $p_i$ by $m_i^{3/2}$ gives the massive Hamiltonian $H_m$; however, this rescaling is not symplectic, and so it won't preserve the flows.
Another Idea
In the case of integer $m_i$, one possibility is to work with $N\times N$ matrices rather than $n\times n$ matrices, where $N=\sum m_i$. Then it is possible to construct a matrix $X$ with eigenvalues $q_i$, each occuring with multiplicity $m_i$, as well as a matrix $Y$ such that $(X,Y)$ defines a point in $\overline{CM}_N$. The Hamiltonian $Tr(Y^2)$ even restricts to the correct 'massive' Hamiltonian $H_m$. However, the flow described by this Hamiltonian on $\overline{CM}_n$ will in almost all cases immediately separate eigenvalues that started together, which we don't want. If we restrict the Hamiltonian to the closed subspace where the eigenvalues are required to stay together, then this gives the desired flow. Unfortunately, restricting to a closed subvariety doesn't preserve a Hamiltonian being integrable.