Bifurcation Theory

Total Page:16

File Type:pdf, Size:1020Kb

Bifurcation Theory Bifurcation Theory Ale Jan Homburg The axiomatization and algebraization of mathematics has led to the illegibility of such a large number of mathematical texts that the threat of complete loss of contact with physics and the natural sciences has been realized. (Vladimir Arnold) Contents 1 Introduction 5 1.1 The Rosenzweig-MacArthur model...............................8 1.2 Exercises............................................. 12 2 Elements of nonlinear analysis 17 2.1 Calculus.............................................. 17 2.1.1 Nemytskii operators................................... 19 2.2 Implicit function theorem.................................... 21 2.2.1 Applications of the implicit function theorem..................... 23 2.3 Lyapunov-Schmidt reduction.................................. 24 2.4 Floquet theory.......................................... 26 2.5 Poincar´ereturn maps...................................... 30 2.6 Center manifolds......................................... 33 2.7 Normal forms........................................... 36 2.8 Manifolds and transversality................................... 41 2.9 Sard's theorem.......................................... 42 2.10 Exercises............................................. 46 3 Local bifurcations 51 3.1 The saddle-node bifurcation................................... 51 3.1.1 One dimensional saddle-node bifurcation........................ 51 3.1.2 Higher dimensional saddle-node bifurcation...................... 53 3.1.3 Saddle-node bifurcation of periodic orbits....................... 54 3.2 The Hopf bifurcation....................................... 56 3.2.1 The planar Hopf bifurcation............................... 56 3.2.2 Higher dimensional Hopf bifurcation.......................... 64 3.3 The period-doubling bifurcation................................. 66 3.4 The Bogdanov-Takens bifurcation................................ 69 3.5 Exercises............................................. 71 3 4 4 Nonlocal bifurcations 77 4.1 Homoclinic loop.......................................... 78 4.2 Shil'nikov homoclinic loop.................................... 84 4.3 Blue sky catastrophe....................................... 85 4.4 Homoclinic tangencies...................................... 86 4.5 Exercises............................................. 94 5 Structural stability and bifurcation 99 5.1 Hyperbolicity and Peixoto's theorem.............................. 101 5.2 The Hartman-Grobman theorem................................ 108 5.3 The Smale horseshoe and the Arnold cat map......................... 109 5.4 Stability theorems........................................ 113 5.5 Global aspects of bifurcations.................................. 117 5.6 Exercises............................................. 119 Chapter 1 Introduction The qualitative or geometric study of dynamical systems originates with Henri Poincar´ewho studied differential equations appearing in problems from celestial mechanics. The study of dynamical systems has been one of the most successful fields of mathematical research with an explosive development in the last fifty years. The research field has gradually developed into a more fragmented field. There is a theoretical core aimed at the description of typical dynamics, i.e., common properties of most dynamical systems, ignoring pathological systems. This part of the research follows a natural development where dynamics of increasing complexity is studied. The focus of the research here has shifted from geometrical techniques to stochastic techniques, always in combination with analysis. Such a combination of techniques has given a special flavor to the study of dynamical systems. On the other hand, with the development of a toolbox for the study of dynamics, the body of applied research using techniques from dynamical systems has expanded dramatically. Whether in physics, engineering, biology, psychology or economics, tools from dynamical systems have been highly successful in understanding phenomena where evolution is the focus. Dynamical systems ideas also appear frequently in other areas of mathematics; there is for instance a lively interaction between number theory and combinatorics with the part of dynamical systems theory called ergodic theory. This syllabus concerns the study of changes of dynamical properties, as the rules defining the dynam- ical system changes. This is the research area of bifurcation theory. It has been powerful in the study of typical dynamics, and is especially powerful in applications in other sciences. My aim was to compose a syllabus on bifurcation theory that serves as a first course on it, assuming standard bachelor courses on topology, analysis and ordinary differential equations. I wanted it to address techniques, albeit sometimes in simplified settings, in a detailed manner and further to give a broad overview of results, and demonstrate its usefulness in applications and the modern study of dynamical systems. I also tried to make visible the scientific heritage of the field; the collection of techniques and results that have been guiding its development. A course on bifurcation theory does not substitute a course on dynamical systems theory. It leaves out many topics, including topological dynamical systems, differentiable chaotic dynamical systems, and ergodic theory. Various aspects and 5 6 CHAPTER 1. INTRODUCTION topics are only skimmed or left out in this course. A first course such as this one leaves out advanced topics, which are interesting from a bifurcation theory point of view. This includes equivariant dynamical systems, Hamiltonian systems, dynamics on networks, infinite dimensional systems (partial differential equations, delay differential equations) and random dynamical systems (stochastic differential equations). The term bifurcation was originally used by Poincar´eto describe the splitting of equilibria in a family of differential equations. In modern use, a bifurcation of a dynamical system is a qualitative change in its dynamics produced by varying parameters. 2 Figure 0.1: The bifurcation diagram for fµ(x) = µ − x . The dashed curve stands for sources, the solid curve for sinks. Example. Consider the differential equationx _ = f(µ, x) with µ, x 2 R, f(µ, x) = µ − x2. Equilibria occur if f(µ, x) = 0 and lie on the curve µ = x2. Thus there are no equilibria if µ < 0, one equilibrium p p at 0 if µ = 0 and two equilibria ± µ if µ > 0. If µ = 0, then Df0(0) = 0. If µ > 0, then Dfµ( µ) = p p p p p −2 µ < 0 and Dfµ(− µ) = 2 µ > 0. This means that µ is a stable equilibrium (a sink), and − µ is unstable (a source). We obtain the qualitative picture of Figure 0.1, in which the branches of equilibria are shown in (µ, x) space. This figure is an example of a bifurcation diagram. More in general, consider a smooth differential equationx _ = f(µ, x) with f(µ, x) = µ − x2 + O(jxj3 + jµxj + jµj2) as x; µ ! 0. Equilibria are solutions of the equation f(µ, x) = 0. Our objective is to compare this with the solution of the equation µ − x2 = 0 obtained by considering only the lowest order terms in the Taylor expansion. The equation f(µ, x) = 0 can be solved near (0; 0) for µ = µ∗(x), µ∗(0) = 0. One computes Dµ∗(0) = 0, D2µ∗(0) > 0. To find the number of solutions on either side of the curve, we first note that 2 Dxf(µ, x) = −2x + O(jxj + jµj) 7 as x; µ ! 0. In a neighborhood of (µ, x) = (0; 0), there is a unique function x∗(µ) near µ = 0, x∗(0) = 0, ∗ ∗ such that Dxf((µ, x (µ)) = 0 since Dxxf(0; 0) 6= 0. The point x (µ) corresponds to the maximum of f(µ, x). If α(µ) = f(µ, x∗(µ)), then α(µ) < 0 implies no solutions near (0; 0) and α(µ) > 0 implies two solutions near (0; 0). It is easy to verify that α(µ) = µ + O(jµj2) as µ ! 0. 3 Figure 0.2: The surface x − λ1x − λ2 = 0 in (λ1; λ2; x) space, with projection of the folds onto the cusp in the (λ1; λ2) plane. 2 Example. As a slightly more complicated example, suppose x 2 R, λ = (λ1; λ2) 2 R and consider the equation 3 4 2 2 2 f(λ, x) = x − λ1x − λ2 + O(jxj + jλ1x j + jλ1xj + jλ2xj + jλ2j + jλ1λ2j) = 0 (0.1) in a small neighborhood of (λ2; λ2; x) = (0; 0; 0). As before, consider the polynomial equation 3 x − λ1x − λ2 = 0 (0.2) and try to compare (0.1) with (0.2). For (0.2) it is easy to determine the bifurcation curves in (λ1; λ2)- space. In fact, these curves are exactly the curves across which the number of solutions changes. That is, multiple solutions occur on these curves. Any multiple solution of (0.2) must satisfy (0.2) as well as 2 3x − λ1 = 0: (0.3) 2 3 Equations (0.2) and (0.3) uniquely define λ1; λ2 as functions of x, namely, λ1 = 3x , λ2 = 2x . This is 3 2 a parametric representation of the cusp 4λ1 − 27λ2 = 0. The basic question is whether the solutions to (0.2) have the same qualitative properties as those to (0.1). The procedure is the same as before. The bifurcation curves solve f(λ, x) = 0 and Dxf(λ, x) = 0. Near λ = 0; x = 0, these equations uniquely determine λ1; λ2 as function of x. One checks λ1 = 2 3 3 4 3x + O(jxj ) and λ2 = 2x + O(jxj ) as x ! 0. These are the parametric representation for a cusp with 3 2 essentially the same form as 4λ1 − 27λ2 = 0. One must show that the number of solutions changes by two as this cusp is crossed. This involves discussing Dxxf(λ, x) and Dλ2 f(λ, x) along this curve. 8 CHAPTER 1. INTRODUCTION We are often interested in the dynamics after long times and less in transients. We will concentrate on attractors, and bifurcations thereof. Definition 0.1. A compact invariant set A of a differential equation x_ = f(x) with flow φt is an attractor if (i) A is Lyapunov stable: for each open neighborhood U of A there is a neighborhood V ⊂ U of A so that φt(x) 2 U for each t ≥ 0; x 2 V . (ii) A is asymptotically stable: there is an open neighborhood U of A so that φt(x) ! A (i.e. the distance between φt(x) and A goes to zero) as t ! 1. (iii) A is minimal: there is no strict and nontrivial subset of A satisfying the first two properties. Other definitions have been proposed to deal with compact invariant sets that attract many points but not necessarily an open neighborhood.
Recommended publications
  • Approximating Stable and Unstable Manifolds in Experiments
    PHYSICAL REVIEW E 67, 037201 ͑2003͒ Approximating stable and unstable manifolds in experiments Ioana Triandaf,1 Erik M. Bollt,2 and Ira B. Schwartz1 1Code 6792, Plasma Physics Division, Naval Research Laboratory, Washington, DC 20375 2Department of Mathematics and Computer Science, Clarkson University, P.O. Box 5815, Potsdam, New York 13699 ͑Received 5 August 2002; revised manuscript received 23 December 2002; published 12 March 2003͒ We introduce a procedure to reveal invariant stable and unstable manifolds, given only experimental data. We assume a model is not available and show how coordinate delay embedding coupled with invariant phase space regions can be used to construct stable and unstable manifolds of an embedded saddle. We show that the method is able to capture the fine structure of the manifold, is independent of dimension, and is efficient relative to previous techniques. DOI: 10.1103/PhysRevE.67.037201 PACS number͑s͒: 05.45.Ac, 42.60.Mi Many nonlinear phenomena can be explained by under- We consider a basin saddle of Eq. ͑1͒ which lies on the standing the behavior of the unstable dynamical objects basin boundary between a chaotic attractor and a period-four present in the dynamics. Dynamical mechanisms underlying periodic attractor. The chosen parameters are in a region in chaos may be described by examining the stable and unstable which the chaotic attractor disappears and only a periodic invariant manifolds corresponding to unstable objects, such attractor persists along with chaotic transients due to inter- as saddles ͓1͔. Applications of the manifold topology have secting stable and unstable manifolds of the basin saddle.
    [Show full text]
  • Synchronization of Chaotic Systems
    Synchronization of chaotic systems Cite as: Chaos 25, 097611 (2015); https://doi.org/10.1063/1.4917383 Submitted: 08 January 2015 . Accepted: 18 March 2015 . Published Online: 16 April 2015 Louis M. Pecora, and Thomas L. Carroll ARTICLES YOU MAY BE INTERESTED IN Nonlinear time-series analysis revisited Chaos: An Interdisciplinary Journal of Nonlinear Science 25, 097610 (2015); https:// doi.org/10.1063/1.4917289 Using machine learning to replicate chaotic attractors and calculate Lyapunov exponents from data Chaos: An Interdisciplinary Journal of Nonlinear Science 27, 121102 (2017); https:// doi.org/10.1063/1.5010300 Attractor reconstruction by machine learning Chaos: An Interdisciplinary Journal of Nonlinear Science 28, 061104 (2018); https:// doi.org/10.1063/1.5039508 Chaos 25, 097611 (2015); https://doi.org/10.1063/1.4917383 25, 097611 © 2015 Author(s). CHAOS 25, 097611 (2015) Synchronization of chaotic systems Louis M. Pecora and Thomas L. Carroll U.S. Naval Research Laboratory, Washington, District of Columbia 20375, USA (Received 8 January 2015; accepted 18 March 2015; published online 16 April 2015) We review some of the history and early work in the area of synchronization in chaotic systems. We start with our own discovery of the phenomenon, but go on to establish the historical timeline of this topic back to the earliest known paper. The topic of synchronization of chaotic systems has always been intriguing, since chaotic systems are known to resist synchronization because of their positive Lyapunov exponents. The convergence of the two systems to identical trajectories is a surprise. We show how people originally thought about this process and how the concept of synchronization changed over the years to a more geometric view using synchronization manifolds.
    [Show full text]
  • Hartman-Grobman and Stable Manifold Theorem
    THE HARTMAN-GROBMAN AND STABLE MANIFOLD THEOREM SLOBODAN N. SIMIC´ This is a summary of some basic results in dynamical systems we discussed in class. Recall that there are two kinds of dynamical systems: discrete and continuous. A discrete dynamical system is map f : M → M of some space M.A continuous dynamical system (or flow) is a collection of maps φt : M → M, with t ∈ R, satisfying φ0(x) = x and φs+t(x) = φs(φt(x)), for all x ∈ M and s, t ∈ R. We call M the phase space. It is usually either a subset of a Euclidean space, a metric space or a smooth manifold (think of surfaces in 3-space). Similarly, f or φt are usually nice maps, e.g., continuous or differentiable. We will focus on smooth (i.e., differentiable any number of times) discrete dynamical systems. Let f : M → M be one such system. Definition 1. The positive or forward orbit of p ∈ M is the set 2 O+(p) = {p, f(p), f (p),...}. If f is invertible, we define the (full) orbit of p by k O(p) = {f (p): k ∈ Z}. We can interpret a point p ∈ M as the state of some physical system at time t = 0 and f k(p) as its state at time t = k. The phase space M is the set of all possible states of the system. The goal of the theory of dynamical systems is to understand the long-term behavior of “most” orbits. That is, what happens to f k(p), as k → ∞, for most initial conditions p ∈ M? The first step towards this goal is to understand orbits which have the simplest possible behavior, namely fixed and periodic ones.
    [Show full text]
  • Chapter 3 Hyperbolicity
    Chapter 3 Hyperbolicity 3.1 The stable manifold theorem In this section we shall mix two typical ideas of hyperbolic dynamical systems: the use of a linear approxi- mation to understand the behavior of nonlinear systems, and the examination of the behavior of the system along a reference orbit. The latter idea is embodied in the notion of stable set. Denition 3.1.1: Let (X, f) be a dynamical system on a metric space X. The stable set of a point x ∈ X is s ∈ k k Wf (x)= y X lim d f (y),f (x) =0 . k→∞ ∈ s In other words, y Wf (x) if and only if the orbit of y is asymptotic to the orbit of x. If no confusion will arise, we drop the index f and write just W s(x) for the stable set of x. Clearly, two stable sets either coincide or are disjoint; therefore X is the disjoint union of its stable sets. Furthermore, f W s(x) W s f(x) for any x ∈ X. The twin notion of unstable set for homeomorphisms is dened in a similar way: Denition 3.1.2: Let f: X → X be a homeomorphism of a metric space X. Then the unstable set of a ∈ point x X is u ∈ k k Wf (x)= y X lim d f (y),f (x) =0 . k→∞ u s In other words, Wf (x)=Wf1(x). Again, unstable sets of homeomorphisms give a partition of the space. On the other hand, stable sets and unstable sets (even of the same point) might intersect, giving rise to interesting dynamical phenomena.
    [Show full text]
  • Role of Nonlinear Dynamics and Chaos in Applied Sciences
    v.;.;.:.:.:.;.;.^ ROLE OF NONLINEAR DYNAMICS AND CHAOS IN APPLIED SCIENCES by Quissan V. Lawande and Nirupam Maiti Theoretical Physics Oivisipn 2000 Please be aware that all of the Missing Pages in this document were originally blank pages BARC/2OOO/E/OO3 GOVERNMENT OF INDIA ATOMIC ENERGY COMMISSION ROLE OF NONLINEAR DYNAMICS AND CHAOS IN APPLIED SCIENCES by Quissan V. Lawande and Nirupam Maiti Theoretical Physics Division BHABHA ATOMIC RESEARCH CENTRE MUMBAI, INDIA 2000 BARC/2000/E/003 BIBLIOGRAPHIC DESCRIPTION SHEET FOR TECHNICAL REPORT (as per IS : 9400 - 1980) 01 Security classification: Unclassified • 02 Distribution: External 03 Report status: New 04 Series: BARC External • 05 Report type: Technical Report 06 Report No. : BARC/2000/E/003 07 Part No. or Volume No. : 08 Contract No.: 10 Title and subtitle: Role of nonlinear dynamics and chaos in applied sciences 11 Collation: 111 p., figs., ills. 13 Project No. : 20 Personal authors): Quissan V. Lawande; Nirupam Maiti 21 Affiliation ofauthor(s): Theoretical Physics Division, Bhabha Atomic Research Centre, Mumbai 22 Corporate authoifs): Bhabha Atomic Research Centre, Mumbai - 400 085 23 Originating unit : Theoretical Physics Division, BARC, Mumbai 24 Sponsors) Name: Department of Atomic Energy Type: Government Contd...(ii) -l- 30 Date of submission: January 2000 31 Publication/Issue date: February 2000 40 Publisher/Distributor: Head, Library and Information Services Division, Bhabha Atomic Research Centre, Mumbai 42 Form of distribution: Hard copy 50 Language of text: English 51 Language of summary: English 52 No. of references: 40 refs. 53 Gives data on: Abstract: Nonlinear dynamics manifests itself in a number of phenomena in both laboratory and day to day dealings.
    [Show full text]
  • Invariant Manifolds and Their Interactions
    Understanding the Geometry of Dynamics: Invariant Manifolds and their Interactions Hinke M. Osinga Department of Mathematics, University of Auckland, Auckland, New Zealand Abstract Dynamical systems theory studies the behaviour of time-dependent processes. For example, simulation of a weather model, starting from weather conditions known today, gives information about the weather tomorrow or a few days in advance. Dynamical sys- tems theory helps to explain the uncertainties in those predictions, or why it is pretty much impossible to forecast the weather more than a week ahead. In fact, dynamical sys- tems theory is a geometric subject: it seeks to identify critical boundaries and appropriate parameter ranges for which certain behaviour can be observed. The beauty of the theory lies in the fact that one can determine and prove many characteristics of such bound- aries called invariant manifolds. These are objects that typically do not have explicit mathematical expressions and must be found with advanced numerical techniques. This paper reviews recent developments and illustrates how numerical methods for dynamical systems can stimulate novel theoretical advances in the field. 1 Introduction Mathematical models in the form of systems of ordinary differential equations arise in numer- ous areas of science and engineering|including laser physics, chemical reaction dynamics, cell modelling and control engineering, to name but a few; see, for example, [12, 14, 35] as en- try points into the extensive literature. Such models are used to help understand how different types of behaviour arise under a variety of external or intrinsic influences, and also when system parameters are changed; examples include finding the bursting threshold for neurons, determin- ing the structural stability limit of a building in an earthquake, or characterising the onset of synchronisation for a pair of coupled lasers.
    [Show full text]
  • Math Morphing Proximate and Evolutionary Mechanisms
    Curriculum Units by Fellows of the Yale-New Haven Teachers Institute 2009 Volume V: Evolutionary Medicine Math Morphing Proximate and Evolutionary Mechanisms Curriculum Unit 09.05.09 by Kenneth William Spinka Introduction Background Essential Questions Lesson Plans Website Student Resources Glossary Of Terms Bibliography Appendix Introduction An important theoretical development was Nikolaas Tinbergen's distinction made originally in ethology between evolutionary and proximate mechanisms; Randolph M. Nesse and George C. Williams summarize its relevance to medicine: All biological traits need two kinds of explanation: proximate and evolutionary. The proximate explanation for a disease describes what is wrong in the bodily mechanism of individuals affected Curriculum Unit 09.05.09 1 of 27 by it. An evolutionary explanation is completely different. Instead of explaining why people are different, it explains why we are all the same in ways that leave us vulnerable to disease. Why do we all have wisdom teeth, an appendix, and cells that if triggered can rampantly multiply out of control? [1] A fractal is generally "a rough or fragmented geometric shape that can be split into parts, each of which is (at least approximately) a reduced-size copy of the whole," a property called self-similarity. The term was coined by Beno?t Mandelbrot in 1975 and was derived from the Latin fractus meaning "broken" or "fractured." A mathematical fractal is based on an equation that undergoes iteration, a form of feedback based on recursion. http://www.kwsi.com/ynhti2009/image01.html A fractal often has the following features: 1. It has a fine structure at arbitrarily small scales.
    [Show full text]
  • An Overview of Bifurcation Theory
    Appendix A An Overview of Bifurcation Theory A.I Introduction In this appendix we want to provide a brief introduction and discussion of the con­ cepts of dynamical systems and bifurcation theory which has been used in the pre­ ceding sections . We refer the reader interested in a more thorough discussion of the mathematical results of dynamical systems and bifurcation theory to the books of Wiggins (1990) and Kuznetsov (1995). A discussion intended to more econom­ ically motivated problems of dynamical systems and chaos can be found in Medio (1992) or Day (1994). It is noteworthy here that chaos is only possible in the non­ autonomous case; that is, the dynamical system depends explicity on time; but not for the autonomous case where time does not enter directly in the dynamical sys­ tem (cf. Wiggins (1990». Since both models we have considered in the preceding sections are autonomous they cannot reveal any chaos . However, as the reader has seen, the dynamical systems we have derived reveal some bifurcations. Roughly spoken, bifurcation theory describes the way in which dynamical sys­ tem changes due to a small perturbation of the system-parameters. A qualitative change in the phase space of the dynamical system occurs at a bifurcation point, that means that the system is structural unstable against a small perturbation in the parameter space and the dynamic structure of the system has changed due to this slight variation in the parameter space. We can distinguish bifurcations into local, global and local/global types, where only the pure types are interesting here. Local bifurcations such as Andronov-Hopfor Fold bifurcations can be detected in a small 182 An Overview of Bifurcation Theory neighborhood ofa single fixed point or stationary point.
    [Show full text]
  • Computation of Stable and Unstable Manifolds of Hyperbolic Trajectories in Two-Dimensional, Aperiodically Time-Dependent Vector fields Ana M
    Physica D 182 (2003) 188–222 Computation of stable and unstable manifolds of hyperbolic trajectories in two-dimensional, aperiodically time-dependent vector fields Ana M. Mancho a, Des Small a, Stephen Wiggins a,∗, Kayo Ide b a School of Mathematics, University of Bristol, University Walk, Bristol, BS8 1TW, UK b Department of Atmospheric Sciences, and Institute of Geophysics and Planetary Physics, UCLA, Los Angeles, CA 90095-1565, USA Received 14 October 2002; received in revised form 29 January 2003; accepted 31 March 2003 Communicated by E. Kostelich Abstract In this paper, we develop two accurate and fast algorithms for the computation of the stable and unstable manifolds of hyperbolic trajectories of two-dimensional, aperiodically time-dependent vector fields. First we develop a benchmark method in which all the trajectories composing the manifold are integrated from the neighborhood of the hyperbolic trajectory. This choice, although very accurate, is not fast and has limited usage. A faster and more powerful algorithm requires the insertion of new points in the manifold as it evolves in time. Its numerical implementation requires a criterion for determining when to insert those points in the manifold, and an interpolation method for determining where to insert them. We compare four different point insertion criteria and four different interpolation methods. We discuss the computational requirements of all of these methods. We find two of the four point insertion criteria to be accurate and robust. One is a variant of a criterion originally proposed by Hobson. The other is a slight variant of a method due to Dritschel and Ambaum arising from their studies of contour dynamics.
    [Show full text]
  • Theoretical-Heuristic Derivation Sommerfeld's Fine Structure
    Theoretical-heuristic derivation Sommerfeld’s fine structure constant by Feigenbaum’s constant (delta): perodic logistic maps of double bifurcation Angel Garcés Doz [email protected] Abstract In an article recently published in Vixra: http://vixra.org/abs/1704.0365. Its author (Mario Hieb) conjectured the possible relationship of Feigen- baum’s constant delta with the fine-structure constant of electromag- netism (Sommerfeld’s Fine-Structure Constant). In this article it demon- strated, that indeed, there is an unequivocal physical-mathematical rela- tionship. The logistic map of double bifurcation is a physical image of the random process of the creation-annihilation of virtual pairs lepton- antilepton with electric charge; Using virtual photons. The probability of emission or absorption of a photon by an electron is precisely the fine structure constant for zero momentum, that is to say: Sommerfeld’s Fine- Structure Constant. This probability is coded as the surface of a sphere, or equivalently: four times the surface of a circle. The original, conjectured calculation of Mario Hieb is corrected or improved by the contribution of the entropies of the virtual pairs of leptons with electric charge: muon, tau and electron. Including a correction factor due to the contributions of virtual bosons W and Z; And its decay in electrically charged leptons and quarks. Introduction The main geometric-mathematical characteristics of the logistic map of double universal bifurcation, which determines the two Feigenbaum’s constants; they are: 1) The first Feigenbaum constant is the limiting ratio of each bifurcation interval to the next between every period doubling, of a one-parameter map.
    [Show full text]
  • INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS of NONAUTONOMOUS DIFFERENTIAL EQUATIONS Bernd Aulbach1 and Martin Rasmussen2 Stefan S
    DISCRETE AND CONTINUOUS Website: http://aimSciences.org DYNAMICAL SYSTEMS Volume 15, Number 2, June 2006 pp. 579-596 INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS OF NONAUTONOMOUS DIFFERENTIAL EQUATIONS Bernd Aulbach1 and Martin Rasmussen2 Department of Mathematics University of Augsburg D-86135 Augsburg, Germany Stefan Siegmund3 Department of Mathematics University of Frankfurt D-60325 Frankfurt am Main, Germany (Communicated by Glenn F. Webb) Abstract. We discuss the relationship between invariant manifolds of nonau- tonomous di®erential equations and pullback attractors. This relationship is essential, e.g., for the numerical approximation of these manifolds. In the ¯rst step, we show that the unstable manifold is the pullback attractor of the di®er- ential equation. The main result says that every (hyperbolic or nonhyperbolic) invariant manifold is the pullback attractor of a related system which we con- struct explicitly using spectral transformations. To illustrate our theorem, we present an application to the Lorenz system and approximate numerically the stable as well as the strong stable manifold of the origin. 1. Introduction. In the theory of dynamical systems (e.g., autonomous di®eren- tial equations) the stable and unstable manifold theorem|¯rst proved by Poincar¶e [13] and Hadamard [10]|plays an important role in the study of the flow near a hyperbolic rest point or a periodic solution. The idea is to show that the rest point x = 0, say, locally inherits the dynamical structure of its hyperbolic linearization xAx_ , i.e., the spectrum s u σ(A) = σ (A) [ σ (A)f¸1; : : : ; ¸rg [ f¸r+1; : : : ; ¸ng consists of eigenvalues ¸1; : : : ; ¸r with negative real parts and ¸r+1; : : : ; ¸n with positive real parts.
    [Show full text]
  • Matrix Algebra I N THIS APPENDIX We Review Important Concepts From
    ApPENDIX A I Matrix Algebra IN THIS APPENDIX we review important concepts from matrix algebra. A.I EIGENVALUES AND EIGENVECTORS Definition A.I LetA be an m X m matrix and x a nonzero m-dimensional real or complex vector. If Ax = Ax for some real or complex number A, then A is called an eigenvalue of A and x the corresponding eigenvector. For example, the matrix A = (~ ~) has an eigenvector (~), and cor­ responding eigenvalue 4. 557 MATRIX ALGEBRA Eigenvalues are found as roots A of the characteristic polynomial det(A - AI). If A is an eigenvalue of A, then any nonzero vector in the nullspace of A - AI is an eigenvector corresponding to A. For this example, det(A - AI) = det e; A 2 ~ A) = (A - l)(A - 2) - 6 = (A - 4)(A + 1), (A1) so the eigenvalues are A = 4, -1. The eigenvectors corresponding to A = 4 are found in the nullspace of = (A2) A- 41 (-32 -23) and so consist of all nonzero multiples of ( ~ ) . Similarly, the eigenvectors corre­ sponding to A = -1 are all nonzero multiples of ( _ ~ ) . If A and Bare m X m matrices, then the set of eigenvalues of the product matrix AB is the same as the set of eigenvalues of the product BA In fact, let A be an eigenvalue of AB, so that ABx = Ax for some x '* O. If A = 0, then 0= det(AB - AI) = det(AB) = (detA)(detB) = det(BA - AI), so A = 0 is also an eigenvalue of BA If A '* 0, then ABx = Ax implies BA(Bx) = ABx.
    [Show full text]