ODE-Project The Hopf Bifurcation

Section 5.4 The Hopf Bifurcation

Objectives

To understand that Hopf bifurcations; that is, bifurcations for systems parameterized by $α,$ occur if an attracting periodic solution encircling an equilibrium solution develops as $α$ passes through the bifurcation value.

If we consider the phase portrait of a system

\begin{aligned} \frac{d x}{d t} & = f (x, y) \\ \frac{d y}{d t} & = g (x, y), \end{aligned}

we may see closed orbits and various types of equilibrium solutions. If we make a slight change in the system, then we might reasonably expect the phase portrait to change slightly. For example, a closed orbit might expand or contract, an equilibrium solution might shift, or a spiral source might tighten up, but we would not expect new equilibrium solutions to appear or a spiral source to turn into a spiral sink. However, this may be exactly what happens. If such a change occurs in the system through an adjustment of a parameter, we say a bifurcation has occurred.

🔗

Subsection 5.4.1 Bifurcations

The simplest bifurcations occur when an equilibrium solution appears, disappears, or splits into two or more equilibrium solutions. For example, the system

🔗

\begin{aligned} \frac{d x}{d t} & = α + x^{2} \\ \frac{d y}{d t} & = - y \end{aligned}

has no equilibrium solution if

α > 0

(Figure 5.4.1). However, we have a single equilibrium solution at the origin if

α = 0

and two equilibrium solutions,

(\pm \sqrt{- α}, 0),

α < 0 .

🔗

Figure 5.4.1. Solution curves $α = - 5 .$
🔗

The Jacobian of the system

🔗

\begin{aligned} \frac{d x}{d t} & = α + x^{2} \\ \frac{d y}{d t} & = - y \end{aligned}

🔗

J (x, y) = (\begin{matrix} 2 x & 0 \\ 0 & - 1 \end{matrix}) .

α = 0,

then

J

has eigenvalues 0 and

- 1 .

This is not an elementary equilibrium point such a saddle, a sink, or a spiral source. We call this type of equilibrium solution a saddle-node. On the right half of the

x y

-plane, the equilibrium solution resembles a saddle, while on the left half it resembles a sink (Figure 5.4.2).

🔗

α < 0,

then

🔗

J = (\begin{matrix} 2 \sqrt{- α} & 0 \\ 0 & - 1 \end{matrix}) or (\begin{matrix} - 2 \sqrt{- α} & 0 \\ 0 & - 1 \end{matrix}) .

In the first case, we have a saddle since the eigenvalues are real and of opposite sign. In the second case, we have a nodal sink, since both eigenvalues are negative. The two equilibrium solutions move in opposite directions as

α

decreases. We can summarize what happens for various values of

α

with a bifurcation diagram such as the one in Figure 5.4.3.

🔗

🔗

Subsection 5.4.2 The Hopf Bifurcation

Now let us consider an entirely different type of bifurcation by examining the system

🔗

\begin{aligned} \frac{d x}{d t} & = α x + 5 y - x (x^{2} + y^{2}) \\ \frac{d y}{d t} & = - 5 x + α y - y (x^{2} + y^{2}) . \end{aligned}

The origin is an equilibrium solution for all values of

α .

The linearization of our system is

🔗

\begin{aligned} \frac{d x}{d t} & = α x + 5 y \\ \frac{d y}{d t} & = - 5 x + α y . \end{aligned}

Since the eigenvalues of this system are

α \pm 5 i,

we can easily determine the nature of of our linearization. If

α < 0,

we have a spiral sink. As

α

moves from negative values to positive values, the origin changes to a center (

α = 0

), and then to a spiral source (

α > 0

🔗

In the case of the nonlinear system, the origin is still a spiral sink for

α > 0

(Figure 5.4.4). If

α > 0,

the origin is still a spiral source (Figure 5.4.5). However, something quite different happens at

α = 0 .

The origin is stable equilibrium solution with solutions spiraling into the origin very slowly. As

α

increases past zero, the equilibrium solution destabilizes and becomes a source. In addition, a closed orbit of radius

\sqrt{α}

develops. Solutions inside of this closed orbit spiral out towards the orbit, while solutions outside of the orbit spiral inward (Figure 5.4.6).

🔗

To see exactly what happens as

α

passes zero and becomes positive, we will rewrite our system in polar coordinates. If we make the substitution

x = r \cos θ

and

y = r \sin θ,

our nonlinear system can be rewritten as

🔗

\begin{aligned} \frac{d r}{d t} & = α r - r^{3} \\ \frac{d θ}{d t} & = - 5. \end{aligned}

α < 0,

the origin is a sink since

α r - r^{3} < 0

for all

r > 0 .

In this case all solutions tend towards the origin as

t \to \infty .

When

α > 0,

the origin is still an equilibrium solution. Moreover,

r^{'} = 0,

when

r = \sqrt{α} .

We also know that

r^{'} > 0

for

0 < r < \sqrt{α}

and

r^{'} < 0

r > \sqrt{α} .

So the circle of radius

\sqrt{α}

is a periodic solution with the trajectory moving clockwise since

θ^{'} = - 5 < 0 .

All nonzero solutions spiral towards this orbit as

t \to \infty .

This type of bifurcation is called a Hopf bifurcation. No new equilibrium solutions arise, but a periodic solution develops as

α

passes through the bifurcation value.

🔗

Theorem 5.4.7. Hopf Bifurcation.

Suppose that

(x_{0} (α), y_{0} (α))

is a equilibrium solution for the family of systems

🔗

\begin{aligned} \frac{d x}{d t} & = f_{α} (x, y) \\ \frac{d y}{d t} & = g_{α} (x, y) \end{aligned}

parameterized by

α,

and the Jacobian matrix for the system evaluated at this equilibrium has eigenvalues

a (α) \pm i b (α) .

Assume that at some

α = α_{0},

we have

a (α_{0}) = 0,

a^{'} (α_{0}) > 0,

b (α_{0}) \neq 0,

and that the equilibrium point is asymptotically stable. Then there exists an

α_{1}

such that the system has a periodic solution encircling the equilibrium solution for

α_{0} < α < α_{1} .

For a proof and description of the Hopf Bifurcation Theorem see C. Chicone. Ordinary Differential Equations with Applications. Springer-Verlag, New York, 1999.

🔗

Recall that our example system

🔗

\begin{aligned} \frac{d x}{d t} & = α x + 5 y - x (x^{2} + y^{2}) \\ \frac{d y}{d t} & = - 5 x + α y - y (x^{2} + y^{2}) \end{aligned}

has an equilibrium at

(0, 0)

with eigenvalues

α \pm 5 i .

If we let

α_{0} = 0,

then

a (0) = 0,

a^{'} (0) = 1,

and

b (0) = 5 .

Thus, the hypothesis of the Hopf Bifurcation Theorem are satisfied and we are guaranteed a period solution surrounding our equilibrium solution.

🔗

Example 5.4.8. Van der Pol’s equation.

Van der Pol’s equation, a simple nonlinear equation having applications in electrical engineering and mathematical biology, is

x^{″} - x (1 - x^{2}) x^{'} + x = 0 .

This equation can be written as the system

🔗

\begin{aligned} (5.4.1) & \frac{d x}{d t} & = y \\ (5.4.2) & \frac{d y}{d t} & = - x + x (1 - x^{2}) y . \end{aligned}

The phase portrait for Van der Pol’s equation is given in Figure 5.4.9. The origin is the only equilibrium solution to Van der Pol’s equation, and one might guess that this solution acts like a spiral source by examining the phase portrait. If we examine the system consisting of only linear terms on the right-hand side of the equation,

🔗

\begin{aligned} \frac{d x}{d t} & = y \\ \frac{d y}{d t} & = - x + y \end{aligned}

we might get a better sense of what is happening. The matrix for this linear system,

🔗

A = (\begin{matrix} 0 & 1 \\ - 1 & 1 \end{matrix}),

has eigenvalues

λ = (1 \pm i \sqrt{3}) / 2 .

This suggests that the origin looks like a spiral source at least locally. Indeed, the Jacobian matrix for (5.4.1)–(5.4.1) is

🔗

J = (\begin{matrix} 0 & 1 \\ - 1 - 2 x_{0} y_{0} & 1 - x_{0}^{2} \end{matrix}),

which agrees with

A

for

(x_{0}, y_{0}) = (0, 0) .

🔗

🔗

Let us examine how a bifurcation might occur in a predator-prey model. In our model, we will assume that a predator’s appetite is satiated when food is abundant. If this is the case, an increase in the prey population has little effect on the interaction terms in our model. We might model this satiable predator-prey scenario by the equations

🔗

\begin{aligned} \frac{d x}{d t} & = (a - b x) x - \frac{c x y}{d + α x} \\ \frac{d y}{d t} & = - e y + \frac{f x y}{d + α x}, \end{aligned}

where

x (t)

is the prey population at time

t

and

y (t)

is the predator population at time

t .

Here,

a, \dots, f

and

α

are constants. Let us choose the constants

a, \dots, f

so that our system becomes

🔗

\begin{aligned} \frac{d x}{d t} & = (1 - x) x - \frac{x y}{0.3 + α x} \\ \frac{d y}{d t} & = - 0.5 y + \frac{x y}{0.3 + α x} . \end{aligned}

The constant

α

is know as the satiation constant. The larger the value of

α,

the more quickly a predator’s appetite is satiated as the prey population increases.

🔗

Let us see examine the phase portrait for several values of

α .

α = 0,

then there is no satiation effect. In this case, all trajectories inside the population quadrant spiral asymptotically towards a single stable equilibrium point (Figure 5.4.10). If

α = 1.35,

we have a very similar phase portrait with solutions spiraling in towards a stable equilibrium solution (Figure 5.4.11). However, if

α = 0.9,

we have a periodic trajectory. All other trajectories are pulled toward this attracting periodic orbit (Figure 5.4.12).

🔗

Figure 5.4.10. Predator-prey model with $α = 0$
🔗

Figure 5.4.11. Predator-prey model with $α = 1.35$
🔗

Figure 5.4.12. Predator-prey model with $α = 0.9$
🔗

Diagram Exploration Keyboard Controls

Key	Action
Enter, A	Activate keyboard driven exploration
B	Activate menu driven exploration
Escape	Leave exploration mode
Cursor down	Explore next lower level
Cursor up	Explore next upper level
Cursor right	Explore next element on level
Cursor left	Explore previous element on level
X	Toggle expert mode
W	Extra details if available
Space	Repeat speech
M	Activate step magnification
Comma	Activate direct magnification
N	Deactivate magnification
Z	Toggle subtitles
C	Cycle contrast settings
T	Monochrome colours
L	Toggle language (if available)
K	Kill current sound
Y	Stop sound output
O	Start and stop sonification
P	Repeat sonification output

Figure 5.4.13. Bifurcation diagram for the satiable predator-prey model

🔗

α

continues to increase a Hopf bifurcation occurs at

α \approx 0.5 .

At this point the equilibrium solution destabilizes and spawns an attracting periodic orbit. As

α

continues to increase, the amplitude of this periodic orbit increases. However, at

α \approx 0.85,

the amplitude of the periodic solution begins to decrease. At

α \approx 1.2,

the equilibrium point re-stabilizes and the periodic solution is absorbed. For

α > 1.2,

we only have a stable equilibrium solution in the first quadrant. The

x

-coordinate of our equilibrium solution is

🔗

x = \frac{0.15}{1 - 0.5 α},

and we can summarize our findings with a bifurcation diagram (Figure 5.4.13).

🔗

Subsection 5.4.3 Important Lessons

If we make a slight change in the system, then we might reasonably expect the phase portrait to change slightly. For example, a closed orbit might expand or contract, an equilibrium solution might shift, or a spiral source might tighten up, but we would not expect new equilibrium solutions to appear or a spiral source to turn into a spiral sink. However, this may be exactly what happens. If such a change occurs in the system through an adjustment of a parameter, we say a bifurcation has occurred.
A Hopf bifurcation if an attracting periodic solution encircling an equilibrium solution develops as $α$ passes through the bifurcation value.
Suppose that $(x_{0} (α), y_{0} (α))$ is a equilibrium solution for the family of systems
$\begin{aligned} \frac{d x}{d t} & = f_{α} (x, y) \\ \frac{d y}{d t} & = g_{α} (x, y) \end{aligned}$
parameterized by $α,$ and the Jacobian matrix for the system evaluated at this equilibrium has eigenvalues $a (α) \pm i b (α) .$ Assume that at some $α = α_{0},$ we have $a (α_{0}) = 0,$ $a^{'} (α_{0}) > 0,$ $b (α_{0}) \neq 0,$ and that the equilibrium point is asymptotically stable. Then there exists an $α_{1}$ such that the system has a periodic solution encircling the equilibrium solution for $α_{0} < α < α_{1} .$
A satiable predator-prey scenario can be modeled by the equations
$\begin{aligned} \frac{d x}{d t} & = (a - b x) x - \frac{c x y}{d + α y} \\ \frac{d y}{d t} & = - e y + \frac{f x y}{d + α y}, \end{aligned}$
where $x (t)$ is the prey population at time $t$ and $y (t)$ is the predator population at time $t .$ Here, $a, \dots, f$ and $α$ are constants. The constant $α$ is know as the satiation constant. The larger the value of $α,$ the more quickly a predator’s appetite is satiated as the prey population increases.

🔗

Reading Questions 5.4.4 Reading Questions

1.

Explain what a Hopf bifurcation is.

🔗

2.

What scenerio does a satiated predator-prey model describe? Explain.

🔗

Exercises 5.4.5 Exercises

1.

Consider the system

🔗

\begin{aligned} x^{'} & = y \\ y^{'} & = - x - \frac{y}{4} + x^{2} . \end{aligned}

Verify that the function
$L (x, y) = \frac{x^{2}}{2} + \frac{y^{2}}{2} - \frac{x^{3}}{3}$
is a Lyapunov function for the system.
Sketch the level sets of $L .$
What can you conclude about the phase portrait of the system from (1) and (2)?

🔗

2.

Let

G (x, y) = x^{3} - 3 x y^{2} .

🔗

What is the gradient system with vector field given by the gradient of $G ?$
Sketch the graph of $G$ and the level sets of $G .$
Sketch the phase portrait of the gradient system in (1).

🔗

3.

Suppose that the smell of a bunch of dead fish in the region

- 2 \leq x \leq 2,

- 2 \leq y \leq 2

is given by the function

🔗

S (x, y) = x^{2} + y^{2} - \frac{x^{4} + y^{4}}{4} - 3 x^{2} y^{2} + 100.

What is the gradient system who vector field is the gradient of $S ?$
Use Sage to graph the phase portrait of the system.
How many dead fish are there and where are they?
Using the results from part (2), sketch the level sets of $S .$
Why is this model not realistic for large values of $x$ or $y ?$

🔗

4.

Consider the system of differential equations

🔗

\begin{aligned} x^{'} & = x (- x - y + 1) \\ y^{'} & = y (- a x - y + b), \end{aligned}

where

a

and

b

are parameters with

a, b > 0 .

Suppose that this system is only defined for

x, y \geq 0 .

🔗

Use nullclines to sketch the phase portrait for this system for several values of $a$ and $b .$
Determine the values of $a$ and $b$ at which a bifurcation occurs.