Stability Analysis and Chaos Diagnosis of a Tumor-Host Immune Cell Interaction Model with Neimark Sacker Bifurcation ()
1. Introduction
Mathematical modeling is one of the tools that can enhance the understanding of physical, chemical and biological phenomena. Physical modeling and chemical modeling are completely different disciplines compared to biological systems modeling. Biological systems are subject to the law of chemistry and physics, yet the specific functions of representing living systems are what distinguish biology from other sciences. In recent years, the use of integer order and fractional order mathematical modeling has become increasingly popular the fixed-point theory has also been used to study the properties of these systems of the study by many researchers [1] [2] [3] [4] [5].
Mathematical models representing the number of cancer cells and their interaction with Immune cells are generally a language that explains complex phenomena; and for more than a century, a lot of resources, whether economic or human, have been expended to fight cancer. The increase in the incidence of cancer among people has proven to be an important area of interdisciplinary research that requires expert medical and biological researchers and recent studies of cancer.
I realized the importance and support of mathematics and computer science in treating cancer and the immune system is the body’s first line of defense against parasites, infections, etc. Although studies of tumor immune dynamics date back to the 1890s, questions about tumor growth and its interaction with immune cells remain today. And diagnosing the tumor early is the most important stage, as the chances of detecting the tumor in the early stages are very less It plays an important role in the development of the tumor immune model and has a vital role in the analysis of parameters that affect the success or failure of the immune response against the tumor [6] [7] [8].
The study of tumor model immune reaction dynamics is complex but has gained great interest and attracted many researchers to research in this field. The dynamic system may undergo changes in the stability zone at various parameters. These changes in the system in the initial state are considered chaotic behavior. Chaotic behavior is not always unexpected and in certain cases, chaos can be predicted the chaos in the tumor can occur it often indicates that the immune reaction model leads to long-term tumor relapse and interpretation of the disorder of tumor growth is the lack of response to treatment. Thus, chaotic behavior somehow provides knowledge of some necessary controls that can be implemented for development in the fight against cancer [9].
Two-dimensional system is dissipative if
, where
is the Jacobi matrix of f at p [10].,A dissipative map satisfies
, so that the two Lyapunov numbers of each trajectory satisfy
[10].
In this paper, the stability of a three-dimensional kinetic system representing the interaction of host cells with cancer cells and immune cells [11] was studied and analyzed, and system chaos was diagnosed. Thus, distributed as follows: system description in Section 2 and the fixed points existence in Section 3, stability analysis (characteristic equation roots test, Jury table, Lyapunov function test) [12] [13] [14] [15] in Section 4, dissipativity [16] in Section 5, numerical analysis and dynamic behavior of system [17], and Neimark-Sacker [18] in Section 6, Lyapunov exponent in Section 7, and adaptive control technique [15] in Section 8.
2. System Description
The nonlinear dynamical system that is three-dimensional and discrete time represent Host-Immune-Tumor Cells Interaction Model given as follows:
(1)
where:
: Host cell population represent.
: Effector Immune cells represent.
: Tumor cells represent.
r: growth rate host cell.
: rate of tumor cells killing host cells.
: rate growth effector immune cells.
: Tumor cells inhibition rate of effector immune cells.
: rate of effector cells natural death.
: tumor cells growth rate.
: rate host cells tumor killing.
: rate tumor killing by effector immune cells.
When
,
,
,
,
,
,
,
.
3. The Fixed Points Existence
The following non-negative points are the fixed points of System (1), obtained by solving the following equations:
(2)
4. Stability Analysis
First, the Jacobi matrix of System (1) must be found before analyzing its stability
(3)
When:
a) Characteristic Equation Roots Test
Proposition (1): let
(4)
Is the characteristic equation of the system, the cases following are true:
1) Sink fixed point if all
,
and is asymptotically stable.
2) Saddle fixed point if at least one of
,
and is unstable.
3) Source if all
(i.e., the fixed point unstable).
Now, to check the system’s fixed points stability by characteristic equation roots substitute the fixed points into the Jacobi matrix we get characteristic polynomials are the form:
(5)
Were
To test the stability of the system, we substituted the fixed points in the Jacobi matrix.
Now we check stability of point
, substitute the point
in to the Jacobi matrix weget:
Find (Det(λI-J) = 0) we get:
From Equation (5) then
We will get the eigenvalues:
Then from Proposition (1) we get
and
Then
is called saddle point and is unstable.
Now we check stability of point
, substitute the point
in to the Jacobi matrix we get:
Find (Det(λI-J) = 0) we get:
From Equation (5) then
We will get the eigenvalues:
Then from Proposition (1) we get
and
Then
is called saddle point and is unstable.
Similarly, we test the point
and we get all fixed points are unstable, so System (1) is unstable.
Table 1 summarizes the stability results for all fixed points of System (1) using the roots of the characteristic equation.
b) Jury Test
Proposition (2): let
Characteristic equation from Jacobi matrix for System (1) and Table 2 is the Jury’s table defined as following:
Such that
,
,
,
,
Table 1. Stability of the system using the roots characteristic test.
The fixed points are said to be stable if they satisfy the following conditions:
1)
.
2)
.
3)
.
4)
.
Otherwise, conditions, the fixed points are unstable.
We check the stability of point
from the equation following we got by substituting
into the Jacobi matrix
We form Table 3.
From Proposition (2) Condition (3) is not achieved therefore the point
is unstable.
and in the same way test all of the fixed points
we get all points are unstable, so System (1) is unstable.
c) Lyapunov Function Test
Proposition (3): Let the Lyapunov function equation of System (1) as following:
When
by using ΔV we get:
(6)
Then called is stable fixed point if they are
, otherwise the point is unstable.
We check stability of System (1) by Lyapunov function method
Now test the stability of point
and
substitute in Equation (6) we get:
Then
and
are unstable points.
In the same way, we tested the rest of the fixed points and got an unstable system.
5. Dissipativity
The Jacobi matrix of the point
is:
when calculating its determinant, we get:
Then is dissipative, and similarly we calculate the rest of the fixed points and we get a dissipative system.
6. Numerical Analysis and Dynamic Behavior
a) Fixed point iteration method
In this section Fixed point iteration method was used to find the solution of System (1) with less error (0.9) written program in MABLE gives a result with (x, y, z) = (0.29315, 0.68038, 0.12214).
b) Time Behavior of System (1)
A time series of System (1) was generated for 1000 iterations with System (1) with the values we got from the fixed point iteration method we get Figure 1: [(1): x(t) and time(t), (2): y(t) and time(t), (3): z(t) and time(t)].
We get phase space of the system in Figure 2.
c) Neimark-Sacker Bifurcation
For a Neimark-Sacker bifurcation to occur the eigenvalues (
) must satisfy the following conditions:
1) At the bifurcation point the true eigenvalue is not equal to ±1.
2) The absolute value of the eigenvalue equals 1.
The bifurcation from the eigenvalues at point
was obtained from the bifurcation taking into account the above from the parameters are
and
(
).
Thus, the system undergoes Neimark-Sacker bifurcation at r = 0.25 in the fixed point
where all the conditions for a Neimark-Sacker bifurcation are satisfied. Figure 3 shows us the bifurcation diagram.
7. Lyapunov Exponent
In this part of the article, we will use the Lyapunov exponential test, to study the chaos of System (1) and if at least one value of the Lyapunov exponent is greater than zero, then this means that the system is chaotic, and the values we obtained are as follows:
,
,
.
We use the following basic formula to calculate the Lyapunov dimension
Figure 2. Time trajectories of the system variables
.
Thus, System (1) is highly chaotic and Figure 4 shows it.
8. Adaptive Control Technique
We will use adaptive control technology to manipulate the chaos of the chaotic System (1) and design a law for the adaptive control with the unknown parameter r, we get the following system after adding the controllers to System (1):
(7)
As they
are the control units for air-conditioned feeding and defined as follows:
(8)
Since
they are positive numbers and the parameter
is an estimation parameter for parameter r and with a substitution (8) in (7)
We get:
(10)
Since
is the error for the estimating parameter and
,
,
and
is the [10]. estimated parameter of r.
Fixed point stability analysis of the System (10).
a) Characteristic equation roots
Before analyzing the stability of the System (10), we find the Jacobi matrix of the system:
(11)
To test the point
, we substitute the fixed point
in Equation (11) we get:
And by finding the determinant (Det(λI-J) = 0), we get:
(12)
,
,
We will get the eigenvalues:
,
,
According to Proposition (1), we obtain:
,
,
, Therefore, the point
is stable.
The same method was tested for the rest of the points, and the results of the test points obtained are shown in Table 4 which turned out to be stable points, that is System (10) is stable.
b) Jury test
To test the point
, from the coefficients of the characteristic Equation (12) for the point, we get:
Accordingly, we form Table 5, which represents a Jury’s table for point
, as follows:
From Table 5, it is clear that all the conditions of the Jury test for point
are fulfilled (Proposition (2))
is stable.
In the same way, the rest of the fixed points were tested, which turned out to be stable as well the System (10) is a stable system.
c) Lyapunov function test
Table 4. Results of the test fixed points using the test for the roots of the characteristic equation.
Here we will analyze the stability of the System (10) using the Lyapunov function test, now test for the fixed point
based on Proposition (3) we get:
Then the fixed point
is stable. In the same way, the stability of the rest of the fixed points was tested, and we obtained that all the fixed points are stable, and this indicates that the system is stable.
d) Lyapunov Exponent of system
We diagnose the chaos in the controlled System (10) using the adaptive control technique of the Lyapunov exponent of System (10) and the following values were obtained:
,
,
Since all the obtained values are negative, this indicates that System (10) is regular, and Figure 5 shows us the behavior of the controlled system.
Figure 5. Lyapunov exponent after adaptive control.
9. Conclusion
The separate three-dimensional immune tumor cell model that was dealt with in this research, and after finding the fixed points of the system, five fixed points were obtained, all of which are real points, so that the system is self-attracting, and the stability of the system was analyzed using (the characteristic equation test and Jury’s table, and the Lyapunov function test), and we got an unstable system for all fixed points. The chaos of the system was diagnosed and the Lyapunov exponent was found where two values were found, which two positives are 1.250000 and 0.890000, and a third negative is −2.55334, and thus the system was very chaotic. As well as finding the Neimark-sacker bifurcation of the system at the complex eigenvalues resulting from the substitution of the fixed point
in the Jacobi matrix, the application of adaptive control technology on the system, thus controlling the system and obtaining a new, stable and regular system with Lyapunov exponent all are negative.