In recent years, the study of dynamic systems has been undertaken in a wide range of fields, which are located at the crossroads of differential geometry, algebraic geometry, number theory, Lie algebra, Intellectual and material means.
A considerable renewed interest appeared for Hamiltonian dynamic systems with two degrees of freedom, one of which is the study of the topological properties of the flow of these systems, their integrability, their chaotic behavior, the study of Periodic solutions and their bifurcation.
A first fundamental task in this field was the search for integrable systems which give rise to non chaotic behavior. For a Hamiltonian system with n degrees of freedom, the most general definition is that of Liouville. In addition to direct analytical methods, various criteria are developed to determine candidates for integrability, namely the Painlevé criterion, the Ziglin criterion and the Poincaré sections.
Integrability is clearly a central issue in understanding the origins and implications of the behaviour of the dynamical systems. Physically interesting integrable systems are rare, and consequently, it stirs up considerable excitement when one is discovered. Moreover, the Painlevé analysis as described in   has been contributing -for some time now-a great deal in this direction.
Most integrable problems characterizing the motion of a rigid body around a fixed point were collected in  , furthermore, the study by Mikhail P. Kharlamov et al. and A.V. Tsiganov et al., as described in   has been contributing for some time already in this direction and others types integrable problems describing the motion of a particle in the Euclidean plane were also collected in the Hietarinta study  nevertheless the problems correspond to the trapping of Ions   , and new integrable problems in this field have been added by many authors.
We examine here an integrable mechanical system which exhibits a great richness of behavior. The proposed system is the system of Paul trap, the Hamiltonian flows are generated by the Hamiltonian:
where is a constant, and , is known to be integrable in the following three cases  : it is also demonstrated that the ion dynamics in a Paul trap can be classified into four dynamical regimes  . This classification seems to be rather universal and shows up in the dynamics of periodically perturbed polar molecules, the hydrogen atom in strong magnetic fields  , and the periodically perturbed hydrogen atom.
The plan of the paper is as follows: In Section 2 we give a detailed description of the real phase space topology of the system (1) in the integrable case, for doing that, we separate the Hamiltonian system from two canonical transformations. This separability implies a description of the topology of the common-level sets of the first integrals (invariant level sets) however, in our study we consider the common-level sets of the first integrals:
where H and F are respectively the Hamiltonian and the second invariant of the system.
According to the classical Liouville theorem, for noncritical values of the first integrals h and f, the regular level sets of a completely integrable Hamiltonian system consists of tori. All generic bifurcations of these tori, corresponding to these critical values will be described by using Fomenko theorem  . We also give in Section 2.3 an explicit periodic solution for singular level sets of bifurcation studied above. Finally, for nonintegrable regimes we carry out a numerical analysis to bring out the order-chaos transition.
In the Hamiltonian (1) there is a singularity at , which necessitates an infinitesimally small step size for numerical integration of the corresponding equation of motion. So, one has to introduce appropriate coordinate transformation to remove this singularity. For this purpose, you can use two canonical transformations, the first is:
here and are the canonical momenta conjugate to the coordinates and respectively.
Then, Equation (1) Can be rewritten as for
Equation (4) is a three degrees of freedom Hamiltonian system in which is a cyclic variable, and so the corresponding canonically conjugate momenta is conserved, or
Then, Equation (4) can be rewritten as
the second canonical transformation is:
Then, Equation (5) can be rewritten as
and the Hamilton’s equations of motion of (7) start as
which is obtained by making .
2. Topological Analysis
In the integrable case, we recall the Hamilton-Jacobi equation corresponding to the system (8) that separates into coordinates defined by
It is easy to check that and can be expressed in terms of and characteristic polynomials in the following way:
where #Math_34# (10)
And F denotes the second integral of motion:
With the rescaled time variable:
Therefore, the differential equations satisfied by and are:
2.1. Topology of Regular Level Set
In order to give a complete description of the topology of , we find first the bifurcation diagram B in the -plane, i.e. the set of the critical values of the energy-momentum mapping
Definition. The bifurcation diagram of an integrable system is defined to be the region of possible motion depicted on the plane of first integrals  .
It turns out (like in the Hénon-Heils  , Gorjatchev-Tchaplygin top  , Fokker-Planck system   and Kolossoff potential   , the phase topology of a special case of Goryachev integrability in rigid body dynamics  ) that B is exactly the discriminant locus of the polynomial whose coefficients are functions of h and f.
The set consists of three connected components (as it is shown in Figure 1). Thus, in each connected component of the set the level set has the same topological type and this latter may be changed only if passes through B.
Figure 1. Bifurcation diagram B.
As θ is acyclical variable
Theorem. The set consists of three connected and nonintersecting with each other domains. The topological type of is a disjoint union of two-dimensional two-tori 2T, two-dimensional tori T and the empty set .
Proof. Consider the complexified system
consider also the elliptic curves
and the corresponding Riemann surfaces and of the same genus. We obtain the explicit solutions of the initial problem (13) by solving the Jacobi inversion problem  .
Define the natural projection
(where Ä is the symmetric product), and the complex conjugation on
Consider also the natural projection on the Riemann surface given in coordinates by .
It induces an involution on the Jacobi variety and hence on by the natural projection . By Equations (10) and (13) imply that this involution coincides with the complex conjugation (19) on . The upshot is that in order to describe it is enough to study the projection:
Definition. A connected component of the set of fixed points of on the curve and is called an oval.
To determine the ovals of and it suffices to study the real roots of the polynomial for different values of h and f as shown in Table 1. Using the Formulae (10), and the condition that we find exactly two admissible ovals whose projections on the and are given by and (see Table 2). The product of the admissible ovals in and the projection of such as,#Math_89#, gives :
1) is a two-dimensional two-tori 2T in domain 3.
2) is a two-dimensional tori T in domain 2.
3) is the empty set in domain 1.
2.2. Topology of Singular Level Sets
Suppose now that the constants h and f are changed in such a way that passes through the bifurcation diagram B. Then the topological type of may change and the bifurcation of Liouville tori takes place. In order to describe all generic bifurcations of Liouville tori, we use Fomenko’s theorem of bifurcation for Liouville tori. We can have in our case two types of bifurcation (see Figure 2).
To prove that, it suffices to look at the bifurcations of roots of the polynomial , the correspondence between bifurcation of roots and Liouville tori is shown in Figure 2.
Figure 2. Correspondence between bifurcations of roots of polynomial. and bifurcations of invariant Liouville tori.
Table 1. Topological type of and the real roots of the polynomials for .
Table 2. Admissible ovals and topological type of for .
2.3. Periodic Solutions
When the bifurcation of Liouville tori takes place, the level set becomes completely degenerate. Then we can have exceptional families of periodic solutions. It is seen from Table 3 that if is on the smooth curves (see Figure 1), contains a isolated circle S which is periodic solution. Consider now a fixed periodic solution belonging to the curve . The parameter takes values in the admissible interval and r is equal to the double root of the polynomial , (see Table 3).
Then we obtain from (9) and (10) the following parameterization of fixed periodic solution:
The expression of can put in the form:
Using the differential equations of Hamilton
It is easy to deduce solution
and the period associated is given by
3. Numerical Illustration
Using a surface of the section map, we give numerical illustrations of the topological analysis studied in Section 2.
For fixed values of energy h and f varies, the Liouville tori contained in the level set and change their topological type. The surfaces of section map shown in Figure 3 gives an illustration of the sequence of bifurcations
(a) (b) (c) (d) (e) (f) (g) (h) (i)
(j) (k) (l) (m) (n) (o) (p)
Figure 3. Surfaces of section map for different values of h, f and , : (a) Domain 2 (h = 3.016, f = 3.345) ; (b) Domain 2 (h = 3.016, f = 4.493) ; (c) Domain 2 (h = 3.016, f = 3.345, 4.493, 5.023, 5.156, 5.377) ; (d) Domain 2 (h = 3.016, f = 3.345, 4.493, 5.023, 5.156, 5.377) ; (e) Domain 2 (h = 3.016,, f = 3.345, 4.493, 5.023, 5.156, 5.377) ; (f) Domain 2 (h = 0.718, f = 0.828) ; (g) Domain 2 (h = 0.718, f = 0.828) ; (h) Domain 2 (h = 1.997, f = 4.077) ; (i) Domain 2 (h = 1.997, f = 4.052) ; (j) View 2D of domain 2 (λ = 1.16); (k) View 3D of domain 2 (λ = 1.16); (l) Domain 2 (λ = 1.35); (m) Domain 3 (h = 4.041, f = 3.017) ; (n) Domain 3 (h = 3.128, f = 3.017) ; (o) Domain 3 (h = 3.0755, f = 3.017) ; (p) View 3D of domain 3 (λ = 1.18).
Table 3. Topological type of for .
of Liouville tori and the order-chaos transition when one of the system parameters is varied. This map is constructed using a clever method introduced by Poincaré and extended by Hénon  .
The Figures 3(a)-(e) represent the sections for five values of the second invariant and . These values correspond to five points of domain 2 on the bifurcation diagram B where is a two-dimensional tori T. Moreover, the Figure 3(f) and Figure 3(g) show the sections for a value of the second invariant and . This value corresponds to a point in domain 2 on the bifurcation diagram B where is a two-dimensional T.
The Figure 3(m) show the sections for a value of the first invariant and . These values correspond to a point of Domain 3 on the bifurcation diagram B where is a two-dimensional two-tori 2T.
For critical values of a control parameter , we observe a fairly random distribution of points which correspond to a dramatic change in the Poincarésections indicating the order-chaos transition, as it is shown respectively on the Figure 3(j), Figure 3(k), Figure 3(l), and Figure 3(p).
In this study we have treated the classical dynamics of an integrable Hamiltonian system with two degrees of freedom. The system is characterized by a polynomial dependent on the invariants of the motion H and F. The different results obtained show the capacity of the method used to provide precise information on this Hamiltonian system. We have shown how this system can be converted by canonical transformations to easily exploitable Hamiltonians.
The very important question that we have studied is the topological analysis of the real invariant manifolds of the system. Fomenko’s theory on surgery and bifurcations of the Liouville tori has been combined with that of the algebraic structure to give a rigorous and detailed description of the topology of the invariant manifolds. For noncritical values of H and F, the variety contains torus or is empty.
In the same way we have shown how the periodic orbits can be found for singular values of first integrals, how the period of solutions is determined, and how explicit formulas can be established.
We have also highlighted numerically the topology of the invariant manifolds, the bifurcations of the Liouville tori and the order-chaos transition when the system control parameter varies.
 Ramani, A., Grammaticos, B. and Bountis, T. (1989) The Painleve Property and Singularity Analysis of Integrable and Non-Integrable Systems. Physics Reports, 180, 159-245.
 Yaremko, Y., Przybylska, M. and Maciejewski, A.J. (2015) Dynamics of a Relativistic Charge in the Penning Trap. Chaos: An Interdisciplinary Journal of Nonlinear Science, 25, 053102-1-16.
 Lanchares, V., Pascual, A.I., Palacián, J., Yanguas, P. and Salas, J.P. (2002) Perturbed Ion Traps: A Generalization of the Three-Dimensional Hénon-Heiles Problem. Chaos: An Interdisciplinary Journal of Nonlinear Science, 12, 87-99.
 Blümel, R., Chen, J.M., Diedrich, F., Peik, E., Quint, W., Schleich, W., Shen, Y.R. and Walther, H. (1988) Phase Transitions of Stored Laser-Cooled Ions. Nature, 334, 309-313.
 Ouazzani, A.T.H., Kharbach, J., Dekkaki, S. and Ouazzani-Jamil, M. (2000) Bifurcations Sets of the Motion of a Heavy Rigid Body around a Fixed Point in Goryatchev-Tchaplygin Case. IL Nuovo Cimento, 115, 1175-1193.
 Kharbach, J., Dekkaki, S., Ouazzani, A.T.H. and Ouazzani-Jamil, M. (2001) Topology and Bifurcations of the Invariant Level Sets of a Fokker-Planck Hamiltonian through Two Coupled Anisotropic Quartic Anharmonic Oscillators. Journal of Physics A: Mathematical and General, 34, 3437-3446.
 Kharbach, J., Dekkaki, S., Ouazzani, A.T.H. and Ouazzani-Jamil, M. (2003) Bifurcations of the Common Level Sets of Atomic Hydrogen in Van Der Waals Potential. International Journal of Bifurcation and Chaos, 13, 107-114.
 Gavrilov, L., Ouazzani-Jamil, M. and Caboz, R. (1993) Bifurcation Diagrams and Fomenko's Surgery on Liouville Tori of the Kolossoff Potential . Annales scientifiques de l'école normale supérieure, 36, 545-564.