*Article* **Integrable Deformations and Dynamical Properties of Systems with Constant Population**

**Cristian L ˘azureanu**

Department of Mathematics, Politehnica University of Timi¸soara, Piat,a Victoriei 2, 300006 Timis,oara, Romania; cristian.lazureanu@upt.ro

**Abstract:** In this paper we consider systems of three autonomous first-order differential equations **x**˙ = **f**(**x**), **x** = (*x*, *y*, *z*),**f** = (*f*1, *f*2, *f*3) such that *x*(*t*) + *y*(*t*) + *z*(*t*) is constant for all *t*. We present some Hamilton–Poisson formulations and integrable deformations. We also analyze the case of Kolmogorov systems. We study from some standard and nonstandard Poisson geometry points of view the three-dimensional Lotka–Volterra system with constant population.

**Keywords:** Hamilton–Poisson systems; integrable deformations; Lotka–Volterra systems; Kolmogorov systems; stability; periodic orbits; heteroclinic orbits

**MSC:** 70K20; 70K42; 70K44

#### **1. Introduction**

The dynamics of the interactions between the individuals of different categories of a population may be modeled by autonomous first-order differential equations. For a period of time, there are situations when the total size of the population can be assumed constant, that is, the model is a dynamical system with constant population. For instance, we recall here the SIR model [1], a model of the opiate-using career [2], and the NERA model of drug consumption in a given population [3]. Additionally, there are cases of high-dimensional Lotka–Volterra systems with constant population [4,5].

A dynamical system with constant population always has a constant of motion. In the three-dimensional case, if the system has a second constant of motion, then it admits a Hamilton–Poisson formulation. In addition, controlled versions of such a system that are also Hamilton–Poisson may be obtained by integrable deformations [6–8].

The paper is organized as follows. In Section 2, we particularize some known results about Hamilton–Poisson formulations and integrable deformations in the case of a threedimensional dynamical system with constant population. In Section 3, we consider threedimensional Kolmogorov systems. We deduce the general form of a Lotka–Volterra system with constant population and we present a Hamilton–Poisson formulation of it. We give integrable deformations that are also with a constant population. Then, we obtain the general form of a polynomial Kolmogorov system of degree 3 with constant population and a particular version of it, which is Hamilton–Poisson. In Section 4, we study the stability of the equilibrium states of the Lotka–Volterra system with a constant population restricted to [0, ∞)3, the existence of periodic orbits, and we point out the heteroclinic orbits. We also study the properties of the energy-Casimir mapping associated with this system and their connections with the dynamics of the system.

#### **2. Hamilton–Poisson Formulations and Integrable Deformations**

In this section, we consider a system of three autonomous first-order differential equations:

$$
\dot{x} = f\_1(\mathbf{x}, y, z), \; \dot{y} = f\_2(\mathbf{x}, y, z), \; \dot{z} = f\_3(\mathbf{x}, y, z), \tag{1}
$$

**Citation:** L ˘azureanu, C. Integrable Deformations and Dynamical Properties of Systems with Constant Population. *Mathematics* **2021**, *9*, 1378. https://doi.org/10.3390/math9121378

Academic Editors: Mihaela Neamt,u, Eva Kaslik and Anca R ˘adulescu

Received: 24 May 2021 Accepted: 11 June 2021 Published: 14 June 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

where *<sup>u</sup>*˙ <sup>=</sup> *du dt* , such that the function *<sup>H</sup>*(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = *<sup>x</sup>* <sup>+</sup> *<sup>y</sup>* <sup>+</sup> *<sup>z</sup>* is a constant of motion. In this case, we say that (1) is a three-dimensional dynamical system with constant population. System (1) has only one constant of motion, namely *H*, or there may exist another functionally independent constant of motion, denoted by *C*.

In the following, we particularize some known results about Hamilton–Poisson formulations and integrable deformations in the case of system (1). For the sake of completeness, we present some proofs. For details about Poisson geometry see, for example, [9–11].

**Proposition 1.** *Let <sup>C</sup>* <sup>∈</sup> *<sup>C</sup>*1(*P*), *where <sup>P</sup>* <sup>⊆</sup> <sup>R</sup><sup>3</sup> *is a <sup>C</sup>*<sup>∞</sup> *manifold, such that <sup>H</sup> and <sup>C</sup> are functionally independent on P, where H*(*x*, *y*, *z*) = *x* + *y* + *z. Then:*

*(i) System*

$$
\dot{\mathbf{x}} = \mathbf{C}\_z - \mathbf{C}\_y \; \prime \; \dot{\mathbf{y}} = \mathbf{C}\_x - \mathbf{C}\_z \; \prime \; \dot{z} = \mathbf{C}\_y - \mathbf{C}\_{x\prime} \tag{2}
$$

*where fu* :<sup>=</sup> *<sup>∂</sup> <sup>f</sup> <sup>∂</sup>u*, *is a dynamical system with constant population.*


$$
\Pi\_{\mathbb{C}} = \begin{bmatrix}
0 & \mathbb{C}\_{z} & -\mathbb{C}\_{y} \\
\mathbb{C}\_{y} & -\mathbb{C}\_{x} & 0
\end{bmatrix},
\tag{3}
$$

*and H is the Hamiltonian. Furthermore, C is a Casimir of the Poisson structure* Π*C.*

Following [12], we have:

**Proposition 2.** *Let <sup>C</sup>*, *<sup>ν</sup>* <sup>∈</sup> *<sup>C</sup>*1(*P*), *where <sup>P</sup>* <sup>⊆</sup> <sup>R</sup><sup>3</sup> *is a <sup>C</sup>*<sup>∞</sup> *manifold. If <sup>H</sup> and <sup>C</sup> are functionally independent on P, where H*(*x*, *y*, *z*) = *x* + *y* + *z, then system*

$$
\dot{\mathbf{x}} = \nu \left( \mathbf{C}\_z - \mathbf{C}\_y \right), \ \dot{\mathbf{y}} = \nu \left( \mathbf{C}\_x - \mathbf{C}\_z \right), \ \dot{z} = \nu \left( \mathbf{C}\_y - \mathbf{C}\_x \right), \tag{4}
$$

*is a dynamical system with constant population. Moreover, C is also a constant of motion and* (*P*, *ν*Π*C*, *H*) *is a Hamilton–Poisson formulation of system* (4)*, where* Π*<sup>C</sup> is given by* (3) *and H is the Hamiltonian. Furthermore, C is a Casimir of the Poisson structure* Π*<sup>ν</sup> <sup>C</sup>* := *ν*Π*C.*

The next result follows by [13,14].

**Proposition 3.** *Let <sup>C</sup>*, *<sup>β</sup>* <sup>∈</sup> *<sup>C</sup>*1(*P*), *where <sup>P</sup>* <sup>⊆</sup> <sup>R</sup><sup>3</sup> *is a <sup>C</sup>*<sup>∞</sup> *manifold, <sup>g</sup>* <sup>∈</sup> <sup>R</sup> *and <sup>H</sup>*(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = *x* + *y* + *z*.


$$\dot{\mathbf{x}} = f\_1 + \mathbf{g}\boldsymbol{\nu}(\boldsymbol{\beta}\_{\mathbf{z}} - \boldsymbol{\beta}\_{\mathbf{y}}) \text{ , } \dot{\mathbf{y}} = f\_2 + \mathbf{g}\boldsymbol{\nu}(\boldsymbol{\beta}\_{\mathbf{x}} - \boldsymbol{\beta}\_{\mathbf{z}}) \text{ , } \dot{\mathbf{z}} = f\_3 + \mathbf{g}\boldsymbol{\nu}(\boldsymbol{\beta}\_{\mathbf{y}} - \boldsymbol{\beta}\_{\mathbf{x}}) \text{ , } \tag{5}$$

*where g is a deformation parameter. Moreover, H and C* + *gβ are constants of motion of system* (5)*.*

*(iii) If H is the only constant of motion of system* (1)*, then a family of deformations of system* (1) *for which H is also a constant of motion is given by* (5)*.*

**Proof.** (*i*) By hypothesis, *H* and *C* are functionally independent on *P*, and thus at least one of the following inequalities

$$\mathsf{C}\_{\mathfrak{X}} \neq \mathsf{C}\_{\mathfrak{Y}} \; \mathsf{C}\_{\mathfrak{X}} \neq \mathsf{C}\_{\mathfrak{Z}} \; \mathsf{C}\_{\mathfrak{Y}} \neq \mathsf{C}\_{\mathfrak{Y}} \neq \mathsf{C}\_{\mathfrak{Z}}$$

holds.

Let us consider *Cx* = *Cy*. Since *H* and *C* are constants of motion of system (1), it follows that

$$f\_1 + f\_2 + f\_3 = 0 \ \text{, } \mathbb{C}\_x f\_1 + \mathbb{C}\_y f\_2 + \mathbb{C}\_z f\_3 = 0 \ \text{, }$$

Therefore *<sup>f</sup>*<sup>1</sup> <sup>=</sup> *Cz* <sup>−</sup> *Cy Cy* − *Cx <sup>f</sup>*3, *<sup>f</sup>*<sup>2</sup> <sup>=</sup> *Cx* <sup>−</sup> *Cz Cy* − *Cx <sup>f</sup>*3. Considering *<sup>ν</sup>* <sup>=</sup> *<sup>f</sup>*<sup>3</sup> *Cy* − *Cx* , the conclu-

sion follows.

(*ii*) From (*i*) we have *ν*∇*H* × ∇*C* = (*f*1, *f*2, *f*3)=(*x*˙, *y*˙, *z*˙). Then *ν*∇*H* × ∇(*C* + *gβ*) = *ν*∇*H* × ∇*C* + *gν*∇*H* × ∇*β* = (*f*1, *f*2, *f*3) + *gν*(*β<sup>z</sup>* − *βy*, *β<sup>x</sup>* − *βz*, *β<sup>y</sup>* − *βx*). Therefore system (5) is obtained. It is easy to see that *H* and *C*˜ = *C* + *gβ* are constants of motion of system (5). Then, by Proposition 2, where *C* is replaced by *C*˜, (5) is a Hamilton–Poisson system. Moreover, if the deformation parameter *g* vanishes then system (5) becomes (1), as required. (*iii*) Using (5), it immediately follows that *H*˙ = 0.

We end this section with some considerations about the dynamics of a three-dimensional dynamical system with constant population.

Let (*x*0, *y*0, *z*0) be an initial condition of the system with constant population (1) and *h* = *x*<sup>0</sup> + *y*<sup>0</sup> + *z*0. Since *H* is a constant of motion, that is *H*(*x*(*t*), *y*(*t*), *z*(*t*)) = *H*(*x*0, *y*0, *z*0) for all *t*, it follows that *x*(*t*) + *y*(*t*) + *z*(*t*) = *h* for all *t*. Therefore, it is clear that the dynamics of such a system takes place in the level surface *<sup>S</sup>* <sup>=</sup> {(*x*, *<sup>y</sup>*, *<sup>z</sup>*) <sup>∈</sup> <sup>R</sup><sup>3</sup> : *<sup>x</sup>* <sup>+</sup> *<sup>y</sup>* <sup>+</sup> *<sup>z</sup>* <sup>=</sup> *h*}. Moreover, if *C* is the second constant of motion, then the trajectories belong to the intersection of the level surfaces *<sup>S</sup>* and {(*x*, *<sup>y</sup>*, *<sup>z</sup>*) <sup>∈</sup> <sup>R</sup><sup>3</sup> : *<sup>C</sup>*(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = *<sup>C</sup>*(*x*0, *<sup>y</sup>*0, *<sup>z</sup>*0)}.

As A. J. Lotka himself highlighted, "Periodic phenomena play an important role in nature, both organic and inorganic" [15]. In our case, if *H* and *C* are functionally independent constants of motion and the intersection of the corresponding level surfaces is a closed curve, then the trajectory may be a periodic orbit. On the other hand, if *H* is the only one constant of motion of system (1) the Bendixon condition div **f** = 0 in a simply connected region of *S* precludes periodic solutions in that region (see, e.g., [16] and references therein).

#### **3. A Particular Case: Polynomial Kolmogorov Systems**

Let *P*1, *P*2, *P*<sup>3</sup> be three continuously differentiable functions on R3. Consider the differential equations

$$
\dot{\mathbf{x}} = \mathbf{x}P\_1(\mathbf{x}, y, z) \; , \; \dot{y} = yP\_2(\mathbf{x}, y, z) \; , \; \dot{z} = zP\_3(\mathbf{x}, y, z) \; , \tag{6}
$$

called the three-dimensional ecological equations or Kolmogorov equations [17,18]. If *P*1, *P*2, *P*<sup>3</sup> are polynomial such that the maximum of their degrees is *k* − 1, then (6) is a three-dimensional polynomial Kolmogorov system of degree *k*. Particularly, if *k* = 2, then (6) is a three-dimensional Lotka–Volterra system. In this section, we deal with Kolmogorov systems of degree *k* ∈ {2, 3} with constant population. More precisely, we deduce the general forms of such systems that admit Hamilton–Poisson formulations. We recall here that the Liouville, Darboux integrability, respectively, of Kolmogorov and Lotka–Volterra systems and their dynamical behavior have been widely investigated (see, e.g., [19–25] and references therein).

**Proposition 4.** *If H*(*x*, *y*, *z*) = *x* + *y* + *z is a constant of motion of the three-dimensional Lotka– Volterra Equation* (6)*, then*

$$P\_1(\mathbf{x}, y, z) = cy - bz \text{ , } P\_2(\mathbf{x}, y, z) = -c\mathbf{x} + az \text{ , } P\_3(\mathbf{x}, y, z) = b\mathbf{x} - ay,\tag{7}$$

*where a*, *b*, *c* ∈ R*.*

**Proof.** Let *Pk*(*x*, *y*, *z*) = *akx* + *bky* + *ckz* + *dk*, *ak*, *bk*, *ck*, *dk* ∈ R, *k* ∈ {1, 2, 3}. The condition *<sup>H</sup>*˙ <sup>=</sup> 0 implies *xP*1(*x*, *<sup>y</sup>*, *<sup>z</sup>*) + *yP*2(*x*, *<sup>y</sup>*, *<sup>z</sup>*) + *zP*3(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = 0, <sup>∀</sup>*x*, *<sup>y</sup>*, *<sup>z</sup>* <sup>∈</sup> <sup>R</sup>, or equivalently,

*<sup>a</sup>*1*x*<sup>2</sup> <sup>+</sup> *<sup>b</sup>*2*y*<sup>2</sup> <sup>+</sup> *<sup>c</sup>*3*z*<sup>2</sup> + (*b*<sup>1</sup> <sup>+</sup> *<sup>a</sup>*2)*xy* + (*c*<sup>1</sup> <sup>+</sup> *<sup>a</sup>*3)*xz* + (*c*<sup>2</sup> <sup>+</sup> *<sup>b</sup>*3)*yz* <sup>+</sup> *<sup>d</sup>*1*<sup>x</sup>* <sup>+</sup> *<sup>d</sup>*2*<sup>y</sup>* <sup>+</sup> *<sup>d</sup>*4*<sup>z</sup>* <sup>=</sup> 0, <sup>∀</sup>*x*, *<sup>y</sup>*, *<sup>z</sup>* <sup>∈</sup> <sup>R</sup>.

Then *a*<sup>1</sup> = *b*<sup>2</sup> = *c*<sup>3</sup> = *d*<sup>1</sup> = *d*<sup>2</sup> = *d*<sup>3</sup> = 0, and denoting *b*<sup>1</sup> = *c*, *a*<sup>3</sup> = *b*, *c*<sup>2</sup> = *a*, the conclusion follows.

**Remark 1.** *We obtained that a three-dimensional Lotka–Volterra system with constant population has the form*

$$\dot{\mathbf{x}} = \mathbf{x}(c\mathbf{y} - b\mathbf{z}) \text{ , } \dot{\mathbf{y}} = \mathbf{y}(-c\mathbf{x} + az) \text{ , } \dot{z} = z(b\mathbf{x} - a\mathbf{y}) \text{ , } a, b, c \in \mathbb{R}. \tag{8}$$

*Moreover, the function C*(*x*, *y*, *z*) = *xaybzc is a constant of motion of this system (see also [26]). In Section 4, we will study this system from some standard and nonstandard Poisson geometry points of view (for details about such approach see, e.g., [27]).*

The next result is a consequence of Propositions 2 and 3.

**Proposition 5.** *Let a*, *b*, *c* > 0.

*(i) The three-dimensional Lotka–Volterra system with constant population* (8) *has the Hamilton– Poisson formulation* (*P*, Π, *H*), *where P* = (0, ∞)3, *H*(*x*, *y*, *z*) = *x* + *y* + *z*, *and the Poisson structure*

$$
\Pi = \begin{bmatrix}
0 & cxy & -bxz \\
bxz & -ayz & 0
\end{bmatrix}.
\tag{9}
$$

*In addition, system* (8) *reads* (*x*˙, *<sup>y</sup>*˙, *<sup>z</sup>*˙) = *<sup>ν</sup>*∇*<sup>H</sup>* × ∇*C, where <sup>C</sup>*(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = *<sup>x</sup>aybzc and the rescaling function ν is given by ν*(*x*, *y*, *z*) = *x*1−*ay*<sup>1</sup>−*bz*<sup>1</sup>−*c*.

*(ii) If <sup>β</sup>* <sup>∈</sup> *<sup>C</sup>*1(*P*) *such that <sup>H</sup> and <sup>C</sup>* <sup>+</sup> *<sup>g</sup><sup>β</sup> are functionally independent on P, then a family of integrable deformations of Lotka–Volterra system* (8) *is given by*

$$\begin{cases} \dot{\mathbf{x}} = \mathbf{x}(cy - bz) + \mathbf{g}\mathbf{x}^{1-a}y^{1-b}z^{1-c}(\boldsymbol{\beta}\_{\bar{z}} - \boldsymbol{\beta}\_{\bar{y}})\\ \dot{y} = y(-cx + az) + \mathbf{g}\mathbf{x}^{1-a}y^{1-b}z^{1-c}(\boldsymbol{\beta}\_{\bar{x}} - \boldsymbol{\beta}\_{\bar{z}})\\ \dot{z} = z(b\mathbf{x} - ay) + \mathbf{g}\mathbf{x}^{1-a}y^{1-b}z^{1-c}(\boldsymbol{\beta}\_{\bar{y}} - \boldsymbol{\beta}\_{\bar{x}}) \end{cases} \tag{10}$$

*where g* ∈ R *is a deformation parameter. Moreover,* (10) *is also a system with constant population and C* + *gβ is a constant of motion.*

Taking into account the form of system (10), it is natural to ask whether there are functions *β* such that (10) is a polynomial Kolmogorov system. The answer is affirmative. Indeed, let us consider *β*(*x*, *y*, *z*) = *xaybzcQ*(*x*, *y*, *z*), where *Q* is an arbitrary polynomial of degree *n*. Then we obtain that system (10) becomes

$$\begin{cases} \dot{x} = x(cy - bz + g(cy - bz)Q + gyz(Q\_z - Q\_y)) \\ \dot{y} = y(-cx + az + g(-cx + az)Q + gxz(Q\_x - Q\_z)) \\ \dot{z} = z(bx - ay + g(bx - ay)Q + gxy(Q\_y - Q\_x)) \end{cases} \tag{11}$$

which is a polynomial Kolmogorov system of degree *n* + 2. Moreover, by Proposition 5, it is a system with constant population for which *C*˜(*x*, *y*, *z*) = *xaybzc*(1 + *gQ*(*x*, *y*, *z*)) is a constant of motion.

In the sequel, we deduce the general form of a three-dimensional polynomial Kolmogorov system of degree 3 with a constant population. In addition, we present a family of such systems that are Hamilton–Poisson.

Proceeding as in the proof of Proposition 4, we obtain the following.

**Proposition 6.** *The family of three-dimensional polynomial Kolmogorov systems of degree 3 with constant population is given by*

$$\begin{cases} \dot{\mathbf{x}} = \mathbf{x}(cy - bz + a\_1y^2 + a\_2z^2 + a\_3xy + a\_4xz + a\_5yz) \\ \dot{\mathbf{y}} = \mathbf{y}(-cx + az - a\_3x^2 + a\_6z^2 - a\_1xy + a\_7xz + a\_8yz) \\ \dot{z} = z(bx - ay - a\_4x^2 - a\_8y^2 - (a\_5 + a\_7)xy - a\_2xz - a\_6yz) \end{cases}, a, b, c, a\_k \in \mathbb{R}. \tag{12}$$

We note that system (12) is a deformation of the Lotka–Volterra system (8). Therefore, if this deformation is an integrable deformation, then system (12) has a second constant of motion and, consequently, it is a Hamilton–Poisson system.

**Proposition 7.** *A family of three-dimensional Hamilton–Poisson polynomial Kolmogorov systems of degree 3 with constant population is given by*

$$\begin{cases} \dot{x} = x(cy - bz + cqy^2 - brz^2 + cpxy - bpxz + ((c+1)r - (b+1)q)yz) \\ \dot{y} = y(-cx + az - cpx^2 + arz^2 - cqxy + ((a+1)p - (c+1)r)xz + aqyz) \\ \dot{z} = z(bx - ay + bpx^2 - aqy^2 + ((b+1)q - (a+1)p)xy + brxz - aryz) \end{cases} \tag{13}$$

*where* , *<sup>a</sup>*, *<sup>b</sup>*, *<sup>c</sup>*, *<sup>p</sup>*, *<sup>q</sup>*,*<sup>r</sup>* <sup>∈</sup> <sup>R</sup>. *Moreover, <sup>C</sup>*˜(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = *<sup>x</sup>aybzc*(<sup>1</sup> <sup>+</sup> *px* <sup>+</sup> *qy* <sup>+</sup> *rz*) *is also a constant of motion of the above system.*

**Proof.** Using Proposition 6, it is obvious that (13) is a three-dimensional Kolmogorov system with constant population. In order to prove that (13) is a Hamilton–Poisson system, we show that it is a particular case of system (10). Let *β* be as in Proposition 5 such that system (10) becomes (12). We obtain

$$\begin{aligned} \beta\_z - \beta\_y &= \mathbf{x}^a y^{b-1} z^{c-1} (a\_1 y^2 + a\_2 z^2 + a\_3 \mathbf{x} y + a\_4 \mathbf{x} z + a\_5 y z), \\ \beta\_x - \beta\_z &= \mathbf{x}^{a-1} y^b z^{c-1} (-a\_3 \mathbf{x}^2 + a\_6 z^2 - a\_1 \mathbf{x} y + a\_7 \mathbf{x} z + a\_8 y z), \\ \beta\_y - \beta\_x &= \mathbf{x}^{a-1} y^{b-1} z^c (-a\_4 \mathbf{x}^2 - a\_8 y^2 - (a\_5 + a\_7) \mathbf{x} y - a\_2 \mathbf{x} z - a\_6 y z), \end{aligned}$$

and

$$\begin{aligned} \beta\_{xx} - \beta\_{xy} &= \mathbf{x}^{a-1} y^{b-1} z^{c-1} (aa\_1 y^2 + aa\_2 z^2 + aa\_3 xy + aa\_4 \mathbf{x} z + aa\_5 yz + a\_3 \mathbf{x} y + a\_4 \mathbf{x} z), \\ \beta\_{xy} - \beta\_{yz} &= \mathbf{x}^{a-1} y^{b-1} z^{c-1} (-ba\_3 \mathbf{x}^2 + ba\_6 \mathbf{z}^2 - ba\_1 \mathbf{x} y + ba\_7 \mathbf{x} z + ba\_8 yz - a\_1 \mathbf{x} y + a\_8 yz), \\ \beta\_{yz} - \beta\_{xz} &= \mathbf{x}^{a-1} y^{b-1} z^{c-1} (-ca\_4 \mathbf{x}^2 - ca\_8 y^2 - c(a\_5 + a\_7) \mathbf{x} y - ca\_2 \mathbf{x} z - ca\_6 yz - a\_2 \mathbf{x} z - a\_6 yz). \end{aligned}$$

Adding the above relations and equating the coefficients of the obtained polynomial to zero, we obtain

$$ca\_4 = -ba\_3,\\ a a\_5 = (\mathfrak{c} + 1)a\_6 - (b + 1)a\_8,\\ ba\_6 = -aa\_2,\\ ba\_7 = (\mathfrak{c} + 1)a\_2 - (\mathfrak{a} + 1)a\_4,\\ ca\_8 = aa\_1. \tag{14}$$

We denote *a*<sup>1</sup> = *cq*, *a*<sup>2</sup> = −*br*, *a*<sup>3</sup> = *cp*, where *p*, *q*,*r* ∈ R. Using (14), we obtain that system (12) becomes (13). We also consider the function *β* given by *β*(*x*, *y*, *z*) = *xaybzcQ*(*x*, *y*, *z*), where *Q*(*x*, *y*, *z*) = *p*1*x* + *q*1*y* + *r*1*z* such that *p* = *gp*1, *q* = *gq*1,*r* = *gr*1. Now, on one the hand, function *β* fulfills the above conditions regarding partial derivatives. On the other hand, we have already seen that for such a function *β*, system (10) becomes (11), which also takes the particular form (13).

By Proposition 5 it follows that *C*˜(*x*, *y*, *z*) = *xaybzc*(1 + *px* + *qy* + *rz*) is a constant of motion and, consequently, (13) is a Hamilton–Poisson system, as required.

#### **4. Dynamical Properties of the Three-Dimensional Lotka–Volterra System with Constant Population**

In this section we study the dynamical properties of the three-dimensional Lotka– Volterra system with constant population (8) and their connections with the corresponding

energy-Casimir mapping. Lotka–Volterra systems [28] are widely investigated from dynamical point of view (see, e.g., [29–31] and references therein).

As we have seen in Section 3, the three-dimensional Lotka–Volterra system with constant population is given by

$$\begin{cases} \dot{\mathbf{x}} = \mathbf{x}(cy - bz) \\ \dot{\mathbf{y}} = \mathbf{y}(-c\mathbf{x} + az) \\ \dot{z} = z(b\mathbf{x} - ay) \end{cases} \tag{15}$$

In addition, the constants of motion of this system is

$$H(\mathbf{x}, y, z) = \mathbf{x} + y + z \,, \, \mathbf{C}(\mathbf{x}, y, z) = \mathbf{x}^a y^b z^c. \tag{16}$$

In the following we consider *a*, *b*, *c* > 0 and we restrict our study on the set [0, ∞)3, which is the region of ecological interest.

System (15) has the Hamilton–Poisson formulation given by Proposition 5. Therefore, the corresponding energy-Casimir mapping is given by

$$\mathcal{EC}: [0, \infty)^3 \to \mathbb{R}^2, \; \mathcal{EC}(\mathbf{x}, \mathbf{y}, z) = (\mathbf{x} + \mathbf{y} + z, \mathbf{x}^a y^b z^c). \tag{17}$$

.

There are some connections between the dynamics of a Hamilton–Poisson system and the properties of the corresponding energy-Casimir mapping (see, e.g., [27,32–38]). Most of the connections are provided by the partition of the image of the energy-Casimir mapping given by the equilibrium states. In the abovementioned papers, the energy-Casimir mapping is defined on R3; thus, we are interested to see which of the properties remain true in our case.

**Lemma 1.** *Let a*, *b*, *c* > 0. *Then*

$$x^a y^b z^c \le \frac{a^a b^b c^c}{(a+b+c)^{a+b+c}} (x+y+z)^{a+b+c} \text{ , } \forall x, y, z \ge 0.$$

*Moreover, the equality holds if and only if <sup>x</sup> <sup>a</sup>* <sup>=</sup> *<sup>y</sup> <sup>b</sup>* <sup>=</sup> *<sup>z</sup> c*

**Proof.** If *xyz* = 0, then the inequality holds. Let *x*, *y*, *z* > 0. We apply Jensen's inequality for a real concave function *ϕ* : *A* ⊂ R → R, namely

$$\varphi\left(\sum\_{i=1}^{n} p\_i \mathbf{x}\_i\right) \ge \sum\_{i=1}^{n} p\_i \varphi(\mathbf{x}\_i), \text{ for all } \mathbf{x}\_i \in A, p\_i > 0, i \in \{1, 2, \dots, n\}, \text{ where } \sum\_{i=1}^{n} p\_i = 1.$$

The equality holds if and only if *x*<sup>1</sup> = *x*<sup>2</sup> = ... = *xn* or *ϕ* is linear.

$$\text{Let } \varphi: (0, \infty) \to \mathbb{R}, \varphi(t) = \ln t \text{ and } \mathbf{x}\_1 = \frac{\mathbf{x}}{a}, \mathbf{x}\_2 = \frac{y}{b}, \mathbf{x}\_3 = \frac{z}{c}, p\_1 = \frac{a}{a+b+c}, p\_2 = \frac{b}{c}$$

$$\frac{b}{a+b+c}, p\_3 = \frac{c}{a+b+c}. \text{ Then the conclusion follows. } \square$$

**Remark 2.** *For H*(*x*, *y*, *z*) = *x* + *y* + *z and C*(*x*, *y*, *z*) = *xaybzc, the inequality from Lemma 1 becomes*

$$\mathbb{C}(x, y, z) \le \frac{a^a b^b c^c}{(a+b+c)^{a+b+c}} H^{a+b+c}(x, y, z) \text{ , } \forall x, y, z \ge 0,\tag{18}$$

*and the equality holds if and only if* (*x*, *y*, *z*)=(*aM*, *bM*, *cM*), *M* ≥ 0.

Now, we can present the image of the energy-Casimir mapping (17).

**Proposition 8.** *The image of the energy-Casimir mapping* (17) *is the set*

$$\operatorname{Im}(\mathcal{E}\mathcal{C}) = \{ (h\_1, h\_2) \in \mathbb{R}^2 \, : \, h\_1 \ge 0, \, 0 \le h\_2 \le \frac{a^a b^b c^c}{(a+b+c)^{a+b+c}} h\_1^{a+b+c} \}. \tag{19}$$

**Proof.** By definition,

$$\begin{array}{rcl} \operatorname{Im}(\mathcal{E}\mathcal{C}) &=& \left\{ (h\_1, h\_2) \in \mathbb{R}^2 : \exists (\mathbf{x}, y, z) \in [0, \infty)^3 \text{ such that } \mathcal{E}\mathcal{C}(\mathbf{x}, y, z) = (h\_1, h\_2) \right\} \\ &=& \left\{ (h\_1, h\_2) \in \mathbb{R}^2 : \exists (\mathbf{x}, y, z) \in [0, \infty)^3 \text{ such that } \mathbf{x} + y + z = h\_1, \mathbf{x}^a y^b z^c = h\_2 \right\}. \end{array}$$

Using Lemma 1, we obtain the conclusion.

**Remark 3.** *The image of the considered energy-Casimir mapping is drawn in Figure 1. In the abovementioned papers, if Im*(EC) *is a proper subset of* R2*, then its boundary is given by the images of the stable equilibrium states through* EC*. Moreover, a partition of Im*(EC) *is given by the images of all equilibrium states through* EC. *If the image of an unstable family of equilibrium states through* EC *is a half-line, then there are heteroclinic orbits that connect pairs of such states. Furthermore, if* Σ*<sup>p</sup> denotes an open subset of the partition of Im*(EC) *that has dimension 2, then each point of* Σ*<sup>p</sup> is the image of at least one periodic orbit through* EC*. In our case, because* EC *is restricted to* [0, ∞)3*, the image of the energy-Casimir mapping is a proper subset of* R2*, but without this restriction, it is* R2*. If the image of the energy-Casimir mapping is* R2*, the dynamics is more complicated, but we can expect to find periodic orbits and heteroclinic orbits.*

**Figure 1.** Image of the energy-Casimir mapping (*a* + *b* + *c* < 1, *a* + *b* + *c* = 1, and *a* + *b* + *c* > 1 respectively).

In the following we study the abovementioned connections.

The equilibrium states of system (15) are given by the families

$$\mathcal{E}\_1 = \{(aM, bM, CM) \mid M \in [0, \infty)\}\,, \mathcal{E}\_2 = \{(M, 0, 0) \mid M \in [0, \infty)\}\,,$$

$$\mathcal{E}\_3 = \{(0, M, 0) \mid M \in [0, \infty)\}\,, \mathcal{E}\_4 = \{(0, 0, M) \mid M \in [0, \infty)\}.\tag{20}$$

**Proposition 9.** *If M* = 0*, the equilibrium states of system* (15) *that belong to the family* E<sup>1</sup> (20) *are Lyapunov stable.*

**Proof.** Let *M* > 0. We use the Arnold method [39]. Consider the function

$$F(\mathbf{x}, y, z) = \mathbb{C}(\mathbf{x}, y, z) - \lambda H(\mathbf{x}, y, z) = \mathbf{x}^a y^b z^c - \lambda (\mathbf{x} + y + z)^c$$

and let *<sup>W</sup>* <sup>=</sup> ker *dH*(*aM*, *bM*, *cM*) = SpanR{(1, 0, 0),(0, 1, 0)}. For *<sup>λ</sup>* <sup>=</sup> *<sup>a</sup>abbccMa*<sup>+</sup>*b*+*c*−1, we obtain:


Therefore the equilibrium state (*aM*, *bM*, *cM*), *M* > 0 is nonlinearly stable.

**Proposition 10.** *For every M* ∈ (0, ∞)*, the equilibrium states of system* (15) *that belong to the families* E*i*, *i* ∈ {2, 3, 4} (20) *are unstable.*

**Proof.** Let *J* be the Jacobian matrix of system (15). Then

$$J(x,y,z) = \begin{bmatrix} cy - bz & cx & -bx \\ -cy & -cx + az & ay \\ bz & -az & bx - ay \end{bmatrix} . \tag{21}$$

The eigenvalues of *J*(*M*, 0, 0) are given by *λ*<sup>1</sup> = 0, *λ*<sup>2</sup> = −*cM*, *λ*<sup>3</sup> = *bM*. Since *λ*<sup>3</sup> > 0, the equilibrium state (*M*, 0, 0) is unstable for all *M* > 0. Analogously, we study the stability of the other equilibrium states.

It remains to study the stability of *O*(0, 0, 0), the common point of all families of equilibrium states.

Consider that system (15) is without the restriction *x*, *y*, *z* ≥ 0. We choose the initial conditions *x*(0) = 0, *y*(0) = −*z*0, *z*(0) = *z*0, where *z*<sup>0</sup> > 0 is near to zero. Then we obtain the particular solution of system (15) given by *<sup>x</sup>*(*t*) = 0, *<sup>y</sup>*(*t*) = <sup>−</sup>*z*<sup>0</sup> 1 − *az*0*t* , *<sup>z</sup>*(*t*) = *<sup>z</sup>*<sup>0</sup> 1 − *az*0*t* , *t* ∈ 0, <sup>1</sup> *az*<sup>0</sup> . Therefore, the equilibrium state *O*(0, 0, 0) is unstable.

On the other hand, our system is restricted to [0, ∞)3. The coordinate axes and the subspaces spanned by them, and also the first octant, are all invariant sets for a Lotka– Volterra system with constant population (15). Hence, populations that start non-negative remain non-negative. If *<sup>V</sup>* is a neighborhood of *<sup>O</sup>*(0, 0, 0), then we say that *<sup>V</sup>* <sup>∩</sup> (0, <sup>∞</sup>)<sup>3</sup> is a positive neighborhood of *O*(0, 0, 0).

The next result shows the behavior of a trajectory that starts in a positive neighborhood of *O*(0, 0, 0).

**Proposition 11.** *For every positive neighborhood U of O*(0, 0, 0) *there is a positive neighborhood V of O*(0, 0, 0) *such that the trajectories* (*x*(*t*), *y*(*t*), *z*(*t*)) *of system* (15) *initially in V never leave U*.

**Proof.** Consider the initial conditions of system (15) given by (*x*0, *<sup>y</sup>*0, *<sup>z</sup>*0) <sup>∈</sup> (0, <sup>∞</sup>)3. Thus (*x*(*t*), *<sup>y</sup>*(*t*), *<sup>z</sup>*(*t*)) <sup>∈</sup> (0, <sup>∞</sup>)<sup>3</sup> for *<sup>t</sup>* <sup>&</sup>gt; 0. We have to prove that for all *<sup>ε</sup>* <sup>&</sup>gt; 0 there is *<sup>δ</sup>* <sup>&</sup>gt; 0 such that if *x*2 <sup>0</sup> + *<sup>y</sup>*<sup>2</sup> <sup>0</sup> + *<sup>z</sup>*<sup>2</sup> <sup>0</sup> <sup>&</sup>lt; *<sup>δ</sup>*, then \$*x*(*t*)<sup>2</sup> <sup>+</sup> *<sup>y</sup>*(*t*)<sup>2</sup> <sup>+</sup> *<sup>z</sup>*(*t*)<sup>2</sup> <sup>&</sup>lt; *<sup>ε</sup>* for *<sup>t</sup>* <sup>&</sup>gt; 0.

For every *<sup>ε</sup>* <sup>&</sup>gt; 0 let *<sup>δ</sup>* <sup>&</sup>gt; 0 such that *<sup>δ</sup>* <sup>&</sup>lt; *<sup>ε</sup>* √3 <sup>3</sup> . Then the following relations hold:

$$\sqrt{\mathbf{x}(t)^2 + y(t)^2 + z(t)^2} \le \mathbf{x}(t) + y(t) + z(t) = \mathbf{x}\_0 + y\_0 + z\_0 \le 3\sqrt{\frac{\mathbf{x}\_0^2 + y\_0^2 + z\_0^2}{3}} \le \sqrt{3}\delta < \varepsilon\_r$$

for every *t* > 0, as required.

**Remark 4.** *In the case when system* (15) *is restricted to* [0, ∞)3*, Proposition 11 tells us that O*(0, 0, 0) *behaves as a Lyapunov stable equilibrium state.*

**Remark 5.** *For every M* ∈ [0, ∞)*, the sets of images of the equilibrium states through* EC *are given by*

$$\begin{aligned} \Sigma\_1^s &= \{ (aM + bM + cM, a^a b^b c^c M^{a+b+c}) : M > 0 \} \\ &= \{ (h\_1, h\_2) \in \mathbb{R}^2 : h\_1 > 0, h\_2 = \frac{a^a b^b c^c}{(a+b+c)^{a+b+c}} h\_1^{a+b+c} \}, \\ \mathcal{EC}(0, 0, 0) &= \{ 0, 0 \}, \\ \Sigma\_2^u &= \{ (M, 0) : M > 0 \} = \{ (h\_1, h\_2) \in \mathbb{R}^2 : h\_1 > 0, h\_2 = 0 \}, \end{aligned}$$

*where "s" and "u" stand for stable and unstable, respectively. We note that in our case the boundary of Im*(EC) *is given by* Σ*<sup>s</sup>* <sup>1</sup> <sup>∪</sup> <sup>Σ</sup>*<sup>u</sup>* <sup>2</sup> ∪ {(0, 0)}*. We also denote by* <sup>Σ</sup>*<sup>p</sup> the interior of Im*(EC) *(see Figure 1). Therefore,*

$$Im(\mathcal{EC}) = \Sigma\_1^s \cup \{ (0,0) \} \cup \Sigma\_2^u \cup \Sigma^p. \tag{22}$$

In the next result we prove that there are periodic orbits of the considered system around the stable equilibrium states. In [40], the existence of periodic orbits of a similar system has been proven by another approach.

**Proposition 12.** *Let E*<sup>1</sup> *<sup>M</sup>* = (*aM*, *bM*, *cM*) ∈ E<sup>1</sup> *be such that M* ∈ (0, ∞)*. Then, for each sufficiently small ε* ∈ R<sup>∗</sup> <sup>+</sup>*, any integral surface*

$$\Sigma\_{\varepsilon}^{E\_{\mathbf{M}}^{\varepsilon}} : \quad a^a b^b c^c M^{a+b+c-1} (\mathbf{x} + \mathbf{y} + \mathbf{z}) \\ - \ge^a y^b z^c - a^a b^b c^c M^{a+b+c} (a+b+c-1) \\ = \varepsilon^2$$

*contains at least one periodic orbit <sup>γ</sup>E*<sup>1</sup> *M <sup>ε</sup> of system (15) whose period is close to* <sup>2</sup>*<sup>π</sup> <sup>M</sup>*√*abc*(*a*+*b*+*c*) *.*

**Proof.** We apply a version of Moser's theorem in the case of a zero eigenvalue [41]. By (21), *<sup>J</sup>*(*aM*, *bM*, *cM*) has the eigenvalues *<sup>λ</sup>*<sup>1</sup> <sup>=</sup> 0, *<sup>λ</sup>*2,3 <sup>=</sup> <sup>±</sup>*iM*\$*abc*(*<sup>a</sup>* <sup>+</sup> *<sup>b</sup>* <sup>+</sup> *<sup>c</sup>*) and the eigenspace corresponding to the eigenvalue zero has dimension 1. Consider the constants of motion *<sup>C</sup>*1(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = *<sup>x</sup>* <sup>+</sup> *<sup>y</sup>* <sup>+</sup> *<sup>z</sup>* and *<sup>I</sup>*(*x*, *<sup>y</sup>*, *<sup>z</sup>*) = *<sup>a</sup>abbccMa*<sup>+</sup>*b*+*c*−1(*<sup>x</sup>* <sup>+</sup> *<sup>y</sup>* <sup>+</sup> *<sup>z</sup>*) <sup>−</sup> *<sup>x</sup>aybzc*. We obtain *W* = ker *dC*1(*aM*, *bM*, *cM*) = SpanR{(1, 0, 0),(0, 1, 0)}, *dI*(*aM*, *bM*, *cM*) = 0, and *d*<sup>2</sup> *I*(*aM*, *bM*, *cM*) *<sup>W</sup>*×*<sup>W</sup>* <sup>=</sup> *<sup>a</sup>a*−1*bb*−1*ccMa*<sup>+</sup>*b*+*c*−2(*b dx*<sup>2</sup> <sup>+</sup> *a dy*2), which is positive definite for all *M* > 0.

By Theorem 2.1 from [41], for each sufficiently small *ε* ∈ R<sup>∗</sup> <sup>+</sup>, any integral surface *<sup>I</sup>*(*x*, *<sup>y</sup>*, *<sup>z</sup>*) <sup>−</sup> *<sup>I</sup>*(*aM*, *bM*, *cM*) = *<sup>ε</sup>*<sup>2</sup> contains at least one periodic solution of system (15) whose period is close to the period of the corresponding linear system around *E*<sup>1</sup> *<sup>M</sup>*, which finishes the proof.

Another special orbit of a Hamilton–Poisson system is the heteroclinic orbit. In our case we have the following.

**Proposition 13.** *For each M* ∈ (0, ∞) *there is a cycle of heteroclinic orbits* HE<sup>1</sup> ∪ HE<sup>2</sup> ∪ HE<sup>3</sup> *that connect the unstable equilibrium states* (*M*, 0, 0),(0, 0, *M*),(0, *M*, 0) *given by*

$$\begin{split} \mathcal{H}\mathcal{E}\_{1}(t) &= \left(\frac{M\mathbf{x}\_{0}}{(M-\mathbf{x}\_{0})e^{bMt}+\mathbf{x}\_{0}}, 0, \frac{M(M-\mathbf{x}\_{0})e^{bMt}}{(M-\mathbf{x}\_{0})e^{bMt}+\mathbf{x}\_{0}}\right), \\ \mathcal{H}\mathcal{E}\_{2}(t) &= \left(0, \frac{M(M-\mathbf{z}\_{0})e^{cMt}}{(M-\mathbf{z}\_{0})e^{aMt}+\mathbf{z}\_{0}}, \frac{M\mathbf{z}\_{0}}{(M-\mathbf{z}\_{0})e^{aMt}+\mathbf{z}\_{0}}\right), \\ \mathcal{H}\mathcal{E}\_{3}(t) &= \left(\frac{M(M-\mathbf{y}\_{0})e^{cMt}}{(M-\mathbf{y}\_{0})e^{cMt}+\mathbf{y}\_{0}}, \frac{M\mathbf{y}\_{0}}{(M-\mathbf{y}\_{0})e^{cMt}+\mathbf{y}\_{0}}, 0\right). \end{split}$$

*where t* ∈ R *and x*0, *y*0, *z*<sup>0</sup> ∈ (0, *M*).

**Proof.** To obtain a parametric form of the heteroclinic orbit HE1, we first reduce our system from three degrees of freedom to one degree of freedom by using the level sets *H*(*x*, *y*, *z*) = *H*(*M*, 0, 0) and *C*(*x*, *y*, *z*) = *C*(*M*, 0, 0) and then integrate the resulting reduced differential equation. We have

$$x + y + z = M \text{ and } x^a y^b z^c = 0, a, b, c > 0.$$

Consider *y* = 0. Thus, *z* = *M* − *x* and system (15) reduces to *x*˙ = *bx*(*x* − *M*). Let *x*<sup>0</sup> ∈ (0, *M*) be the initial condition of the above equation. We obtain the solution *x*(*t*) = *Mx*<sup>0</sup> (*M* − *x*0)*ebMt* + *x*<sup>0</sup> and then HE1(*t*). We similarly obtain the other parametric forms of the heteroclinic orbits.

**Remark 6.** *By geometrical point of view, the heteroclinic orbits given in Proposition 13 are open line segments with endpoints given by the unstable equilibrium states. Another approach in the study of the existence of heteroclinic orbits of a similar system has been used in [40].*

The dynamics of system (15) takes place at the intersection of the level sets *H*(*x*, *y*, *z*) = *h*1, *C*(*x*, *y*, *z*) = *h*2, where (*h*1, *h*2) ∈ Im(EC). The fiber of the energy-Casimir mapping EC corresponding to (*h*1, *h*2) is the set

$$\mathcal{F}\_{\left(h\_{1},h\_{2}\right)} = \left\{ \left(\mathbf{x}, \mathbf{y}, z\right) \in \left[0, \infty\right)^{3} \mid \mathcal{EC}(\mathbf{x}, \mathbf{y}, z) = \left(h\_{1}, h\_{2}\right) \right\}. \tag{23}$$

Then, in other words, the dynamics of system (15) takes place on the fibers of the corresponding energy-Casimir mapping. In the following we give a topological classification of these fibers and, as a consequence, we point out the connections between the partition of Im(EC) and the dynamics of the considered system. This classification is provided by the partition of the image of the energy-Casimir mapping (22).

**Proposition 14.** *Let* (*h*1, *h*2) ∈ *Im*(EC) (19) *and* F(*h*1,*h*2) (23) *be the corresponding fiber of the energy-Casimir mapping.*


**Proof.** (*i*) Let (*x*, *<sup>y</sup>*, *<sup>z</sup>*) <sup>∈</sup> <sup>F</sup>(*h*1,*h*2) (23), (*h*1, *<sup>h</sup>*2) <sup>∈</sup> <sup>Σ</sup>*<sup>s</sup>* <sup>1</sup>, where <sup>Σ</sup>*<sup>s</sup>* <sup>1</sup> is given in Remark 5. Using Lemma 1, we obtain *<sup>x</sup> <sup>a</sup>* <sup>=</sup> *<sup>y</sup> <sup>b</sup>* <sup>=</sup> *<sup>z</sup> c* . Since *x* + *y* + *z* = *h*1, we obatin the conclusion.

(*ii*) We have *xaybzc* = 0 and *x* + *y* + *z* = *M*. Using Proposition 13 and Remark 6, the conclusion follows.

(*iii*) The dynamics of system (15) takes place at the intersection of the level sets *H*(*x*, *y*, *z*) = *h*<sup>1</sup> and *C*(*x*, *y*, *z*) = *h*2, which is the closed curve drawn in Figure 2. On the other hand, by Proposition 12 we know that there are periodic orbits around the stable equilibrium (*aM*, *bM*, *cM*), for every *M* > 0. Hence, in this case, we conclude that the fiber is a periodic orbit.

We notice that the level set *Sh* <sup>=</sup> {(*x*, *<sup>y</sup>*, *<sup>z</sup>*) <sup>∈</sup> <sup>R</sup><sup>3</sup> : *<sup>x</sup>* <sup>+</sup> *<sup>y</sup>* <sup>+</sup> *<sup>z</sup>* <sup>=</sup> *<sup>h</sup>*, *<sup>x</sup>*, *<sup>y</sup>*, *<sup>z</sup>* <sup>≥</sup> <sup>0</sup>}, *<sup>h</sup>* <sup>&</sup>gt; 0 is locally foliated by the periodic orbits that expand from the stable equilibrium *ah a* + *b* + *c* , *bh a* + *b* + *c* , *ch a* + *b* + *c* . Moreover, these orbits tend to the heteroclinic cycle given by the boundary of *Sh* (see Figure 2, right).

**Figure 2. Left**: The intersection of level sets: a periodic orbit. **Right**: A local foliation of the level set *Sh* <sup>=</sup> {(*x*, *<sup>y</sup>*, *<sup>z</sup>*) <sup>∈</sup> <sup>R</sup><sup>3</sup> : *<sup>x</sup>* <sup>+</sup> *<sup>y</sup>* <sup>+</sup> *<sup>z</sup>* <sup>=</sup> *<sup>h</sup>*, *<sup>x</sup>*, *<sup>y</sup>*, *<sup>z</sup>* <sup>≥</sup> <sup>0</sup>} by periodic orbits around the stable equilibrium ( *ah <sup>a</sup>*+*b*+*<sup>c</sup>* , *bh <sup>a</sup>*+*b*+*<sup>c</sup>* , *ch <sup>a</sup>*+*b*+*<sup>c</sup>* ).

**Remark 7.** *Consider that system* (15) *describes the dynamics of the interactions between three categories of a constant population of total size h. After an interaction between individuals from distinct categories, they can change their category. At a moment t, we denote by x*(*t*), *y*(*t*), *z*(*t*) *the number of individuals of these categories. At t* = 0*, the initial conditions are x*(0) = *x*0, *y*(0) = *y*0, *z*(0) = *z*<sup>0</sup> *such that x*0, *y*0, *z*<sup>0</sup> ∈ [0, *h*) *and x*<sup>0</sup> + *y*<sup>0</sup> + *z*<sup>0</sup> = *h*. *The above analysis shows the following types of behaviors. First, if x*0, *y*0, *z*<sup>0</sup> ∈ (0, *h*), *then we obtain a periodic behavior around the stable equilibrium ah a* + *b* + *c* , *bh a* + *b* + *c* , *ch a* + *b* + *c* . *Second, if z*<sup>0</sup> = 0*, then z*(*t*) = 0 *for all t* > 0 *since the plane of equation z* = 0 *is an invariant set of our system. In this case, system* (15) *reduces to x*˙ = *cxy*, *y*˙ = −*cxy, and hence x increases in time and y decreases, that is* (*x*(*t*), *y*(*t*), 0) *approaches* (*h*, 0, 0) *in infinite time (the heteroclinic orbit* HE3*). Analogously, if y*<sup>0</sup> = 0, *then* (*x*(*t*), 0, *z*(*t*)) *approaches* (0, 0, *h*)*, and if x*<sup>0</sup> = 0, *then* (0, *y*(*t*), *z*(*t*)) *approaches* (0, *h*, 0). *Finally, if* (*x*0, *y*0, *z*0) *is an equilibrium state, then the number of individuals from each category stays unchanged.*

#### **5. Conclusions**

In this paper, Hamilton–Poisson formulations and integrable deformations of a threedimensional dynamical system with constant population were presented. Particularly, general forms of some Kolmogorov systems with constant population were deduced. Using an integrable deformation of a Lotka–Volterra system with a constant population, a Hamilton–Poisson version of a Kolmogorov system of degree threes with a constant population were constructed.

The dynamical properties of the three-dimensional Lotka–Volterra system with a constant population were analyzed in connection with the associated energy-Casimir mapping. Such connections were reported in the case of other Hamilton–Poisson systems and, even if there is no a general result, the expected properties hold.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** We would like to thank the referees very much for their valuable comments and suggestions.

**Conflicts of Interest:** The author declares no conflict of interest.

#### **References**


MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Mathematics* Editorial Office E-mail: mathematics@mdpi.com www.mdpi.com/journal/mathematics

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34

www.mdpi.com

ISBN 978-3-0365-5866-0