Next Article in Journal
Higher Order Methods for Nonlinear Equations and Their Basins of Attraction
Previous Article in Journal
Solution of Differential Equations with Polynomial Coefficients with the Aid of an Analytic Continuation of Laplace Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Mathematics, University of Rhode Island, Kingston, RI 02881-0816, USA
2
Department of Mathematics, University of Sarajevo, 71000 Sarajevo, Bosnia and Herzegovina
*
Author to whom correspondence should be addressed.
Mathematics 2016, 4(1), 20; https://doi.org/10.3390/math4010020
Submission received: 6 February 2016 / Accepted: 14 March 2016 / Published: 18 March 2016

Abstract

:
By using the KAM(Kolmogorov-Arnold-Moser) theory and time reversal symmetries, we investigate the stability of the equilibrium solutions of the system:
x n + 1 = 1 y n , y n + 1 = β x n 1 + y n , n = 0 , 1 , 2 , ,
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 prove the existence of periodic points with arbitrarily large periods in every neighborhood of the unique positive equilibrium. We use invariants to find a Lyapunov function and Morse’s lemma to prove closedness of invariants. We also use the time reversal symmetry method to effectively find some feasible periods and the corresponding periodic orbits.

1. Introduction

The following rational system of difference equations:
x n + 1 = 1 y n y n + 1 = β x n 1 + y n , n = 0 , 1 ,
and the corresponding equation:
y n + 1 = β y n - 1 ( 1 + y n ) n = 0 , 1 , ,
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:
I ( x n , y n ) = β x n + y n + 1 x n + β y n + y n x n .
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:
y n + 1 = B y n y n - 1 + E y n - 1 + F b y n y n - 1 + e y n - 1 + f n = 0 , 1 ,
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:
y n + 1 = y n - 2 y n - 3 ( 1 + y n - 1 ) ( 1 + y n - 2 ) β + y n - 2 ( 1 + y n - 1 ) , 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.,
T x y = 1 y β x 1 + y .
The map (5) has the unique fixed point ( 1 / y ¯ , y ¯ ) in the positive quadrant, where
y ¯ 2 ( 1 + y ¯ ) = β .
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 ( y ¯ x ) , and v = ln ( y / 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
J T ( x , y ) = 0 - 1 y 2 β y + 1 - β x ( y + 1 ) 2
with
d e t J T ( x , y ) = β y 2 ( y + 1 ) .
We substitute u = ln ( y ¯ x ) , v = ln ( y / y ¯ ) and rewrite the map in ( u , v ) coordinates to obtain the map
u v - v ln β + u - ln e v y ¯ + 1 - 2 ln y ¯ .
The Jacobian matrix of this map is
J T ( u , v ) = 0 - 1 1 1 e v y ¯ + 1 - 1 ,
and so d e t J T ( u , v ) = 1 . ☐
A fixed point ( x ¯ , y ¯ ) is an elliptic point of an area preserving map if the eigenvalues of J T ( x ¯ , y ¯ ) 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 / y ¯ , y ¯ ) . In the logarithmic coordinates, the corresponding fixed point is ( 0 , 0 ) .
Proof
For the fixed points in ( x , y ) coordinates, solving 1 / y = x and β x / ( 1 + y ) = y yields the fixed point ( 1 / y ¯ , y ¯ ) where y ¯ is the unique positive solution of y ¯ 2 ( 1 + y ¯ ) = β . Evaluating the Jacobian matrix (6) of T at ( 1 / y ¯ , y ¯ ) gives
J T ( 1 / y ¯ , y ¯ ) = 0 - 1 y ¯ 2 β y ¯ + 1 - β y ¯ y ¯ + 1 2 .
By using β = y ¯ 3 + y ¯ 2 , we obtain that the eigenvalues of J T ( 1 / y ¯ , y ¯ ) are λ , λ ¯ where
λ = - y ¯ + i y ¯ + 2 3 y ¯ + 2 2 y ¯ + 1 .
Since | λ | = 1 , we have that ( 1 / y ¯ , y ¯ ) is an elliptic fixed point.
Under the logarithmic coordinate change ( x , y ) ( u , v ) , the fixed point ( 1 / y ¯ , y ¯ ) becomes ( 0 , 0 ) . Evaluating the Jacobian matrix (7) of T at ( 0 , 0 ) gives
J T ( 0 , 0 ) = 0 - 1 1 1 y ¯ + 1 - 1
with eigenvalues which are given by Equation (8). ☐
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.

2. 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
z λ z + g ( z , z ¯ ) ,
where λ is the eigenvalue of the elliptic fixed point, z = x + i y and z ¯ = x - i y 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
ζ λ ζ e i τ ( ζ ζ ¯ ) + h ( ζ , ζ ¯ )
in a neighborhood of the elliptic fixed point, where τ ( ζ ζ ¯ ) = τ 1 | ζ | 2 + + τ s | ζ | 2 s is a real polynomial, s = [ ( q - 2 ) / 2 ] 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 = c o n s t . 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
ζ λ ζ + c 1 ζ 2 ζ ¯ + O ( | ζ | 4 )
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 I m ( λ ) 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
T x 1 x 2 = R e ( λ ) - I m ( λ ) I m ( λ ) R e ( λ ) x 1 x 2 + g 1 ( x 1 , x 2 ) g 2 ( x 1 , x 2 ) .
One can now pass to the complex coordinates z = x 1 + i x 2 to obtain the complex form of the system
z λ z + ξ 20 z 2 + ξ 11 z z ¯ + ξ 02 z ¯ 2 + ξ 30 z 3 + ξ 21 z 2 z ¯ + ξ 12 z z ¯ 2 + ξ 03 z ¯ 3 + O ( | z | 4 ) .
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:
c 1 = ξ 20 ξ 11 ( λ ¯ + 2 λ - 3 ) ( λ 2 - λ ) ( λ ¯ - 1 ) + | ξ 11 | 2 1 - λ ¯ + 2 | ξ 02 | 2 λ 2 - λ ¯ + ξ 21 ,
where
ξ 20 = 1 8 ( g 1 ) x 1 x 1 - ( g 1 ) x 2 x 2 + 2 ( g 2 ) x 1 x 2 + i ( g 2 ) x 1 x 1 - ( g 2 ) x 2 x 2 - 2 ( g 1 ) x 1 x 2 ,
ξ 11 = 1 4 ( g 1 ) x 1 x 1 + ( g 1 ) x 2 x 2 + i ( g 2 ) x 1 x 1 + ( g 2 ) x 2 x 2 ,
ξ 02 = 1 8 ( g 1 ) x 1 x 1 - ( g 1 ) x 2 x 2 - 2 ( g 2 ) x 1 x 2 + i ( g 2 ) x 1 x 1 - ( g 2 ) x 2 x 2 + 2 ( g 1 ) x 1 x 2 ,
ξ 21 = 1 16 ( g 1 ) x 1 x 1 x 1 + ( g 1 ) x 1 x 2 x 2 + ( g 2 ) x 1 x 1 x 2 + ( g 2 ) x 2 x 2 x 2 ,
+ i 16 ( g 2 ) x 1 x 1 x 1 + ( g 2 ) x 1 x 2 x 2 - ( g 1 ) x 1 x 1 x 2 - ( g 1 ) x 2 x 2 x 2 .
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
F u v = - v ln β + u - ln e v y ¯ + 1 - 2 ln y ¯ .
Then, F has the unique elliptic fixed point ( 0 , 0 ) . The Jacobian matrix of F at ( u , v ) is given by
J F ( u , v ) = 0 - 1 1 1 e v y ¯ + 1 - 1 .
At ( 0 , 0 ) , J F ( u , v ) has the form
J 0 = J F ( 0 , 0 ) = 0 - 1 1 1 y ¯ + 1 - 1 .
The eigenvalues of (9) are λ and λ ¯ where
λ = - y ¯ + i y ¯ + 2 3 y ¯ + 2 2 y ¯ + 1 .
One can prove that
| λ | = 1
λ 2 = - y ¯ y ¯ + 4 + 2 2 y ¯ + 1 2 - i y ¯ + 2 3 y ¯ + 2 y ¯ y ¯ + 1 2 y ¯ + 2 λ 3 = y ¯ 2 y ¯ y ¯ + 3 + 3 2 y ¯ + 1 3 - i 2 y ¯ + 1 y ¯ + 2 3 y ¯ + 2 2 y ¯ + 1 3 λ 4 = y ¯ - y ¯ 3 + 8 y ¯ + 8 + 2 2 y ¯ + 1 4 + i y ¯ y ¯ + 2 3 y ¯ + 2 y ¯ y ¯ + 4 + 2 2 y ¯ + 1 4 ,
from which follows that λ k 1 for k = 1 , 2 , 3 , 4 and β > 0 .
Now, we have that
F u v = 0 - 1 1 1 y ¯ + 1 - 1 u v + f 1 ( δ , u , v ) f 2 ( δ , u , v ) ,
where
f 1 ( δ , u , v ) = 0 , f 2 ( δ , u , v ) = v y ¯ y ¯ + 1 - ln e v y ¯ + 1 - 2 ln y ¯ + ln β .
Then, the system ( u n + 1 , v n + 1 ) = F ( u n , v n ) is equivalent to
u n + 1 v n + 1 = 0 - 1 1 1 y ¯ + 1 - 1 u n v n + f 1 ( u n , v n ) f 2 ( u n , v n ) .
Let
u n v n = P u ˜ n v ˜ n ,
where
P = 1 D y ¯ 2 y ¯ + 1 - y ¯ + 2 3 y ¯ + 2 2 y ¯ + 1 1 0 , P - 1 = D 0 1 - 2 y ¯ + 1 y ¯ + 2 3 y ¯ + 2 y ¯ y ¯ + 2 3 y ¯ + 2
and
D = y ¯ + 2 3 y ¯ + 2 2 y ¯ + 1 .
Thus, the system ( u n + 1 , v n + 1 ) = F ( u n , v n ) is transformed into its Birkhoff normal form
u ˜ n + 1 v ˜ n + 1 = - y ¯ 2 y ¯ + 2 - y ¯ + 2 3 y ¯ + 2 2 y ¯ + 1 y ¯ + 2 3 y ¯ + 2 2 y ¯ + 2 - y ¯ 2 y ¯ + 2 u ˜ n v ˜ n + P - 1 H P u ˜ n v ˜ n ,
where
H u v : = f 1 ( u , v ) f 2 ( u , v ) .
Let
G u v = g 1 ( u , v ) g 2 ( u , v ) = P - 1 H P u v .
By a straightforward calculation, we obtain that
g 1 ( u , v ) = D u y ¯ D y ¯ + 1 - ln y ¯ e u D + 1 - 2 ln y ¯ + ln β g 2 ( u , v ) = y ¯ D y ¯ + 1 - ln y ¯ e u D + 1 - 2 ln y ¯ + ln β + u y ¯ y ¯ + 1 y ¯ + 2 3 y ¯ + 2 .
Another calculation gives
ξ 20 | u = v = 0 = - y ¯ y ¯ + 2 3 y ¯ + 2 + i y ¯ 8 D y ¯ + 1 2 y ¯ + 2 3 y ¯ + 2 ξ 11 | u = v = 0 = - y ¯ 1 + i y ¯ y ¯ + 2 3 y ¯ + 2 4 D y ¯ + 1 2 ξ 02 | u = v = 0 = - y ¯ y ¯ + 2 3 y ¯ + 2 + i y ¯ 8 D y ¯ + 1 2 y ¯ + 2 3 y ¯ + 2 ξ 21 | u = v = 0 = y ¯ - 1 y ¯ y ¯ + 2 3 y ¯ + 2 + i y ¯ 16 D y ¯ + 1 3 y ¯ + 2 3 y ¯ + 2 .
By using
ξ 20 ξ 11 = y ¯ 2 y ¯ + 2 3 y ¯ + 2 + i y ¯ 2 32 D y ¯ + 1 4 y ¯ + 2 3 y ¯ + 2 ξ 11 ξ 11 ¯ = y ¯ 2 4 D y ¯ + 1 2 3 y ¯ 2 + 8 y ¯ + 4 ξ 02 ξ 02 ¯ = y ¯ 2 16 D y ¯ + 1 2 3 y ¯ 2 + 8 y ¯ + 4 ,
a straightforward calculation yields
c 1 = ξ 20 ξ 11 ( λ ¯ + 2 λ - 3 ) ( λ 2 - λ ) ( λ ¯ - 1 ) + | ξ 11 | 2 1 - λ ¯ + 2 | ξ 02 | 2 λ 2 - λ ¯ + ξ 21 = y ¯ - 1 y ¯ y ¯ + 1 3 y ¯ + 2 y ¯ + 2 3 y ¯ + 2 4 + y ¯ 3 y ¯ - i y ¯ + 2 3 y ¯ + 2 + 8 .
It can be proved that
τ 1 = - i λ ¯ c 1 = - y ¯ - 1 y ¯ 2 y ¯ + 2 3 y ¯ + 2 2 ,
which implies that τ 1 0 for β 2 since y ¯ 2 ( 1 + y ¯ ) = β .  ☐
 ☐
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 ¯ , y ¯ ) be a non-degenerate elliptic fixed point. If β 2 , then there exist periodic points with arbitrarily large periods in every neighborhood of ( x ¯ , y ¯ ) . In addition, ( x ¯ , y ¯ ) is a stable fixed point.

3. 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).
Proof
The system (1) possesses an invariant given by Equation (3). The function I ( x , y ) associated with Equation (3) has partial derivatives
I x = β - 1 x 2 - y x 2 , I y = 1 - β y 2 + 1 x .
The unique equilibrium of Equation (1) satisfies that x ¯ = 1 y ¯ and y ¯ 2 ( 1 + y ¯ ) = β . Hence x ¯ is the unique positive solution of the equation β x 3 - x - 1 = 0 . Equation (10) implies that any critical point ( x , y ) of Equation (3) satisfies y = β x 2 - 1 and ( x + 1 ) y 2 - β x = 0 . Substitution yields
β 2 x 5 + β 2 x 4 - 2 β x 3 - 2 β x 2 + ( 1 - β ) x + 1 = 0 .
Equation (11) can be rewritten as
( β x 3 - x - 1 ) ( β x 2 + β x - 1 ) = 0 ,
which has x ¯ as a solution. Ⅰ
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 ¯ , y ¯ ) is stable.
Proof
The Hessian matrix associated with I ( x , y ) is
H ( x , y ) = 2 x 3 + 2 y x 3 - 1 x 2 - 1 x 2 2 β y 3
with determinant
det ( H ( x , y ) ) = 4 β ( 1 + y ) x 3 y 3 - 1 x 4 .
For ( x , y ) a critical point of I ( x , y ) ,
det ( H ( x , y ) ) = 4 β 2 x y 3 - 1 x 4 .
For ( x ¯ , y ¯ ) the unique equilibrium of Equation (1), we can further reduce Equation (12) to
det ( H ( x ¯ , y ¯ ) ) = 4 β 2 y ¯ y ¯ 3 - y ¯ 4 = 4 β 2 - y ¯ 6 y ¯ 2 .
Note that since the equation β x 3 - x - 1 = 0 has x ¯ as its unique positive solution, the equation β / y 3 - 1 / y - 1 = β - y 2 - y 3 y 3 = 0 has y ¯ as its unique positive solution. Let us define f ( y ) = y 3 + y 2 - β . We observe that f ( 0 ) = - β < 0 and f ( 4 β 2 6 ) = 2 β + 4 β 2 3 - β > 0 , which guarantees that 0 < y ¯ < 4 β 2 6 . Now the Morse’s lemma [18] guarantees the result provided det ( H ( x ¯ , y ¯ ) ) > 0 . However, 4 β 2 - y ¯ 6 y ¯ 2 > 0 if and only if y ¯ < 4 β 2 6 , 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 ¯ , y ¯ )
is a Lyapunov function, and so the equilibrium point ( x ¯ , y ¯ ) 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.

4. 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 = i d , 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 , f o r n e v e n ; T ( n ± 1 ) / 2 ( x , y ) S 1 , 0 , f o r n o d d .
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 2 j + i = T j ( S i ) , S 2 j - i = I j ( S i ) , for all i , j .
The inverse of the map (5) is the map T - 1 ( x , y ) = y ( 1 + 1 / x ) β , 1 x . The involution R ( x , y ) = 1 y , 1 x is a reversor for (3). Indeed,
( R T R ) ( x , y ) = ( R T ) 1 y , 1 x = R x , β / y 1 + 1 / x = y 1 + 1 / x β , 1 x = T - 1 ( x , y ) .
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 ) : x y = 1 } , S 1 = { ( x , y ) : β x = y 2 ( x + 1 ) } .
Periodic orbits on the symmetry line S 0 with even period n are searched for by starting with points ( x 0 , 1 / x 0 ) S 0 and imposing that ( x n / 2 , y n / 2 ) S 0 , where ( x n / 2 , y n / 2 ) = T n / 2 ( x 0 , 1 / x 0 ) . This reduces to a one-dimensional root finding for the equation x n / 2 y n / 2 = 1 , where the unknown is x 0 . Furthermore, periodic orbits on S 0 with odd period n are obtained by solving for x 0 the equation β x ( n + 1 ) / 2 = y ( n + 1 ) / 2 2 ( 1 + x ( n + 1 ) / 2 ) , where ( x ( n + 1 ) / 2 , y ( n + 1 ) / 2 ) = T ( n + 1 ) / 2 ( x 0 , 1 / x 0 ) .
For example, for β = 1 . 8 , in Figure 3, we have an intersection between the symmetry lines S 0 and S 22 = T 11 ( S 0 ) , S 4 = T 2 ( S 0 ) and S 22 = T 11 ( S 0 ) , and S 1 and S 23 = T 11 ( S 1 ) of the map T. The intersection points of this lines correspond to the periodic orbits of period 22 , 18 and 22, respectively.
See Figure 3 for some examples of the periodic orbits of periods 18 and 22.

5. 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:
x n + 1 = 1 y n , y n + 1 = β x n 1 + y n , 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.

Author Contributions

All three authors have written this paper and the final form of this paper is approved by all three authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Amleh, A.M.; Camouzis, E.; Ladas, G. On the Dynamics of a Rational Difference Equation, Part I. Int. J. Differ. Equ. 2008, 3, 1–35. [Google Scholar]
  2. Drymmonis, E.; Camouzis, E.; Ladas, G.; Tikjha, G.W. Patterns of boundedness of the rational system x n + 1 = α 1 A 1 + B 1 x n + C 1 y n and y n + 1 = α 2 + α 2 x n + γ 2 y n A 2 + B 2 x n + C 2 y n . J. Differ. Equ. Appl. 2012, 18, 89–110. [Google Scholar]
  3. Feuer, J.; Janowski, E.; Ladas, G. Invariants for some rational recursive sequences with periodic coefficients. J. Differ. Equ. Appl. 1996, 2, 167–174. [Google Scholar] [CrossRef]
  4. Grove, E.A.; Janowski, E.J.; Kent, C.M.; Ladas, G. On the rational recursive sequence x n + 1 = α x n + β ( γ x n + δ ) x n 1 . Comm. Appl. Nonlinear Anal. 1994, 1, 61–72. [Google Scholar]
  5. Kulenović, M.R.S. Invariants and related Liapunov functions for difference equations. Appl. Math. Lett. 2000, 13, 1–8. [Google Scholar] [CrossRef]
  6. Kulenović, M.R.S.; Merino, O. Discrete Dynamical Systems and Difference Equations with Mathematica; Chapman and Hall/CRC: Boca Raton, FL, USA; London, UK, 2002. [Google Scholar]
  7. Chan, D.M.; Kent, C.M.; Ortiz-Robinson, N.L. Convergence results on a second-order rational difference equation with quadratic terms. Adv. Differ. Equ. 2009. [Google Scholar] [CrossRef]
  8. Dehghan, M.; Kent, C.M.; Mazrooei-Sebdani, R.; Ortiz, N.L.; Sedaghat, H. Monotone and oscillatory solutions of a rational difference equation containing quadratic terms. J. Differ. Equ. Appl. 2008, 14, 1045–1058. [Google Scholar] [CrossRef]
  9. Dehghan, M.; Mazrooei-Sebdani, R.; Sedaghat, H. Global behaviour of the Riccati difference equation of order two. J. Differ. Equ. Appl. 2011, 17, 467–477. [Google Scholar] [CrossRef]
  10. Kulenović, M.R.S.; Moranjkić, S.; Nurkanović, Z. Naimark-Sacker bifurcation of second order rational difference equation with quadratic terms. J. Nonlinear Sci. Appl. 2016, in press. [Google Scholar]
  11. Cima, A.; Gasull, A.; Manosa, V. Non-integrability of measure preserving maps via Lie symmetries. J. Differ. Equ. 2015, 259, 5115–5136. [Google Scholar] [CrossRef] [Green Version]
  12. Kulenović, M.R.S.; Nurkanović, Z.; Pilav, E. Birkhoff Normal Forms and KAM theory for Gumowski-Mira Equation. Sci. World J. Math. Anal. 2014. [Google Scholar] [CrossRef] [PubMed]
  13. Gidea, M.; Meiss, J.D.; Ugarcovici, I.; Weiss, H. Applications of KAM Theory to Population Dynamics. J. Biol. Dyn. 2011, 5, 44–63. [Google Scholar] [CrossRef] [PubMed]
  14. Kocic, V.L.; Ladas, G.; Tzanetopoulos, G.; Thomas, E. On the stability of Lyness’ equation. Dynam. Contin. Discrete Impuls. Syst. 1995, 1, 245–254. [Google Scholar]
  15. Kulenović, M.R.S.; Nurkanović, Z. Stability of Lyness’ Equation with Period-Two Coefficient via KAM Theory. J. Concr. Appl. Math. 2008, 6, 229–245. [Google Scholar]
  16. Ladas, G.; Tzanetopoulos, G.; Tovbis, A. On May’s host parasitoid model. J. Differ. Equ. Appl. 1996, 2, 195–204. [Google Scholar] [CrossRef]
  17. Hrustić, S.J.; Kulenović, M.R.S.; Nurkanović, Z.; Pilav, E. Birkhoff Normal Forms, KAM theory and Symmetries for Certain Second Order Rational Difference Equation with Quadratic Term. Int. J. Differ. Equ. 2015, 10, 181–199. [Google Scholar]
  18. Hale, J.K.; Kocak, H. Dynamics and Bifurcation; Springer-Verlag: New York, NY, USA, 1991. [Google Scholar]
  19. Tabor, M. Chaos and Integrability in Nonlinear Dynamics: An introduction; A Wiley-Interscience Publication; John Wiley & Sons, Inc.: New York, NY, USA, 1989. [Google Scholar]
  20. Del-Castillo-Negrete, D.; Greene, J.M.; Morrison, E.J. Area preserving nontwist maps: Periodic orbits and transition to chaos. Phys. D 1996, 91, 1–23. [Google Scholar] [CrossRef]
  21. Wan, Y.H. Computation of the stability condition for the Hopf bifurcation of diffeomorphisms on 2 . SIAM J. Appl. Math. 1978, 34, 167–175. [Google Scholar] [CrossRef]
  22. Moeckel, R. Generic bifurcations of the twist coefficient. Ergodic Theory Dyn. Syst. 1990, 10, 185–195. [Google Scholar] [CrossRef]
  23. MacKay, R.S. Renormalization in Area-Preserving Maps; World Scientific: River Edge, NJ, USA, 1993. [Google Scholar]
  24. Siegel, C.; Moser, J. Lectures on Celestial Mechanics; Springer-Varlag: New York, NY, USA, 1971. [Google Scholar]
Figure 1. Some orbits of the map T for (a) a = 0 . 5 and (b) a = 1 . 5 . The plots are generated by Dynamica 3 [6].
Figure 1. Some orbits of the map T for (a) a = 0 . 5 and (b) a = 1 . 5 . The plots are generated by Dynamica 3 [6].
Mathematics 04 00020 g001
Figure 2. A bifurcation diagram in ( β - x )-plane. The plots are generated by Dynamica 3 [6].
Figure 2. A bifurcation diagram in ( β - x )-plane. The plots are generated by Dynamica 3 [6].
Mathematics 04 00020 g002
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).
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).
Mathematics 04 00020 g003

Share and Cite

MDPI and ACS Style

Denette, E.; Kulenović, M.R.S.; Pilav, E. Birkhoff Normal Forms, KAM Theory and Time Reversal Symmetry for Certain Rational Map. Mathematics 2016, 4, 20. https://doi.org/10.3390/math4010020

AMA Style

Denette E, Kulenović MRS, Pilav E. Birkhoff Normal Forms, KAM Theory and Time Reversal Symmetry for Certain Rational Map. Mathematics. 2016; 4(1):20. https://doi.org/10.3390/math4010020

Chicago/Turabian Style

Denette, Erin, Mustafa R. S. Kulenović, and Esmir Pilav. 2016. "Birkhoff Normal Forms, KAM Theory and Time Reversal Symmetry for Certain Rational Map" Mathematics 4, no. 1: 20. https://doi.org/10.3390/math4010020

APA Style

Denette, E., Kulenović, M. R. S., & Pilav, E. (2016). Birkhoff Normal Forms, KAM Theory and Time Reversal Symmetry for Certain Rational Map. Mathematics, 4(1), 20. https://doi.org/10.3390/math4010020

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop