Integrable Deformations and Dynamical Properties of Systems with Constant Population

: In this paper we consider systems of three autonomous ﬁrst-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. 70K20; 70K42; 70K44


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][7][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.

Hamilton-Poisson Formulations and Integrable Deformations
In this section, we consider a system of three autonomous first-order differential equations:ẋ = f 1 (x, y, z) ,ẏ = f 2 (x, y, z) ,ż = f 3 (x, y, z), whereu = du dt , such that the function H(x, y, z) = x + y + z 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][10][11]. Proposition 1. Let C ∈ C 1 (P), where P ⊆ R 3 is a C ∞ manifold, such that H and C are functionally independent on P, where H(x, y, z) = x + y + z. Then: where f u := ∂ f ∂u , is a dynamical system with constant population.
(ii) C is also a constant of motion of system (2).
(iii) (P, Π C , H) is a Hamilton-Poisson formulation of system (2), where the Poisson structure Π C is given by and H is the Hamiltonian. Furthermore, C is a Casimir of the Poisson structure Π C .
Following [12], we have: If H and C are functionally independent on P, where H(x, y, z) = x + y + z, then systeṁ 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 Π C is given by (3) and H is the Hamiltonian. Furthermore, C is a Casimir of the Poisson structure Π ν C := νΠ C .
If H and C are constants of motion of system (1) and they are functionally independent on P, then there is a differentiable function ν on an open and dense subset of P such that system (1) takes the form (4). (ii) If H and C + gβ are functionally independent on P, then a family of integrable deformations of system (1) is given by: 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 Let us consider C x = C y . Since H and C are constants of motion of system (1), it follows that It is easy to see that H andC = C + gβ are constants of motion of system (5). Then, by Proposition 2, where C is replaced byC, (5) is a Hamilton-Poisson system. Moreover, if the deformation parameter g vanishes then system (5) becomes (1), as required.
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 0 + y 0 + 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 S = {(x, y, z) ∈ R 3 : x + y + z = h}. Moreover, if C is the second constant of motion, then the trajectories belong to the intersection of the level surfaces S and {(x, y, z) ∈ R 3 : C(x, y, z) = C(x 0 , y 0 , z 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).

A Particular Case: Polynomial Kolmogorov Systems
Let P 1 , P 2 , P 3 be three continuously differentiable functions on R 3 . Consider the differential equationsẋ called the three-dimensional ecological equations or Kolmogorov equations [17,18]. If P 1 , P 2 , P 3 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][20][21][22][23][24][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 where a, b, c ∈ R.

Remark 1.
We obtained that a three-dimensional Lotka-Volterra system with constant population has the formẋ Moreover, the function C(x, y, z) = x a y b z c 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.
In addition, system (8) reads (ẋ,ẏ,ż) = ν∇H × ∇C, where C(x, y, z) = x a y b z c and the rescaling function ν is given by ν(x, y, z) = x 1−a y 1−b z 1−c . (ii) If β ∈ C 1 (P) such that H and C + gβ are functionally independent on P, then a family of integrable deformations of Lotka-Volterra system (8) is given by 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) = x a y b z c Q(x, y, z), where Q is an arbitrary polynomial of degree n. Then we obtain that system (10) becomes which is a polynomial Kolmogorov system of degree n + 2. Moreover, by Proposition 5, it is a system with constant population for whichC(x, y, z) = x a y b z c (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 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 where , a, b, c, p, q, r ∈ R. Moreover,C(x, y, z) = x a y b z c (1 + px + qy + rz) is also a constant of motion of the above system.
We denote a 1 = cq, a 2 = −br, a 3 = 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) = x a y b z c Q(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 thatC(x, y, z) = x a y b z c (1 + px + qy + rz) is a constant of motion and, consequently, (13) is a Hamilton-Poisson system, as required.

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][30][31] and references therein).
As we have seen in Section 3, the three-dimensional Lotka-Volterra system with constant population is given by In addition, the constants of motion of this system is 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 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][33][34][35][36][37][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 R 3 ; 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 ≤ a a b b c c (a + b + c) a+b+c (x + y + z) a+b+c , ∀x, y, z ≥ 0.

Moreover, the equality holds if and only if
x Proof. If xyz = 0, then the inequality holds. Let x, y, z > 0. We apply Jensen's inequality for a real concave function ϕ : The equality holds if and only if x 1 = x 2 = . . . = x n or ϕ is linear.
. Then the conclusion follows.

Remark 2.
For H(x, y, z) = x + y + z and C(x, y, z) = x a y b z c , the inequality from Lemma 1 becomes C(x, y, z) ≤ a a b b c c (a + b + c) a+b+c H a+b+c (x, y, z) , ∀x, y, z ≥ 0, (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). (17) is the set

Proposition 8. The image of the energy-Casimir mapping
Proof. By definition, 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 R 2 , 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 Σ p denotes an open subset of the partition of Im(EC) that has dimension 2, then each point of Σ p 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 R 2 , but without this restriction, it is R 2 . If the image of the energy-Casimir mapping is R 2 , the dynamics is more complicated, but we can expect to find periodic orbits and heteroclinic orbits. In the following we study the abovementioned connections. The equilibrium states of system (15) are given by the families Proposition 9. If M = 0, the equilibrium states of system (15) that belong to the family E 1 (20) are Lyapunov stable.
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 0 > 0 is near to zero. Then we obtain the particular solution of system (15) given by 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 V is a neighborhood of O(0, 0, 0), then we say that V ∩ (0, ∞) 3 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.

Remark 5.
For every M ∈ [0, ∞), the sets of images of the equilibrium states through EC are given by where "s" and "u" stand for stable and unstable, respectively. We note that in our case the boundary of Im(EC) is given by Σ s 1 ∪ Σ u 2 ∪ {(0, 0)}. We also denote by Σ p the interior of Im(EC) (see Figure 1). Therefore, 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 1 M = (aM, bM, cM) ∈ E 1 be such that M ∈ (0, ∞). Then, for each sufficiently small ε ∈ R * + , any integral surface contains at least one periodic orbit γ Proof. We apply a version of Moser's theorem in the case of a zero eigenvalue [41]. By (21), J(aM, bM, cM) has the eigenvalues λ 1 = 0, λ 2,3 = ±iM abc(a + b + c) and the eigenspace corresponding to the eigenvalue zero has dimension 1. Consider the constants of motion C 1 (x, y, z) = x + y + z and I(x, y, z) = a a b b c c M a+b+c−1 (x + y + z) − x a y b z c . We obtain W = ker dC 1 (aM, bM, cM) = Span R {(1, 0, 0), (0, 1, 0)}, dI(aM, bM, cM) = 0, and d 2 I(aM, bM, cM) W×W = a a−1 b b−1 c c M a+b+c−2 (b dx 2 + a dy 2 ), which is positive definite for all M > 0.
By Theorem 2.1 from [41], for each sufficiently small ε ∈ R * + , any integral surface I(x, y, z) − I(aM, bM, cM) = ε 2 contains at least one periodic solution of system (15) whose period is close to the period of the corresponding linear system around E 1 M , which finishes the proof.
Another special orbit of a Hamilton-Poisson system is the heteroclinic orbit. In our case we have the following.

Proof.
To obtain a parametric form of the heteroclinic orbit HE 1 , 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 and x a y b z c = 0, a, b, c > 0.
Consider y = 0. Thus, z = M − x and system (15) reduces toẋ = bx(x − M). Let x 0 ∈ (0, M) be the initial condition of the above equation. We obtain the solution x(t) = Mx 0 (M − x 0 )e bMt + x 0 and then HE 1 (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 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).
, that is, a stable equilibrium state. In addition, , where M = h 1 , that is, three unstable equilibrium states and the cycle of heteroclinic orbits that connect them (see Proposition 13).
is, a periodic orbit.
(ii) We have x a y b z c = 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 1 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 S h = {(x, y, z) ∈ R 3 : x + y + z = h, x, y, z ≥ 0}, h > 0 is locally foliated by the periodic orbits that expand from the stable equilibrium . Moreover, these orbits tend to the heteroclinic cycle given by the boundary of S h (see Figure 2, right). x + y + z = h, x, y, z ≥ 0} by periodic orbits around the stable equilibrium ( ah a+b+c , bh a+b+c , ch a+b+c ).

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 0 such that x 0 , y 0 , z 0 ∈ [0, h) and x 0 + y 0 + z 0 = h. The above analysis shows the following types of behaviors. First, if x 0 , y 0 , z 0 ∈ (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 0 = 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ẋ = cxy,ẏ = −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 HE 3 ). Analogously, if y 0 = 0, then (x(t), 0, z(t)) approaches (0, 0, h), and if x 0 = 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.

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.

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.