1. Introduction
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 [1] [2] 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 [3] , furthermore, the study by Mikhail P. Kharlamov et al. and A.V. Tsiganov et al., as described in [4] [5] 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 [6] nevertheless the problems correspond to the trapping of Ions [7] [8] , 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:
(1)
where
is a constant, and
, is known to be integrable in the following three cases [9] :
it is also demonstrated that the ion dynamics in a Paul trap can be classified into four dynamical regimes [10] . 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 [11] , 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:
(2)
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 [12] . 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:
(3)
here
and
are the canonical momenta conjugate to the coordinates
and
respectively.
Then, Equation (1) Can be rewritten as for
(4)
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
(5)
the second canonical transformation is:
(6)
Then, Equation (5) can be rewritten as
(7)
and the Hamilton’s equations of motion of (7) start as
(8)
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
(9)
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:
(11)
(12)
With the rescaled time variable:
Therefore, the differential equations satisfied by
and
are:
(13)
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
(14)
Definition. The bifurcation diagram of an integrable system is defined to be the region of possible motion depicted on the plane of first integrals
[13] .
It turns out (like in the Hénon-Heils [14] , Gorjatchev-Tchaplygin top [15] , Fokker-Planck system [16] [17] and Kolossoff potential [18] [19] , the phase topology of a special case of Goryachev integrability in rigid body dynamics [20] ) that B is exactly the discriminant locus of the polynomial
whose coefficients are functions of h and f.
(15)
where
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.
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
(16)
consider also the elliptic curves
and
(17)
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 [21] .
Define the natural projection
(18)
(where Ä is the symmetric product), and the complex conjugation on
(19)
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:
(20)
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:
(21)
Knowing that
(22)
The expression of
can put in the form:
(23)
Using the differential equations of Hamilton
(24)
It is easy to deduce solution
(25)
and the period associated is given by
(26)
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 [22] .
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 fixed points in Figure 3(h) and Figure 3(i) show the sections representing the periodic solution where
is an isolated circle S for
and
on the curve
of B.
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.
The Figure 3(n) and Figure 3(o) show the sections corresponding to the bifurcation on the curve
of B for
and
where
is a
.
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).
4. Conclusions
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.