Birkhoff Normal Forms, KAM Theory and Time Reversal Symmetry Birkhoff Normal Forms, KAM Theory and Time Reversal Symmetry for Certain Rational Map for Certain Rational Map

: By using the KAM(Kolmogorov-Arnold-Moser) theory and time reversal symmetries, we investigate the stability of the equilibrium solutions of the system:


Introduction
The following rational system of difference equations: and the corresponding equation: where the parameter β > 0 and initial conditions x 0 , y 0 are positive numbers were considered in [1] and [2].The authors established the boundedness of all solutions of system (1) by using the invariant: Equation ( 2) and its invariant (3), where x n = 1/y n−1 were obtained in [3,4] and the stability of the equilibrium by means of Lyapunov function generated by invariant (3) was derived in [5,6], (pp.247-250).Equation ( 2) is also a special case of equation: By n y n−1 + Ey n−1 + F by n y n−1 + ey n−1 + f n = 0, 1, . . .(4) with all nonnegative coefficients and initial conditions.Equation ( 4) is a rational difference equation with quadratic terms which is a subject of recent research, see [7][8][9][10].
In this paper, we will show that the corresponding map can be transformed into an area preserving map for which we will find the Birkhoff Normal form, and, using it, we will apply the KAM theorem to prove the stability of the unique positive equilibrium and the existence of periodic points with an arbitrarily large period in every neighborhood of the unique positive equilibrium.In addition, we prove that the corresponding map is conjugate to its inverse map through the involution map.Then, we will use this conjugacy to find some feasible periods of this map.The KAM theory will be enough to prove the stability of the equilibrium for β = 2, and then we use the invariant (3) and Morse's lemma to prove the stability in the remaining case β = 2, see [5].In addition, Morse's lemma implies that all invariants are locally simple closed curves.A very recent paper [11] gives some effective tests for difference equation to have a continuous invariant.The method of invariants for the construction of a Lyapunov function and proving stability of the equilibrium points was used successfully in [5,6,12], and the KAM theory was used for the same objective in [12][13][14][15][16].The class of difference equations which admit an invariant is not a large class even in the case of rational difference equations, see [11].In the case when a difference equation's corresponding map is area preserving and does not possess an invariant, the only tool left seems to be KAM theory, see [17] for such an example.Furthermore, the corresponding Equation ( 4) can be embedded by iteration into a fourth order difference equation: , n = 0, 1, . . ., which is increasing in all its arguments and yet exhibits the chaos.
Let T be the map associated to the system (1), i.e., The map (5) has the unique fixed point (1/ ȳ, ȳ) in the positive quadrant, where An invertible map T : R 2 → R 2 is area preserving if the area of T(A) equals the area of A for all measurable subsets A [6,18,19].As is known, a differentiable map T is area preserving if the determinant of its Jacobian matrix is equal ±1, that is det J T = ±1 at every point of domain of T, see [18,19].We claim that in logarithmic coordinates (u, v) where u = ln ( ȳx), and v = ln (y/ ȳ), the map (5) is area preserving.Lemma 1.The map (5) is an area preserving map in the logarithmic coordinates.
Proof.The Jacobian matrix of the map T is .
We substitute u = ln ( ȳx), v = ln (y/ ȳ) and rewrite the map in (u, v) coordinates to obtain the map The Jacobian matrix of this map is and so detJ T (u, v) = 1.
A fixed point ( x, ȳ) is an elliptic point of an area preserving map if the eigenvalues of J T ( x, ȳ) form a purely imaginary, complex conjugate pair λ, λ, see [6,18].Lemma 2. The map T in (x, y) coordinates has an elliptic fixed point (1/ ȳ, ȳ).In the logarithmic coordinates, the corresponding fixed point is (0, 0).
The rest of the paper is organized into three sections.The second section contains a derivation of the Birkhoff normal form for map T and an application of the KAM theory, which proves stability of the equilibrium and the existence of an infinite number of periodic solutions for β = 2.The third section makes use of the invariant (3) in proving stability for β = 2 and the construction of a Lyapunov function.The fourth section uses the symmetries for the map T showing that this map is conjugate to its inverse through an involution.Then, we use time reversal symmetry method [13,20] based on the symmetries to effectively find some feasible periods and corresponding orbits of the map T.

The KAM Theory and Birkhoff Normal Form
The KAM Theorem asserts that, in any sufficiently small neighborhood of a non degenerate elliptic fixed point of a smooth area-preserving map, there exists many invariant closed curves.We explain this theorem in some detail.Consider a smooth, area-preserving map (x, y) → T(x, y) of the plane that has (0, 0) as an elliptic fixed point.After a linear transformation, one can represent the map in the form where λ is the eigenvalue of the elliptic fixed point, z = x + iy and z = x − iy are complex variables, and g vanishes with its derivative at z = 0. Assume that the eigenvalue λ of the elliptic fixed point satisfies the non-resonance condition λ k = 1 for k = 1, . . ., q, for some q ≥ 4.Then, Birkhoff showed that there exists new, canonical complex coordinates (ζ, ζ) relative to which the mapping takes the normal form and h vanishes with its derivatives up to order q − 1.The numbers τ 1 , . . ., τ s are called twist coefficients.Consider an invariant annulus < |ζ| < 2 in a neighborhood of an elliptic fixed point, for , a very small positive number.Note that, if we neglect the remainder h, the normal form approximation ζ → λζe iτ(ζ ζ) leaves invariant all circles |ζ| 2 = const.The motion restricted to each of these circles is a rotation by some angle.In addition, please note that, if at least one of the twist coefficients τ j is nonzero, the angle of rotation will vary from circle to circle.A radial line through the fixed point will undergo twisting under the map.The KAM theorem (Moser's twist theorem) says that, under the addition of the remainder term, most of these invariant circles will survive as invariant closed curves under the full map.Theorem 3. Assume that τ(ζ ζ) is not identically zero and is sufficiently small, then the map T has a set of invariant closed curves of positive Lebesque measure close to the original invariant circles.Moreover, the relative measure of the set of surviving invariant curves approaches full measure as approaches 0. The surviving invariant closed curves are filled with dense irrational orbits.
The KAM theorem requires that the elliptic fixed point be non-resonant and non degenerate.Note that for q = 4 the non-resonance condition λ k = 1 requires that λ = ±1 or ±i.The above normal form yields the approximation with c 1 = iλτ 1 and τ 1 being the first twist coefficient.We will call an elliptic fixed point non-degenerate if τ 1 = 0.
Consider a general map T that has a fixed point at the origin with complex eigenvalues λ and λ satisfying |λ| = 1 and Im(λ) = 0.By putting the linear part of such a map into Jordan Normal form, we may assume that T has the following form near the origin One can now pass to the complex coordinates z = x 1 + ix 2 to obtain the complex form of the system The coefficient c 1 can be computed directly using the formula below derived by Wan in the context of Hopf bifurcation theory [21].In [22], it is shown that, when one uses area-preserving coordinate changes, Wan's formula yields the twist coefficient τ 1 that is used to verify the non-degeneracy condition necessary to apply the KAM theorem.We use the formula: where Theorem 4. The elliptic fixed point (0, 0), in the (u, v) coordinates, is non-degenerate for β = 2 and non-resonant for β > 0.
Proof.Let F be the function defined by Then, F has the unique elliptic fixed point (0, 0).The Jacobian matrix of F at (u, v) is given by At (0, 0), J F (u, v) has the form The eigenvalues of ( 9) are λ and λ where One can prove that |λ| = 1 , from which follows that λ k = 1 for k = 1, 2, 3, 4 and β > 0. Now, we have that where where where By a straightforward calculation, we obtain that .
Another calculation gives . By using , a straightforward calculation yields It can be proved that The following result is a consequence of Moser's twist map theorem [13,19,23,24].
Theorem 5. Let T be a map (5) associated to the system (1), and ( x, ȳ) be a non-degenerate elliptic fixed point.If β = 2, then there exist periodic points with arbitrarily large periods in every neighborhood of ( x, ȳ).
In addition, ( x, ȳ) is a stable fixed point.

Invariant
In this section, we use the invariant to find a Lyapunov function and prove stability of the equilibrium for all values of parameter β > 0, see [5,6] for similar results.Lemma 6.The unique equilibrium ( 1 y , y) of Equation ( 1) is a critical point of the invariant (3).
Lemma 7. The graph of the function I(x, y) associated with Equation (3) is a simple closed curve in a neighborhood of the equilibrium of Equation (1).The equilibrium point ( x, ȳ) is stable.
Proof.The Hessian matrix associated with I(x, y) is For (x, y) a critical point of I(x, y), For (x, y) the unique equilibrium of Equation (1), we can further reduce Equation (12) to Note that since the equation βx 3 − x − 1 = 0 has x as its unique positive solution, the equation = 0 has ȳ as its unique positive solution.Let us define f (y) = y 3 + y 2 − β.We observe that f (0) = −β < 0 and f ( 6 4β 2 ) = 2β + 3 4β 2 − β > 0, which guarantees that 0 < y < 6 4β 2 .Now the Morse's lemma [18] guarantees the result provided det(H(x, y)) > 0. However, > 0 if and only if y < 6 4β 2 , which is indeed the case.In view of the Morse's lemma [18], the level sets of the function I(x, y) are diffeomorphic to circles in the neighborhood of (x, y).In addition, the function V(x, y) = I(x, y) − I( x, ȳ) is a Lyapunov function, and so the equilibrium point ( x, ȳ) is stable, see [5].
See Figure 1 for the family of invariant curves around the equilibrium.See Figure 2 for the bifurcation diagrams which indicate the appearance of chaos. .

Symmetries
In the study of area-preserving maps, symmetries play an important role since they yield special dynamic behavior.A transformation R of the plane is said to be a time reversal symmetry for T if R −1 • T • R = T −1 , meaning that applying the transformation R to the map T is equivalent to iterating the map backwards in time, see [13,20].If the time reversal symmetry R is an involution, i.e., R 2 = id, then the time reversal symmetry condition is equivalent to R • T • R = T −1 , and T can be written as the composition of two involutions T = I 1 • I 0 , with I 0 = R and I 1 = T • R. Note that if I 0 = R is a reversor, then so is I 1 = T • R. In addition, the jth involution, defined as I j := T j • R, is also a reversor.
The invariant sets of the involution maps, S 0,1 = {(x, y)|I 0,1 (x, y) = (x, y)} are one-dimensional sets called the symmetry lines of the map.Once the sets S 0,1 are known, the search for periodic orbits can be reduced to a one-dimensional root finding problem using the following result, see [13,20]: Theorem 8.If (x, y) ∈ S 0,1 then T n (x, y) = (x, y) if and only if T n/2 (x, y) ∈ S 0,1 , for n even; T (n±1)/2 (x, y) ∈ S 1,0 , for n odd.
That is, according to this result, periodic orbits can be found by searching in the one-dimensional sets S 0,1 , rather than in the whole domain.Periodic orbits of different orders can then be found at the intersection of the symmetry lines S j j = 1, 2, ... associated to the jth involution; for example, if (x, y) ∈ S j ∩ S k , then T j−k (x, y) = (x, y).In addition, the symmetry lines are related to each other by the following relations: S 2j+i = T j (S i ), S 2j−i = I j (S i ), for all i, j.The inverse of the map (5) is the map T x is a reversor for (3).Indeed, Thus, T = I 1 • I 0 where I 0 (x, y) = R(x, y) and I 1 (x, y) = T • R = x, β y(1+1/x) .The symmetry lines corresponding to I 0 and I 1 are S 0 = {(x, y) : xy = 1}, S 1 = {(x, y) : βx = y 2 (x + 1)}.
See Figure 3 for some examples of the periodic orbits of periods 18 and 22.

Conclusions
By using the KAM (Kolmogorov-Arnold-Moser) theory, invariants and corresponding Lyapunov function and time reversal symmetries, we proved the stability of the equilibrium solution of the system: , n = 0, 1, . . ., where the parameter β > 0, and initial conditions x 0 and y 0 are positive numbers.We obtain the Birkhoff normal form for this system and used them to prove the existence of periodic points with arbitrarily large periods in every neighborhood of the unique positive equilibrium.

Figure 3 .
Figure 3. (a) The first twelve iterations of symmetry line S 0 of the map T for β = 0.18; (b) the first eleven iterations of symmetry line S 1 of the map T for β = 0.18; (c) the periodic orbits of period 22 (red) and 18 (blue).