
Citation: Chao Hu, Qing-Hua Qin. Bone remodeling and biological effects of mechanical stimulus[J]. AIMS Bioengineering, 2020, 7(1): 12-28. doi: 10.3934/bioeng.2020002
[1] | Massoumeh Nazari, Mahmoud Dehghan Nayeri, Kiamars Fathi Hafshjani . Developing mathematical models and intelligent sustainable supply chains by uncertain parameters and algorithms. AIMS Mathematics, 2024, 9(3): 5204-5233. doi: 10.3934/math.2024252 |
[2] | T. K. Buvaneshwari, D. Anuradha . Solving bi-objective bi-item solid transportation problem with fuzzy stochastic constraints involving normal distribution. AIMS Mathematics, 2023, 8(9): 21700-21731. doi: 10.3934/math.20231107 |
[3] | Zhimin Liu, Ripeng Huang, Songtao Shao . Data-driven two-stage fuzzy random mixed integer optimization model for facility location problems under uncertain environment. AIMS Mathematics, 2022, 7(7): 13292-13312. doi: 10.3934/math.2022734 |
[4] | Zhimin Liu . Data-driven two-stage sparse distributionally robust risk optimization model for location allocation problems under uncertain environment. AIMS Mathematics, 2023, 8(2): 2910-2939. doi: 10.3934/math.2023152 |
[5] | Vishwas Deep Joshi, Medha Sharma, Huda Alsaud . Solving a multi-choice solid fractional multi objective transportation problem: involving the Newton divided difference interpolation approach. AIMS Mathematics, 2024, 9(6): 16031-16060. doi: 10.3934/math.2024777 |
[6] | Jun Wang, Dan Wang, Yuan Yuan . Research on low-carbon closed-loop supply chain strategy based on differential games-dynamic optimization analysis of new and remanufactured products. AIMS Mathematics, 2024, 9(11): 32076-32101. doi: 10.3934/math.20241540 |
[7] | Xiaodie Lv, Yi Liu, Yihua Zhong . A novel method to solve the optimization problem of uncertain network system based on uncertainty theory. AIMS Mathematics, 2023, 8(3): 5445-5461. doi: 10.3934/math.2023274 |
[8] | Chantapish Zamart, Thongchai Botmart, Wajaree Weera, Prem Junsawang . Finite-time decentralized event-triggered feedback control for generalized neural networks with mixed interval time-varying delays and cyber-attacks. AIMS Mathematics, 2023, 8(9): 22274-22300. doi: 10.3934/math.20231136 |
[9] | Zifu Fan, Wanyan Xu, Wei Zhang . Research on digital transformation strategy and subsidy mechanism of manufacturing supply chain based on differential game. AIMS Mathematics, 2023, 8(10): 23850-23870. doi: 10.3934/math.20231216 |
[10] | Biwen Li, Qiaoping Huang . Synchronization issue of uncertain time-delay systems based on flexible impulsive control. AIMS Mathematics, 2024, 9(10): 26538-26556. doi: 10.3934/math.20241291 |
Chronic Hepatitis C virus (HCV) infection is a significant global health concern, affecting an estimated 71 million individuals worldwide. This bloodborne infection, primarily targeting the liver, led to approximately 399,000 deaths in 2016 due to severe complications such as liver cirrhosis and hepatocellular carcinoma. A large proportion (60%–80%) of those infected with HCV develop chronic infection, with 15%–30% of these individuals at risk of progressing to liver cirrhosis over time. The primary modes of transmission include unscreened blood transfusions, injection drug use, and the use of improperly sterilized medical equipment, making it a preventable yet persistent health issue in many parts of the world. Despite the recent advancements in antiviral treatments, particularly with the introduction of direct-acting antivirals (DAAs), which have significantly improved cure rates, interferon-alpha IFN-α therapy has long played a central role in managing HCV infection. Traditionally administered either alone or in combination with ribavirin, IFN-α has been instrumental in achieving sustained virologic response (SVR) in chronic HCV patients, thereby preventing disease progression. However, the precise mechanisms through which IFN-α exerts its therapeutic effects, especially in terms of HCV dynamics within the host, remain not fully understood. HCV infection triggers a complex interplay of biological processes within the liver. Once the virus enters hepatocytes, it initiates viral replication, leading to an immune response from the host. At the same time, the liver's natural regenerative capacity attempts to maintain homeostasis by proliferating healthy hepatocytes. This dynamic interaction between viral replication, immune activation, and hepatocellular proliferation presents a challenging scenario for effective treatment. The availability of uninfected hepatocytes, the rate of hepatocyte proliferation, and the overall homeostatic balance of the liver become critical factors in determining the progression of the infection and the host's ability to control it. To better understand these intricate processes, mathematical modeling has become a crucial tool in studying HCV dynamics. Mathematical models allow researchers to simulate the spread of the virus within the liver, analyze how the infection interacts with treatment protocols like IFN-α and DAAs, and explore the role of immune response and hepatocyte regeneration. These models not only provide insights into the transmission dynamics of HCV but also help to refine treatment strategies by identifying key factors influencing viral clearance or persistence. Through such modeling, it becomes possible to investigate the intricate balance between viral proliferation, host immune defense, and the liver's regenerative capabilities, offering valuable perspectives on both disease progression and therapeutic intervention. To illustrate, Pan and Chakrabarty [1] explored the Hopf bifurcation of the HCV model
˙T=λ−β1TV−d1T−β2TI+αI,˙I=β1TV−αI−d2I+β2TI,˙V=kI−pVZ−d3V,˙Z=cZ(t−τ)V(t−τ)−d4Z, | (1.1) |
where T, V, I, and Z respectively denote the densities of uninfected hepatocytes, virions, actively infected hepatocytes, and humoral immunity. Mojaver and Kheiri [2] explored the behavior of a class of HCV infection models with application of optimal control. Gümüş and Türk [3] examined the dynamical analysis of a HBV model. Kong et al. [4] examined the Hopf bifurcation of the HCV model
˙T=s−βTV−d1T,˙I=βTV−d2I,˙V=kI−uV, | (1.2) |
where T and I respectively denote uninfected and infected cells, while V is free virus. Furthermore, cells T are generated, die, and are infected at rates s, d1T, and βTV, cells I emerge and die at rates βTV and bI, and finally, free virus are released and decay at rates kI and uV. Jiang et al. [5] examined the bifurcation, chaos, and circuit realisation of a four-dimensional memristor system. Chen et al. [6] explored the bifurcation, chaos and fixed-time synchronization of memristor cellular neural networks. On the other hand, many authors have performed theoretical and mathematical analysis of the mathematical models [7,8,9]. Furthermore, different mathematicians [10,11,12] have first analyzed the continuous-time HCV model
˙T=s+rT(1−T+ITmax)−dT−(1−η)βVT,˙I=(1−η)βVT+rT(1−T+ITmax)−δI,˙V=(1−ϵ)pI−cV, | (1.3) |
and then showed that it takes the following form after decomposing the viral dynamics to two distinct time scales
˙u=u(1−u−v)−gu−(1−α)buv,˙v=v(1−u−v)−ζu+(1−α)buv. | (1.4) |
It is worthwhile to mention here that discrete-time models for Hepatitis C virus (HCV) dynamics are often preferred over continuous-time models due to several practical advantages. One of the primary reasons is their compatibility with the nature of clinical data collection. In practice, data such as viral load measurements, liver enzyme levels, and patient responses to treatments are typically recorded at discrete time intervals, such as daily, weekly, or monthly. Discrete-time models naturally align with these sampling intervals, allowing for direct integration of empirical data without the need for interpolation or extrapolation, which can introduce errors and complicate model validation. Furthermore, discrete-time models tend to be computationally more efficient. Solving continuous-time models involves complex numerical methods and often requires fine time discretization to ensure accuracy, which increases computational costs. In contrast, discrete-time models can achieve accurate results with simpler algorithms and larger time steps, making them more practical for large-scale simulations or real-time applications. This efficiency is particularly beneficial when modeling the long-term progression of HCV or evaluating various treatment strategies over extended periods. The flexibility of discrete-time models is another significant advantage. They can more easily incorporate complex biological processes and treatment protocols. For instance, antiviral drug effects, immune responses, and patient adherence to therapy, which often change at specific time points (e.g., weekly dosing schedules), can be modeled more intuitively in a discrete framework. This capability is crucial for accurately reflecting the real-world dynamics of HCV and its treatment, where interventions are frequently time-dependent. From an implementation perspective, discrete-time models are generally easier to understand and use. Researchers and clinicians can work with models that update state variables at each time step based on straightforward recurrence relations. This simplicity facilitates broader adoption of the models in both clinical and research settings, promoting better communication and collaboration among multidisciplinary teams. The ease of implementation also lowers the barrier for developing and testing new hypotheses or treatment strategies, accelerating progress in HCV research. Lastly, discrete-time models offer greater robustness to uncertainties and variabilities inherent to biological systems and clinical data. They can handle irregular sampling intervals and missing data more effectively, which is critical in real-world scenarios where ideal, continuous data collection is rarely feasible. This robustness ensures that the models remain reliable and useful even when the data is imperfect, enhancing their applicability and trustworthiness in guiding clinical decisions and policy-making. In summary, while continuous-time models have their theoretical advantages, the practical benefits of discrete-time models such as data compatibility, computational efficiency, flexibility, ease of implementation, and robustness to uncertainties make them a preferable choice for modeling HCV dynamics in many contexts. Due to the aforementioned facts, our aim in this paper is to explore the dynamics of a discrete-time HCV model instead of continuous-time HCV model, which is depicted in (1.4), where the discrete version is obtained by a non-standard finite difference scheme. The dynamics of a HCV infection model can vary depending on whether it is a discrete-time model or a continuous-time model. Below, we shall highlight the key differences in their dynamical behaviors, focusing on stability, bifurcations, and chaotic behaviors. In a continuous-time model, the system evolves continuously over time, with changes happening smoothly at every instant; in a discrete-time model, the system evolves at specific time intervals, and changes happen at discrete points. Discrete systems are prone to phenomena like Neimark-Sacker and period-doubling bifurcations, which can lead to quasi-periodic behavior and chaos more readily than in continuous models. Numerically, continuous models are often simulated using methods like Euler or Runge-Kutta, requiring fine time steps to approximate smooth dynamics. Discrete models, however, are simpler to implement and can reveal complex behaviors with fewer time steps. The contrast between the smooth progression of continuous models and the stepwise, often abrupt behavior of discrete models highlights how the mathematical formulation influences the biological insights derived from HCV infection dynamics. In order to derive the discrete version, we replace ˙u by ut+1−uth, ˙v by vt+1−vth, u2 by utut+1, v2 by vtvt+1, u by ut, and v by vt, and so continuous HCV (1.4) becomes
ut+1−uth=ut−utut+1−utvt−gut−(1−α)butvt,vt+1−vth=vt−utvt−vtvt+1−ζut+(1−α)butvt, | (1.5) |
with h being the integral step size. Furthermore, after simplification, (1.5) takes the form
ut+1=(1+h(1−g))ut−h(1+(1−α)b)utvt1+hut,vt+1=(1+h(1−ζ))vt−h(1−(1−α)b)utvt1+hvt. | (1.6) |
More precisely, our main findings in this paper include:
● Existence of fixed points of a discrete model (1.6).
● Dynamics at fixed points of a discrete model (1.6).
● Identification of bifurcation sets for fixed points.
● Bifurcation analysis of a discrete model (1.6).
● Study of chaos.
● Numerical validation of theoretical results.
The rest of the paper is organized as follows: Existence of fixed points are given in Section 2. Section 3 regards the local behavior at fixed points of the model (1.6). Bifurcation analysis is given in Section 4, whereas in Section 5 chaos control is studied. The numerical simulation of obtained results is given in Section 6, while the conclusions are given in Section 7.
In this section, we will study the existence of fixed points of the discrete model (1.6). If (ˆu,ˆv) is a fixed point of model (1.6), then
ˆu=(1+h(1−g))ˆu−h(1+(1−α)b)ˆuˆv1+hˆu,ˆv=(1+h(1−ζ))ˆv−h(1−(1−α)b)ˆuˆv1+hˆv. | (2.1) |
Since system (2.1) is obviously satisfied if (ˆu,ˆv)=(0,0),(1−g,0),(0,1−ζ). Therefore, ∀ b, α, g, ζ, h, discrete HCV model (1.6) has liver-free fixed point (LFFP), disease-free fixed point (DFFP) (1−g,0) if 0<g<1, and total infection fixed point (TIFP) (0,1−ζ) if 0<ζ<1. In order to find the partial-infection fixed point (PIFP), from (2.1), one has
−ˆu+(−1−(1−α)b)ˆv+1−g=0,−ˆv+(−1+(1−α)b)ˆu+1−ζ=0. | (2.2) |
From the first equation of (2.2), one gets
ˆv=−ˆu+1−g1+(1−α)b. | (2.3) |
Using (2.3) into the second equation of (2.2), we get
ˆu=ζ−g−(1−α)(1−ζ)b(1−α)2b2. | (2.4) |
Finally, using (2.4) into (2.3), one gets
ˆv=g−ζ+(1−α)(1−g)b(1−α)2b2. | (2.5) |
Equations (2.4) and (2.5) imply that if ζ−g(1−α)(1−g)<b<ζ−g(1−α)(1−ζ), then HCV model (1.6) has PIFP (ζ−g−(1−α)(1−ζ)b(1−α)2b2,g−ζ+(1−α)(1−g)b(1−α)2b2).
In this section, the local stability at fixed points of the HCV model (1.6) is analyzed using linear stability theory [13,14]. Specifically, we explore local stability by evaluating the eigenvalues of the Jacobian matrix at fixed points. This approach allows us to assess whether perturbations around the fixed points are stable or unstable. For this, we first compute the Jacobian matrix at fixed points, and then eigenvalues will be determined in order to apply linear stability theory [15,16]. So, Jacobian J|(ˆu,ˆv) evaluated at (ˆu,ˆv) under the map (g1,g2)↦(ut+1,vt+1) is
J|(ˆu,ˆv):=(1−h(−1+g+(1+b(1−α))ˆv)(1+hˆu)2−h(1+b(1−α))ˆu1+hˆu−h(1−b(1−α))ˆv1+hˆv1−h(−1+ζ+(1−b(1−α))ˆu)(1+hˆv)2), | (3.1) |
where
g1=(1+h(1−g))ut−h(1+(1−α)b)utvt1+hut,g2=(1+h(1−ζ))vt−h(1−(1−α)b)utvt1+hvt. | (3.2) |
Now, at LFFP, (3.1) becomes
J|LFFP:=(1+h−gh001+h−hζ), | (3.3) |
with
λ1=1+h−gh, λ2=1+h−hζ. | (3.4) |
From (3.4), it is easy to see that LFFP is a saddle if 2ζ−1<h<2g−1 or 2g−1<h<2ζ−1, a source if 0<h<min{2g−1,2ζ−1}, a sink if h>max{2g−1,2ζ−1}, and finally, non-hyperbolic if
h=2g−1, | (3.5) |
or
h=2ζ−1. | (3.6) |
At DFFP, (3.1) becomes
J|DFFP:=(11+h(1−g)h(−1−(1−α)b)(1−g)1+h(1−g)01+h(g−ζ+(1−α)(1−g)b)), | (3.7) |
with
λ1=11+h−hg, λ2=1+h(g−ζ+(1−α)(1−g)b). | (3.8) |
From (3.8), one can also obtain that DFFP is a saddle if 2g−1<h<2−g+ζ−(1−g)(1−α)b or 2−g+ζ−(1−g)(1−α)b<h<2g−1, a source if 0<h<min{2g−1,2−g+ζ−(1−g)(1−α)b}, a sink if h>max{2g−1,2−g+ζ−(1−g)(1−α)b}, and non-hyperbolic if
h=2g−1, | (3.9) |
or
h=2−g+ζ−(1−α)(1−g)b. | (3.10) |
Furthermore, at TIFP, (3.1) becomes
J|TIFP:=(1+h(ζ−g−(1−α)(1−ζ)b)0h(−1+(1−α)b)(1−ζ)1+h(1−ζ)11+h(1−ζ)), | (3.11) |
with
λ1=1+h(ζ−g)−h(1−α)(1−ζ)b, λ2=11+h(1−ζ), | (3.12) |
from which it is easy to prove that TIFP is a saddle if 2ζ−1<h<2−ζ+g+(1−α)(1−ζ)b or 2−ζ+g+(1−α)(1−ζ)b<h<2ζ−1, a source if 0<h<min{2ζ−1,2−ζ+g+(1−α)(1−ζ)b}, a sink if h>max{2ζ−1,2−ζ+g+(1−α)(1−ζ)b}, and non-hyperbolic if
h=2ζ−1, | (3.13) |
or
h=2−ζ+g+(1−α)(1−ζ)b. | (3.14) |
Finally, at PIFP, (3.1) becomes
J|PIFP:=((1−α)2b2b2(1−α)2−h(g−ζ)−hb(1−α)(1−ζ)h(b(1−α)+1)(g+b(1−α)(1−ζ)−ζ)(1−α)2b2−h(g−ζ)−hb(1−α)(1−ζ)h(b(1−α)−1)(g+b(1−α)(1−g)−ζ)(1−α)2b2+h(g−ζ)+hb(1−α)(1−g)(1−α)2b2(1−α)2b2+h(g−ζ)+hb(1−α)(1−g)), | (3.15) |
with
λ2−ϱ1λ+ϱ2=0, | (3.16) |
where
ϱ1=b3(1−α)3(2b(1−α)−h(g−ζ))((1−α)2b2−h(g−ζ)−hb(1−α)(1−ζ))((1−α)2b2+h(g−ζ)+hb(1−α)(1−g)),ϱ2=b4(1−α)4+h2(1−b2(1−α)2)(g+b(1−g)(1−α)−ζ)(g+b(1−ζ)(1−α)−ζ)((1−α)2b2−h(g−ζ)−hb(1−α)(1−ζ))((1−α)2b2+h(g−ζ)+hb(1−α)(1−g)). | (3.17) |
The J|PIFP has characteristics roots
λ1,2=ϱ1±√Δ2, | (3.18) |
where
Δ=ϱ21−4ϱ2=(b3(1−α)3(2b(1−α)−h(g−ζ)){((1−α)2b2−h(g−ζ)−hb(1−ζ)(1−α))×((1−α)2b2+h(g−ζ)+hb(1−α)(1−g))})2−4({b4(1−α)4+h2(1−b2(1−α)2)(g+b(1−g)(1−α)−ζ)×(g+b(1−α)(1−ζ)−ζ)}{((1−α)2b2−h(g−ζ)−hb(1−α)(1−ζ))×(b2(1−α)2+h(g−ζ)+hb(1−α)(1−g))}). | (3.19) |
If Δ<0, then PIFP is an unstable focus if h>(1−α)3(ζ−g)b3(2−(1−α)2b2)(g+b(1−g)(1−α)−ζ)(g+b(1−ζ)(1−α)−ζ), a stable focus if 0<h<(1−α)3(ζ−g)b3(2−(1−α)2b2)(g+b(1−g)(1−α)−ζ)(g+b(1−ζ)(1−α)−ζ), and non-hyperbolic if
h=(1−α)3(ζ−g)b3(2−(1−α)2b2)(g+b(1−g)(1−α)−ζ)(g+b(1−ζ)(1−α)−ζ). | (3.20) |
In addition, if Δ>0, then PIFP is an unstable node if
h>max{4b(1−α)g−ζ−√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)},4b(1−α)g−ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)}}, |
a stable node if
0<h<min{4b(1−α)g−ζ−√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)2(g−ζ)(−2+g+ζ)},4b(1−α)g−ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)}}, |
and non-hyperbolic if
h=4b(1−α)g−ζ−√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)}, | (3.21) |
or
h=4b(1−α)g−ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)}. | (3.22) |
Example 1. If b=12.3,g=0.0006,ζ=4.5, then at PIFP, the dynamics of the discrete HCV model (1.6) are given in Figure 1.
Example 2. If b=0.5,g=0.09,ζ=0.55, then at PIFP, the dynamics of the discrete HCV model (1.6) are given in Figure 2.
In this section, detailed bifurcation analysis is given based on bifurcation theory [17,18,19]. Before detailed bifurcation analysis, we first give the existence of possible bifurcation sets for fixed points as follows:
(a) If (3.5) holds, then roots of J|LFFP satisfying λ1|(3.5)=−1 but λ2|(3.5)=g−1+2(1−ζ)g−1≠1 or −1. So, period-doubling bifurcation may take place if (b,α,g,ζ,h)=L crosses through the curve
F|LFFP={L: h=2g−1}; | (4.1) |
(b) if (3.6) holds, then roots of J|LFFP satisfying λ1|(3.6)=1−2g+ζζ−1≠1 or −1 but λ2|(3.6)=−1. So, period-doubling bifurcation may take place if L crosses through the curve
F|LFSS={L: h=2ζ−1}; | (4.2) |
(c) if (3.9) holds, then roots of J|DFFP satisfying λ1|(3.9)=−1 but λ2|(3.9)=1+2(g−ζ+b(1−α)(1−g))g−1≠1 or −1. So, period-doubling bifurcation may take place if L crosses through the curve
F|DFFP={L: h=2g−1}; | (4.3) |
(d) if (3.10) holds, then roots of J|DFFP satisfying λ1|(3.10)=−g+ζ−b(1−α)(1−g)2−3g−b(1−α)(1−g)+ζ≠1 or −1 but λ2|(3.10)=−1. So, period-doubling bifurcation may take place if L crosses through the curve
F|DFFP={L: h=2−g+ζ−(1−α)(1−g)b}; | (4.4) |
(e) if (3.10) holds, then roots of J|TIFP satisfying λ1|(3.10)=1+2ζ−1(ζ−g−b(1−α)(1−ζ))≠1 or −1 but λ2|(3.10)=−1. So, period-doubling bifurcation may take place if L crosses through the curve
F|TIFP={L: h=2ζ−1}; | (4.5) |
(f) if (3.14) holds, then roots of J|TIFP satisfying λ1|(3.14)=−1 but λ2|(3.14)=g−ζ+b(1−α)(1−ζ)2−3ζ+b(1−α)(1−ζ)+g. So, period-doubling bifurcation may take place if L crosses through the curve
F|TIFP={L: h=2−ζ+g+(1−ζ)(1−α)b}; | (4.6) |
(g) if (1.3) holds, then a pair of complex roots of J|PIFP satisfying |λ1,2|(1.3)=1. So, Neimark-Sacker bifurcation may take place if L crosses through the curve
N|PIFP={L:Δ<0 and h=(1−α)3b3(ζ−g){(2−b2(1−α)2)(g+b(1−g)(1−α)−ζ)(g+b(1−ζ)(1−α)−ζ)}}; | (4.7) |
(h) if (3.21) is a true, then the value of λ1,2 at (3.21) has been calculated and is depicted in (A.1). This implies that no period-doubling bifurcation takes place if L crosses through the curve
F|PIFP={L:Δ>0 and h=4b(1−α)g−ζ−√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)}}; | (4.8) |
(i) if (3.22) holds, then the value of λ1,2 at (3.22) has been calculated and is depicted in (A.2). This implies that no period-doubling bifurcation takes place if L crosses through the curve
F|PIFP={L:Δ>0 and h=4b(1−α)g−ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)}}. | (4.9) |
Now we will study detailed bifurcation analysis at fixed points based on the bifurcation sets. It is mentioned that no period-doubling bifurcation occurs at LFFP if (4.1) and (4.2) hold, and therefore, LFFP is degenerate. Now, we examine bifurcation analysis at DFFP, TIFP, and PIFP of HCV model (1.6) if (4.3)–(4.9) hold.
Lemma 4.1. The discrete HCV model (1.6) does not undergo period-doubling at respective fixed points if (4.3)–(4.6) hold.
Proof. The HCV model (1.6), takes the following form after being restricted to v=0,
ut+1=(1+h(1−g))ut1+hut. | (4.10) |
From (4.10), we define
g1:=(1+h(1−g))u1+hu. | (4.11) |
Now, if h=ˆh=2g−1 and u=ˆu=1−g, then from (4.11), we get
∂g1∂u|u=ˆu=1−g, h=ˆh=2g−1=−1, | (4.12) |
∂2g1∂u2|u=ˆu=1−g, h=ˆh=2g−1=−4g−1, | (4.13) |
and
∂g1∂h|u=ˆu=1−g, h=ˆh=2g−1=0. | (4.14) |
So, from (4.14), no period-doubling bifurcation occurs at DFFP if (4.3) holds. Similarly, it is also easy to obtain that no period-doubling exists at the respective fixed points if (4.4)–(4.6) hold.
Remark 4.2. In order to study the Neimark-Sacker bifurcation at PIFP of the HCV model (1.6), first it is imperative to highlight the seminal work of Rahman and Blackmore [20]. This paper presented significant advancements in understanding the complex behavior exhibited by discrete dynamical systems. The study meticulously examined the onset of chaos through Neimark-Sacker bifurcations, providing crucial insights into the mechanisms that drive such phenomena. Rahman and Blackmore [20] contributions to this research are particularly noteworthy for their depth and clarity in elucidating the mathematical underpinnings of chaotic behavior in these systems. In the broader context of modern dynamical systems analysis, Rahman and Blackmore's work has consistently demonstrated a pioneering spirit. Their investigations into the stability and bifurcation of nonlinear systems have laid a foundational framework that subsequent researchers have built upon. By systematically exploring the conditions under which chaotic dynamics emerge, Rahman and Blackmore have advanced our understanding of the intricate balance between order and disorder in mathematical models of physical phenomena. The meticulous approach adopted in their study underscores the critical role that rigorous mathematical analysis plays in unveiling the complexities of dynamical systems. Moreover, the implications of Rahman and Blackmore's research extend beyond theoretical exploration, influencing practical applications in various fields where understanding chaotic behavior is essential. From engineering to biological systems, the insights derived from their work on Neimark-Sacker bifurcations have provided valuable tools for predicting and managing chaotic dynamics. This integration of theoretical rigor with practical relevance is a hallmark of Rahman and Blackmore's legacy in the field, reflecting the profound impact on both the academic community and applied sciences. In conclusion, the work of Rahman and Blackmore on the study of Neimark-Sacker bifurcations and chaotic behavior stands as a testament to their pioneering contributions to dynamical systems analysis. By bridging the gap between abstract mathematical theory and real-world applications, Rahman and Blackmore left an indelible mark on the field, inspiring future generations of researchers to explore the rich and intricate world of nonlinear dynamics. This legacy continues to influence contemporary studies and will undoubtedly remain a cornerstone of dynamical systems research for years to come.
Theorem 4.3. If (4.7) holds, then the HCV model (1.6) undergoes a Neimark-Sacker bifurcation with h being a bifurcation parameter.
Proof. By considering h as a bifurcation parameter, the discrete HCV model (1.6) becomes
ut+1=(1+(h∗+ϵ)(1−g))ut−(h∗+ϵ)(1+(1−α)b)utvt1+(h∗+ϵ)ut,vt+1=(1+(h∗+ϵ)(1−ζ))vt−(h∗+ϵ)(1−(1−α)b)utvt1+(h∗+ϵ)vt, | (4.15) |
where ϵ<<1, and the characteristic roots of J|PIFP at PIFP of ϵ-dependence HCV model (4.15) are
λ1,2=ϱ1(ϵ)±ι√4ϱ2(ϵ)−ϱ21(ϵ)2, | (4.16) |
where
ϱ1(ϵ)=b3(1−α)3(2b(1−α)−(h∗+ϵ)(g−ζ)){((1−α)2b2−(h∗+ϵ)(g−ζ)−(h∗+ϵ)b(1−α)(1−ζ))((1−α)2b2+(h∗+ϵ)(g−ζ)+(h∗+ϵ)b(1−α)(1−g))},ϱ2(ϵ)={b4(1−α)4+(h∗+ϵ)2(1−b2(1−α)2)(g+b(1−g)(1−α)−ζ)(g+b(1−ζ)(1−α)−ζ)}{((1−α)2b2−(h∗+ϵ)(g−ζ)−(h∗+ϵ)b(1−α)(1−ζ))((1−α)2b2+(h∗+ϵ)(g−ζ)+(h∗+ϵ)b(1−α)(1−g))}. | (4.17) |
From (4.16) and (4.17), one has
d|λ1,2|dϵ|ϵ=0={−b4h(1−α)4(g+b(1−α)(1−ζ)−ζ)(−2g−b(1−α)×(2−g−b(1−α)(g+b(1−g)(1−α)−ζ)−ζ)+2ζ)}{(b4(1−α)4−h2(g−ζ)2−2bh2(1−α)(g−ζ)(1−ζ)−b2h2(1−α)2(1−ζ)2)3/2×√{b4(1−α)4−h2(−1+(1−α)2b2)(g+b(1−α)(1−g)−ζ)(g+(1−α)b(1−ζ)−ζ)}}≠0. | (4.18) |
Now, (4.15) becomes
It+1=(1+(h∗+ϵ)(1−g))(It+ˆu)−(h∗+ϵ)(1+(1−α)b)(It+ˆu)(St+ˆv)1+(h∗+ϵ)(It+ˆu)−ˆu,St+1=(1+(h∗+ϵ)(1−ζ))(St+ˆv)−(h∗+ϵ)(1−(1−α)b)(It+ˆu)(St+ˆv)1+(h∗+ϵ)(St+ˆv)−ˆv, | (4.19) |
by
It=ut−ˆu,St=vt−ˆv, | (4.20) |
where ˆu=ζ−g−(1−α)(1−ζ)b(1−α)2b2 and ˆv=g−ζ+(1−α)(1−g)b(1−α)2b2. Hereafter, following Kuznetsov's argument for a Neimark-Sacker bifurcation involves a detailed analysis of the model's dynamics as it undergoes this type of bifurcation, which is characterized by the emergence of a closed invariant curve from a fixed point as parameters are varied. The first step in this analysis is to bring the system into a normal form that simplifies the bifurcation process. This normal form representation allows for the computation of the first Lyapunov coefficients, which are essential in determining the stability and nature of the bifurcating solutions. The first Lyapunov coefficient, in particular, indicates whether the bifurcation will lead to stable or unstable quasiperiodic orbits. By computing these coefficients, one can predict the onset of chaotic dynamics and the subsequent behavior of the system. Kuznetsov's methodology provides a robust framework for analyzing the local behavior near the bifurcation point, ensuring a thorough understanding of the system's transition to chaos. This approach is crucial for identifying the critical parameters and conditions that govern the complex behavior in nonlinear dynamical systems. For more detail, we also refer the reader to the existing theory [17,18,20]. Now, in order to study the normal form of the HCV model (4.19) if ϵ=0, the Taylor's series expansion at (0,0) yields
It+1=l11It+l12St+l13I2t+l14ItSt,St+1=l21It+l22St+l23ItSt+l24S2t, | (4.21) |
where
l11=1+h(1−g)−h(1+(1−α)b)ˆv(1+hˆu)2,l12=−h(1+(1−α)b)ˆu1+hˆu,l13=−2h(1+h(1−g)−h(1+(1−α)b)ˆv)(1+hˆu)3,l14=h(−1−(1−α)b)(1+hˆu)2,l21=h(−1+(1−α)b)ˆv1+hˆv,l22=1+h(1−ζ)−h(1−(1−α)b)ˆu(1+hˆv)2,l23=h(−1+(1−α)b)(1+hˆv)2,l24=−2h(1+h(1−ζ)−h(1−(1−α)b)ˆu)(1+hˆv)3. | (4.22) |
Now, if
(ItSt):=(l120η−l11−ρ)(utvt), | (4.23) |
then (4.21) becomes
(ut+1vt+1):=(η−ρρη)(utvt)+(Υ1(ut,vt)Υ2(ut,vt)), | (4.24) |
where
η=b3(1−α)3(2b(1−α)−h(g−ζ))2((1−α)2b2−h(g−ζ)−hb(1−α)(1−ζ))((1−α)2b2+h(g−ζ)+hb(1−α)(1−g)),ρ=√{−b6(1−α)6(−2(1−α)b+h(g−ζ))2+4((1−α)2b2+h(g+b(1−α)(1−g)−ζ))(b4(1−α)4−h2(−1+b2(1−α)2)×(g+b(1−g)(1−α)−ζ)(g+b(1−α)(1−ζ)−ζ))(b2(1−α)2+ζ))}2((1−α)2b2−h(g−ζ)−hb(1−ζ)(1−α))((1−α)2b2+h(g−ζ)+hb(1−α)(1−g)), | (4.25) |
and
Υ1=σ11u2t+σ12utvt,Υ2=σ21u2t+σ22utvt+σ23v2t, | (4.26) |
with
σ11=l13l12+l14η−l14l11,σ12=−ρl14,σ21=(η−l11)l12ρ(l13−l23)+(η−l11)2ρ(l14−l24),σ22=l23l12+(η−l11)(2l24−l14),σ23=−l24ρ. | (4.27) |
From (4.26), we get
∂2Υ1∂u2t|(0,0)=2σ11, ∂2Υ1∂ut∂vt|(0,0)=σ12,∂2Υ1∂v2t|(0,0)=∂3Υ1∂u3t|(0,0)=∂3Υ1∂u2t∂vt|(0,0)=∂3Υ1∂ut∂v2t|(0,0)=∂3Υ∂v3t|(0,0)=0,∂2Υ2∂u2t|(0,0)=2σ21, ∂2Υ2∂ut∂vt|(0,0)=σ22,∂2Υ2∂v2t|(0,0)=2σ23,∂3Υ2∂u3t|(0,0)=∂3Υ2∂u2t∂vt|(0,0)=∂3Υ2∂ut∂v2t|(0,0)=∂3Υ2∂v3t|(0,0)=0. | (4.28) |
Now, by the existing theory [21,22,23]. the map (4.24) undergoes Neimark-Sacker bifurcation if ψ≠0, that is,
ψ=−ℜ((1−2ˉλ)ˉλ21−λρ11ρ20)−12‖ρ11‖2−‖ρ02‖2+ℜ(ˉλρ21), | (4.29) |
where
ρ02=18(∂2Υ1∂u2t−∂2Υ1∂v2t+2∂2Υ2∂ut∂vt+ι(∂2Υ2∂u2t−∂2Υ2∂v2t+2∂2Υ1∂ut∂vt))|(0,0),ρ11=14(∂2Υ1∂u2t+∂2Υ1∂v2t+ι(∂2Υ2∂u2t+∂2Υ2∂v2t))|(0,0),ρ20=18(∂2Υ1∂u2t−∂2Υ1∂v2t+2∂2Υ2∂ut∂vt+ι(∂2Υ2∂u2t−∂2Υ2∂v2t−2∂2Υ1∂ut∂vt))|(0,0),ρ21=116(∂3Υ1∂u3t+∂3Υ1∂v3t+∂3Υ2∂u2t∂vt+∂3Υ2∂v3t+ι(∂3Υ2∂u3t+∂3Υ2∂ut∂v2t−∂3Υ1∂u2t∂vt−∂3Υ1∂v3t))|(0,0). | (4.30) |
After calculation, one gets
ρ02=14(σ11+σ22+ι(σ21−σ23+σ12)),ρ11=12(σ11+ι(σ21+σ23)),ρ20=14(σ11+σ22+ι(σ21−σ12−σ23)),ρ21=0. | (4.31) |
Using (4.31) along with (4.25) in (4.29) if one obtains ψ≠0, then at PIFP, the HCV model (1.6) undergoes N-S bifurcation. Additionally, it also requires that λm1,2≠1, m=1,⋯,4 if ϵ=0 that corresponds to ϱ1(0)≠−2,0,1,2. But if h=(1−α)3b3(ζ−g)(2−b2(1−α)2)(g+b(1−g)(1−α)−ζ)(g+b(1−ζ)(1−α)−ζ) holds, then from (4.17), we get ϱ2(0)=1, and so ϱ1(0)≠−2,2. Therefore, it is only required that ϱ1(0)≠0,1, i.e.,
g≠{8b(−1+α)+8b3(1−α)3−2b5(1−α)5−b6(1−α)6(−1+ζ)+8ζ+b4(1−α)4(−4+5ζ)−2b2(−1+α)2(−2+5ζ)±√b6(−2+b2(1−α)2)3(1−α)6(1−ζ)2}8−8b(1−α)−6b2(1−α)2+8b3(1−α)3+b4(1−α)4−2b5(1−α)5. | (4.32) |
On the other hand, also by existing theory [24,25,26,27], the Neimark-Sacker bifurcation at PIFP is also supercritical (resp. subcritical) if ψ<0 (resp. ψ>0).
Remark 4.4. It is important to note that no period-doubling bifurcation occurs if (4.8) and (4.9) hold. In conclusion, at PIFP, no other bifurcations occur except the Neimark-sacker bifurcation.
In this section, the OGY method, which was proposed by Ott et al. [28], is utilized for the HCV model (1.6). For this, one writes the HCV model (1.6) as follows
ut+1=(1+h(1−g))ut−h(1+(1−α)b)utvt1+hut=g1(ut,vt,ζ),vt+1=(1+h(1−ζ))vt−h(1−(1−α)b)utvt1+hvt=g2(ut,vt,ζ), | (5.1) |
where the control parameter is ζ to acquire the desired chaos control by small perturbations. Due to this, we restrict it where ζ∈(ζ0−σ,ζ0+σ) with σ>0, and ζ0 indicates small value corresponding to chaotic region. The trajectory is moved in the direction of the target orbit using the stabilizing feedback control strategy. Assuming that PIFP is an unstable point in the chaotic region created by the occurrence of N-S bifurcation, then the following map can be used to approximate the model (5.1)
(ut+1−ˆuvt+1−ˆv)≈J(ˆu,ˆv,ζ0)(ut−ˆuvt−ˆv)+P(ζ−ζ0), | (5.2) |
where
J(ˆu,ˆv,ζ0)=((1−α)2b2(1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1)h(1+b(1−α))(g+b(1−α)(1−ζ0)−ζ0)(1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1)h(−1+b(1−α))(g+b(1−α)(1−g)−ζ0)(1−α)2b2+h(g−ζ0)+hb(1−α)(1−g)(1−α)2b2(1−α)2b2+h(g−ζ0)+hb(1−α)(1−g)), | (5.3) |
and
P=(∂g1(ˆu,ˆv,ζ0)∂ζ∂g2(ˆu,ˆv,ζ0)∂ζ)=(0h(−g+ζ0−(1−g)(1−α)b)(1−α)2b2+h(g−ζ0+(1−g)(1−α)b)). | (5.4) |
Moreover, the HCV model (5.1) is controlled provided that the following matrix
C=(P:J(ˆu,ˆv,ζ0)P)= (0h(−g+ζ0−(1−α)(1−g)b)b2(1−α)2+h(g−ζ0+(1−g)(1−α)b)h2(g−ζ0+b(1−ζ0)(1−α))(g−ζ0+(1−α)(1−g)b)(−1−b(1−α))((1−α)2b2+h(ζ0−g)+hb(1−ζ0)(1−α))((1−α)2b2+h(g−ζ0)+hb(1−α)(1−g))hb2(1−α)2(−g+ζ0−(1−α)(1−g)b)((1−α)2b2+h(g−ζ0+(1−α)(1−g)b))2), | (5.5) |
is of rank 2. So, to apply the OGY method for the HCV model (1.6), we choose ζ as a control parameter. Now, if ζ−ζ0=−T(ut−ˆuvt−ˆv), where T=(^k1 ^k2), then model (5.2) can be written as
(ut+1−ˆuvt+1−ˆv)≈(J(ˆu,ˆv,ζ0)−PT)(ut−ˆuvt−ˆv). | (5.6) |
Thus, the corresponding controlled model of (1.6) is given by
ut+1=(1+h(1−g))ut−h(1+(1−α)b)utvt1+hut,vt+1=(1+h(1−(ζ0−^k1(ut−ˆu)−^k2(vt−ˆv))))vt−h(1−(1−α)b)utvt1+hvt. | (5.7) |
By the theory of stability, PIFP is a sink iff roots of J(ˆu,ˆv,ζ0)−PT satisfying |λ1,2|<1 where its Jacobian is
J(ˆu,ˆv,ζ0)−PT= ((1−α)2b2(1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1)−h(g+b(1−α)(1−g)−ζ0)(−^k1+1−b(1−α))(1−α)2b2+h(g−ζ0)+hb(1−α)(1−g)h(1+b(1−α))(g+b(1−α)(1−ζ0)−ζ0)(1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1)(1−α)2b2+h(g−ζ0+b(1−α)(1−g))^k2(1−α)2b2+h(g−ζ0)+hb(1−α)(1−g)). | (5.8) |
J(ˆu,ˆv,ζ0)−PT has the following characteristic equation
λ2−{(1−α)3b3(2b(1−α)−h(g−ζ))+h(g+(1−α)b(1−g)−ζ0)((1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1))^k2}{((1−α)2b2+h(ζ0−g)+hb(ζ0−1)(1−α))((1−α)2b2+h(g−ζ0)+hb(1−g)(1−α))}λ+{(1−α)4b4−h2((1−α)b+1)(g+b(1−ζ0)(1−α)−ζ0)(g+b(1−α)×(1−g)−ζ0)(−1+b(1−α)+^k1)+h(1−α)2b2(g−ζ0+(1−g)(1−α)b)^k2}((1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1))((1−α)2b2+h(g−ζ0)+hb(1−α)(1−g))=0. | (5.9) |
If one denotes roots of (5.9) are λ1,2, then
λ1+λ2={(1−α)3b3(2b(1−α)−h(g−ζ))+h(g+b(1−α)(1−g)−ζ0)×((1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1))^k2}{((1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1))×((1−α)2b2+h(g−ζ0)+hb(1−α)(1−g))}, | (5.10) |
and
λ1λ2={(1−α)4b4−h2(b(1−α)+1)(g+b(1−ζ0)(1−α)−ζ0)(g+b(1−α)×(1−g)−ζ0)(−1+b(1−α)+^k1)+h(1−α)2b2(g−ζ0+(1−g)(1−α)b)^k2}{((1−α)2b2+h(ζ0−g)+hb(1−α)(ζ0−1))×((1−α)2b2+h(g−ζ0)+hb(1−α)(1−g))}. | (5.11) |
Now, to get lines of marginal stability for (5.7), one requires λ1=±1 and λ1λ2=1, where these restrictions give characteristics roots satisfying |λ1,2|<1. If λ1λ2=1, then from (5.11), we have
L1: h2((1−α)b+1)(g+b(1−ζ0)(1−α)−ζ0)(g+b(1−g)(1−α)−ζ0)^k1−h(1−α)2b2(g−ζ0+(1−g)(1−α)b)^k2+((1−α)2b2+h(ζ0−g)+hb(ζ0−1)(1−α))×((1−α)2b2+h(g−ζ0)+hb(1−g)(1−α))+b4(1−α)4−h2(1−(1−α)2b2)×(g+b(1−ζ0)(1−α)−ζ0)(g+b(1−g)(1−α)−ζ0)=0. | (5.12) |
If λ1=1, then from (5.10) and (5.11), we have
L2: h2((1−α)b+1)(g+(1−α)b(1−ζ0)−ζ0)(g+(1−α)b(1−g)−ζ0)^k1−h2(g−ζ0+(1−g)(1−α)b)(g+b(1−α)(1−ζ0)−ζ0)^k2+h2(1−α)2b2(g−ζ0+(1−g)(1−α)b)(g+b(1−α)(1−ζ0)−ζ0)=0. | (5.13) |
Finally, if λ1=−1, then from (5.10) and (5.11), we have
L3: h2((1−α)b+1)(g+(1−α)b(1−ζ0)−ζ0)(g+(1−α)b(1−g)−ζ0)^k1−h(g−ζ0+(1−g)(1−α)b)(2(1−α)2b2−h(g+b(1−α)(1−ζ0)−ζ0))^k2+b2(1−α)2(−4b2(1−α)2+hb2(1−α)2(1−g)(1−ζ0)+h2(g−ζ0)2−bh(1−α)(g−ζ0)(−2+h(−2+g+ζ0)))=0. | (5.14) |
Therefore, the stable eigenvalues are located within the triangular region (^k1^k2)-plane bounded by Li(i=1,2,3) for parametric values b, α, g, ζ, and h.
The following two cases are to be considered for the completeness of this section.
Case Ⅰ. If b=12.01,α=0.665,g=0.1232,ζ=2.002,h=[0.55,1] with the initial condition (u0,v0)=(0.012352,0.232), then at h=0.8850211364486024, the HCV model (1.6) undergoes the Neimark-Sacker bifurcation, and Figure 3 shows the maximum Lyapunov exponent with bifurcation diagrams. Further, if b=12.01, α=0.665, g=0.1232, ζ=2.002, h=0.8850211364486024, then model (1.6) has PIFP=(0.36511216748407743,0.1018618715629854) and, moreover, from (3.15), one gets
J|PIFP=(0.36511216748407743,0.1018618715629854)=(0.7557825001858909−1.22678997769120480.25001582446986680.9173050343262055), | (6.1) |
with
λ2−1.6730875345120961λ+1=0. | (6.2) |
The roots are λ1,2=0.8365437672560481±0.5479001053705491ι with |λ1,2|=1, which implies from the respective bifurcation curve (4.7) that (b,α,g,ζ,h) = (12.01,0.665,0.1232,2.002,0.8850211364486024)∈N|PIFP=(0.36511216748407743,0.1018618715629854). Furthermore, by fixing b=12.01, α=0.665, g=0.1232, ζ=2.002 and varying h=0.8850211364486024, then one can conclude that the corresponding PIFP is a stable focus. To verify this, if h=0.795<0.8850211364486024, then Figure 4(a) indicates that PIFP is a stable focus. Similarly, if h=0.805, 0.835, 0.855, 0.865, 0.875<0.8850211364486024, then Figures 4(b)–(f) also show that the corresponding PIFP is stable focus. Conversely, if h>0.8850211364486024, then respective PIFP is an unstable focus, and as a result, supercritical N-S bifurcation takes place for the HCV model (1.6). For instance, if b=12.01,α=0.665,g=0.1232, and ζ=2.002, then from (4.18), one gets d|λ1,2|dϵ|ϵ=0=0.6070024854770324>0, and if h=0.895>0.8850211364486024, then from (4.22) and (4.25), one gets
l11=0.7537070762796483,l12=−1.2372155583706284,l13=−1.0168530987328883,l14=−2.5540045067599264,l21=0.252599297693408,l22=0.916450527496516,l23=2.272634068707745,l24=−1.50338800913407, | (6.3) |
η=0.8350788018880821,ζ=0.39028624560846026. | (6.4) |
From (6.3), (6.4), and (4.27), one gets
σ11=1.0502427204029419,σ12=0.9967928302104191,σ21=0.8306993329208056,σ22=−2.848581027486128,σ23=0.5867516617777137. | (6.5) |
From (6.5) and (4.31), one gets
ρ02=−0.4495845767707965+0.6035609562272346ι,ρ11=0.5251213602014709+0.7087254973492596ι,ρ20=−0.4495845767707965−0.1882112897668318ι,ρ21=0. | (6.6) |
Now, utilizing (6.6) with λ,ˉλ=0.8350788018880821±0.39028624560846026ι into (4.29), one obtains ψ=−0.6734309764490927<0, indicating the existence of a closed invariant curve. Thus, at the PIFP, the HCV model undergoes a supercritical N-S bifurcation (see Figure 5(a). Similarly, for h=0.905, 0.915, 0.925, 0.935, 0.945>0.8850211364486024, stable invariant curves also emerge, as depicted in Figures 5(b)–(f). It is worth highlighting that similar bifurcation behavior, as depicted in Figures 4 and 5, has been observed in another system [29].
Case Ⅱ. If one reconsiders the data from Case Ⅰ with ζ0=0.0001, then (5.7) becomes
ut+1=2.652824992ut−9.46934126775utvt1+0.885065ut,vt+1={5.699211267749999utvt+(1+0.885065(0.9999+^k1(0.25612894566782324+ut)+^k2(−0.22553255211518677+vt)))vt}1+0.885065vt, | (6.7) |
where
JPIFP−PT= (1.9335617021802910.9019145397657196+0.29831628483824885^k14.6896071766473661.+0.4251435203530145^k2), | (6.8) |
with lines of marginal stability being
L1:1272.0324493907451+270.18749966926106^k1−111.40041843372863^k2=0, | (6.9) |
L2:870.6576949546483+270.18749966926106^k1−53.786317829588036^k2=0, | (6.10) |
and
L3:−530.511530819259+270.18749966926106^k1−169.01451903786923^k2=0. | (6.11) |
Therefore, the lines (6.9)–(6.11) determine a triangular region that gives |λ1,2|<1 (see Figure 6).
This work investigates the dynamical characteristics of a discrete two-dimensional Hepatitis C virus (HCV) model (1.6), specifically focusing on the local behavior at the fixed points of the discrete model. By examining these fixed points, we gain insights into the stability and potential bifurcations that may arise within the system. The study particularly emphasizes the absence of period-doubling bifurcations at various fixed points under specific conditions. For example, no period-doubling bifurcation occurs at liver-free fixed point (LFFP), disease-free fixed point (DFFP), and total infection fixed point (TIFP). Furthermore, the study proves that at the PIFP, the HCV model undergoes a N-S bifurcation. Additionally, no period-doubling bifurcation exists at the PIFP. To manage the chaos induced by the Neimark-Sacker bifurcation, the OGY control strategy is applied, showcasing the practical implications of controlling chaotic behavior in the HCV model (1.6). Numerical simulations confirm the theoretical results, reinforcing the validity and applicability of the discrete HCV model (1.6) in understanding the dynamics of Hepatitis C virus. The discrete HCV model is crucial for comprehending the complexities of HCV dynamics. Local stability analysis at various fixed points, such as the liver-free fixed point, disease-free fixed point, total infection fixed point, and partial infection fixed point, provides a comprehensive understanding of the virus's behavior under different conditions. The study's exploration of Neimark-Sacker bifurcation at the partial infection fixed point offers significant insights into the potential for chaotic dynamics in the virus's spread. Understanding the existence of chaos in the model underscores the importance of this research in devising effective control strategies for Hepatitis C. It is noted here that the existence of a Neimark-Sacker bifurcation at a partial infection fixed point in a discrete HCV model (1.6) indicates a transition from stable periodic behavior to quasi-periodic or chaotic dynamics as a key parameter crosses a critical threshold. Biologically, this implies that the HCV infection dynamics can shift from predictable oscillations of viral load and immune response to more complex and potentially erratic patterns. Such a bifurcation may correspond to scenarios where the immune system's response becomes less synchronized with the viral replication cycle, leading to irregular flare-ups and remissions of the infection. Understanding this transition is crucial for effective treatment strategies, as it suggests that minor changes in treatment or the patient's immune status could lead to significant changes in infection dynamics. Chaos control techniques, such as parameter adjustments or therapeutic interventions, are essential to stabilize the system, preventing chaotic outbreaks and maintaining the infection at manageable levels. This understanding provides insights into optimizing treatment regimens to avoid triggering complex and less controllable infection states. Future work can build upon this study by exploring more complex models and control strategies. Previous research, such as Blackmore et al.'s [29,30] work on discrete dynamical modeling and analysis of the R-S flip-flop circuit and their study on sigma map dynamics and bifurcations, provide a foundation for further investigations. Additionally, Rahman and Blackmore's [31] examination of bifurcations in walking droplet dynamics and Joshi et al.'s [32] exploration of generalized attracting horseshoes and chaotic strange attractors, suggest avenues for future research. Studies like Murthy et al. [33] on generalized attracting horseshoes in the Rössler attractor also offer valuable insights that can be integrated into future models of HCV dynamics. These directions can enhance our understanding of the complex behavior of the HCV model and improve strategies for managing and controlling the spread of Hepatitis C.
Abdul Qadeer Khan: Conceptualization, formal Analysis, investigation, methodology, resources, software, supervision, validation, visualization, writing-original draft, writing-review and editing; Ayesha Yaqoob: Formal Analysis, investigation, writing-original draft, writing-review and editing; Ateq Alsaadi: Conceptualization, funding acquisition, resources, writing-original draft, writing-review and editing. All authors have read and approved the final version of the manuscript for publication.
The authors extend their appreciation to Taif University, Saudi Arabia, for supporting this work through project number (TU-DSPP-2024-259).
The authors declare that they have no conflict of interest.
λ1,2|(3.21)=b3(1−α)3(2b(1−α)−h(g−ζ))2((1−α)2b2−h(g−ζ)−hb(1−α)(1−ζ))((1−α)2b2+h(g−ζ)+hb(1−α)(1−g))±2(√64(1−b2(1−α)2)(g−ζ)2(g+b(1−g)(1−α)−ζ)2(g+b(1−α)(1−ζ)−ζ)2−16b2(−1+b2(1−α)2)(1−α)2(g−ζ)(g+b(1−g)(1−α)−ζ)(g+b(1−α)(1−ζ)−ζ)(−g+ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)})+b4(1−α)4(9(g−ζ)2+8b2(1−g)(1−α)2(1−ζ)+8b(−1+α)(g−ζ)(−2+g+ζ))(−g+ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)})2+2b4(1−α)4(g−ζ)(−g+ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)})3{(−g+ζ+√{5(g−ζ)2+4b2(1−g)(1−α)2(1−ζ)−4b(g−ζ)(1−α)(−2+g+ζ)})4(16(g−ζ)2+16b2(1−g)(1−α)2(1−ζ)+16b(−1+α)(g−ζ)(−2+g+ζ)−4b2(1−α)2(g−ζ)(g−ζ−√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)2(g−ζ)(−2+g+ζ)})−b2(1−α)2(−g+ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)}})≠1 or −1 | (A.1) |
λ1,2|(3.22)=(1−α)3b3(2b(1−α)−h(g−ζ))2((1−α)2b2−h(g−ζ)−hb(1−α)(1−ζ))((1−α)2b2+h(g−ζ)+hb(1−α)(1−g))±2(√64(1−b2(1−α)2)(g−ζ)2(g+b(1−g)(1−α)−ζ)2(g+b(1−α)(1−ζ)−ζ)2+16b2(−1+b2(1−α)2)(1−α)2(g−ζ)(g+b(1−g)(1−α)−ζ)(g+b(1−α)(1−ζ)−ζ)(−g+ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)})+b4(1−α)4(9(g−ζ)2+8b2(1−g)(1−α)2(1−ζ)+8b(−1+α)(g−ζ)(−2+g+ζ))(−g+ζ+√{5(g−ζ)2+4b2(1−g)(1−α)2(1−ζ)−4b(g−ζ)(1−α)(−2+g+ζ)})2−2b4(1−α)4(g−ζ)(−g+ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)})3{(−g+ζ+√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)})4(16(g−ζ)2+16b2(1−g)(1−α)2(1−ζ)+16b(−1+α)(g−ζ)(−2+g+ζ)−4b2(1−α)2(g−ζ)(g−ζ−√{5(g−ζ)2+4b2(1−α)2(1−g)(1−ζ)−4b(1−α)(g−ζ)(−2+g+ζ)})−b2(1−α)2(−g+ζ+√{5(g−ζ)2+4b2(1−g)(1−α)2(1−ζ)−4b(g−ζ)(1−α)(−2+g+ζ)}})≠1 or −1 | (A.2) |
[1] |
Hu C, Ashok D, Nisbet DR, et al. (2019) Bioinspired surface modification of orthopedic implants for bone tissue engineering. Biomaterials 119366. doi: 10.1016/j.biomaterials.2019.119366
![]() |
[2] |
Karsenty G, Olson EN (2016) Bone and muscle endocrine functions: unexpected paradigms of inter-organ communication. Cell 164: 1248-1256. doi: 10.1016/j.cell.2016.02.043
![]() |
[3] |
Rossi M, Battafarano G, Pepe J, et al. (2019) The endocrine function of osteocalcin regulated by bone resorption: A lesson from reduced and increased bone mass diseases. Int J Mol Sci 20: 4502. doi: 10.3390/ijms20184502
![]() |
[4] |
Loebel C, Burdick JA (2018) Engineering stem and stromal cell therapies for musculoskeletal tissue repair. Cell Stem Cell 22: 325-339. doi: 10.1016/j.stem.2018.01.014
![]() |
[5] |
Dimitriou R, Jones E, McGonagle D, et al. (2011) Bone regeneration: current concepts and future directions. BMC Med 9: 66. doi: 10.1186/1741-7015-9-66
![]() |
[6] | Nordin M, Frankel VH (2001) Basic Biomechanics of the Musculoskeletal System, 3 Eds USA: Lippincott Williams & Wilkins. |
[7] |
Kobayashi S, Takahashi HE, Ito A, et al. (2003) Trabecular minimodeling in human iliac bone. Bone 32: 163-169. doi: 10.1016/S8756-3282(02)00947-X
![]() |
[8] |
Bartl R, Bartl C (2019) Control and regulation of bone remodelling. The Osteoporosis Manual Cham: Springer, 31-39. doi: 10.1007/978-3-030-00731-7_4
![]() |
[9] |
Kenkre JS, Bassett JHD (2018) The bone remodelling cycle. Ann Clin Biochem 55: 308-327. doi: 10.1177/0004563218759371
![]() |
[10] |
Prendergast PJ, Huiskes R (1995) The biomechanics of Wolff's law: recent advances. Irish J Med Sci 164: 152-154. doi: 10.1007/BF02973285
![]() |
[11] |
Wegst UGK, Bai H, Saiz E, et al. (2015) Bioinspired structural materials. Nat Mater 14: 23-36. doi: 10.1038/nmat4089
![]() |
[12] |
Reznikov N, Shahar R, Weiner S (2014) Bone hierarchical structure in three dimensions. Acta Biomater 10: 3815-3826. doi: 10.1016/j.actbio.2014.05.024
![]() |
[13] |
Weiner S, Wagner HD (1998) The material bone: structure-mechanical function relations. Annu Rev Mater Sci 28: 271-298. doi: 10.1146/annurev.matsci.28.1.271
![]() |
[14] |
Recker RR, Kimmel DB, Dempster D, et al. (2011) Issues in modern bone histomorphometry. Bone 49: 955-964. doi: 10.1016/j.bone.2011.07.017
![]() |
[15] |
Eriksen EF, Vesterby A, Kassem M, et al. (1993) Bone remodeling and bone structure. Physiology and Pharmacology of Bone Heidelberg: Springer, 67-109. doi: 10.1007/978-3-642-77991-6_2
![]() |
[16] |
Augat P, Schorlemmer S (2006) The role of cortical bone and its microstructure in bone strength. Age Ageing 35: ii27-ii31. doi: 10.1093/ageing/afl081
![]() |
[17] |
Kozielski M, Buchwald T, Szybowicz M, et al. (2011) Determination of composition and structure of spongy bone tissue in human head of femur by Raman spectral mapping. J Mater Sci: Mater Med 22: 1653-1661. doi: 10.1007/s10856-011-4353-0
![]() |
[18] |
Cross LM, Thakur A, Jalili NA, et al. (2016) Nanoengineered biomaterials for repair and regeneration of orthopedic tissue interfaces. Acta Biomater 42: 2-17. doi: 10.1016/j.actbio.2016.06.023
![]() |
[19] |
Zebaze R, Seeman E (2015) Cortical bone: a challenging geography. J Bone Miner Res 30: 24-29. doi: 10.1002/jbmr.2419
![]() |
[20] |
Liu Y, Luo D, Wang T (2016) Hierarchical structures of bone and bioinspired bone tissue engineering. Small 12: 4611-4632. doi: 10.1002/smll.201600626
![]() |
[21] |
Brodsky B, Persikov AV (2005) Molecular structure of the collagen triple helix. Adv Protein Chem 70: 301-339. doi: 10.1016/S0065-3233(05)70009-7
![]() |
[22] |
Cui FZ, Li Y, Ge J (2007) Self-assembly of mineralized collagen composites. Mater Sci Eng R Rep 57: 1-27. doi: 10.1016/j.mser.2007.04.001
![]() |
[23] |
Wang Y, Azaïs T, Robin M, et al. (2012) The predominant role of collagen in the nucleation, growth, structure and orientation of bone apatite. Nat Mater 11: 724-733. doi: 10.1038/nmat3362
![]() |
[24] | Bentmann A, Kawelke N, Moss D, et al. (2010) Circulating fibronectin affects bone matrix, whereas osteoblast fibronectin modulates osteoblast function. J Bone Miner Res 25: 706-715. |
[25] |
Szweras M, Liu D, Partridge EA, et al. (2002) α2-HS glycoprotein/fetuin, a transforming growth factor-β/bone morphogenetic protein antagonist, regulates postnatal bone growth and remodeling. J Biol Chem 277: 19991-19997. doi: 10.1074/jbc.M112234200
![]() |
[26] |
Boskey AL, Robey PG (2013) The regulatory role of matrix proteins in mineralization of bone. Osteoporosis, 4 Eds Academic Press, 235-255. doi: 10.1016/B978-0-12-415853-5.00011-X
![]() |
[27] |
Boskey AL (2013) Bone composition: relationship to bone fragility and antiosteoporotic drug effects. Bonekey Rep 2: 447. doi: 10.1038/bonekey.2013.181
![]() |
[28] |
Stock SR (2015) The mineral–collagen interface in bone. Calcified Tissue Int 97: 262-280. doi: 10.1007/s00223-015-9984-6
![]() |
[29] |
Nikel O, Laurencin D, McCallum SA, et al. (2013) NMR investigation of the role of osteocalcin and osteopontin at the organic–inorganic interface in bone. Langmuir 29: 13873-13882. doi: 10.1021/la403203w
![]() |
[30] |
He G, Dahl T, Veis A, et al. (2003) Nucleation of apatite crystals in vitro by self-assembled dentin matrix protein 1. Nat Mater 2: 552-558. doi: 10.1038/nmat945
![]() |
[31] |
Clarke B (2008) Normal bone anatomy and physiology. Clin J Am Soc Nephro 3: S131-S139. doi: 10.2215/CJN.04151206
![]() |
[32] |
Olszta MJ, Cheng X, Jee SS, et al. (2007) Bone structure and formation: A new perspective. Mater Sci Eng R Rep 58: 77-116. doi: 10.1016/j.mser.2007.05.001
![]() |
[33] |
Nair AK, Gautieri A, Chang SW, et al. (2013) Molecular mechanics of mineralized collagen fibrils in bone. Nature Commun 4: 1724. doi: 10.1038/ncomms2720
![]() |
[34] |
Landis WJ (1995) The strength of a calcified tissue depends in part on the molecular structure and organization of its constituent mineral crystals in their organic matrix. Bone 16: 533-544. doi: 10.1016/8756-3282(95)00076-P
![]() |
[35] |
Hunter GK, Hauschka PV, POOLE RA, et al. (1996) Nucleation and inhibition of hydroxyapatite formation by mineralized tissue proteins. Biochem J 317: 59-64. doi: 10.1042/bj3170059
![]() |
[36] |
Oikeh I, Sakkas P, Blake D P, et al. (2019) Interactions between dietary calcium and phosphorus level, and vitamin D source on bone mineralization, performance, and intestinal morphology of coccidia-infected broilers. Poult Sci 11: 5679-5690. doi: 10.3382/ps/pez350
![]() |
[37] |
Boyce BF, Rosenberg E, de Papp AE, et al. (2012) The osteoclast, bone remodelling and treatment of metabolic bone disease. Eur J Clin Invest 42: 1332-1341. doi: 10.1111/j.1365-2362.2012.02717.x
![]() |
[38] |
Teitelbaum SL (2000) Bone resorption by osteoclasts. Science 289: 1504-1508. doi: 10.1126/science.289.5484.1504
![]() |
[39] |
Yoshida H, Hayashi SI, Kunisada T, et al. (1990) The murine mutation osteopetrosis is in the coding region of the macrophage colony stimulating factor gene. Nature 345: 442-444. doi: 10.1038/345442a0
![]() |
[40] |
Roodman GD (2006) Regulation of osteoclast differentiation. Ann NY Acad Sci 1068: 100-109. doi: 10.1196/annals.1346.013
![]() |
[41] |
Martin TJ (2013) Historically significant events in the discovery of RANK/RANKL/OPG. World J Orthop 4: 186-197. doi: 10.5312/wjo.v4.i4.186
![]() |
[42] |
Coetzee M, Haag M, Kruger MC (2007) Effects of arachidonic acid, docosahexaenoic acid, prostaglandin E2 and parathyroid hormone on osteoprotegerin and RANKL secretion by MC3T3-E1 osteoblast-like cells. J Nutr Biochem 18: 54-63. doi: 10.1016/j.jnutbio.2006.03.002
![]() |
[43] |
Steeve KT, Marc P, Sandrine T, et al. (2004) IL-6, RANKL, TNF-alpha/IL-1: interrelations in bone resorption pathophysiology. Cytokine Growth F R 15: 49-60. doi: 10.1016/j.cytogfr.2003.10.005
![]() |
[44] |
Mellis DJ, Itzstein C, Helfrich M, et al. (2011) The skeleton: a multi-functional complex organ. The role of key signalling pathways in osteoclast differentiation and in bone resorption. J Endocrinol 211: 131-143. doi: 10.1530/JOE-11-0212
![]() |
[45] | Silva I, Branco J (2011) Rank/Rankl/opg: literature review. Acta Reumatol Port 36: 209-218. |
[46] |
Martin TJ, Sims NA (2015) RANKL/OPG; Critical role in bone physiology. Rev Endocr Metab Dis 16: 131-139. doi: 10.1007/s11154-014-9308-6
![]() |
[47] |
Wang Y, Qin QH (2012) A theoretical study of bone remodelling under PEMF at cellular level. Comput Method Biomec 15: 885-897. doi: 10.1080/10255842.2011.565752
![]() |
[48] |
Weitzmann MN, Pacifici R (2007) T cells: unexpected players in the bone loss induced by estrogen deficiency and in basal bone homeostasis. Ann NY Acad Sci 1116: 360-375. doi: 10.1196/annals.1402.068
![]() |
[49] |
Duong LT, Lakkakorpi P, Nakamura I, et al. (2000) Integrins and signaling in osteoclast function. Matrix Biol 19: 97-105. doi: 10.1016/S0945-053X(00)00051-2
![]() |
[50] |
Stenbeck G (2002) Formation and function of the ruffled border in osteoclasts. Semin Cell Dev Biol 13: 285-292. doi: 10.1016/S1084952102000587
![]() |
[51] |
Jurdic P, Saltel F, Chabadel A, et al. (2006) Podosome and sealing zone: specificity of the osteoclast model. Eur J Cell Biol 85: 195-202. doi: 10.1016/j.ejcb.2005.09.008
![]() |
[52] |
Väänänen HK, Laitala-Leinonen T (2008) Osteoclast lineage and function. Arch Biochem Biophys 473: 132-138. doi: 10.1016/j.abb.2008.03.037
![]() |
[53] | Vaananen HK, Zhao H, Mulari M, et al. (2000) The cell biology of osteoclast function. J cell Sci 113: 377-381. |
[54] |
Sabolová V, Brinek A, Sládek V (2018) The effect of hydrochloric acid on microstructure of porcine (Sus scrofa domesticus) cortical bone tissue. Forensic Sci Int 291: 260-271. doi: 10.1016/j.forsciint.2018.08.030
![]() |
[55] |
Delaissé JM, Engsig MT, Everts V, et al. (2000) Proteinases in bone resorption: obvious and less obvious roles. Clin Chim Acta 291: 223-234. doi: 10.1016/S0009-8981(99)00230-2
![]() |
[56] |
Logar DB, Komadina R, Preželj J, et al. (2007) Expression of bone resorption genes in osteoarthritis and in osteoporosis. J Bone Miner Metab 25: 219-225. doi: 10.1007/s00774-007-0753-0
![]() |
[57] |
Lorget F, Kamel S, Mentaverri R, et al. (2000) High extracellular calcium concentrations directly stimulate osteoclast apoptosis. Biochem Bioph Res Co 268: 899-903. doi: 10.1006/bbrc.2000.2229
![]() |
[58] |
Nesbitt SA, Horton MA (1997) Trafficking of matrix collagens through bone-resorbing osteoclasts. Science 276: 266-269. doi: 10.1126/science.276.5310.266
![]() |
[59] |
Xing L, Boyce BF (2005) Regulation of apoptosis in osteoclasts and osteoblastic cells. Biochem Bioph Res Co 328: 709-720. doi: 10.1016/j.bbrc.2004.11.072
![]() |
[60] |
Hughes DE, Wright KR, Uy HL, et al. (1995) Bisphosphonates promote apoptosis in murine osteoclasts in vitro and in vivo. J Bone Miner Res 10: 1478-1487. doi: 10.1002/jbmr.5650101008
![]() |
[61] |
Choi Y, Arron JR, Townsend MJ (2009) Promising bone-related therapeutic targets for rheumatoid arthritis. Nat Rev Rheumatol 5: 543-548. doi: 10.1038/nrrheum.2009.175
![]() |
[62] |
Harvey NC, McCloskey E, Kanis JA, et al. (2017) Bisphosphonates in osteoporosis: NICE and easy? Lancet 390: 2243-2244. doi: 10.1016/S0140-6736(17)32850-7
![]() |
[63] |
Ducy P, Schinke T, Karsenty G (2000) The osteoblast: a sophisticated fibroblast under central surveillance. Science 289: 1501-1504. doi: 10.1126/science.289.5484.1501
![]() |
[64] |
Katagiri T, Takahashi N (2002) Regulatory mechanisms of osteoblast and osteoclast differentiation. Oral dis 8: 147-159. doi: 10.1034/j.1601-0825.2002.01829.x
![]() |
[65] |
Kretzschmar M, Liu F, Hata A, et al. (1997) The TGF-beta family mediator Smad1 is phosphorylated directly and activated functionally by the BMP receptor kinase. Gene Dev 11: 984-995. doi: 10.1101/gad.11.8.984
![]() |
[66] |
Bennett CN, Longo KA, Wright WS, et al. (2005) Regulation of osteoblastogenesis and bone mass by Wnt10b. P Natl A Sci 102: 3324-3329. doi: 10.1073/pnas.0408742102
![]() |
[67] | Wang Y, Qin QH, Kalyanasundaram S (2009) A theoretical model for simulating effect of parathyroid hormone on bone metabolism at cellular level. Mol Cell Biomech 6: 101-112. |
[68] |
Elefteriou F, Ahn JD, Takeda S, et al. (2005) Leptin regulation of bone resorption by the sympathetic nervous system and CART. Nature 434: 514-520. doi: 10.1038/nature03398
![]() |
[69] |
Proff P, Römer P (2009) The molecular mechanism behind bone remodelling: a review. Clin Oral Invest 13: 355-362. doi: 10.1007/s00784-009-0268-2
![]() |
[70] |
Katsimbri P (2017) The biology of normal bone remodelling. Eur J Cancer Care 26: e12740. doi: 10.1111/ecc.12740
![]() |
[71] |
Fratzl P, Weinkamer R (2007) Nature's hierarchical materials. Prog Mater Sci 52: 1263-1334. doi: 10.1016/j.pmatsci.2007.06.001
![]() |
[72] |
Athanasiou KA, Zhu CF, Lanctot DR, et al. (2000) Fundamentals of biomechanics in tissue engineering of bone. Tissue Eng 6: 361-381. doi: 10.1089/107632700418083
![]() |
[73] |
Takahashi N, Udagawa N, Suda T (1999) A new member of tumor necrosis factor ligand family, ODF/OPGL/TRANCE/RANKL, regulates osteoclast differentiation and function. Biocheml Bioph Res Co 256: 449-455. doi: 10.1006/bbrc.1999.0252
![]() |
[74] |
Nakashima T, Hayashi M, Fukunaga T, et al. (2011) Evidence for osteocyte regulation of bone homeostasis through RANKL expression. Nat Med 17: 1231-1234. doi: 10.1038/nm.2452
![]() |
[75] |
Prideaux M, Findlay DM, Atkins GJ (2016) Osteocytes: the master cells in bone remodelling. Curr Opin Pharmacol 28: 24-30. doi: 10.1016/j.coph.2016.02.003
![]() |
[76] |
Dallas SL, Prideaux M, Bonewald LF (2013) The osteocyte: an endocrine cell… and more. Endocr Rev 34: 658-690. doi: 10.1210/er.2012-1026
![]() |
[77] |
Rochefort GY, Pallu S, Benhamou CL (2010) Osteocyte: the unrecognized side of bone tissue. Osteoporosis Int 21: 1457-1469. doi: 10.1007/s00198-010-1194-5
![]() |
[78] |
Rowe PSN (2012) Regulation of bone–renal mineral and energy metabolism: The PHEX, FGF23, DMP1, MEPE ASARM pathway. Crit Rev Eukaryot Gene Expr 22: 61-86. doi: 10.1615/CritRevEukarGeneExpr.v22.i1.50
![]() |
[79] |
Pajevic PD, Krause DS (2019) Osteocyte regulation of bone and blood. Bone 119: 13-18. doi: 10.1016/j.bone.2018.02.012
![]() |
[80] | Frost HM (1987) The mechanostat: a proposed pathogenic mechanism of osteoporoses and the bone mass effects of mechanical and nonmechanical agents. Bone Miner 2: 73-85. |
[81] |
Tate MLK, Adamson JR, Tami AE, et al. (2004) The osteocyte. Int J Biochem Cell Biol 36: 1-8. doi: 10.1016/S1357-2725(03)00241-3
![]() |
[82] |
Bonewald LF, Johnson ML (2008) Osteocytes, mechanosensing and Wnt signaling. Bone 42: 606-615. doi: 10.1016/j.bone.2007.12.224
![]() |
[83] |
Manolagas SC, Parfitt AM (2010) What old means to bone. Trends Endocrinol Metab 21: 369-374. doi: 10.1016/j.tem.2010.01.010
![]() |
[84] |
Wang Y, Qin QH (2010) Parametric study of control mechanism of cortical bone remodeling under mechanical stimulus. Acta Mech Sinica 26: 37-44. doi: 10.1007/s10409-009-0313-z
![]() |
[85] |
Qu C, Qin QH, Kang Y (2006) A hypothetical mechanism of bone remodeling and modeling under electromagnetic loads. Biomaterials 27: 4050-4057. doi: 10.1016/j.biomaterials.2006.03.015
![]() |
[86] |
Parfitt AM (2002) Targeted and nontargeted bone remodeling: relationship to basic multicellular unit origination and progression. Bone 1: 5-7. doi: 10.1016/S8756-3282(01)00642-1
![]() |
[87] |
Hadjidakis DJ, Androulakis II (2006) Bone remodeling. Ann NYAcad Sci 1092: 385-396. doi: 10.1196/annals.1365.035
![]() |
[88] | Vaananen HK, Zhao H, Mulari M, et al. (2000) The cell biology of osteoclast function. J cell Sci 113: 377-381. |
[89] |
Goldring SR (2015) The osteocyte: key player in regulating bone turnover. RMD Open 1: e000049. doi: 10.1136/rmdopen-2015-000049
![]() |
[90] |
Silver IA, Murrills RJ, Etherington DJ (1988) Microelectrode studies on the acid microenvironment beneath adherent macrophages and osteoclasts. Exp Cell Res 175: 266-276. doi: 10.1016/0014-4827(88)90191-7
![]() |
[91] |
Delaissé JM, Andersen TL, Engsig MT, et al. (2003) Matrix metalloproteinases (MMP) and cathepsin K contribute differently to osteoclastic activities. Microsc Res Techniq 61: 504-513. doi: 10.1002/jemt.10374
![]() |
[92] |
Delaisse JM (2014) The reversal phase of the bone-remodeling cycle: cellular prerequisites for coupling resorption and formation. Bonekey Rep 3: 561. doi: 10.1038/bonekey.2014.56
![]() |
[93] | Bonewald LF, Mundy GR (1990) Role of transforming growth factor-beta in bone remodeling. Clin Orthop Relat R 250: 261-276. |
[94] |
Locklin RM, Oreffo ROC, Triffitt JT (1999) Effects of TGFβ and bFGF on the differentiation of human bone marrow stromal fibroblasts. Cell Biol Int 23: 185-194. doi: 10.1006/cbir.1998.0338
![]() |
[95] |
Lee B, Oh Y, Jo S, et al. (2019) A dual role of TGF-β in human osteoclast differentiation mediated by Smad1 versus Smad3 signaling. Immunol Lett 206: 33-40. doi: 10.1016/j.imlet.2018.12.003
![]() |
[96] |
Koseki T, Gao Y, Okahashi N, et al. (2002) Role of TGF-β family in osteoclastogenesis induced by RANKL. Cell Signal 14: 31-36. doi: 10.1016/S0898-6568(01)00221-2
![]() |
[97] |
Anderson HC (2003) Matrix vesicles and calcification. Curr Rheumatol Rep 5: 222-226. doi: 10.1007/s11926-003-0071-z
![]() |
[98] |
Bellido T, Plotkin LI, Bruzzaniti A (2019) Bone cells. Basic and Applied Bone Biology, 2 Eds Elsevier, 37-55. doi: 10.1016/B978-0-12-813259-3.00003-8
![]() |
[99] |
Weinstein RS, Jilka RL, Parfitt AM, et al. (1998) Inhibition of osteoblastogenesis and promotion of apoptosis of osteoblasts and osteocytes by glucocorticoids. Potential mechanisms of their deleterious effects on bone. J Clin Invest 102: 274-282. doi: 10.1172/JCI2799
![]() |
[100] |
Vezeridis PS, Semeins CM, Chen Q, et al. (2006) Osteocytes subjected to pulsating fluid flow regulate osteoblast proliferation and differentiation. Biochem Bioph Res Co 348: 1082-1088. doi: 10.1016/j.bbrc.2006.07.146
![]() |
[101] |
Lind M, Deleuran B, Thestrup-Pedersen K, et al. (1995) Chemotaxis of human osteoblasts: Effects of osteotropic growth factors. Apmis 103: 140-146. doi: 10.1111/j.1699-0463.1995.tb01089.x
![]() |
[102] |
Russo CR, Lauretani F, Seeman E, et al. (2006) Structural adaptations to bone loss in aging men and women. Bone 38: 112-118. doi: 10.1016/j.bone.2005.07.025
![]() |
[103] |
Ozcivici E, Luu YK, Adler B, et al. (2010) Mechanical signals as anabolic agents in bone. Nat Rev Rheumatol 6: 50-59. doi: 10.1038/nrrheum.2009.239
![]() |
[104] |
Rosa N, Simoes R, Magalhães FD, et al. (2015) From mechanical stimulus to bone formation: a review. Med Eng Phys 37: 719-728. doi: 10.1016/j.medengphy.2015.05.015
![]() |
[105] |
Noble BS, Peet N, Stevens HY, et al. (2003) Mechanical loading: biphasic osteocyte survival and targeting of osteoclasts for bone destruction in rat cortical bone. Am J Physiol-Cell Ph 284: C934-C943. doi: 10.1152/ajpcell.00234.2002
![]() |
[106] |
Robling AG, Castillo AB, Turner CH (2006) Biomechanical and molecular regulation of bone remodeling. Annu Rev Biomed Eng 8: 455-498. doi: 10.1146/annurev.bioeng.8.061505.095721
![]() |
[107] |
Qin QH, Mai YW (1999) A closed crack tip model for interface cracks inthermopiezoelectric materials. Int J Solids Struct 36: 2463-2479. doi: 10.1016/S0020-7683(98)00115-2
![]() |
[108] |
Yu SW, Qin QH (1996) Damage analysis of thermopiezoelectric properties: Part I—crack tip singularities. Theor Appl Fract Mec 25: 263-277. doi: 10.1016/S0167-8442(96)00026-2
![]() |
[109] |
Qin QH, Mai YW, Yu SW (1998) Effective moduli for thermopiezoelectric materials with microcracks. Int J Fracture 91: 359-371. doi: 10.1023/A:1007423508650
![]() |
[110] |
Jirousek J, Qin QH (1996) Application of hybrid-Trefftz element approach to transient heat conduction analysis. Comput Struct 58: 195-201. doi: 10.1016/0045-7949(95)00115-W
![]() |
[111] |
Qin QH (1995) Hybrid-Trefftz finite element method for Reissner plates on an elastic foundation. Comput Method Appl M 122: 379-392. doi: 10.1016/0045-7825(94)00730-B
![]() |
[112] |
Qin QH (1994) Hybrid Trefftz finite-element approach for plate bending on an elastic foundation. Appl Math Model 18: 334-339. doi: 10.1016/0307-904X(94)90357-3
![]() |
[113] |
Qin QH (2013) Mechanics of Cellular Bone Remodeling: Coupled Thermal, Electrical, and Mechanical Field Effects CRC Press. doi: 10.1201/b13728
![]() |
[114] | Wang H, Qin QH (2010) FE approach with Green's function as internal trial function for simulating bioheat transfer in the human eye. Arch Mech 62: 493-510. |
[115] | Qin QH (2003) Fracture analysis of cracked thermopiezoelectric materials by BEM. Electronic J Boundary Elem 1: 283-301. |
[116] |
Qin QH, Ye JQ (2004) Thermoelectroelastic solutions for internal bone remodeling under axial and transverse loads. Int J Solids Struct 41: 2447-2460. doi: 10.1016/j.ijsolstr.2003.12.026
![]() |
[117] |
Qin QH, Qu C, Ye J (2005) Thermoelectroelastic solutions for surface bone remodeling under axial and transverse loads. Biomaterials 26: 6798-6810. doi: 10.1016/j.biomaterials.2005.03.042
![]() |
[118] |
Ducher G, Jaffré C, Arlettaz A, et al. (2005) Effects of long-term tennis playing on the muscle-bone relationship in the dominant and nondominant forearms. Can J Appl Physiol 30: 3-17. doi: 10.1139/h05-101
![]() |
[119] |
Robling AG, Hinant FM, Burr DB, et al. (2002) Improved bone structure and strength after long-term mechanical loading is greatest if loading is separated into short bouts. J Bone Miner Res 17: 1545-1554. doi: 10.1359/jbmr.2002.17.8.1545
![]() |
[120] |
Rubin J, Rubin C, Jacobs CR (2006) Molecular pathways mediating mechanical signaling in bone. Gene 367: 1-16. doi: 10.1016/j.gene.2005.10.028
![]() |
[121] |
Tatsumi S, Ishii K, Amizuka N, et al. (2007) Targeted ablation of osteocytes induces osteoporosis with defective mechanotransduction. Cell Metab 5: 464-475. doi: 10.1016/j.cmet.2007.05.001
![]() |
[122] |
Robling AG, Turner CH (2009) Mechanical signaling for bone modeling and remodeling. Crit Rev Eukar Gene 19: 319-338. doi: 10.1615/CritRevEukarGeneExpr.v19.i4.50
![]() |
[123] |
Galli C, Passeri G, Macaluso GM (2010) Osteocytes and WNT: the mechanical control of bone formation. J Dent Res 89: 331-343. doi: 10.1177/0022034510363963
![]() |
[124] |
Robling AG, Duijvelaar KM, Geevers JV, et al. (2001) Modulation of appositional and longitudinal bone growth in the rat ulna by applied static and dynamic force. Bone 29: 105-113. doi: 10.1016/S8756-3282(01)00488-4
![]() |
[125] |
Burr DB, Milgrom C, Fyhrie D, et al. (1996) In vivo measurement of human tibial strains during vigorous activity. Bone 18: 405-410. doi: 10.1016/8756-3282(96)00028-2
![]() |
[126] |
Sun W, Chi S, Li Y, et al. (2019) The mechanosensitive Piezo1 channel is required for bone formation. Elife 8: e47454. doi: 10.7554/eLife.47454
![]() |
[127] |
Goda I, Ganghoffer JF, Czarnecki S, et al. (2019) Topology optimization of bone using cubic material design and evolutionary methods based on internal remodeling. Mech Res Commun 95: 52-60. doi: 10.1016/j.mechrescom.2018.12.003
![]() |
[128] |
Goda I, Ganghoffer JF (2018) Modeling of anisotropic remodeling of trabecular bone coupled to fracture. Arch Appl Mech 88: 2101-2121. doi: 10.1007/s00419-018-1438-y
![]() |
[129] |
Louna Z, Goda I, Ganghoffer JF, et al. (2017) Formulation of an effective growth response of trabecular bone based on micromechanical analyses at the trabecular level. Arch Appl Mech 87: 457-477. doi: 10.1007/s00419-016-1204-y
![]() |
[130] |
Goda I, Ganghoffer JF (2017) Construction of the effective plastic yield surfaces of vertebral trabecular bone under twisting and bending moments stresses using a 3D microstructural model. ZAMM Z Angew Math Mech 97: 254-272. doi: 10.1002/zamm.201600141
![]() |
[131] |
Qin QH, Wang YN (2012) A mathematical model of cortical bone remodeling at cellular level under mechanical stimulus. Acta Mech Sinica-Prc 28: 1678-1692. doi: 10.1007/s10409-012-0154-z
![]() |
1. | Noppasorn Sutthibutr, Navee Chiadamrong, Kunihiko Hiraishi, Suttipong Thajchayapong, A five-phase combinatorial approach for solving a fuzzy linear programming supply chain production planning problem, 2024, 11, 2331-1916, 10.1080/23311916.2024.2334566 | |
2. | Marjia Haque, Sanjoy Kumar Paul, Ruhul Sarker, Daryl Essam, A novel heuristic approach for planning decentralised supply chain under uncertainties, 2023, 10, 2330-2674, 10.1080/23302674.2023.2258778 | |
3. | Jérémie Schutz, Christophe Sauvey, Teaching of system reliability based on challenging practical works using a spreadsheet software, 2023, 8, 2473-6988, 24764, 10.3934/math.20231263 | |
4. | Noppasorn Sutthibutr, Kunihiko Hiraishi, Navee Chiadamrong, A fuzzy multi-criteria decision-making for optimizing supply chain aggregate production planning based on cost reduction and risk mitigation, 2024, 10, 21998531, 100377, 10.1016/j.joitmc.2024.100377 |