# Stability and Energy-Casimir Mapping for Integrable Deformations of the Kermack-McKendrick System

Authors: Cristian Lăzureanu; Camelia Petrişor
Journal: Advances in Mathematical Physics
 (2018)

---

## Abstract

Integrable deformations of a Hamilton-Poisson system can be obtained altering its constants of motion. These deformations are integrable systems that can have various dynamical properties. In this paper, we give integrable deformations of the Kermack-McKendrick model for epidemics, and we analyze a particular integrable deformation. More precisely, we point out two Poisson structures that lead to infinitely many Hamilton-Poisson realizations of the considered system. Furthermore, we study the stability of the equilibrium points, we give the image of the energy-Casimir mapping, and we point out some of its properties.


---

## Body

## 1. Introduction

Evolution equations represent models for describing phenomena that appear in physics, biology, chemistry, economy, and engineering. In many situations these evolution equations can be analyzed in the frame of the Lagrangian mechanics or the Hamiltonian mechanics. Furthermore, there are phenomena that are modeled by three-dimensional systems of differential equations, particularly Hamilton-Poisson systems. Such systems can be perturbed in order to obtain a desired behavior. A way to perturb a three-dimensional Hamilton-Poisson systems consists in alteration of its constants of motion. This method leads to integrable deformations of the initial system.In recent papers, integrable deformations of some particular Hamilton-Poisson systems were analyzed. In [1], observing that the constants of motion of the Euler top determine its equations, integrable deformations of the Euler top were given. In [2], integrable deformations of the three-dimensional real valued Maxwell-Bloch equations were obtained by altering the constants of motion of the considered system. In the same manner, in [3], integrable deformations of the Rikitake system were constructed. These integrable deformations can be viewed as controlled systems and, in consequence, a study of modifications in their dynamics can be performed. Moreover, the integrable deformations of the above systems are also Hamilton-Poisson systems. Consequently, they can be analyzed from some standard and nonstandard Poisson geometry points of view [4].The study of a three-dimensional Hamilton-Poisson system from some standard and nonstandard Poisson geometry points of view tries to answer the following open problem formulated by Tudoran et al. [4]: “Is there any connection between the dynamical properties of a given dynamical system and the geometry of the image of the energy-Casimir mapping, and if yes, how can one detect as many as possible dynamical elements (e.g., equilibria, periodic orbits, homoclinic and heteroclinic connections) and dynamical behavior (e.g., stability, bifurcation phenomena for equilibria, periodic orbits, homoclinic and heteroclinic connections) by just looking at the image of this mapping?” Affirmative answers were given for some particular systems [5–9]. In these cases the image of the energy-Casimir mapping E
C
=
(
H
,
C
), where H is the Hamiltonian and C is a Casimir function, is a closed subset of R
2, namely, the convex hull of the images of the stable equilibrium points through E
C. Furthermore, the images of the equilibrium points through the energy-Casimir mapping give an algebraic partition of the set I
m
(
E
C
), and the orbits of these systems are bounded. On the other hand, the image of the energy-Casimir mapping can be R
2 ([10, 11]), and other connections were observed (for example, there are unbounded orbits). In [12], taking into account these facts, some questions regarding the connections between the dynamics of a Hamilton-Poisson systems and the associated energy-Casimir mapping were asked. We recall some of them: “are the observed properties in every case when I
m
(
E
C
)
⊊
R
2 true?” or “can I
m
(
E
C
) be a nonconvex set? If yes, do the observed properties remain true?” One of the goals of this paper is to give some answers to these questions.The finding of some counterexamples assumes the study of many Hamilton-Poisson systems, and how we have already seen such systems can be obtained using integrable deformations of known integrable systems. Moreover, because the constants of motionH and C of the above-mentioned systems are polynomials, it is a good idea to analyze systems that have nonpolynomials constants of motion. Such a system is the well-known system introduced in 1927 by Kermack and McKendrick [13] and brought back in attention by Anderson and May, in 1979 [14]. The Kermack-McKendrick system and its generalizations were widely investigated. We mention a very short list of works [15–17]. We also notice the applications of such type of systems in health, networks, informatics, economics, and finance (see, for example, [18] and references therein).The paper is organized as follows. In Section2, we recall the Kermack-McKendrick model and we give integrable deformations of this system. In Section 3, we analyze a particular integrable deformation of the Kermack-McKendrick system. More precisely, we point out two Poisson structures and, in consequence, we obtain two Hamilton-Poisson realizations of the considered system. In addition, using these structures, we construct infinitely many Hamilton-Poisson realizations of our system. Furthermore, we study the stability of the equilibrium points, we give the image of the energy-Casimir mapping, and we point out some of its properties. The conclusions are presented in the last section.


## 2. Integrable Deformations of the Kermack-McKendrick System

Following [1], in this section, we give integrable deformations of the Kermack-McKendrick system. First, we recall the epidemic model introduced by Kermack and McKendrick [13] (for details, see also [19]).In the mathematical theory of epidemics, a basic model is given by the Kermack-McKendrick system. This model intends to describe the spread of the infection within the population as a function of time. It is considered that the total population is constant, and it is divided into three distinct groups. First group is formed by individuals who can catch the disease, named the susceptibles. At a momentt their number is S
(
t
). The second group, the infected population, consists in individuals who have the disease and can transmit it. Their number is I
(
t
). Finally, the group of the removed subjects, in number of R
(
t
), formed by those who had the disease, cannot be reinfected and cannot infect other individuals. In order to obtain the evolution equations, some assumptions were made. Firstly, the gain in the infective group is at a rate proportional to the number of infected subjects and susceptibles, that is, a
S
(
t
)
I
(
t
), where a
>
0 is the infection rate. The susceptibles are lost at the same rate. Furthermore, the rate of removal of the infected subjects to the removed group is proportional to the number of the infected subjects, that is, b
I
(
t
), where b
>
0 is the removal rate of the infected subjects. In addition, the incubation period is negligible, and every pair of individuals has equal probability of coming into contact with one another. Therefore, the following equations were deduced:(1)
d
S
d
t
=
-
a
S
I
d
I
d
t
=
a
S
I
-
b
I
d
R
d
t
=
b
I
.We denote S
≔
x
,
I
≔
y
,
R
≔
z. Then the Kermack-McKendrick system is written as follows:(2)
x
˙
=
-
a
x
y
y
˙
=
a
x
y
-
b
y
z
˙
=
b
y
,where a
,
b are positive constants.It is obvious that a constant of motion is given by the total number of individuals; namely,(3)
I
1
x
,
y
,
z
=
x
+
y
+
z
.We recall the second constant of motion(4)
I
2
x
,
y
,
z
=
ln
⁡
x
+
a
b
z
.Differentiating the above constants of motion we obtain (5)
x
˙
+
y
˙
=
-
z
˙
,
x
˙
x
=
-
a
b
z
˙
,and considering z
˙
=
b
y, we get system (2). Therefore the constants of motion (3), (4) generate system (2). This property allows us to obtain integrable deformations of system (2) by alteration of its constants of motion [1].Consider as constants of motion the functions(6)
I
~
1
x
,
y
,
z
=
x
+
y
+
z
+
α
x
,
y
,
z
,
I
~
2
x
,
y
,
z
=
ln
⁡
x
+
a
b
z
+
β
x
,
y
,
z
,where α
,
β are arbitrary differentiable functions. As above, we obtain that these functions generate the following system:(7)
x
˙
=
-
a
x
y
+
b
x
y
∂
y
β
-
∂
z
β
-
a
x
y
∂
y
α
-
b
x
y
∂
y
α
∂
z
β
-
∂
z
α
∂
y
β
y
˙
=
a
x
y
-
b
y
+
a
x
y
∂
x
α
-
b
y
∂
z
α
+
b
x
y
∂
z
β
-
∂
x
β
+
b
x
y
∂
x
α
∂
z
β
-
∂
z
α
∂
x
β
z
˙
=
b
y
+
b
y
∂
y
α
-
b
x
y
∂
y
β
-
∂
x
β
-
b
x
y
∂
x
α
∂
y
β
-
∂
y
α
∂
x
β
,where ∂
x
f
≔
∂
f
/
∂
x. If α and β are constant functions, then (7) reduces to (2). Therefore, for any differential functions α and β, system (7) is an integrable deformation of the Kermack-McKendrick system.Remark 1.
In order to maintain constant the total population, the functionα vanishes. In this case system (7) becomes(8)
x
˙
=
-
a
x
y
+
b
x
y
∂
y
β
-
∂
z
β
y
˙
=
a
x
y
-
b
y
+
b
x
y
∂
z
β
-
∂
x
β
z
˙
=
b
y
-
b
x
y
∂
y
β
-
∂
x
β
.


## 3. Dynamical Properties of an Integrable Deformation of the Kermack-McKendrick System

In this section, we consider some particular deformation functions, and we give some dynamical properties of the corresponding integrable deformation of the Kermack-McKendrick system. First, we give Hamilton-Poisson realizations of this system that provides the geometric framework of our study. Furthermore, we study the stability of the equilibrium points. We also give some properties of the energy-Casimir mapping associated with the considered system.We consider the following deformation functions:(9)
α
x
,
y
,
z
=
g
2
y
2
+
g
2
z
2
,
β
x
,
y
,
z
=
0
,where g
∈
R is the deformation parameter. Then system (7) becomes(10)
x
˙
=
-
a
x
y
-
a
g
x
y
2
y
˙
=
a
x
y
-
b
y
-
b
g
y
z
z
˙
=
b
y
+
b
g
y
2
,and its constants of motion are given by (6); namely,(11)
C
1
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
,
I
2
x
,
y
,
z
=
ln
⁡
x
+
a
b
z
.In what follows we need constants of motion defined on R
3. We immediately get that the function C
2,(12)
C
2
x
,
y
,
z
=
x
·
e
a
/
b
z
,is a constant of motion of system (10).


### 3.1. Hamilton-Poisson Realizations

We recall that the dynamical system generated by theC
1 vector field f
=
(
f
1
,
f
2
,
f
3
) on a C
∞ manifold P
⊆
R
3,(13)
x
˙
=
f
x
,
x
=
x
1
,
x
2
,
x
3
,has the Hamilton-Poisson realization (
P
,
{
·
,
·
}
,
H
), if it can be put in the form(14)
x
˙
i
=
x
i
,
H
,
i
∈
1,2
,
3
,where H is the Hamiltonian function, and {
·
,
·
} is a Poisson bracket on P.Considering a smooth functionC on P, a Poisson structure on P is generated by the Poisson bracket(15)
f
,
g
C
≔
∂
C
,
f
,
g
∂
x
1
,
x
2
,
x
3
=
∂
x
1
C
∂
x
2
C
∂
x
3
C
∂
x
1
f
∂
x
2
f
∂
x
3
f
∂
x
1
g
∂
x
2
g
∂
x
3
g
,for every f
,
g
∈
C
∞
(
R
3
,
R
). In addition, see, for example, [20], using any smooth function ν on P, a new Poisson bracket is given by(16)
f
,
g
C
ν
≔
ν
·
f
,
g
C
.In both cases, the function C is a Casimir of the Poisson structure; that is, {
C
,
f
}
=
0 for every f. In coordinates, the Poisson structure is given in matrix notation(17)
Π
C
ν
=
ν
Π
C
=
ν
0
∂
z
C
-
∂
y
C
-
∂
z
C
0
∂
x
C
∂
y
C
-
∂
x
C
0
.Moreover, if there is a smooth function H such that system (13) takes the form x
˙
t
=
Π
C
ν
·
∇
H, then (13) is a Hamilton-Poisson system.In our case, letg
≠
0 and let(18)
C
x
,
y
,
z
=
C
1
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
.We obtain the rescaling function ν
(
x
,
y
,
z
)
=
b
y
·
e
-
(
a
/
b
)
z and the Poisson structure generated by C
1:(19)
Π
C
1
=
0
g
z
+
1
-
g
y
-
1
-
g
z
-
1
0
1
g
y
+
1
-
1
0
=
0
g
z
-
g
y
-
g
z
0
0
g
y
0
0
+
0
1
-
1
-
1
0
1
1
-
1
0
.Due to its linearity, the above Poisson structure is a Lie-Poisson structure on the dual vector space of a Lie algebra, namely, s
e
(
2
)
∗. Indeed, consider the special Euclidean Lie group S
E
(
2
) of all orientation-preserving isometries (see, for example, [21]), given by(20)
S
E
2
=
A
∈
G
L
3
,
R
∣
A
=
1
0
0
u
cos
⁡
φ
-
sin
⁡
φ
v
sin
⁡
φ
cos
⁡
φ
,
u
,
v
∈
R
,
φ
∈
0,2
π
.The corresponding Lie algebra of S
E
(
2
) is(21)
s
e
2
=
X
∈
g
l
3
,
R
∣
X
=
0
0
0
β
0
-
α
γ
α
0
,
α
,
β
,
γ
∈
R
,with the commutator bracket [
X
,
Y
]
=
X
Y
-
Y
X
.As vector space,s
e
(
2
) has the basis B
s
e
(
2
)
=
{
e
1
,
e
2
,
e
3
}, where(22)
e
1
=
0
0
0
0
0
-
g
0
g
0
,
e
2
=
0
0
0
1
0
0
0
0
0
,
e
3
=
0
0
0
0
0
0
1
0
0
.We obtain the following bracket relations: (23)
e
1
,
e
2
=
g
e
3
,
e
1
,
e
3
=
-
g
e
2
,
e
2
,
e
3
=
0
.We consider the bilinear form Θ
:
s
e
(
2
)
×
s
e
(
2
)
→
R given by the matrix (
Θ
i
j
)
1
≤
i
,
j
≤
3, Θ
12
=
-
Θ
21
=
1, Θ
13
=
-
Θ
31
=
-
1, and Θ
23
=
-
Θ
32
=
1. By straightforward computations we obtain that Θ satisfies condition (24)
Θ
X
,
Y
,
Z
+
Θ
Y
,
Z
,
X
+
Θ
Z
,
X
,
Y
=
0
,for every triplet (
X
,
Y
,
Z
) of elements in s
e
(
2
)
. Therefore Θ is a symplectic cocycle of the Lie algebra s
e
(
2
). Moreover, it is not a coboundary since Θ
(
e
2
,
e
3
)
=
1
≠
0
=
f
(
[
e
2
,
e
3
]
), for every linear map f
,
f
:
s
e
(
2
)
→
R
. Following [22], the modified Lie-Poisson structure Π
C
1 is defined on the dual space s
e
(
2
)
∗
≃
R
3.The following result gives a Hamilton-Poisson realization of the considered system.Proposition 2.
Letν be the rescaling function given by (25)
ν
x
,
y
,
z
=
b
y
·
e
-
a
/
b
z
.If Π
C
1 is the Poisson structure generated by the Casimir function C
1, (26)
C
1
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
,given by (27)
x
,
y
C
1
=
g
z
+
1
,
x
,
z
C
1
=
-
g
y
-
1
,
y
,
z
C
1
=
1
,then system (10) has the Hamilton-Poisson realization (28)
R
3
,
·
,
·
C
1
ν
,
C
2
,where {
f
,
g
}
C
1
ν
=
ν
·
{
f
,
g
}
C
1 for any f
,
g
∈
C
∞
(
R
3
,
R
), and C
2 is the Hamiltonian function given by (29)
C
2
x
,
y
,
z
=
x
·
e
a
/
b
z
.Proof.
An easy computation shows that system (10) takes the form Π
C
1
ν
·
∇
C
2
=
(
x
˙
,
y
˙
,
z
˙
)
t; hence the conclusion follows.Using the same notations as in Proposition2, we similarly obtain the following results.Proposition 3.
Letμ be the rescaling function given by μ
(
x
,
y
,
z
)
=
-
ν
(
x
,
y
,
z
) and let C
2 be a Casimir function. If Π
C
2
μ is the Poisson structure generated by C
2 and μ, given in matrix notation by (30)
Π
C
2
μ
=
0
-
a
x
y
0
a
x
y
0
-
b
y
0
b
y
0
,then system (10) has the Hamilton-Poisson realization (31)
R
3
,
·
,
·
C
2
μ
,
C
1
,where C
1 is the Hamiltonian function.Remark 4.
The above Poisson structures are compatible and(32)
Π
C
1
ν
∇
C
2
=
Π
C
2
μ
∇
C
1
,and hence (10) is a bi-Hamiltonian system. Moreover, this pair of Hamilton-Poisson realizations gives rise to infinitely many Hamilton-Poisson realizations of system (10) (see Proposition 5).Proposition 5.
Letp
,
q
,
r
,
s
∈
R such that p
s
-
q
r
=
1, and let ν
(
x
,
y
,
z
)
=
b
y
·
e
-
(
a
/
b
)
z be the rescaling function. There exists infinitely many Hamilton-Poisson realizations of system (10) given by (
R
3
,
Π
p
,
q
,
H
r
,
s
), where the Hamiltonian H
r
,
s is given by(33)
H
r
,
s
=
r
C
1
+
s
C
2
=
r
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
+
s
x
·
e
a
/
b
z
,the Poisson structure is defined by(34)
Π
p
,
q
=
p
Π
C
1
ν
+
q
Π
C
2
ν
,and a Casimir function of the above Poisson structure is given by(35)
C
p
,
q
=
p
C
1
+
q
C
2
=
p
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
+
q
x
·
e
a
/
b
z
.Proof.
It is clear thatΠ
p
,
q
=
ν
(
p
Π
C
1
+
q
Π
C
2
). Therefore Π
p
,
q is the Poisson structure generated by the Casimir function C
p
,
q and the rescaling function ν. Using the Hamilton-Poisson realizations given in Propositions 2 and 3, we immediately obtain Π
p
,
q
·
∇
H
r
,
s
=
(
x
˙
,
y
˙
,
z
˙
)
t, which finishes the proof.


### 3.2. Stability of the Equilibrium Points

The equilibrium points of system (10) are given by the following families:(36)
E
1
=
M
,
-
1
g
,
a
M
-
b
b
g
:
M
∈
R
,
E
2
=
N
,
0
,
P
:
N
,
P
∈
R
.We remark that the second family represents the set of all equilibrium points of Kermack-McKendrick system (2). Therefore, the existence of another family of equilibrium points produces changes in the dynamics of initial system (2). We are concerned with the study of these new equilibrium points.Proposition 6.
Letg
∈
R
∖
{
0
} and let e
1
M
=
M
,
-
1
/
g
,
(
a
M
-
b
)
/
b
g
∈
E
1 be an equilibrium point of system (10).
(i) Ifg
(
a
2
M
+
b
2
g
)
<
0, then e
1
M is an unstable equilibrium point.
(ii) Ifg
(
a
2
M
+
b
2
g
)
>
0, then e
1
M is a nonlinearly stable equilibrium point.Proof.
(i) LetA be the matrix of the linear part of our system; that is, (37)
A
=
-
a
y
1
+
g
y
-
a
x
-
2
a
g
x
0
a
y
a
x
-
b
-
b
g
z
-
b
g
y
0
b
+
2
b
g
y
0
.The characteristic roots of A
(
e
1
M
) are given by (38)
λ
1
=
0
,
λ
2,3
=
±
-
b
2
-
a
2
M
g
.Considering M such that g
(
a
2
M
+
b
2
g
)
<
0, we conclude that the equilibrium point e
1
M is unstable.
(ii) Now, letM be such that g
(
a
2
M
+
b
2
g
)
>
0. We use the energy-Casimir method [23]. Let H
φ be the energy-Casimir function: (39)
H
φ
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
+
φ
x
·
e
a
/
b
z
,where φ
:
R
→
R is a smooth real valued function defined on R.
The first variation ofH
φ is given by (40)
δ
H
φ
x
,
y
,
z
=
1
+
e
a
/
b
z
·
φ
˙
x
·
e
a
/
b
z
δ
x
+
1
+
g
y
δ
y
+
1
+
g
z
+
a
x
b
·
e
a
/
b
z
·
φ
˙
x
·
e
a
/
b
z
δ
z
,where φ
˙
(
u
)
≔
d
φ
/
d
u. We have (41)
δ
H
φ
e
1
M
=
1
+
e
a
a
M
-
b
/
b
2
g
·
φ
˙
M
·
e
a
a
M
-
b
/
b
2
g
δ
x
+
a
M
b
δ
z
,which vanishes if and only if(42)
φ
˙
M
·
e
a
a
M
-
b
/
b
2
g
=
-
e
-
a
a
M
-
b
/
b
2
g
.The second variation of H
φ is given by (43)
δ
2
H
φ
x
,
y
,
z
=
e
2
a
/
b
z
·
φ
¨
x
·
e
a
/
b
z
δ
x
2
+
g
δ
y
2
+
x
·
a
2
b
2
·
e
a
/
b
z
g
+
φ
˙
x
·
e
a
/
b
z
+
x
·
e
a
/
b
z
φ
¨
x
·
e
a
/
b
z
δ
z
2
+
2
a
b
·
e
a
/
b
z
φ
˙
x
·
e
a
/
b
z
+
x
·
e
a
/
b
z
φ
¨
x
·
e
a
/
b
z
δ
x
δ
z
.Taking into account relation (42), we obtain (44)
δ
2
H
φ
e
1
M
=
e
2
a
a
M
-
b
/
b
2
g
·
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
δ
x
2
+
g
δ
y
2
+
g
-
M
·
a
2
b
2
+
M
2
·
a
2
b
2
·
e
2
a
a
M
-
b
/
b
2
g
·
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
δ
z
2
+
2
a
b
M
·
e
2
a
a
M
-
b
/
b
2
g
·
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
-
1
δ
x
δ
z
.If g
>
0, then we choose a function φ such that relation (42) holds and (45)
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
>
a
2
b
2
g
+
a
2
M
·
e
-
2
a
a
M
-
b
/
b
2
g
>
0
.For example, let (46)
φ
u
=
1
2
g
+
a
2
b
2
g
+
a
2
M
·
e
-
2
a
a
M
-
b
/
b
2
g
u
-
M
·
e
a
a
M
-
b
/
b
2
g
2
-
e
-
a
a
M
-
b
/
b
2
g
u
.We get that δ
2
H
φ
(
e
1
M
) is positive definite. Therefore the equilibrium point e
1
M is nonlinearly stable.
Ifg
<
0, then the same function φ has the property (47)
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
<
a
2
b
2
g
+
a
2
M
·
e
-
2
a
a
M
-
b
/
b
2
g
<
0
.We obtain that δ
2
H
φ
(
e
1
M
) is negative definite, and, in consequence, the equilibrium point e
1
M is nonlinearly stable.Remark 7.
Lete
2
N
,
P
=
(
N
,
0
,
P
)
∈
E
2. The eigenvalues of the characteristic polynomial associated with the linearization of system (10) at e
2
N
,
P are given by λ
1,2
=
0
,
λ
3
=
a
N
-
b
-
b
g
P
. Therefore the equilibrium point e
2
N
,
P is unstable in the case a
N
-
b
-
b
g
P
>
0.


### 3.3. Energy-Casimir Mapping

We consider the Hamilton-Poisson realization of system (10) given in Proposition 3. The corresponding energy-Casimir mapping E
C
:
R
3
→
R
2 is given by(48)
E
C
x
,
y
,
z
=
H
x
,
y
,
z
,
C
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
,
x
·
e
a
/
b
z
.The set(49)
I
m
E
C
=
h
,
c
∈
R
2
∣
∃
x
,
y
,
z
∈
R
3
:
H
x
,
y
,
z
=
h
,
C
x
,
y
,
z
=
cis called the image of the energy-Casimir mapping.We denote byΓ
M the set of images of the equilibrium points (
M
,
-
1
/
g
,
(
a
M
-
b
)
/
b
g
) through the energy-Casimir mapping; namely,(50)
Γ
M
≔
h
,
c
∈
R
2
:
h
=
M
-
1
g
+
a
2
2
b
2
g
M
2
,
c
=
M
·
e
a
2
/
b
2
g
M
-
a
/
b
g
,
M
∈
R
.We also consider the subsets(51)
Γ
M
+
≔
h
,
c
∈
Γ
M
:
M
>
0
,
(52)
Γ
M
-
≔
h
,
c
∈
Γ
M
:
M
<
0
.For g
>
0, we deduce that Γ
M
+ is the graph of a function c
=
φ
(
h
), h
>
-
1
/
g. Also, for g
<
0, we have that Γ
M
- is the graph of a function c
=
ψ
(
h
), h
<
-
1
/
g. We define the sets(53)
Σ
g
+
≔
h
,
c
∈
R
2
:
h
>
-
1
g
,
0
<
c
<
φ
h
for
g
>
0
,
(54)
Σ
g
-
≔
h
,
c
∈
R
2
:
h
<
-
1
g
,
ψ
h
<
c
<
0
for
g
<
0
.The setI
m
(
E
C
) is described in the next result.Proposition 8.
(i) Letg
>
0. The image of the energy-Casimir mapping E
C (48) is given by(55)
I
m
E
C
=
h
,
c
:
c
<
0
∪
h
,
c
:
h
≥
-
1
g
,
c
=
0
∪
Σ
g
+
∪
Γ
M
+
,where Γ
M
+ and Σ
g
+ are given by (51) and (53), respectively (Figure 1).
(ii) Letg
<
0. The image of the energy-Casimir mapping E
C (48) is given by(56)
I
m
E
C
=
h
,
c
:
c
>
0
∪
h
,
c
:
h
≤
-
1
g
,
c
=
0
∪
Σ
g
-
∪
Γ
M
-
,where Γ
M
- and Σ
g
- are given by (52) and (54), respectively (Figure 2).Figure 1
The image of the energy-Casimir mapping (g
=
10
,
a
=
2
,
b
=
1).Figure 2
The image of the energy-Casimir mapping (g
=
-
10
,
a
=
2
,
b
=
1).Proof.
(i) Leth
∈
R
,
c
∈
(
-
∞
,
0
), arbitrary. Consider y
=
-
1
/
g. The conditions H
(
x
,
y
,
z
)
=
h
,
C
(
x
,
y
,
z
)
=
c from (49) become (57)
x
=
h
+
1
g
-
g
2
z
+
1
g
2
≔
f
1
z
x
=
c
·
e
-
a
/
b
z
≔
f
2
z
.We deduce that the image of the function f
=
f
1
-
f
2 is R; hence the above system has solution for every h and c. Therefore {
(
h
,
c
)
:
c
<
0
}
⊂
I
m
(
E
C
).
The conditionc
=
0 implies x
=
0; hence the equation H
(
x
,
y
,
z
)
=
h has solution if and only if h
≥
-
1
/
g; that is, {
(
h
,
c
)
:
h
≥
-
1
/
g
,
c
=
0
}
⊂
I
m
(
E
C
).
It is clear thatΓ
M
+
⊂
I
m
(
E
C
). It only remains to prove Σ
g
+
⊂
I
m
(
E
C
) and other pairs (
h
,
c
) do not belong to I
m
(
E
C
).
Consider the functions(58)
u
z
=
M
+
a
2
M
2
2
b
2
g
-
g
2
z
+
1
g
2
,
v
z
=
M
·
e
a
/
b
a
M
/
b
g
-
1
/
g
-
z
,
f
z
=
u
z
-
v
z
.We deduce that there is z
0
∈
R such that f
(
z
0
)
=
0 and f
(
z
)
<
0 for all z
∈
R
∖
{
z
0
}. Moreover, I
m
(
f
)
=
(
-
∞
,
0
].
Now, we fix an arbitrary pair(
h
M
e
,
c
M
e
)
∈
Γ
M
+, where h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2, c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g
,
M
>
0
. We show that (
h
,
c
M
e
)
∉
I
m
(
E
C
) for every h
<
h
M
e and (
h
,
c
M
e
)
∈
I
m
(
E
C
) for every h
>
h
M
e.
With the above notations, for a pair(
h
,
c
M
e
) the system given by (49) becomes(59)
x
=
h
+
1
g
-
g
2
y
+
1
g
2
-
g
2
z
+
1
g
2
=
u
z
+
h
-
h
M
e
-
g
2
y
+
1
g
2
x
=
c
M
e
·
e
-
a
/
b
z
=
v
z
,and hence(60)
f
z
+
h
-
h
M
e
=
g
2
y
+
1
g
2
.Because I
m
(
f
)
=
(
-
∞
,
0
] we obtain that there is z such that f
(
z
)
+
h
-
h
M
e
≥
0 for any h
>
h
M
e. Consequently, (60) has solution; that is, (
h
,
c
M
e
)
∈
I
m
(
E
C
) for every h
>
h
M
e. Therefore Σ
g
+
⊂
I
m
(
E
C
).
On the other hand, we get that (60) does not have solutions for h
<
h
M
e. Therefore (
h
,
c
e
)
∉
I
m
(
E
C
) for every h
<
h
M
e, which finishes the proof of (i).
(ii) The conclusion follows using the same arguments as in the first case.Remark 9.
The image of the energy-Casimir mapping is a nonconvex subset ofR
2. Moreover, it is not a closed set, and, clearly, it is not the convex hull of the set of the images of the stable equilibrium points of the system through the map E
C.
Taking into account the results that have been reported in the papers [5–9], we notice that our example shows there is no a general result regarding the properties of the image of the energy-Casimir mapping. Furthermore, the answer to the question “can I
m
(
E
C
) be a nonconvex set?” is affirmative.
Because one of the constants of motion is not a polynomial function, it remains an open problem to establish that the results observed in the above-mentioned papers are true in the cases when the constants of motion are polynomials.Remark 10.
Another property that has been reported is the following. As a closed set, the setI
m
(
E
C
) has the boundary given by images of some stable equilibrium points of the system through the energy-Casimir mapping. In our case, this property is partially true, in the sense that only a part of the boundary of I
m
(
E
C
), namely, the set Γ
M
+ (51), is formed by the images of stable equilibrium points through E
C (g
>
0, see Figure 1). If g
<
0, a similar result is obtained for the set Γ
M
- (52) (see Figure 2).Remark 11.
It is easy to see that the image through the energy-Casimir mapping of a family of equilibrium points that has the formE
(
M
)
=
(
x
(
M
)
,
y
(
M
)
,
z
(
M
)
), M
∈
R, is a curve included in Im(
E
C
). In our case, for g
>
0, we have E
C
(
E
1
)
=
Γ
M
+
∪
Γ
M
-
,
s
∪
Γ
M
-
,
u (Figure 1), where the superscripts s and u mean stable and unstable, respectively. On the other hand, the second family of equilibrium points depends of two parameters. It is natural to ask about the image of this family through the energy-Casimir mapping.In the next result we give the set of all images of the equilibrium points that belong toE
2 through the energy-Casimir mapping.Proposition 12.
(i) Letg
>
0.
(a) For all(
h
,
c
)
∈
I
m
(
E
C
), c
<
0, there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
).
(b) For everyc
≥
0, let h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2 such that c
=
c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g
,
M
≥
0 (see Figure 3). Then for every h
≥
h
M
e
+
1
/
2
g there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
M
e
).
(ii) Letg
<
0.
(a) For all(
h
,
c
)
∈
I
m
(
E
C
), c
>
0, there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
).
(b) For everyc
≤
0, let h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2 such that c
=
c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g
,
M
≤
0. Then for every h
≤
h
M
e
+
1
/
2
g there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
M
e
).Figure 3
The images of the equilibrium points(
N
,
0
,
P
) through the energy-Casimir mapping: Σ
2
∪
Σ
3 (g
=
1
,
a
=
2
,
b
=
1).Proof.
(i) We have(61)
E
C
N
,
0
,
P
=
N
+
P
+
g
2
P
2
,
N
·
e
a
/
b
P
=
h
,
c
,and hence (62)
N
=
c
·
e
-
a
/
b
P
,
c
·
e
-
a
/
b
P
+
P
+
g
2
P
2
=
h
.We denote(63)
f
P
=
c
·
e
-
a
/
b
P
+
P
+
g
2
P
2
-
h
.(a) If c
<
0, then the image of the function f is R. Therefore there is P
∈
R such that f
(
P
)
=
0. Consequently, there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
).
(b) For eachM
≥
0, let (
h
M
e
,
c
M
e
)
∈
Γ
M
+; that is, h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2, c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g. If c
=
c
M
e, then the function f (63) becomes (64)
f
P
=
M
·
e
a
/
b
a
M
/
b
g
-
1
/
g
-
P
+
P
+
g
2
P
2
-
h
.We obtain Imf
=
[
h
M
e
+
1
/
2
g
-
h
,
∞
), where h
M
e
+
1
/
2
g
-
h
=
f
(
a
M
/
b
g
-
1
/
g
)
. Therefore the equation f
(
P
)
=
0 has solutions if and only if h
≥
h
M
e
+
1
/
2
g, and conclusion follows.
(ii) It is analogous.Remark 13.
BecauseE
C
(
E
2
)
⊂
I
m
(
E
C
), from the above Proposition we deduce that if g
>
0, then E
C
(
E
2
)
=
Σ
2
∪
Σ
3, where Σ
2 contains all the points (
h
,
c
) situated within E
C
(
0,0
,
P
), that is, c
=
0, h
≥
-
1
/
2
g, and E
C
(
M
,
0
,
a
M
/
b
g
-
1
/
g
), M
≥
0, that is, Γ
2 (see Figure 3). A similar result is obtained in the case g
<
0.


## 3.1. Hamilton-Poisson Realizations

We recall that the dynamical system generated by theC
1 vector field f
=
(
f
1
,
f
2
,
f
3
) on a C
∞ manifold P
⊆
R
3,(13)
x
˙
=
f
x
,
x
=
x
1
,
x
2
,
x
3
,has the Hamilton-Poisson realization (
P
,
{
·
,
·
}
,
H
), if it can be put in the form(14)
x
˙
i
=
x
i
,
H
,
i
∈
1,2
,
3
,where H is the Hamiltonian function, and {
·
,
·
} is a Poisson bracket on P.Considering a smooth functionC on P, a Poisson structure on P is generated by the Poisson bracket(15)
f
,
g
C
≔
∂
C
,
f
,
g
∂
x
1
,
x
2
,
x
3
=
∂
x
1
C
∂
x
2
C
∂
x
3
C
∂
x
1
f
∂
x
2
f
∂
x
3
f
∂
x
1
g
∂
x
2
g
∂
x
3
g
,for every f
,
g
∈
C
∞
(
R
3
,
R
). In addition, see, for example, [20], using any smooth function ν on P, a new Poisson bracket is given by(16)
f
,
g
C
ν
≔
ν
·
f
,
g
C
.In both cases, the function C is a Casimir of the Poisson structure; that is, {
C
,
f
}
=
0 for every f. In coordinates, the Poisson structure is given in matrix notation(17)
Π
C
ν
=
ν
Π
C
=
ν
0
∂
z
C
-
∂
y
C
-
∂
z
C
0
∂
x
C
∂
y
C
-
∂
x
C
0
.Moreover, if there is a smooth function H such that system (13) takes the form x
˙
t
=
Π
C
ν
·
∇
H, then (13) is a Hamilton-Poisson system.In our case, letg
≠
0 and let(18)
C
x
,
y
,
z
=
C
1
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
.We obtain the rescaling function ν
(
x
,
y
,
z
)
=
b
y
·
e
-
(
a
/
b
)
z and the Poisson structure generated by C
1:(19)
Π
C
1
=
0
g
z
+
1
-
g
y
-
1
-
g
z
-
1
0
1
g
y
+
1
-
1
0
=
0
g
z
-
g
y
-
g
z
0
0
g
y
0
0
+
0
1
-
1
-
1
0
1
1
-
1
0
.Due to its linearity, the above Poisson structure is a Lie-Poisson structure on the dual vector space of a Lie algebra, namely, s
e
(
2
)
∗. Indeed, consider the special Euclidean Lie group S
E
(
2
) of all orientation-preserving isometries (see, for example, [21]), given by(20)
S
E
2
=
A
∈
G
L
3
,
R
∣
A
=
1
0
0
u
cos
⁡
φ
-
sin
⁡
φ
v
sin
⁡
φ
cos
⁡
φ
,
u
,
v
∈
R
,
φ
∈
0,2
π
.The corresponding Lie algebra of S
E
(
2
) is(21)
s
e
2
=
X
∈
g
l
3
,
R
∣
X
=
0
0
0
β
0
-
α
γ
α
0
,
α
,
β
,
γ
∈
R
,with the commutator bracket [
X
,
Y
]
=
X
Y
-
Y
X
.As vector space,s
e
(
2
) has the basis B
s
e
(
2
)
=
{
e
1
,
e
2
,
e
3
}, where(22)
e
1
=
0
0
0
0
0
-
g
0
g
0
,
e
2
=
0
0
0
1
0
0
0
0
0
,
e
3
=
0
0
0
0
0
0
1
0
0
.We obtain the following bracket relations: (23)
e
1
,
e
2
=
g
e
3
,
e
1
,
e
3
=
-
g
e
2
,
e
2
,
e
3
=
0
.We consider the bilinear form Θ
:
s
e
(
2
)
×
s
e
(
2
)
→
R given by the matrix (
Θ
i
j
)
1
≤
i
,
j
≤
3, Θ
12
=
-
Θ
21
=
1, Θ
13
=
-
Θ
31
=
-
1, and Θ
23
=
-
Θ
32
=
1. By straightforward computations we obtain that Θ satisfies condition (24)
Θ
X
,
Y
,
Z
+
Θ
Y
,
Z
,
X
+
Θ
Z
,
X
,
Y
=
0
,for every triplet (
X
,
Y
,
Z
) of elements in s
e
(
2
)
. Therefore Θ is a symplectic cocycle of the Lie algebra s
e
(
2
). Moreover, it is not a coboundary since Θ
(
e
2
,
e
3
)
=
1
≠
0
=
f
(
[
e
2
,
e
3
]
), for every linear map f
,
f
:
s
e
(
2
)
→
R
. Following [22], the modified Lie-Poisson structure Π
C
1 is defined on the dual space s
e
(
2
)
∗
≃
R
3.The following result gives a Hamilton-Poisson realization of the considered system.Proposition 2.
Letν be the rescaling function given by (25)
ν
x
,
y
,
z
=
b
y
·
e
-
a
/
b
z
.If Π
C
1 is the Poisson structure generated by the Casimir function C
1, (26)
C
1
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
,given by (27)
x
,
y
C
1
=
g
z
+
1
,
x
,
z
C
1
=
-
g
y
-
1
,
y
,
z
C
1
=
1
,then system (10) has the Hamilton-Poisson realization (28)
R
3
,
·
,
·
C
1
ν
,
C
2
,where {
f
,
g
}
C
1
ν
=
ν
·
{
f
,
g
}
C
1 for any f
,
g
∈
C
∞
(
R
3
,
R
), and C
2 is the Hamiltonian function given by (29)
C
2
x
,
y
,
z
=
x
·
e
a
/
b
z
.Proof.
An easy computation shows that system (10) takes the form Π
C
1
ν
·
∇
C
2
=
(
x
˙
,
y
˙
,
z
˙
)
t; hence the conclusion follows.Using the same notations as in Proposition2, we similarly obtain the following results.Proposition 3.
Letμ be the rescaling function given by μ
(
x
,
y
,
z
)
=
-
ν
(
x
,
y
,
z
) and let C
2 be a Casimir function. If Π
C
2
μ is the Poisson structure generated by C
2 and μ, given in matrix notation by (30)
Π
C
2
μ
=
0
-
a
x
y
0
a
x
y
0
-
b
y
0
b
y
0
,then system (10) has the Hamilton-Poisson realization (31)
R
3
,
·
,
·
C
2
μ
,
C
1
,where C
1 is the Hamiltonian function.Remark 4.
The above Poisson structures are compatible and(32)
Π
C
1
ν
∇
C
2
=
Π
C
2
μ
∇
C
1
,and hence (10) is a bi-Hamiltonian system. Moreover, this pair of Hamilton-Poisson realizations gives rise to infinitely many Hamilton-Poisson realizations of system (10) (see Proposition 5).Proposition 5.
Letp
,
q
,
r
,
s
∈
R such that p
s
-
q
r
=
1, and let ν
(
x
,
y
,
z
)
=
b
y
·
e
-
(
a
/
b
)
z be the rescaling function. There exists infinitely many Hamilton-Poisson realizations of system (10) given by (
R
3
,
Π
p
,
q
,
H
r
,
s
), where the Hamiltonian H
r
,
s is given by(33)
H
r
,
s
=
r
C
1
+
s
C
2
=
r
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
+
s
x
·
e
a
/
b
z
,the Poisson structure is defined by(34)
Π
p
,
q
=
p
Π
C
1
ν
+
q
Π
C
2
ν
,and a Casimir function of the above Poisson structure is given by(35)
C
p
,
q
=
p
C
1
+
q
C
2
=
p
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
+
q
x
·
e
a
/
b
z
.Proof.
It is clear thatΠ
p
,
q
=
ν
(
p
Π
C
1
+
q
Π
C
2
). Therefore Π
p
,
q is the Poisson structure generated by the Casimir function C
p
,
q and the rescaling function ν. Using the Hamilton-Poisson realizations given in Propositions 2 and 3, we immediately obtain Π
p
,
q
·
∇
H
r
,
s
=
(
x
˙
,
y
˙
,
z
˙
)
t, which finishes the proof.


## 3.2. Stability of the Equilibrium Points

The equilibrium points of system (10) are given by the following families:(36)
E
1
=
M
,
-
1
g
,
a
M
-
b
b
g
:
M
∈
R
,
E
2
=
N
,
0
,
P
:
N
,
P
∈
R
.We remark that the second family represents the set of all equilibrium points of Kermack-McKendrick system (2). Therefore, the existence of another family of equilibrium points produces changes in the dynamics of initial system (2). We are concerned with the study of these new equilibrium points.Proposition 6.
Letg
∈
R
∖
{
0
} and let e
1
M
=
M
,
-
1
/
g
,
(
a
M
-
b
)
/
b
g
∈
E
1 be an equilibrium point of system (10).
(i) Ifg
(
a
2
M
+
b
2
g
)
<
0, then e
1
M is an unstable equilibrium point.
(ii) Ifg
(
a
2
M
+
b
2
g
)
>
0, then e
1
M is a nonlinearly stable equilibrium point.Proof.
(i) LetA be the matrix of the linear part of our system; that is, (37)
A
=
-
a
y
1
+
g
y
-
a
x
-
2
a
g
x
0
a
y
a
x
-
b
-
b
g
z
-
b
g
y
0
b
+
2
b
g
y
0
.The characteristic roots of A
(
e
1
M
) are given by (38)
λ
1
=
0
,
λ
2,3
=
±
-
b
2
-
a
2
M
g
.Considering M such that g
(
a
2
M
+
b
2
g
)
<
0, we conclude that the equilibrium point e
1
M is unstable.
(ii) Now, letM be such that g
(
a
2
M
+
b
2
g
)
>
0. We use the energy-Casimir method [23]. Let H
φ be the energy-Casimir function: (39)
H
φ
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
+
φ
x
·
e
a
/
b
z
,where φ
:
R
→
R is a smooth real valued function defined on R.
The first variation ofH
φ is given by (40)
δ
H
φ
x
,
y
,
z
=
1
+
e
a
/
b
z
·
φ
˙
x
·
e
a
/
b
z
δ
x
+
1
+
g
y
δ
y
+
1
+
g
z
+
a
x
b
·
e
a
/
b
z
·
φ
˙
x
·
e
a
/
b
z
δ
z
,where φ
˙
(
u
)
≔
d
φ
/
d
u. We have (41)
δ
H
φ
e
1
M
=
1
+
e
a
a
M
-
b
/
b
2
g
·
φ
˙
M
·
e
a
a
M
-
b
/
b
2
g
δ
x
+
a
M
b
δ
z
,which vanishes if and only if(42)
φ
˙
M
·
e
a
a
M
-
b
/
b
2
g
=
-
e
-
a
a
M
-
b
/
b
2
g
.The second variation of H
φ is given by (43)
δ
2
H
φ
x
,
y
,
z
=
e
2
a
/
b
z
·
φ
¨
x
·
e
a
/
b
z
δ
x
2
+
g
δ
y
2
+
x
·
a
2
b
2
·
e
a
/
b
z
g
+
φ
˙
x
·
e
a
/
b
z
+
x
·
e
a
/
b
z
φ
¨
x
·
e
a
/
b
z
δ
z
2
+
2
a
b
·
e
a
/
b
z
φ
˙
x
·
e
a
/
b
z
+
x
·
e
a
/
b
z
φ
¨
x
·
e
a
/
b
z
δ
x
δ
z
.Taking into account relation (42), we obtain (44)
δ
2
H
φ
e
1
M
=
e
2
a
a
M
-
b
/
b
2
g
·
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
δ
x
2
+
g
δ
y
2
+
g
-
M
·
a
2
b
2
+
M
2
·
a
2
b
2
·
e
2
a
a
M
-
b
/
b
2
g
·
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
δ
z
2
+
2
a
b
M
·
e
2
a
a
M
-
b
/
b
2
g
·
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
-
1
δ
x
δ
z
.If g
>
0, then we choose a function φ such that relation (42) holds and (45)
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
>
a
2
b
2
g
+
a
2
M
·
e
-
2
a
a
M
-
b
/
b
2
g
>
0
.For example, let (46)
φ
u
=
1
2
g
+
a
2
b
2
g
+
a
2
M
·
e
-
2
a
a
M
-
b
/
b
2
g
u
-
M
·
e
a
a
M
-
b
/
b
2
g
2
-
e
-
a
a
M
-
b
/
b
2
g
u
.We get that δ
2
H
φ
(
e
1
M
) is positive definite. Therefore the equilibrium point e
1
M is nonlinearly stable.
Ifg
<
0, then the same function φ has the property (47)
φ
¨
M
·
e
a
a
M
-
b
/
b
2
g
<
a
2
b
2
g
+
a
2
M
·
e
-
2
a
a
M
-
b
/
b
2
g
<
0
.We obtain that δ
2
H
φ
(
e
1
M
) is negative definite, and, in consequence, the equilibrium point e
1
M is nonlinearly stable.Remark 7.
Lete
2
N
,
P
=
(
N
,
0
,
P
)
∈
E
2. The eigenvalues of the characteristic polynomial associated with the linearization of system (10) at e
2
N
,
P are given by λ
1,2
=
0
,
λ
3
=
a
N
-
b
-
b
g
P
. Therefore the equilibrium point e
2
N
,
P is unstable in the case a
N
-
b
-
b
g
P
>
0.


## 3.3. Energy-Casimir Mapping

We consider the Hamilton-Poisson realization of system (10) given in Proposition 3. The corresponding energy-Casimir mapping E
C
:
R
3
→
R
2 is given by(48)
E
C
x
,
y
,
z
=
H
x
,
y
,
z
,
C
x
,
y
,
z
=
x
+
y
+
z
+
g
2
y
2
+
g
2
z
2
,
x
·
e
a
/
b
z
.The set(49)
I
m
E
C
=
h
,
c
∈
R
2
∣
∃
x
,
y
,
z
∈
R
3
:
H
x
,
y
,
z
=
h
,
C
x
,
y
,
z
=
cis called the image of the energy-Casimir mapping.We denote byΓ
M the set of images of the equilibrium points (
M
,
-
1
/
g
,
(
a
M
-
b
)
/
b
g
) through the energy-Casimir mapping; namely,(50)
Γ
M
≔
h
,
c
∈
R
2
:
h
=
M
-
1
g
+
a
2
2
b
2
g
M
2
,
c
=
M
·
e
a
2
/
b
2
g
M
-
a
/
b
g
,
M
∈
R
.We also consider the subsets(51)
Γ
M
+
≔
h
,
c
∈
Γ
M
:
M
>
0
,
(52)
Γ
M
-
≔
h
,
c
∈
Γ
M
:
M
<
0
.For g
>
0, we deduce that Γ
M
+ is the graph of a function c
=
φ
(
h
), h
>
-
1
/
g. Also, for g
<
0, we have that Γ
M
- is the graph of a function c
=
ψ
(
h
), h
<
-
1
/
g. We define the sets(53)
Σ
g
+
≔
h
,
c
∈
R
2
:
h
>
-
1
g
,
0
<
c
<
φ
h
for
g
>
0
,
(54)
Σ
g
-
≔
h
,
c
∈
R
2
:
h
<
-
1
g
,
ψ
h
<
c
<
0
for
g
<
0
.The setI
m
(
E
C
) is described in the next result.Proposition 8.
(i) Letg
>
0. The image of the energy-Casimir mapping E
C (48) is given by(55)
I
m
E
C
=
h
,
c
:
c
<
0
∪
h
,
c
:
h
≥
-
1
g
,
c
=
0
∪
Σ
g
+
∪
Γ
M
+
,where Γ
M
+ and Σ
g
+ are given by (51) and (53), respectively (Figure 1).
(ii) Letg
<
0. The image of the energy-Casimir mapping E
C (48) is given by(56)
I
m
E
C
=
h
,
c
:
c
>
0
∪
h
,
c
:
h
≤
-
1
g
,
c
=
0
∪
Σ
g
-
∪
Γ
M
-
,where Γ
M
- and Σ
g
- are given by (52) and (54), respectively (Figure 2).Figure 1
The image of the energy-Casimir mapping (g
=
10
,
a
=
2
,
b
=
1).Figure 2
The image of the energy-Casimir mapping (g
=
-
10
,
a
=
2
,
b
=
1).Proof.
(i) Leth
∈
R
,
c
∈
(
-
∞
,
0
), arbitrary. Consider y
=
-
1
/
g. The conditions H
(
x
,
y
,
z
)
=
h
,
C
(
x
,
y
,
z
)
=
c from (49) become (57)
x
=
h
+
1
g
-
g
2
z
+
1
g
2
≔
f
1
z
x
=
c
·
e
-
a
/
b
z
≔
f
2
z
.We deduce that the image of the function f
=
f
1
-
f
2 is R; hence the above system has solution for every h and c. Therefore {
(
h
,
c
)
:
c
<
0
}
⊂
I
m
(
E
C
).
The conditionc
=
0 implies x
=
0; hence the equation H
(
x
,
y
,
z
)
=
h has solution if and only if h
≥
-
1
/
g; that is, {
(
h
,
c
)
:
h
≥
-
1
/
g
,
c
=
0
}
⊂
I
m
(
E
C
).
It is clear thatΓ
M
+
⊂
I
m
(
E
C
). It only remains to prove Σ
g
+
⊂
I
m
(
E
C
) and other pairs (
h
,
c
) do not belong to I
m
(
E
C
).
Consider the functions(58)
u
z
=
M
+
a
2
M
2
2
b
2
g
-
g
2
z
+
1
g
2
,
v
z
=
M
·
e
a
/
b
a
M
/
b
g
-
1
/
g
-
z
,
f
z
=
u
z
-
v
z
.We deduce that there is z
0
∈
R such that f
(
z
0
)
=
0 and f
(
z
)
<
0 for all z
∈
R
∖
{
z
0
}. Moreover, I
m
(
f
)
=
(
-
∞
,
0
].
Now, we fix an arbitrary pair(
h
M
e
,
c
M
e
)
∈
Γ
M
+, where h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2, c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g
,
M
>
0
. We show that (
h
,
c
M
e
)
∉
I
m
(
E
C
) for every h
<
h
M
e and (
h
,
c
M
e
)
∈
I
m
(
E
C
) for every h
>
h
M
e.
With the above notations, for a pair(
h
,
c
M
e
) the system given by (49) becomes(59)
x
=
h
+
1
g
-
g
2
y
+
1
g
2
-
g
2
z
+
1
g
2
=
u
z
+
h
-
h
M
e
-
g
2
y
+
1
g
2
x
=
c
M
e
·
e
-
a
/
b
z
=
v
z
,and hence(60)
f
z
+
h
-
h
M
e
=
g
2
y
+
1
g
2
.Because I
m
(
f
)
=
(
-
∞
,
0
] we obtain that there is z such that f
(
z
)
+
h
-
h
M
e
≥
0 for any h
>
h
M
e. Consequently, (60) has solution; that is, (
h
,
c
M
e
)
∈
I
m
(
E
C
) for every h
>
h
M
e. Therefore Σ
g
+
⊂
I
m
(
E
C
).
On the other hand, we get that (60) does not have solutions for h
<
h
M
e. Therefore (
h
,
c
e
)
∉
I
m
(
E
C
) for every h
<
h
M
e, which finishes the proof of (i).
(ii) The conclusion follows using the same arguments as in the first case.Remark 9.
The image of the energy-Casimir mapping is a nonconvex subset ofR
2. Moreover, it is not a closed set, and, clearly, it is not the convex hull of the set of the images of the stable equilibrium points of the system through the map E
C.
Taking into account the results that have been reported in the papers [5–9], we notice that our example shows there is no a general result regarding the properties of the image of the energy-Casimir mapping. Furthermore, the answer to the question “can I
m
(
E
C
) be a nonconvex set?” is affirmative.
Because one of the constants of motion is not a polynomial function, it remains an open problem to establish that the results observed in the above-mentioned papers are true in the cases when the constants of motion are polynomials.Remark 10.
Another property that has been reported is the following. As a closed set, the setI
m
(
E
C
) has the boundary given by images of some stable equilibrium points of the system through the energy-Casimir mapping. In our case, this property is partially true, in the sense that only a part of the boundary of I
m
(
E
C
), namely, the set Γ
M
+ (51), is formed by the images of stable equilibrium points through E
C (g
>
0, see Figure 1). If g
<
0, a similar result is obtained for the set Γ
M
- (52) (see Figure 2).Remark 11.
It is easy to see that the image through the energy-Casimir mapping of a family of equilibrium points that has the formE
(
M
)
=
(
x
(
M
)
,
y
(
M
)
,
z
(
M
)
), M
∈
R, is a curve included in Im(
E
C
). In our case, for g
>
0, we have E
C
(
E
1
)
=
Γ
M
+
∪
Γ
M
-
,
s
∪
Γ
M
-
,
u (Figure 1), where the superscripts s and u mean stable and unstable, respectively. On the other hand, the second family of equilibrium points depends of two parameters. It is natural to ask about the image of this family through the energy-Casimir mapping.In the next result we give the set of all images of the equilibrium points that belong toE
2 through the energy-Casimir mapping.Proposition 12.
(i) Letg
>
0.
(a) For all(
h
,
c
)
∈
I
m
(
E
C
), c
<
0, there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
).
(b) For everyc
≥
0, let h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2 such that c
=
c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g
,
M
≥
0 (see Figure 3). Then for every h
≥
h
M
e
+
1
/
2
g there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
M
e
).
(ii) Letg
<
0.
(a) For all(
h
,
c
)
∈
I
m
(
E
C
), c
>
0, there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
).
(b) For everyc
≤
0, let h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2 such that c
=
c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g
,
M
≤
0. Then for every h
≤
h
M
e
+
1
/
2
g there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
M
e
).Figure 3
The images of the equilibrium points(
N
,
0
,
P
) through the energy-Casimir mapping: Σ
2
∪
Σ
3 (g
=
1
,
a
=
2
,
b
=
1).Proof.
(i) We have(61)
E
C
N
,
0
,
P
=
N
+
P
+
g
2
P
2
,
N
·
e
a
/
b
P
=
h
,
c
,and hence (62)
N
=
c
·
e
-
a
/
b
P
,
c
·
e
-
a
/
b
P
+
P
+
g
2
P
2
=
h
.We denote(63)
f
P
=
c
·
e
-
a
/
b
P
+
P
+
g
2
P
2
-
h
.(a) If c
<
0, then the image of the function f is R. Therefore there is P
∈
R such that f
(
P
)
=
0. Consequently, there is an equilibrium point (
N
,
0
,
P
) such that E
C
(
N
,
0
,
P
)
=
(
h
,
c
).
(b) For eachM
≥
0, let (
h
M
e
,
c
M
e
)
∈
Γ
M
+; that is, h
M
e
=
M
-
1
/
g
+
(
a
2
/
2
b
2
g
)
M
2, c
M
e
=
M
·
e
(
a
2
/
b
2
g
)
M
-
a
/
b
g. If c
=
c
M
e, then the function f (63) becomes (64)
f
P
=
M
·
e
a
/
b
a
M
/
b
g
-
1
/
g
-
P
+
P
+
g
2
P
2
-
h
.We obtain Imf
=
[
h
M
e
+
1
/
2
g
-
h
,
∞
), where h
M
e
+
1
/
2
g
-
h
=
f
(
a
M
/
b
g
-
1
/
g
)
. Therefore the equation f
(
P
)
=
0 has solutions if and only if h
≥
h
M
e
+
1
/
2
g, and conclusion follows.
(ii) It is analogous.Remark 13.
BecauseE
C
(
E
2
)
⊂
I
m
(
E
C
), from the above Proposition we deduce that if g
>
0, then E
C
(
E
2
)
=
Σ
2
∪
Σ
3, where Σ
2 contains all the points (
h
,
c
) situated within E
C
(
0,0
,
P
), that is, c
=
0, h
≥
-
1
/
2
g, and E
C
(
M
,
0
,
a
M
/
b
g
-
1
/
g
), M
≥
0, that is, Γ
2 (see Figure 3). A similar result is obtained in the case g
<
0.


## 4. Conclusions

In [4], Tudoran et al. have considered the energy-Casimir mapping associated with a Hamilton-Poisson system and have proposed an open problem regarding the connections between dynamical properties of a Hamilton-Poisson system and the corresponding energy-Casimir mapping. The observed properties remain true for some particular systems [5–9]. It was natural to ask if there are other cases [12]. In our paper, we have considered such a case, obtained by using integrable deformations of the Kermack-McKendrick model. We have given Hamilton-Poisson realizations of the considered system. We have also studied the stability of the new family of equilibrium points that has developed in the considered dynamics. Furthermore, we have pointed out some properties of the energy-Casimir mapping associated with the considered system.In our case, the image of the energy-Casimir mapping has other properties than those reported for other systems, which leaves room for further studies such as the existence of the periodic orbits of the considered system around some nonlinearly stable equilibrium points that belong to the first family, as well as homoclinic and heteroclinic orbits.



---

*Source: 5398768-2018-06-03.xml*