= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :

Similar documents
8.1 Bifurcations of Equilibria

7 Two-dimensional bifurcations

Dynamical Systems in Neuroscience: Elementary Bifurcations

MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation

MATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation.

Computational Neuroscience. Session 4-2

CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION

10 Back to planar nonlinear systems

One Dimensional Dynamical Systems

Nonlinear dynamics & chaos BECS

(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,

EE222 - Spring 16 - Lecture 2 Notes 1

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems

Stability of Dynamical systems

ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

Application demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems

Example of a Blue Sky Catastrophe

Chapter 8 Equilibria in Nonlinear Systems

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

Summary of topics relevant for the final. p. 1

Mathematical Modeling I

Figure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., 2002.

Local Phase Portrait of Nonlinear Systems Near Equilibria

Lecture 3 : Bifurcation Analysis

2 Discrete growth models, logistic map (Murray, Chapter 2)

Introduction to Applied Nonlinear Dynamical Systems and Chaos

Solution to Homework #5 Roy Malka 1. Questions 2,3,4 of Homework #5 of M. Cross class. dv (x) dx

College of Natural Science Department of Mathematics

Edward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology

Dynamical systems tutorial. Gregor Schöner, INI, RUB

Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps

Mixed mode oscillations as a mechanism for pseudo-plateau bursting

Chapter 24 BIFURCATIONS

A Producer-Consumer Model With Stoichiometry

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

Phase Oscillators. and at r, Hence, the limit cycle at r = r is stable if and only if Λ (r ) < 0.

STABILITY. Phase portraits and local stability

Problem Set Number 2, j/2.036j MIT (Fall 2014)

APPPHYS217 Tuesday 25 May 2010

TWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations

arxiv: v2 [physics.optics] 15 Dec 2016

Bifurcation of Fixed Points

Bifurcation Analysis of Neuronal Bursters Models

Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting

Numerical techniques: Deterministic Dynamical Systems

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

DYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS

Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells

Boulder School for Condensed Matter and Materials Physics. Laurette Tuckerman PMMH-ESPCI-CNRS

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

Modelling in Biology

Part II. Dynamical Systems. Year

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d

Feedback-mediated oscillatory coexistence in the chemostat

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

Hello everyone, Best, Josh

Lotka Volterra Predator-Prey Model with a Predating Scavenger

Classification of Phase Portraits at Equilibria for u (t) = f( u(t))

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

A Mathematical Study of Electrical Bursting in Pituitary Cells

Outline. Learning Objectives. References. Lecture 2: Second-order Systems

THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A 2 2 SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS

Chapter 7. Nonlinear Systems. 7.1 Introduction

Continuous Threshold Policy Harvesting in Predator-Prey Models

Nonlinear Autonomous Dynamical systems of two dimensions. Part A

Problem Set 5 Solutions

Chimera States in Nonlocal Phase-coupled Oscillators. Jianbo Xie. A dissertation submitted in partial satisfaction of the

Reliable Computation of Equilibrium States and Bifurcations in Food Chain Models

MathQuest: Differential Equations

11 Chaos in Continuous Dynamical Systems.

Synchronization and Phase Oscillators

Phase Portraits of Nonlinear Differential Equations

Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University

THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL

8 Vector Field Topology

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting

Math Notes on sections 7.8,9.1, and 9.3. Derivation of a solution in the repeated roots case: 3 4 A = 1 1. x =e t : + e t w 2.

Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting

Models Involving Interactions between Predator and Prey Populations

8 Autonomous first order ODE

arxiv: v1 [nlin.ao] 19 May 2017

NBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011

Math 5490 November 5, 2014

Copyright (c) 2006 Warren Weckesser

2 Lecture 2: Amplitude equations and Hopf bifurcations

ANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM. China

NONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis

BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs

BLUE SKY CATASTROPHE IN SINGULARLY-PERTURBED SYSTEMS

A showcase of torus canards in neuronal bursters

Nonlinear Dynamics & Chaos MAT D J B Lloyd. Department of Mathematics University of Surrey

+ i. cos(t) + 2 sin(t) + c 2.

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

Department of Mathematics IIT Guwahati

Yulin, A. V., & Champneys, A. R. (2009). Discrete snaking: multiple cavity solitons in saturable media.

6.3. Nonlinear Systems of Equations

arxiv: v1 [nlin.cd] 28 Sep 2017

Transcription:

1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306

A bifurcation is a qualitative change in the asymptotic structure of a dynamical system. Stationary bifurcations involve changes in the number and/or stability of steady states. Periodic bifurcations involve changes in the number and/or 2 stability of periodic solutions. All of these bifurcations can occur in continuous dynamical systems of two dimensions or higher, but for simplicity we will discuss planar systems (the generalization to higher dimensions is straight forward). Consider a typical nonlinear planar system d dt = F ( ; ) (1) where is a 2-dimensional vector, is a parameter, and F : R 2 R 2 is a nonlinear function. The roots of F are the steady states and their stability can be determined by linearizing the system about these points. Let p be a steady state of the planar system, then the linearization about p is d dt = J (2) where J is the 2 2 Jacobian matri evaluated at p.

3 Stationary bifurcations A steady state p is a node if the two eigenvalues λ 1 and λ 2 of J are real and have the same sign. If the eignenvalues are negative, then p is stable; otherwise it is unstable. The steady state is a saddle point if the eigenvalues are real with opposite signs. If the eigenvalues form a comple conjugate pair, then p is a focus (also called a spiral). If both eigenvalues of p have non-zero real parts, then the steady state is hyperbolic. This is the typical case, and by the Hartman- Grobman theorem the linearization is a good description of the flow in the neighborhood of a hyperbolic steady state. If, in contrast, p has an eigenvalue with zero real part it is called non-hyperbolic and the flow of the linearlized system may not be a good description of the flow of the nonlinear system near p. Suppose that Eq. 1 has a non-hyperbolic equilibrium when parameter =. Then there is a bifurcation at this value of the parameter. Bifurcations occur at non-hyperbolic steady states. At a stationary bifurcation one of the eigenvalues is zero. At

a saddle-node bifurcation a saddle point and a node coalesce, creating a single steady state with a zero eigenvalue. For on one side of this bifurcation there are no steady states; for on the other side there are two steady states. 4 Figure 1: Saddle node bifurcation at =. Solid=stable, dashed=unstable. Figure 2: Pitchfork bifurcation at =. Solid=stable, dashed=unstable. A pitchfork bifurcation can occur in symmetric systems. The number of equilibria goes from 1 to 3 as this bifurcation

is crossed and the stability of the original equilibrium changes. At a transcritical bifurcation the stability is transferred from one equilibrium point to another. There are two equilibria for on either side of, and one equilibrium when =. The stable equilibrium for < becomes unstable for > 5 and vice versa. Figure 3: Transcritical bifurcation at =. Solid=stable, dashed=unstable. Periodic bifurcations A bifurcation involving the change in the number or stability of a periodic solution, or limit cycle, is a periodic bifurcation. A Hopf bifurcation is a periodic bifurcation in which a new limit cycle is born from a stationary solution. The stationary solution changes stability as the limit cycle is born. This oc-

curs when the eigenvalues of the stationary solution cross the imaginary ais. Suppose that p is a steady state and the eigenvalues of the linearization are purely imaginary when =, but for some neighborhood to the left of and some neighborhood to the right of the eigenvalues are comple with non-zero real part. Then there is a Hopf bifurcation at =. Im(λ) 6 Re( λ) Figure 4: The eigenvalues cross through the imaginary ais at a Hopf bifurcation. The red points indicate the Hopf bifurcation. According to the Hopf Bifurcation Theorem, the following things are true at a Hopf bifurcation: The limit cycle born at the bifurcation has 0 amplitude. As is moved away from the amplitude A grows as A for in some neighborhood of.

The period T of the limit cycle born at the bifurcation is T = 2π Im(λ). y 7 Figure 5: The new limit cycle (red) and the steady state (black) have opposite stability properties. At a supercritical Hopf bifurcation the limit cycle that is born is stable. At a subcritical Hopf bifurcation the limit cycle is born unstable. In the associated bifurcation diagrams the periodic branch that emerges from the stationary branch at the bifurcation is stable (supercritical) or unstable (subcritical). Another type of periodic bifurcation is a saddle-node of periodics (SNP) bifurcation. This occurs when a stable limit cycle coalesces with an unstable limit cycle, creating a limit cycle that attracts from one direction and repels from another (it might attract phase points inside the cycle and repel those

8 ma() min() Figure 6: Typical bifurcation diagram near a supercritical Hopf bifurcation. Black=stationary branch, Red=min and ma of stable periodic branch. ma() min() Figure 7: Typical bifurcation diagram near a subcritical Hopf bifurcation. Black=stationary branch, Red=min and ma of unstable periodic branch.

outside, or vice versa). As the parameter is moved past the bifurcation point the limit cycle disappears. y 9 Figure 8: Two limit cycles near a saddle-node of periodics (SNP) bifurcation. ma() min() Figure 9: At an (SNP) bifurcation two periodic branches, one stable and the other unstable, join together at a turning point. In the diagram there appears to be two turning points, but there is only one since the upper curve represents the maimum of the limit cycles and the lower curve represents the minimum. In terms of the bifurcation diagram, a periodic branch has a turning point at an SNP bifurcation. One incoming periodic

10 branch is stable, while the other is unstable. There are several other types of periodic bifurcations, but we will only discuss a homoclinic bifurcation. This bifurcation terminates a periodic solution branch and the oscillation period approaches infinity as the bifurcation is approached (i.e., as ). Beyond the bifurcation point the limit cycle no longer eists. We will discuss this more in the net chapter. Both the SNP and homoclinic bifurcations are known as global bifurcations, since they cannot be studied through a local analysis of a steady state (as is done with the local bifurcations such as the Hopf or the stationary bifurcations). To study their stability one must form Poincaré sections and eamine the eigenvalues of the linearization of the Poincaré map. These are called Floquet multipliers. An SNP bifurcation occurs when of these multipliers crosses through the unit comple circle at 1.