BIFURCATIONS AND ATTRACTORS IN BOGDANOV MAP
Ilhem Djellit, Ibtissem Boukemara
Abstract: In this paper, we study bifurcation space and the phase plane of the Bogdanov map. Specific bifurcation structures can be observed in parameter space, related for instance to embedded boxes structure and configurations of bifurcation curves of periodic points near the cusp. This model is a diffeomorphism. The dynamics is extremely rich, involving periodicity, quasiperiodicity and chaos. The method of the study is a numerical iteration to an attractor in which the guesses are inspired by the theory. Bifurcation diagrams obtained in different parameter planes are given and a sketch showing the cusp bifurcations, for the versal unfolding of Bogdanov map, is given. Phase plane is also studied, different attractors are shown, their evolution giving rise to chaotic attractors is explained. Basins of attraction are considered and fuzzy boundaries of basins are put in evidence. The study of such kind of diffeomorphisms can give an interesting contribution to nonlinear systems. Key words: bifurcation, basin, Bogdanov map, invertible, chaos. Classification: 34D4537D1034C37
Introduction We study a diffeomorphism. The dynamics involves various transitions by bifurcations. A planar quadratic map which is conjugate to the Hénon map in its conservative limit (see [H1]). Perturbation theory here reveals periodicity and quasiperiodicity, parametrically organized by embedded boxes structure. Outside these parameter regimes the perturbation becomes larger and bifurcations complicate the dynamics, often involving chaos. On the one hand, this concerns local bifurcations, e.g. cascades of flips and cusp bifurcations. On the other hand, we have to deal with homoclinic bifurcations. It also undergoes a Hopf bifurcation which creates a stable closed curve growing out, at its destruction appears a saddle connection. Here the dynamical features are explored by numerical methods. Also, we often iterate to a periodic attractor. In several cases we find interesting dynamical objects predicted by the theory and global phenomena in the parameter plane. This kind of scanning has been made for giving a first idea about bifurcation organization. That is, the dynamic behaviors may be, as expected, regular (equilibrium , stable cycles and stable invariant curves) or chaotic. Object of the present study is to investigate the main routes to complex behaviors in attracting sets and their basins of attraction associated with the dynamics of invertible nonlinear systems which are still of great interest. The Bogdanov map provides a good approximation to the dynamics of the Poincaré map of periodically forced oscillators, first considered by Bogdanov (see [B1]), Takens and Arnold [A3] in their study of the double zero eigenvalue singularity, they proved the existence of a codimensiontwo fixed point at the origin, called a BogdanovTakens cusp with a nonzero Jordan canonical form but the best contribution come from Arrowsmith and al. [A1, A2]. We emphasize that our aim is to study simply the map without emulate the Bogdanov vector field. This diffeomorphism can be considered as a model giving rise to a complete set of bifurcations with a fractal structure, and with birth of cyclic invariant curves (impossible to obtain with Hénon map). From this point of view it gives the fundamental mechanisms occurring in the nonautonomous differential equation with quasi periodic solutions, or more complex solutions. The considered diffeomorphism is the following termed by Bogdanov map: T: x_{n}_{+1} = x_{n}+y_{n}_{+1, }y_{n}_{+1} = _{ }y_{n}+ey_{n}+hx_{n}(x_{n}1)+mx_{n}.y_{n} (1) Where x, y are real variables, e, m are real parameters . The parameter h is the step length in the discretization, it takes small values. T has not a constant jacobian determinant J = 1+e+mx, then it results from this situation that a cycle can undergo the Hopf bifurcation and for e = m = 0, we observe the hamiltonian nature and we have a globally area preserving map. Eq. (1) has only two fixed points (0,0) and (1,0). These both fixed points can undergo identical sets of bifurcations in parameter space, this is due to the symmetry in the map (see Fig.1 (a,b,c)) under the transformation ([`x] ,[`y]) = (1x,y) which switches the fixed points. We study the bifurcation diagram which consists to show
bifurcation curves in the parameter (e,h),(e,m)
planes together with representative phase portraits. The Figures 1 (a,b,c)
give the parameter value for which at least one fixed point is attractive
(blue domain corresponding to the value 1). More generally, the Figures
1 (a,b,c) give the regions of parameter ( e,h)
or (e,m) plane for
which at least a periodic orbit of order k exists (k = 1,2,...,14).
Each part coloured corresponds to the existence of a periodic orbit of
order k (period) of which is given by the upper coloured squares.
The black regions (k = 15) corresponds to the existence of bounded
iterated sequences.
Figure (1.a): bifurcation
structure in (
eh) plane, m
= 2.
Figure (1.b): bifurcation structure in ( e,h) plane, m = 2. These figures are typical of maps with dominating quadratic
nonlinear terms, analogous bifurcation diagrams have already been obtained
for systems with non linearity given by hyperbolic tangency [M1,
M2,
T]. We can recognize on the diagram period doubling bifurcation,
saddle area and spring area [M1]. The bifurcation structure
is a quadratic "boxwithinabox" type (see [M3,
F2]),
as is well known infinitely many periodic are opened by fold bifurcations
and are closed by homoclinic bifurcations by the intriguing "boxwithinabox"
bifurcation structure. There is a fold bifurcation on the line h
= 0, a flip bifurcation on the line
e = h/22,
and a Hopf bifurcation on the line e = 0 for
0 < h < 4.
The eigenvalues of the fixed point at the origin depend on the parameters
(e,h), in the interior of the curve e^{2}+h^{2}2eh4h
= 0 we have complex eigenvalues.
Figure (1.c): bifurcation structure in (e,m) plane, h = 1.44.
Definitions and some properties This section essentially concerns a twodimensional map T , considered here as a discrete dynamical system in R^{2} (x_{n+1}, y_{n+1}) = T (x_{n}, y_{n}) = (f (x_{n}, y_{n}, l), g (x_{n}, y_{n}, l)) where f (x, y, l) and g (x, y, l) are continuous and differentiable functions with respect to real variables x, y and continuous with respect to the real parameter l. It is recalled that a closed and invariant set A is called an attracting set if some neighborhood U of A may exist such that T(U) ÌU, and T^{n}(x, y) ® A as n®¥, "(x, y) Î U. An attracting set A may contain one or several attractors coexisting with sets of repulsive points ( strange repulsors) to either chaotic transients towards these attractors, or fuzzy boundaries of their basin [M4]. The set D = È_{n³ 0}T ^{n}(U) is called the total basin (or influence domain) of A. That is, D is the open set of points x whose forward trajectories (set of images of x) converge towards A. D is invariant under backward iteration T^{1} of T, but not necessarily invariant by T: T ^{1}(D) = D, T (D) ÍD. The relations above hold also for the closure D of D. The boundary of D is denoted by ¶D, is defined by the geometrical equality ¶D = D Ç[`C(D)], where C(D) denotes the complementary set of D. This boundary satisfies : T ^{1}(¶D) = ¶D, T (¶D) Í ¶D. The basin D may be connected or nonconnected. A connected basin may be simply connected or multiply connected (which means connected with holes). An attracting set is called periodic of period k if it is made up of k disjoint sets, A = È_{i = 1}^{k}A_{i}, where each A_{i} is an attracting set of the map T^{k}. When A is an attracting set of period k = 1 , then its total basin is given by D = D_{0} if it is connected, by D = È_{n³ 0}T ^{n}(D_{0}) if it is nonconnected. When A is an attracting set of period k, the immediate basin D_{0} of A is the open set D_{0} = È_{i = 1}^{n}D_{0, i }, the D_{0,i} being open disjoints basins of A_{i}. If A is connected attractor , the immediate basin D_{0} of A is defined as the widest connected component of D containing A. When A is the widest attracting set of a map T , its basin D is the total basin of bounded iterates. That is, the open set D contains A such that D is the locus of points of the plane having bounded trajectories. Definition 2.1: A chaotic area A is an invariant absorbing area (T(A) = A)), the points of which give rise to iterated sequences having the property of sensitivity to initial conditions. About chaotic areas, it is important to emphasize that the study of such area has only the purpose to obtain properties giving rise to fractal basin boundaries. Definition 2.2: We say that l = l^{*} is a bifurcation of contact of A , if a contact between the frontier of A and the frontier its basin of attraction takes place. Proposition 2.1: When a bifurcation of contact of a chaotic area A arises for a value l = l^{*}, the crossing of this value leads to the destruction of A, either to a qualitative modification of properties of A (i.e. a sudden and important modification of the size of such an area or its basin of attraction). The destruction of A, after the crossing of the bifurcation value has been demonstrated by Gumowski & Mira [GM]. The qualitative change of properties of A has been described by Mira & al. [M4]. Definition 2.3: Let S be a saddle fixed point of T; a q point is called homocline to S, so q Î W^{s}(S)ÇW^{u}(S) and q ¹ S. q is a transversal homocline point, so W^{s}(S) intersects transversally W^{u}(S). We say the basin boundary is fractal if it contains a transversal homoclinic point. Definition 2.4: One calls homoclinic orbit O_{o}(q) associated with q, q belonging to a U(S) of S, a set constituted of successive iterates of q, and its infinite sequence of preimages obtained by application of the local inverse map T _{l}^{1}of T in U(S). Definition 2.5: Let T be a diffeomorphism of IR^{2} depending on a parameter l and let S be a saddle point of T. A homoclinic bifurcation takes place, if for a value l = l^{*}, there is apparition (or disappearance) of an infinity of homoclinic orbits. Overlapping of global stable and unstable manifolds usually
leads to nonlocal bifurcations. The existence of a transversal point of
a planar map leads to very complicated behavior of orbits nearby. Such
dynamical complexity is often dubbed as chaos, the onset of chaos typically
occurs at the parameters values for which the stable and unstable manifolds
of a saddle point come into contact tangentially.
Basins and attractors In this section we give a sketch on bifurcations which give rise to fractalization of attractor and its basin boundary. The main purpose of this section is to give some results on basin structure of invertible maps and their bifurcations. Especially, we explain basin bifurcations. These basic bifurcations result from the contact of a basin boundary with a chaotic attractor, such a bifurcation leads either to the chaotic area destruction, or a sudden and important modification of this area [GM, M5]. These bifurcations are called contact bifurcations. Different types of fractalization of basin boundary occur, a part of them resulting from a phenomenon of islands aggregation for nonconnected basins [M1]. It is worth noting that the results presented in this section were essentially obtained via a numerical way and using the bifurcation theory tool and the terminology of [M1]. Unfortunately taking into account the complexity of the matter and its nature, it seems difficult to carry out the study with success from another process. The following figures show the corresponding basin structure.
For fixed parameter values, it is possible to define in the phase plane
the basin of an attractor, that means the set of initial conditions giving
to iterated sequences converging towards the considered attractor. Figures
2. (a(1,2), b(1,2), c(1,2), d(1,2)) represent the existing attractors (order
k
cycles, invariant curves or chaotic attractors) and their basins.
Figure 2(a.1): A closed curve
and its basin e=0.15, m= 1.70, h=0.30.
Figure 2(a.2): Enlargement of a part of Fig. 2(a.1).
Figure 2(b.1): Chaotic attractor,
tongues appear on
Figure 2(b.2): Enlargement
of a part of Fig. 2(b.1).
Figure 2(c.1): Contact bifurcation
of the chaotic
Figure 2(c.2): Enlargement of a part of Fig. 2(c.1).
Figure 2(d.1): More and more
loops are created on the boundary
Figure 2(d.2): Situation
just after the disappearance of the
The evolution of attractors and their basins is given
directly in figures, the parameters (e,m)
have been chosen constant and h varies. There exist several ways
in which a dynamical system can become chaotic, of which the perioddoubling
route to chaos is the best known. Another scenario leading to chaos is
the destruction of a closed invariant curve, identified here in this model,
the closed curve can lose its smoothness and turn into an irregular invariant
set near a homoclinic structure formed by the intersection of the stable
and unstable manifolds of the saddle point. When h decreases, a
process gives rise to new appearance of tongues on basin boundary, and
there is appearance of oscillations of the closed invariant curve and then
we see a cyclic chaotic behavior. We can see also that the bifurcation
which is put in evidence, corresponds to a limit contact between the basin
boundary ¶D and the chaotic attractor.
A further decrease of the parameter h causes a contact between these
two boundaries which marks the destruction of the basin of attraction.
As argued in Ref [GM], we can see the ghost of this area
destroyed at the contact.
Multistability and Hopf bifurcation This diffeomorphism can be considered as a model giving
rise to a complete and interesting set of bifurcations with birth of cyclic
invariant curves (impossible to obtain with Hénon map) in relation
with quasi periodic solutions, or more complex solutions occuring in the
nonautonomous differential equation . Invariant closed curves are issued
from Hopf bifurcations (see Fig.3 (a, b, c, d, e, f, g, h, i, j, k, l)),
these complicated attracting sets exist for very small parameter ranges.
Arrowsmith & al. [A1] have explained the relationship
between the tongues in Hamiltonian structure and the invariant closed curves
(circles) and clarified various points. Eventual destruction of the Hopf
circles occurs as they grow outwards and their basin boundaries jump inwards
in a series of steps, accumulating on the stable manifolds of subharmonic
saddles.
Figure 3: A sequence of phase
portraits, we can see the basin of period 5,
Properties of W^{s}(S) and W^{u}(S) This section is devoted to the study of evolution of attractors when parameters vary. This evolution can involve relations with the attractor basin or homoclinic and heteroclinic bifurcations. In the case of a diffeomorphism, the stable invariant
manifold W^{s}(S) of
a saddle point S constitutes in general the frontier of the basin
of attraction of an associated attractor A. When W^{s}(S)
becomes nonconnected each of its connected components approaches again
itself to constitute the frontier of an island. It transforms the connected
attraction basin in nonconnected basin. Holes appear inside the basin,
which becomes multiply connected. All these phenomena generate a very big
sensitivity to initial conditions, due principally to multistability, interconnection
between basins and fractalization of basins.
Some properties of invariant curves For fixed parameter values, we plot the attraction basin of an attractor. When there exist several attracting sets, it is possible to define a global basin, that means the set of initial conditions giving rise to bounded iterated sequences, independently of the fact that they converge to one attractor or another. We shall see that particular homoclinic bifurcations which occur are related with the structure of the foliation of the invertible map. Another phenomena can also occur and generate sensitivity to initial conditions. This is the case where there exists two or more different attractors and where their own basins are considered. in this case, one basin constitute holes in the other. This situation is shown in Fig.5, where three different cycles coexist, two basins give rise to holes inside the third one. Proposition 4.1:Consider the invertible map T depending on a parameter l. If the number of connected components of D changes when l crosses through a bifurcation value l^{*}, then the basin D may undergo a qualitative change of one of the following types: (a) simply connected basin « multiply connected basin (when the number of connected components of D changes). (b) modification of the number of holes (lakes) in D, or new arborescent sequence of holes. The numeric studies give to the following results, for
values: e = 0.0, m
= 0.10, h = 1.44.
Figure 4: Plot of invariant manifolds of the saddle point (1,0).
Numerical plot of stable and unstable manifolds of the saddle point at (1,0) shows that the orbit sweeps around (0,0) twice, we can see homoclinic tangle. One has four attractors, the first is stable focus f (0,0) of which the basin of attraction has for external boundary the stable invariant manifold of the saddle point S (1,0). The second is the stable focus cycle of order 5 (represented in five small pink domains in the Figure 5).
Figure 5: f
coexists with an attractive focus of order
5,
These attractors coexist with a cycle of order 6 represented
by grey domains in the Figure 5, and an invariant stable invariant circle
(in yellow color).
Conclusion In this paper, specific bifurcations arising in invertible
maps of quadratic type have been studied. These bifurcations concern the
evolution of chaotic attractors, which appear or disappear by contact bifurcations,
with their own basin boundary. In such maps, chaotic attractors can be
issued of Hopf bifurcation giving rise to invariant closed curves which
evolve into chaotic attractor or of a succession of perioddoubling bifurcations
[F1, B2, B3]. As
usually in many nonlinear systems, important changes can occur with small
changes of parameter values.
Acknowledgment This work was supported by grants from ANDRU under N^{°}:
CU39904.
References [A1] D.K. Arrowsmith, The Bogdanov Map: Bifurcation, Mode Locking, and Chaos in Dissipative System, Int. J. Bifurcation and Chaos, 3, 803842 [1993]. [A2] D.K Arrowsmith & C.M. Place, An Introduction to dynamical systems, Cambridge University Press [1990]. [A3] V.I. Arnold, Lecture on bifurcations and versal systems, Russ Maths. Surveys 27, 54 [1972]. [B1] R.I. Bogdanov, Bifurcation of the limit cycle of a family of plane vector fields, Trudy Sem. Petrovsk. 2, 23. In Russian: the english translation is Sel. Math.Sov. 1(4), 37388 [1981]. [B2] H. Broer, C. Simo & J.C. Tatjer, Towards global models near homoclinic tangencies of dissipative diffeomorphisms, Nonlinearity, 11, 667770 [1998]. [B3] M. Bosch, J.C., Tatjer, On the codimensiontwo bifurcation in families of one and twodimensional maps, ECIT Batschuns (Australia) [1989]. [F1] D. Fournier, H. Kawakami & C. Mira, Feuilletage du plan des bifurcations d'un difféomorphisme bidimensionnel. Doublement de l'ordre des zones sources et des zones échangeurs. C.R. Acad. Sc. Paris, sér 1, 301(5), 223320 [1986]. [F2] D. Fournier, Contribution à l'étude de structures de bifurcations rencontrées en dynamique chaotique, Thèse de doctorat d'état (Toulouse) [1987]. [H1] M. Hénon, A twodimensional mapping with a strange attractor, Commun. Math. Phys. 50, 6978 [1976]. [H2] P. Holmes & D. Whitley, Bifurcation of one and two dimensional maps, Phil. Trans. R. Lond. A311, 43102 [1984]. [M1] C. Mira, Chaotic dynamics, World Scientific Publishing. Co. Pte. Ltd [1987]. [M2] C. Mira & C. Rauzy, Fractal Aggregation of Bassin Island in TwoDimensional Quadratic NonInversible Maps, Int. J. Bifurcation and Chaos, Vol. 5, No. 4 [1995] 9911019. [M3] C. Mira, Accumulations de bifurcations et structures boitesemboitées dans les récurrences et transformations ponctuelles, Proceeding of the VII^{e} International Conference on Non Linear Oscillations, Band I 2, 8193 [1975]. [M4] C. Mira, D. Fournier, L. Gardini, H. Kawakami & J.C. Cathala, Bassin bifurcations of twodimensional noninvertible maps Fractalization of basins, Int. J. Bifurcation and Chaos, vol. 4, no. 2, pp. 343382 [1994]. [M5] G. Millérioux, Modèles par Transformations Non Inversibles et Application à la Sécurité des Communications Synchronisation du Chaos, Thèse de Doctorat de L'INSA ( Toulouse) [1997]. [GM] I. Gumowski & C. Mira, Bifurcations déstabilisant une solution chaotique d'un endomorphisme du 2 ^{nd}ordre, Comptes Rendus Acad.SC; Paris A286, 427431(1978). [T] J.C. Tatjer, Invariant manifolds
and bifurcations for onedimensional and two dimensional dissipative maps,
PHD. Thesis University of Barcelona [1990].
