The problem of finding solutions to quantum integrable finite-dimensional systems in many cases still remains unsolved. The empirical constructions of such solutions were important at early stages of the development of the theory of these systems and lead to many important results being applied to trigonometric Calogero-Sutherland-Moser systems with the Hamiltonian of the form
for N particles in one dimension with the two-body potential given by
The coupling constant g is supposed to be real and chosen as .
It turned out   that the ground-state wave function of the trigonometric model is of factorized form and the wave functions of all the excitations can be written as products of this function and multivariable Jack polynomials  . These results were also extended to more complicated cases with interaction terms modified by the introduction of more general root systems  .
The quantum elliptic many-particle problem which also has been proven to be integrable   is till now quite far from being solved completely. It has two-particle interaction potential of the form
where is the Weierstrass elliptic function with two periods which do not lie on the line in complex plane. The hermiticity of the Hamiltonian implies . The trigonometric case (2) corresponds to infinite complex period.
In the simplest case of , the eigenvalue problem for the Hamiltonian (1), (3) is just the well-known Lamé equation. For general coupling constants g, its solutions have a branch point at the origin and their expression through the known transcendental functions is not known. However, there is an exception: for integer values of g, Hermite found that the solutions are expressed in terms of an exponent and quasiperiodic Weierstrass functions  . The reason for the existence of such a solution is in fact based on its “good” analytic properties in a complex plane of the variable q: at integer g there is no branch points and the only singularity is a pole at up to the quasiperiodicity.
This fact inspired the authors in the paper  to consider the case of general and . It has been proved that the double quasiperiodic solutions for many-particle wave functions are still expressed in terms of the Weierstrass functions but the procedure of finding them is rather complicated. They were able to find it explicitly only for , . In  , these solutions have been presented analytically for arbitrary and , also in overcomplicated form requiring many nontrivial operations to their explicit writing. As for arbitrary real the solution of the eigenproblem for the elliptic case was constructed by the perturbation theory in the form of infinite series  .
However, there is another approach to finding the solutions for the dynamics of integrable systems, namely separation of variables. It is well known in its simple form using purely co-ordinate transformations. As for the elliptic Calogero-Moser systems, simple forms do not work but separation still takes place as it was proposed in  for 3-particle case at arbitrary values of coupling constant g. The separation of variables occurs after transformation corresponding to a classical canonical transformation of phase space variables mixing coordinates and momenta. The transformation is realized as an integral transform of the wave function in the quantum case. The original two-dimensional problem has been reduced to one-dimensional one and the process of finding the eigenfunction contains investigation of the solution to a third order ordinary differential equation  ,
where are constants (the values of the integrals of motion).
The aim of this paper is to find the explicit solutions to equation (4). We shall show below that for integer values of g, they may be obtained via the solution to the system of g usual transcendental equations.
2. Finding the solution
It should be noted at first that the coefficients in (4) are double periodic functions of x. Hence one try to seek the particular solution as double quasiperiodic function. The second observation concerns possible singularities of this solution. Since has double pole at and it is analytic at the other points of torus , all the singularities of in must be also located at . The assumption
at gives the result
as the leading singularity, taking the most singular possibility only. For non-integer g, this is a branch point and there is no simple ansatz to the solution of (4).
The situation is changed drastically if , . In this case, the leading singularity of is a pole of the order and there are no branch points. Combining this property with double quasiperiodicity allows one to write down the Hermite-like ansatz for the possible solution to (4)
where A is inessential normalization constant, and are parameters which have to be determined, and is the Weierstrass sigma function. It is connected with by the relations
where is the Weierstrass function with the property
We assume that all are mutually different for and different from 0 in .
By consecutive differentiations of (5), one finds
Note that the all right-hand sides of these equalities are elliptic functions of the argument x. Substitution of these expressions into (4) yields
The function is elliptic and might have poles up to third order at the points . However, the direct inspection of the Laurent decompositions near these points shows that all the coefficients at the terms vanish identically for arbitrary and . Hence this function can be written in the form
where the constant coefficients should obey the relation
(statement (III) of par. 20.12 in  , e.g.). The Laurent decomposition of (7) near the points with the use of (6) allows one to find the coefficients explicitly. Due to (7), all of them should vanish. This results in the system of transcendental equations to the parameters :
It remains only to calculate the constant term in (8). Equivalently, we calculate
using the Laurent decomposition of (7) near the point . After long but straightforward calculations (performed by the MATHEMATICA® program), one finds the condition
The algebraic system (10)-(11) allows one to determine the parameters under which the elliptic function B has no poles and equals zero at one point. Then due to the Liouville theorem (statement (IV) of par. 20.12 in  ). The last equation is cubic in γ. This corresponds to three linearly independent solutions to the original Equation (4).
Let us summarize our results. We obtained the explicit solutions of the separated equation (4) at integer couplings g which, in its turn, gives the solution to the three-particle quantum Calogero-Moser problem via the procedure described in  . We conjecture that g equations (10)-(11) determine the g parameters in the generic case at the least. However, it is not clear whether the solution to the above problem in the forms known before   can be transformed into the forms with separated variables. As    , we consider in general singular solutions to the differential Equation (4) leaving aside the right physical boundary conditions which are even not known here  .
The work was supported by the Votruba-Blokhintsev CR-JINR cooperation program in theoretical physics, Czech Science Foundation project 17-01706S and NPI CAS institutional support RVO 61389005.
 Dittrich, J. and Inozemtsev, V.I. (1993) On the Structure of Eigenvectors of Multidimensional Lamé Operator. Journal of Physics A: Mathematical and General, 26, L753-L756.
 Felder, G. and Varchenko, A. (1995) Integral Representation of Solutions of the Elliptic Knizhnik-Zamolodchikov-Bernard Equations. International Mathematics Research Notices, 5, 221-233.
 Komori, Y. and Takemura, K. (2002) The Perturbation of the Quantum Calogero-Moser-Sutherland System and Related Results. Communications in Mathematical Physics, 227, 93-118.