Multiplicity of positive periodic solutions for a discrete impulsive blood cell production model

: In this paper, we investigate the multiplicity of positive periodic solutions of a discrete blood cell production model with impulse e ff ects. This model is described by periodic coe ffi cients and time delays, as well as nonlinear feedback with exponential terms. By employing the Krasnosel’skii fixed point theorem, we establish a su ffi cient condition for the existence of at least two positive periodic solutions. To this end, we construct solution transformation between an impulsive delay di ff erence equation and the corresponding nonimpulsive delay di ff erence equation. Aditionally, a solution representation of the positive periodic solution of the blood cell production model is presented. Moreover, a numerical example and its simulations are given to illustrate the main result.


Introduction
Hematopoiesis is a continuous physiological process of blood cell formation that occurs primarily in the bone marrow [1].In this procedure, as a common origin, hematopoietic stem cells are responsible for the lifetime generation of the three major cell lines: erythrocytes (red blood cells), leukocytes (white blood cells), and platelets.These stem cells have the extraordinary ability to renew, self-replicate, and differentiate into other functional cells.Blood cell numbers in the hematopoietic system maintain a dynamic balance through complex regulatory mechanisms [2].The concentration of blood cells can be measured in venous blood samples [3].Under normal physiological conditions, cell concentrations are maintained within the determinable normal range.However, blood cell counts appear to have wide abnormal fluctuations under pathological conditions, leading to a variety of blood diseases [4,5].For instance, an increased red blood cell count indicates polycythemia vera, cardiovascular disease and stress.A decreased red blood cell count is an indicator of anemia, chronic renal failure and acute hemorrhaging.Therefore, exploring the dynamic changes in blood cell counts is of crucial importance for human life and health.
Conceptually simple differential equations are believed to be fundamental for modeling a large number of different physiological processes [6].Delay is typically incorporated into such models to provide a more accurate description for the dynamical systems [7,8].Delay differential equations are used to describe the development systems that depend on not only the current state but also the historical state and have been widely applied into many fields such as biology, engineering and epidemiology [9,10].To understand the dynamics of blood cell production, the behaviors of blood cell number have been widely employed by delay differential models.In 1977, first-order nonlinear autonomous differential equations, called by Mackey-Glass equations, with time delay were proposed in [11].Blood cells change over time with the increase or decrease in number.Hence, different feedback control mechanisms, responsive to the body's demand for cells, are operated in these models.They showed the dynamic behaviors including limit cycle, oscillations and chaotic solutions of models.Lasota [12] studied a delay feedback model in which the nonlinear feedback has an exponential term.By using the approach of ergodic theory, chaotic behaviour of this biological model was obtained.
Considering the important effects of environmentally driven cyclical factors on organisms, Berezansky et al. [13] studied a continuous hematopoiesis model with periodic coefficient and time delay.The feedback control mechanism was implemented by a monotonically decreasing function.
In their other work [14], a unimodal feedback control mechanism was taken into account to show the different behaviors of blood cell number.More comprehensive hematopoiesis models were presented by considering the different maturation times required for different blood cells in hematopoietic system.Continuous hematopoiesis models with multiple feedback loops dominated by different time delays (maturation times) were used in Liu et al. [15] and Wu et al. [16].On the other hand, Yao [17] and the present author [18] explored such models from a discrete perspective, wherein the discrete models are more effective than continuous ones.All these results suggest the existence or global attractivity of positive periodic solutions for hematopoiesis models under some suitable conditions.The research of positive periodic solutions is of great significance both in theoretical interest and real applications [19][20][21][22].
Impulses maybe generated in evolutionary system due to transient changes in state experienced at certain moments.In the hematopoietic process, these instantaneous and abrupt changes are caused by external interference such as received radiation, medication or other forms of blood cell pressure [23,24].The blood cell model driven by impulse effects allows discrete rapid changes in cell number and is more accurate and effective in capturing external disturbances [25].For relevant literatures about impulse in other research areas, we refer to [26][27][28][29][30][31][32][33][34].
The mechanism by which blood cells increase and decrease have gradually elucidated by clinical studies.On the other hand, it is also important to identify the essence of hematopoietic phenomena by studying relatively simple mathematical models based on various facts.From that point of view, models of blood cells production have been investigated by various studies.Existing studies on continuous models of blood cell dynamics conducted with impulsive effects can be referred to [23,[35][36][37][38].By contrast, relatively few studies on discrete dynamical models have considered the impulse terms.Therefore, in this paper, we formulate a new discrete blood cell production model considering the key impacts of impulse to study the dynamics of blood cell numbers in hematopoietic system.
The present paper is organized as follows: In Section 2, we investigate the multiplicity of positive periodic solutions of discrete impulsive blood cell production model.Section 3 constructs a solution transformation and a solution representation for the positive periodic solution of the model.Section 4 discusses the multiple positive periodic solutions, which are found by using the Krasnosel'skii fixed point theorem.Section 5 provides a numerical example and its simulations to illustrate our main result.

Blood cell production model
In this paper, we focus on the positive periodic solutions of a hematopoietic system with impulse effects.We begin with the Lasota equation [12], which describes the dynamics of the formation of blood cells based on experimental results In this model, x denotes the number of blood cells in blood circulation.The first term on the right-hand side represents the destruction of blood cells and α is the death rate of blood cells.The second term on the right-hand side represents the influx of blood cells into blood circulation, in which β and n are positive constants.The production process of blood cells is inherently nonlinear in nature, and the parameter n is intended to capture this nonlinearity that has nonmonotonic character.The delay τ is the time required for blood cells to attain maturity from immature cells made in bone marrow to mature ones.
Cyclical fluctuations in the environment have an important effect on many living systems and they give the organisms inherent periodicity.However, constant coefficients and constant time delay could not reflect such periodic environmental influence on blood cells.Hence, it is reasonable and realistic to consider the environmental cyclicity by assuming that mortality, productivity, and maturation during hematopoietic processes change periodically.
The purpose of this study is to present a sufficient condition that ensures that system (2.1) has at least two positive ω-periodic solutions.The main result is the following: holds, then (2.1) has at least two positive ω-periodic solutions, where (1 − α(r)).

Krasnosel'skii fixed point theorem and solution transformation
In this section, we start with the well-known Krasnosel'skii fixed point theorem [39] which is an efficient method for searching periodic solutions of nonlinear differential or difference equations.Lemma 3.1.(Krasnosel'skii fixed point theorem) Let (X, || • ||) be a Banach space, and let P ⊂ X be a cone in X. Suppose that Ω 1 and Ω 2 are open bounded subsets of X with θ ∈ Ω 1 ⊂ Ω 1 ⊂ Ω 2 .Let Φ : P → P be a completely continuous operator on P such that, either Then, Φ has a fixed point in P ∩ Ω 2 \Ω 1 .
Under the initial condition y(s) = ϕ(s) > 0 for s ∈ Z[− τ, 0], we consider the following nonimpulsive difference equation with time-delays: and we perform the transformation between the solutions of (2.1) and (3.1).
Proof.We first prove part (i).For any k k j and j ∈ N, we obtain = y(k + 1) It follows from (3.1) that y(k j + 1) = y(k j ).Then, Now, we prove part (ii).For any k k j and j ∈ N, one has .
Moreover, from (2.1), we have Hence, it follows the definition of y that Then, the following representation of a positive ω-periodic solution of (2.1) can be obtained.
Lemma 3.3.Assume that (H) holds.Then a positive ω-periodic sequence x is the solution x(•; φ) of (2.1) if and only if the original sequence x ∈ X satisfies Proof.(Necessity) Suppose that the positive ω-periodic sequence x is the solution x(•; φ) of (2.1).Then it is obvious that x(k) > 0 for k ∈ Z.We rewrite (3.1) as .
Summing up both sides from k to k + ω − 1 results in Using the periodicity of α and γ j , we have (Sufficiency) Direct computation shows that Hence, H is ω-periodic with respect to both variables k and s.Then, we can prove that an element x of X having the expression (3.2) is ω-periodic.In fact, Now, we construct x and φ from the element x of X given by (3.2).To prove that x is the positive ω-periodic solution x(•; φ) of (2.1), it suffices to show that x satisfies (2.1).For k k j , since βi = β i (i = 1, 2, . . ., m) for k ∈ Z, we have Moreover, .
Then, we obtain .
For k = k j , one has Because of the periodicity of γ j , we see that x(k j + 1) = (1 + γ j )x(k j ).Hence, x satisfies system (2.1).The proof is complete.

□ □
We define a cone in Banach space X by P = x ∈ X : Then, it is easy to obtain the following property of the operator Φ.
Proof.In view of the proof of sufficiency of Lemma 3.3, it follows from (3.3) that Φx(k + ω) = Φx(k).Hence, Φx ∈ X.Since 0 < α < 1 for k ∈ Z, we can estimate that Moreover, The above two inequalities lead to Therefore, Thus, Φx ∈ P. The straightforward calculations can show that Φ is a completely continuous operator.The proof is complete.□

Multiplicity of positive ω-periodic solutions
In this section, we investigate the multiplicity of positive ω-periodic solutions.The main result Theorem 2.1 is proved by the Krasnosel'skii fixed point theorem below.
Choose a constant ε 0 satisfying , Then one can find an u 0 with 0 < u 0 < n and an U with n < U such that Let u 1 = U/ρ > U.The condition (2.3) implies that there exists a sufficiently small σ * ∈ (0, Then we define four open bounded subsets Ω i (1 ≤ i ≤ 4) of X by According to the characteristics of the open bounded sets Ω i (1 ≤ i ≤ 4), the discussion is divided into four cases.
Case 1. Suppose that x is an element of P ∩ ∂Ω 1 ⊂ X.For k ∈ Z, we see that 0 ≤ x(k) < u 0 and x is ω-periodic.By (3.4) and (4.1), one has Case 2. Suppose that x is an element of P ∩ ∂Ω 2 ⊂ X.For k ∈ Z, we see that nρ = ρ∥x∥ ≤ x(k) < n and x is ω-periodic.The unimodal property of u n /e u shows that min nρ≤u≤n u n /e u = n n ρ n /e nρ .Hence, from (2.3) and (3.4), it leads to Case 3. Suppose that x is an element of P ∩ ∂Ω 3 ⊂ X.For k ∈ Z, we see that ρ(n It follows from (3.4) and (4.2) that, for k ∈ Z, Hence, ||Φx|| > n + σ * = ||x|| for x ∈ P ∩ ∂Ω 3 .Case 4. Suppose that x is an element of P ∩ ∂Ω 4 ⊂ X.For k ∈ Z, we see that ρu 1 ≤ x(k) ≤ u 1 and x is ω-periodic.By (3.4) and (4.1), one has Now we are in the position to conclude that the conditions of the Krasnosel'skii fixed point theorem are satisfied by Ω 1 , Ω 2 and by Ω 3 , Ω 4 , respectively.Therefore, the operator given by (3.3) has two fixed points.Among them, one is x 1 ∈ P ∩ Ω 2 \Ω 1 which satisfies x 1 (k) ≥ ρ||x 1 || for k ∈ Z and 0 < u 0 ≤ ||x|| ≤ n, and the other is Then, Lemma 3.3 implies that x1 is a positive ω-periodic solution of (2.1) with the initial function φ1 and x2 is another positive ω-periodic solution of (2.1) with the initial function φ2 .Thus, (2.1) has two different positive ω-periodic solutions x1 and x2 satisfying ∥ x1 ∥ ≤ n < n + σ * ≤ ∥ x2 ∥.The proof is complete.□

Numerical example
Let us consider the difference equation Here, all coefficients are positive 4-periodic discrete functions defined as follows: The impulse values are γ 1 = −0.5 for k 1 = 1 and γ 2 = 1.5 for k 2 = 3.We can verify that there are at least two positive 4-periodic solutions of (5.Therefore, condition (2.3) is satisfied.Thus, from Theorem 2.1, we see that (5.1) has at least two positive 4-periodic solutions.
Figure 1a shows that as k gradually increases, the four solutions separately determined by arbitrarily given four initial functions display the same periodic changes.In other words, there exists a positive 4-periodic solution of (5.1).The details of this positive 4-periodic solution are revealed in Figure 1b.
(a) (b) Figure 1.Graphs of four arbitrary positive solutions of system (5.1).The numerical simulations show that there is a positive 4-periodic solution of (5.1) and this positive 4periodic solution is locally asymptotically stable.
Theorem 2.1 ensures the existence of two positive 4-periodic solutions.In addition to the locally asymptotically stable positive 4-periodic solution described above, system (5.1) has another positive 4-periodic solution.In general, this solution is considered unstable and very difficult to find.

Conclusions
A discrete blood cell production model subjected to impulse effects is studied in this paper.This model considers sudden changes in the number of blood cells at certain times that cannot be ignored.The sufficient condition that guarantees the multiplicity of positive periodic solutions is established by Krasnoseli-skii fixed point theorem.To be precise, there exist at least two positive periodic solutions of the model under this condition.By using coefficients and impulse values, this sufficient condition can be easily verified.

Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.