On the Integrable Deformations of the Maximally Superintegrable Systems

: In this paper, we present the integrable deformations method for a maximally superintegrable system. We alter the constants of motion, and using these new functions, we construct a new system which is an integrable deformation of the initial system. In this manner, new maximally superintegrable systems are obtained. We also consider the particular case of Hamiltonian mechanical systems. In addition, we use this method to construct some deformations of an arbitrary system of ﬁrst-order autonomous differential equations.


Introduction
The construction of integrable deformations of a given integrable system was the subject of some recent articles. In [1], taking into account the fact that Poisson algebras endowed with a co-product map give rise to a systematic way of constructing integrable systems, integrable deformations of a class of three-dimensional Lotka-Volterra equations were given. In [2], considering Poisson-Lie groups as deformations of Lie-Poisson (co)algebras, integrable deformations of certain integrable types of Rössler and Lorenz systems were presented. Based on the construction of a family of compatible Poisson structures, in [3], a family of integrable deformations of the Bogoyavlenskij-Itoh systems was constructed. In [4], altering the constants of motion, integrable deformations of the Euler top were constructed. In the same manner, in [5][6][7][8][9], integrable deformations of some three-dimensional systems were obtained. Moreover, in [10], the integrable deformations method for a three-dimensional system of differential equations was presented. The property that allows to construct integrable deformations of the above-mentioned three-dimensional systems is their maximal superintegrability.
Both in classical mechanics and in quantum mechanics, the superintegrable systems have been widely investigated (see, e.g, [11][12][13] and references therein). In classical mechanics, a superintegrable system on a 2n-dimensional phase space is a completely integrable Hamiltonian system which possesses more functionally independent first integrals than degrees of freedom. Moreover, such a system is called maximally superintegrable if the number of the independent first integrals is 2n − 1 (see, e.g., [14,15]). Similarly, for an arbitrary natural number n, a system of first-order differential equations on R n which has n − 1 functionally independent constants of motion is called a maximally superintegrable system.
In this paper, we present some integrable deformations of maximally superintegrable systems. Firstly, using n − 1 functionally independent functions, we construct a family of maximally superintegrable systems for which the considered functions are constants of motion, and we point out a Hamilton-Poisson realization of such a system. Secondly, the Hamilton-Poisson realization allows us to give integrable deformations of the obtained system, which are also maximally superintegrable Hamilton-Poisson systems. Moreover, under certain conditions, the obtained system and its integrable deformations are generalized Nambu-Hamilton systems [16]. In Section 3 we analyze the particular case of Hamiltonian mechanical systems. In Section 4 we extend the method given in [10] to the case of n−dimensional systems of first-order differential equations.

The Integrable Deformations Method for Maximally Superintegrable Systems
In this section, using the idea given in [4], we give integrable deformations of a maximally superintegrable system. For this purpose, we present a construction of a maximally superintegrable system considering n − 1 functionally independent functions as its constants of motion.
Let Ω ⊆ R n be on an open set. Consider the functions C 1 , C 2 , . . . , with . . , C n−1 are n − 1 functionally independent functions on Ω.
Considering that these functions are constants of motion, that is,Ċ k = dC k dt = 0 for every k ∈ {1, 2, . . . , n − 1}, it follows that: (2) This linear system has a unique solution given bẏ where From (1) and (4), it is easy to see that ∆ n = ∆. Setẋ n = ∆. Therefore, the functions C 1 , . . . , C n−1 give rise to the following differential system on R nẋ Considering the Poisson bracket on R n generated by C 1 , C 2 , . . . , C n−2 (see, e.g., [17]), namely, for every f , g ∈ C ∞ (R n , R), system (5) readṡ Thus, it is a Hamilton-Poisson system with the Hamiltonian H := C n−1 and Casimirs C 1 , C 2 , . . . , C n−2 . We notice that {C i , C j } C 1 ,...,C n−2 = 0 for every i, j ∈ {1, 2, . . . , n − 1}. Now, we setẋ n = ν∆ = ν∆ n , where ν ∈ C 1 (Ω, R) is an arbitrary function. Then system (3) becomesẋ and it is obvious that it has the same constants of motion C 1 , . . . , C n−1 . Furthermore, System (8) is the Hamilton-Poisson systeṁ with the Hamiltonian H := C n−1 and the Poisson bracket (see, e.g., [18]) given by In this manner, a family of maximally superintegrable systems with the same constants of motion was constructed. On the other hand, for a given maximally superintegrable system, there is a unique function ν such that the system takes the form (9) [17,18] (see the next Remark). The function ν is called the rescaling function [18] and it is usually defined on an open and dense subset of Ω.

Remark 1.
In the papers [17,18], a systemẋ = X(x) on R n is considered. For this system, it is assumed there are n − 1 functionally independent constants of motion H 1 , . . . , H n−1 . In [17] it was proven that X = (X 1 , . . . , X n ) is orthogonal to the vectors ∇H k , k ∈ {1, 2, . . . , n − 1}, and consequently there is a function µ such that In [18], using the orthogonality of the vector fields X and ∇H k , k ∈ {1, 2, . . . , n − 1}, it was obtained that X is given as the vector field (∇H 1 ∧ . . . ∧ ∇H n−1 ) multiplied by a C 1 real function, where denotes the Hodge star operator for multivector fields.
In the final part of this section, we give integrable deformations of system (9).
. . , f n ), we understand a systemẋ = f(x) + u(x), where u = (u 1 , u 2 , . . . , u n ) and each function u i , i ∈ {1, 2, ..., n} depends on some real parameters such that when these parameters vanish, then u vanishes. When the initial system and its deformation are maximally superintegrable, then we talk about an integrable deformation.
In the next result, we give a family of integrable deformations of a maximally superintegrable system. Theorem 1. Consider system (9) generated by the constants of motion C 1 , . . . , C n , a rescaling function ν, and the Poisson bracket (10).

Remark 2.
In order to apply the above Theorem to a maximally superintegrable . . , f n ), first we determine the rescaling function ν, as in Remark 1. Secondly, we alter the constants of motion and then apply (12). The new system (12) can be considered as the controlled initial system, where deformation parameters become control parameters. By choosing appropriate deformation functions α 1 , . . . , α n−1 , the dynamics of the initial system may be controlled in a desired way. In addition, the methods of geometric mechanics are applicable to such systems.

Remark 3.
Every three-dimensional Hamilton-Poisson system is a maximally superintegrable system; therefore, Theorem 1 gives integrable deformations for such systems. On the other hand, System (5) may be expressed in the Nambu mechanics formalism [16,19], namely,ẋ where the canonical Nambu bracket is given by Under some supplementary conditions, for exampleν = 0 (see, e.g., [20]), the bracket is also a Nambu bracket, therefore system (9) is a generalized Nambu-Hamilton system. Consequently, we have the following integrable deformation of this Nambu system: Remark 5. In order to preserve the Poisson bracket (10) generated by C 1 , . . . , C n−2 , we alter only the Hamiltonian C n−1 , namely, I n−1 = C n−1 + pα n−1 . In this case, the integrable deformation (12) becomeṡ where p is a deformation parameter.

Example 1.
Consider the periodic Volterra system (see, e.g., [21]) that describes an interaction between only the close neighboring variables. This system has the Hamiltonian function H = x 1 + x 2 + x 3 + x 4 and the constants of motion C 1 = x 2 x 4 and C 2 = x 1 x 3 . Therefore, it is maximally superintegrable. Furthermore, this system is the Nambu systeṁ where the Nambu bracket is given by (17) { Consider the deformation parameters g i ∈ R and deformation functions α i , i ∈ {1, 2, 3}. Denote I 1 = C 1 + g 1 α 1 , I 2 = C 2 + g 2 α 2 ,H = H + g 3 α 3 . Therefore, systeṁ gives an integrable deformation of the considered system.

Integrable Deformations of a Maximally Superintegrable Hamiltonian System
In this section, we discuss some integrable deformations of a Hamiltonian system on a 2n-dimensional phase space in the case when it possesses 2n − 1 functionally independent integrals of motion.
Let M ⊂ R 2n be a smooth manifold and let (M, ω, H) be a Hamiltonian mechanical system, where ω = n ∑ i=1 dp i ∧ dq i and H is the Hamiltonian. The equations of motion are where {., .} ω is the Poisson bracket on R 2n , that is, e.g., [22,23]). Consider that C 1 , C 2 , . . . , C 2n−2 , H = C 2n−1 are functionally independent integrals of motion of system (20), thus,Ċ j = {H, C j } ω = 0, j ∈ {1, 2, . . . , 2n − 2}. System (20) has the formẋ = f(x), where x = (x 1 , x 2 , . . . , x 2n ), f = ( f 1 , f 2 , . . . , f 2n ), and the procedure used in the previous section gives integrable deformations of system (20). It is natural to ask whether such a deformation is a Hamiltonian mechanical system. In addition, if the Hamiltonian mechanical system (20) is completely integrable in Arnold-Liouville's sense, then it is maximally superintegrable [15] and, consequently, the following question arises: is the above-mentioned deformation a maximally superintegrable Hamiltonian mechanical system? To give a partial answer to this questions, we consider the case when only the Hamiltonian is deformed.
First, we perturb the Hamiltonian function H; that is, we consider the new Hamiltonian functionH given byH = H + kα, where k ∈ R and α is an arbitrary differentiable function on M such that C 1 , C 2 , . . . , C 2n−2 ,H are functionally independent. Then we obtain the following Hamiltonian mechanical system which is a deformation of system (20). In general, we notice that the equality {H, C j } ω = 0 does not hold. On the other hand, taking into account Remark 1, there is a function ν such that ν · ∂(C 1 , . . . , C 2n−2 , H, q i ) ∂(p 1 , . . . , p n , q 1 , . . . , q n ) Using the Poisson bracket (10) Therefore, via Theorem 1, an integrable deformation of system (22) in the casẽ H = H + kα is given by (21) and (23) represent two deformations of the same Hamiltonian system (20) obtained by a perturbation of the Hamiltonian H. In the first case, system (21) is still a Hamiltonian mechanical system, but it does not preserve the other constants of motion. In the second case, system (23) preserves the integrals of motion, and it is a Hamilton-Poisson system; however, generally, it is not a Hamiltonian mechanical system with the Hamiltonian functionH. We notice that if the deformation function α is chosen such that

Example 2.
Consider the Hamiltonian function H(q 1 , q 2 , p 1 , where U, V are smooth functions. The equations of motion are given bẏ It immediately follows that the functions C 1 = q 1 + q 2 2 and C 2 = p 1 + p 2 are first integrals of the above system. This Hamiltonian mechanical system is maximally superintegrable. In addition, it can be written in the form (22) with ν = 1.
Particularly, if we alter only the Hamiltonian function, that is, g 1 = g 2 = 0, then the above integrable deformation becomeṡ By Remark 6, we deduce that this particular integrable deformation is in fact a Hamiltonian mechanical system with the Hamiltonian functionH = I 3 = H + g 3 α 3 if, and only if α 3 (q 1 , q 2 , p 1 , p 2 ) = W(q 1 − q 2 , p 1 − p 2 ), where W is a smooth function. In this case, the above integrable deformation is a maximally superintegrable Hamiltonian mechanical system.

Remark 7.
While the deformation (21) of system (20) can be obtained in one way, namely by a perturbation of the Hamiltonian function, integrable deformations of the same system, written in the form (22), can be obtained in many more ways by altering one or more of the functions C 1 , C 2 , . . . , C 2n−1 .

Some Deformations of an Arbitrary System of First-Order Autonomous Differential Equations
In [10], some deformations of three-dimensional systems of differential equations were obtained. First, a Hamilton-Poisson part of the considered system is identified. Then an integrable deformations of this part is constructed by the alteration of its constants of motion. Lastly, the non-Hamilton-Poisson part and the obtained integrable deformation are put together to obtain a new system, that is, a deformation of the initial system. In this section we extend this method to systems of differential equations on R n , n > 3. More precisely, the same method works if we can identify a maximally superintegrable part of the considered system. Consider a dynamical system of the forṁ where x = (x 1 , x 2 , . . . , x n ), f = ( f 1 , f 2 , . . . , f n ). We assume that system (24) can be written in the formẋ where the systemẋ = s(x) is maximally superintegrable. Then, using the notations from Section 2, system (26) takes the form (9) and one of its integrable deformations is given by Theorem 1. We denote this deformation byẋ = s(x) +s(x).
Adding h(x) to the right side, we obtain the following deformation of system (24) Note that in the case of maximally superintegrable systems, some constants of motions can be preserved. If a system has a constant of motion, but it is not maximally superintegrable, its above-mentioned deformation does not generally preserve that constant of motion. It is natural to ask whether there are deformations that preserve a constant of motion.

Example 3.
Another Lotka-Volterra type system (see, e.g., [21]) is given bẏ The function H = x 1 + x 2 + x 3 + x 4 is a constant of motion of this system. Apparently, this system is not maximally superintegrable. However, we can consider a maximally superintegrable part of it, namely,ẋ Indeed, the functions C 1 = x 1 , C 2 = x 2 + x 3 , and C 3 = x 2 x 4 are constants of motion. Moreover, the above system takes the forṁ where ν = x 3 and the bracket is given by (10). If we consider I 1 = C 1 + g 1 α 1 , I 2 = C 2 + g 2 α 2 , , then an integrable deformation of the maximally superintegrable part of the considered system is given by Theorem 1 Now we bring back the unused terms into the system and obtain the following deformation of the considered systeṁ A particular deformation is obtained if we consider particular deformation functions. For example, if we take g 1 = g 3 = 0, and α 2 (x 1 , x 2 , x 3 , x 4 ) = β(x 2 ), then we geṫ It is easy to see that H is not a constant of motion of the above system. Nevertheless, we can find some deformation functions such that H will be a constant of motion. For example, if we choose g 1 = g 2 = g, g 3 = 0 and α 1 (x 1 , x 2 , x 3 , x 4 ) = α(x 4 ), α 2 (x 1 , x 2 , x 3 , x 4 ) = β(x 4 ), then we obtaiṅ We obtainḢ = gx 3 x 4 (α (x 4 ) + β (x 4 )).
Therefore, if α(x 4 ) + β(x 4 ) = constant, then H is a constant of motion of the deformed system.

Conclusions
Integrable deformations of a maximally superintegrable system were obtained by the alteration of their constants of motion. The new systems are also maximally superintegrable. Moreover, this method can be used to perturb a maximally superintegrable part of an arbitrary system of first-order differential equations, which gives a deformation of the considered system.
Note that the systems obtained by deformations depend on more parameters than the initial system. Therefore, they can exhibit rich dynamics and can display some types of bifurcations, such as Hopf, zero-Hopf, Bogdanov-Takens, Bautin, or even three codimension bifurcations. Moreover, the functions added to a system by deformation can be viewed as control functions, and as a consequence, some natural queries arise, such as: how the stability changes, whether they can be stabilized by some states/orbits, how the periodic motion is affected, and whether new chaotic systems can be obtained or whether some chaotic trajectories can be stabilized, and so forth.