# Bifurcation

A **bifurcation** of a dynamical system is a qualitative change in its dynamics produced
by varying parameters.

## Definition

Consider an autonomous system of ordinary differential equations (ODEs) \[\tag{1} \dot{x}=f(x,\lambda),\ \ \ x \in {\mathbb R}^n, \ \ \ \lambda \in {\mathbb R}^p \]

where \(f\)
is smooth.
A **bifurcation** occurs at parameter \(\lambda = \lambda_0\) if
there are parameter values \(\lambda_1\) arbitrarily close to \(\lambda_0\) with dynamics
topologically inequivalent from those at \(\lambda_0\ .\) For example, the number
or stability of equilibria or periodic orbits of \(f\)
may change with perturbations of \(\lambda\)
from \(\lambda_0\ .\) One goal of bifurcation theory is to produce *parameter space maps*
or *bifurcation diagrams*
that divide the \(\lambda\) parameter space into regions of topologically equivalent systems.
Bifurcations occur at points that do not lie in the interior of one of these regions.

## Bifurcation theory

Bifurcation theory provides a strategy for investigating the bifurcations that occur within a
family. It does so by identifying ubiquitous patterns of bifurcations. Each *bifurcation type*
or *singularity*
is given a name; for example, Andronov-Hopf bifurcation. No distinction
has been made in the literature between "bifurcation" and "bifurcation type," both being called
"bifurcations."

Associated with each bifurcation type are

**defining equations**that locate bifurcations of that type in a family \(\dot{x} = f(x,\lambda)\)- normal forms that give model systems exemplifying the bifurcation type

Inequalities called **non-degeneracy conditions** are part of the specification of a bifurcation type.
The bifurcation types and their normal forms serve as templates that facilitate construction of parameter
space maps.
Bifurcation theory analyzes the bifurcations within the normal forms and investigates the similarity of
the dynamics within systems having a given bifurcation type. The "gold standard" for
similarity of systems used by the theory is topological equivalence. In some cases, bifurcation
theory proves structural stability of a family.
One of the principal objectives of bifurcation theory is to prove the structural
stability of normal forms. Note, however, that there are bifurcation types for which structurally stable normal forms do not
exist. An important aspect of the definition of structural stability in the context of bifurcation theory
is the specification of which
perturbations of a family are allowed. For example, bifurcation types of systems possessing specified
symmetries have been studied extensively (Equivariant Bifurcation Theory).

## Classification of bifurcations

One can view bifurcations as a failure of structural stability within a family. A starting point for classifying bifurcation types is the Kupka-Smale theorem that lists three generic properties of vector fields:

- hyperbolic equilibrium points
- hyperbolic periodic orbits
- transversal intersections of stable and unstable manifolds of equilibrium points and periodic orbits.

Different ways that these Kupka-Smale conditions fail lead to different bifurcation types. Bifurcation
theory constructs a layered graph of bifurcation types in which successive layers consist of types whose
defining equations specify more failure modes. These layers can be organized by the **codimension** of the
bifurcation types, defined as the minimal number of parameters of families in which that bifurcation type
occurs. Equivalently, the codimension is the number of equality conditions that characterize a bifurcation.

Codimension one bifurcations comprise the top level of bifurcation types. Single failures of the Kupka-Smale properties yield the following types of codimension one bifurcations:

- Equilibria
- Periodic Orbits
- Fold Limit Cycle Bifurcation
- Flip Bifurcation (aka Period Doubling bifurcation)
- Neimark-Sacker Bifurcation (aka Torus bifurcation)

- Global Bifurcations
- Homoclinic Bifurcation of equilibria
- Homoclinic tangencies of stable and unstable manifolds of periodic orbits
- Heteroclinic Bifurcation of equilibria and periodic orbits

This is not a comprehensive list of codimension one bifurcations. Additional types can be found in systems with quasiperiodic oscillations or chaotic dynamics. Moreover, there are subcases in the list above that deal with such issues as whether an Andronov-Hopf bifurcation is sub-critical or super-critical, and the implications of eigenvalue magnitudes for homoclinic bifurcation.

The classification of bifurcation types becomes more complex as their codimension increases. There are five types of "local" codimension two bifurcations of equilibria:

- Bautin Bifurcation
- Bogdanov-Takens Bifurcation
- Cusp Bifurcation
- Fold-Hopf Bifurcation
- Hopf-Hopf Bifurcation

## Numerical Methods

One of the principal uses of bifurcation theory is to analyze the bifurcations that occur in specific families of dynamical systems. Investigations commonly identify the types of bifurcations in parameter space maps either by comparison of simulation results with normal forms or by solving defining equations for those bifurcation types in the systems under investigation and computing coefficients of the normal forms. Several software packages (AUTO, CONTENT, MATCONT, XPPAUT, PyDSTool) give implementations of algorithms that perform the latter type of analysis. The numerical core of these packages consist of

*Regular*implementations of defining equations for the bifurcation types- equation solvers such as Newton's method
- Numerical continuation methods for differential equations
- Computation of normal forms
- initial and
- boundary value solvers for differential equations.

The continuation methods compute curves of solutions to regular systems of \( N \) equations in \( N+1 \) variables. The bifurcation analysis of a system implemented to varying degrees in the packages listed above is based upon the following strategy:

- An initial equilibrium or periodic orbit is located.
- Numerical continuation is used to follow this special orbit as a single
*active*parameter varies. - Defining equations for codimension one bifurcations detect and locate bifurcations that occur on this branch of solutions.
- Starting at one of the located codimension one bifurcations,

two parameters are designated to be active and the continuation methods are used to compute a curve of codimension one bifurcations.

- Defining equations for codimension two bifurcations detect and locate bifurcations that occur on this branch of solutions.
- Starting at one of the located codimension two bifurcations,

three parameters are designated to be active and the continuation methods are used to compute a curve of codimension two bifurcations.

This process can be continued as long as one has regular defining equations for bifurcations of increasing codimension, but these hardly exist beyond codimension three. Moreover, the dynamic behaviour near bifurcations with codimension higher than three is usually so poorly understood that the computation of such points is hardly worthwhile. In many cases, bifurcation analysis identifies additional curves of codimension k bifurcations that meet at a codimension k+1 bifurcation. Continuation methods can be started at one of these codimension k bifurcations to find curves of this type of bifurcation with k+1 active parameters. Switching to the continuation of a periodic orbit at an Andronov-Hopf bifurcation or to the continuation of a saddle homoclinic bifurcation curve from the Bogdanov-Takens bifurcation are examples of such starting techniques based on normal form computations.

## Bifurcation Theory of Chaotic and Quasiperiodic Systems

Bifurcation theory has intensively investigated varied topics that bear on chaotic and quasiperiodic dynamics. Much of this theory has been developed in the context of discrete time dynamical systems defined by iteration of mappings. The bifurcation theory described above has analogous results for this setting. In some areas, bifurcation theory of discrete systems goes farther than that for continuous time systems. In particular, an extensive, deep theory describing the properties of iterations of one dimensional mappings was developed over the last quarter of the twentieth century. This theory characterizes universal sequences of bifurcations and the existence of chaotic attractors. Some of this theory carries over to the setting of invertible mappings in higher dimensions and to continuous time dynamical systems via Poincar\'e maps. There are also results that are specific to continuous time systems, especially those that apply to homoclinic orbits of equilibrium points. Early results in this area include the theory of the Lorenz Attractor and Silnikov's analysis of systems with a homoclinic orbit of a saddle-focus in three dimensional systems. Methods originating in KAM (Kolmogorov-Arnold-Moser) theory describe how quasiperiodic invariant sets arise naturally in families of vector fields. Sophisticated numerical methods have been developed based upon this theory to compute invariant tori with (quasi)periodic motion in families of vector fields.

## References

- W. De Melo and S. Van Strien (1993) One Dimensional Dynamics, Springer.
- J. Guckenheimer and P. Holmes (1983) Nonlinear Oscillations, Dynamical systems and Bifurcations of Vector Fields. Springer
- Yu.A. Kuznetsov (2004) Elements of Applied Bifurcation Theory, Springer, 3rd edition.

**Internal references**

- Yuri A. Kuznetsov (2006) Andronov-Hopf bifurcation. Scholarpedia, 1(10):1858.
- John W. Milnor (2006) Attractor. Scholarpedia, 1(11):1815.
- John Guckenheimer and Yuri A. Kuznetsov (2007) Bautin bifurcation. Scholarpedia, 2(5):1853.
- John Guckenheimer and Yuri A. Kuznetsov (2007) Bogdanov-Takens bifurcation. Scholarpedia, 2(1):1854.
- Yuri A. Kuznetsov (2007) Conjugate maps. Scholarpedia, 2(12):5420.
- John Guckenheimer and Yuri A. Kuznetsov (2007) Cusp bifurcation. Scholarpedia, 2(4):1852.
- James Meiss (2007) Dynamical systems. Scholarpedia, 2(2):1629.
- Eugene M. Izhikevich (2007) Equilibrium. Scholarpedia, 2(10):2014.
- Jeff Moehlis and Edgar Knobloch (2007) Equivariant bifurcation theory. Scholarpedia, 2(9):2511.
- John Guckenheimer and Yuri A. Kuznetsov (2007) Fold-Hopf bifurcation. Scholarpedia, 2(10):1855.
- Lawrence F. Shampine and Skip Thompson (2007) Initial value problems. Scholarpedia, 2(3):2861.
- Willy Govaerts, Yuri A. Kuznetsov, Bart Sautois (2006) MATCONT. Scholarpedia, 1(9):1375.
- James Murdock (2006) Normal forms. Scholarpedia, 1(10):1902.
- Kendall E. Atkinson (2007) Numerical analysis. Scholarpedia, 2(8):3163.
- Jeff Moehlis, Kresimir Josic, Eric T. Shea-Brown (2006) Periodic orbit. Scholarpedia, 1(7):1358.
- Anatoly M. Samoilenko (2007) Quasiperiodic oscillations. Scholarpedia, 2(5):1783.
- Yuri A. Kuznetsov (2006) Saddle-node bifurcation. Scholarpedia, 1(10):1859.
- Leonid Pavlovich Shilnikov and Andrey Shilnikov (2007) Shilnikov bifurcation. Scholarpedia, 2(8):1891.
- Philip Holmes and Eric T. Shea-Brown (2006) Stability. Scholarpedia, 1(10):1838.
- James Murdock (2006) Unfoldings. Scholarpedia, 1(12):1904.
- Bard Ermentrout (2007) XPPAUT. Scholarpedia, 2(1):1399.

## See Also

Catastrophe Theory, Dynamical Systems, Normal Forms, Phase Space, Singularity Theory, Structural Stability, Unfoldings