Abstract
In this paper, we propose a degree-based mean-field SIS epidemic model with a saturated function on complex networks. First, we adopt an edge-compartmental approach to lower the dimensions of such a proposed system. Then we give the existence of the feasible equilibria and completely study their stability by a geometric approach. We show that the proposed system exhibits a backward bifurcation, whose stabilities are determined by signs of the tangent slopes of the epidemic curve at the associated equilibria. Our results suggest that increasing the management and the allocation of medical resources effectively mitigate the lag effect of the treatment and then reduce the risk of an outbreak. Moreover, we show that decreasing the average of a network sufficiently eradicates the disease in a region or a country.
1 Introduction
Mathematical modeling plays a crucial role in fighting against large scale infectious disease such as Tuberculosis, HIV, COVID-19, etc., Compartment models have been used to anticipate the progression of diseases and evaluate the effect of interventions on disease spread [1]. One of such models separates the total population into two distinct categories with respect to disease status. People who have not gotten the disease are labeled “susceptibles”; while those who have been infected by a certain disease are called “infectives”. This kind of compartment model is denoted by “an SIS epidemic model” [2–7], which has been extensively used to address the dynamics of those diseases, describing an individual infected by a disease as having no immunity, thus becoming a susceptible again.
Most of the existing models assume that all the individuals are well-mixed and they have homogeneous mixing of surfaces, which implies that each individual has the same probability to contact other individuals and ignores the degree of social heterogeneity induced by age, household, spatial structures, and social spheres, etc. Generally, the social interactions of individuals generate a certain pattern based on social preferences, which contributes to transmission heterogeneity. Indeed, such factors may play a decisive role in the disease transmission and they also may help health policymakers to take more effective control measures for curbing the disease spread [7, 8]. Epidemic models on complex networks incorporate such contact heterogeneity and take account for how the structures of the networks affect the disease prevalence. A popular degree-based SIS epidemic model has been built [9] and it exhibits threshold dynamics [4, 5]. Since then, many factors including vector-borne [10, 11], infective media [12], awareness reaction [13], and the gene diversity of pathogens [14], etc., have been incorporated into the study of the co-evolution of networks and epidemics.
In epidemiology, the basic reproduction number, , denotes the average number of secondary cases produced by a typical disease carrier during their infectious period [15–17]. Generally, if the infection will be eradicated as each infected individual suffers a disease, on average, to less than one other individual; if it is greater than one, then the outbreak will grow exponentially. For a degree-based SIS epidemic model, the expression of the basic reproduction number is given bywith transmission rate β and treatment rate which serves a vital role in mitigating a disease transmission. However, when people face an emerging infectious disease, such as COVID-19, SARS, H1N1 etc., it may lead to a sharp increase in the number of patients and may cause serious runs of medical resources. A delay phenomena occurs since some infected individuals can not be instantly treated in hospitals. Such properties can be addressed by a Hill function bywhere α denotes a hysteretic effect due to medical limitations. Apparently, when I is small enough and it tends towards to when I is large enough. The constant represents the capacity for treatment. Some scholars proposed a type of step function [18] and a Verhulst-type function [19], which have similar properties to .
In view of such epidemiological models incorporating a saturated treatment function, it is not hard to find that most of them enable such models, essentially changing their dynamics. Once a saturated function has been introduced, there always exists a backward bifurcation, which implies that even if some certain control measures make , it does not guarantee the eradication of a disease in a region or a county. The presence of a backward bifurcation indicates that there exists a bistable phenomena, i.e, a stable endemic equilibria and a stable disease-free equilibrium coexist when The epidemic curve converging to zero or a positive constant is definitely determined by initial sizes of an infection. Some other epidemiological mechanisms, including partial immunity by vaccination [20], death caused by the disease [21], and susceptible heterogeneity [22] etc., have been identified to produce a backward bifurcation. However, the existing results are usually shown in homogeneous models. Up to date, there are few results to study backward bifurcation phenomena in epidemic models on complex networks from theoretical view of points. Recently, Li and Yousef proposed two degree-based SIR and SIS models with a saturated function to study a backward bifurcation phenomena. They showed a sufficient condition for the occurrence of a backward bifurcation [23, 24]. Indeed, since network models have higher dimensions than those of well-mixed models, studying the complex dynamics of such models is a challenging issue from a theoretical view of points.
There are three main contributions in this paper. First, we propose a degree-based SIS epidemic model with a saturated function to study its long-term behaviors. Second, to overcome the difficulty of high dimensions for a network, we adopt an edge-based compartmental approach to lower the dimensions of an SIS epidemic model. Such an approach changes a complete degree-based model to a degree-edge-mixed model, and hence, it lowers the dimension of such a model from to Third, we proposed a geometric method to completely characterize the local stability of each equilibrium.
The organization of this paper is as follows: In Section 2, a degree-based SIS epidemic model on complex networks with a saturated function is proposed. Furthermore, we adopt an edge-compartmental approach to rewrite it as a degree-edge-mixed model. Section 3 gives a geometric approach to study the local stability of each equilibrium. In Section 4, we conducted some numerical simulations to illustrate our theoretical results. We give a brief discussion in the last section.
2 Model Formulation
In this paper, we focus on the complete stability of each equilibrium by a novel approach. Let us assume that the maximum contact number of an individual is n and then the degree set is The total subpopulation is categorized into two statuses: susceptible and infective. Let and represent the densities of susceptible nodes and infected nodes at time t and degree , respectively. Inspired by the modeling approach [9], we account for the limitation of medical resources and propose a mean-field SIS epidemic model as follows:where β denotes the transmission rate and γ represents the treatment rate; α stands for the lag effect of the treatment due to the limitation of medical resources. From epidemiological view of points, the termdenotes the probability of a given node connecting to an infected node at time t. Hence, it can be considered as a density of an edge at time t. Following the steps [10], model (1) can be rewritten as a degree-edge-mixed mode.with initial condition
L emma
roof
Observing the equation of system (2), we solve it and obtainwhereObviously, if Now, we will prove the positivity of for and By the continuity of the solution, if we claim that there exists a such that for some . Then for and However, from the first n equation of (2), it follows thatThis leads to a contradiction with the claim. Therefore, for and .Remark 2.1. Lemma 2.1 ensures that the solution of system (2) is strictly positive if, which supports thatfor allandThe strong connectivity of a network guarantees the positivity of the solution as long asor there exists at leastsuch thatNoting that we have from Lemma 2.1 that and for all and Employing the expression of , we assert thatTherefore, for all and , the setis positively invariant associate with system (2) and In what follows, we focus on the long term behaviors of system (2) taking the initial data from .
3 Stability of Equilibria
In this section, we will consider the local stability of system (2) by a novel-geometric approach, which resolves such a matter once and for all. First, we try to give the basic reproduction number in an explicit form by the approach in [25]. Obviously, system (2) always has a disease-free equilibrium . Linearising the equation in system (2) yields toSolving Eq. 5 by a constant variation method, one drives a renew equationSo that the basic reproduction number is calculated in form ofThe epidemiological meaning of is the average number of secondary infected edges produced by a typical edge during its infectious period. Now, let us account for the local stability of disease-free equilibrium .
. If then the disease-free equilibrium is locally asymptotically stable for any
Proof Linearizing system (2) around the disease-free equilibrium results inLet us assume that system (8) has the solution with exponential forms, i.e, and Plugging them into (8), we have thatIf , we cancel on both sides of system (8) to derive that Otherwise, From the first equation of (8), it follows that From what has been discussed, we conclude that if then the disease-free equilibrium is locally stable for any
. If the disease-free equilibrium is globally asymptotically stable if
Proof Let us candidate a Lyapunov function byDifferentiating V along the solution of system (2) leads tohere we have used the fact that and for all Hence, if then While if thenThe equality holds if and only if Hence, the largest invariant set only contains a singleton point Consequently, the Lyapunov-Kasovskii-LaSalle theorem ensures that is globally asymptotically stable if
Now, we are in a position to study the existence of the endemic steady state, whose components of that feasible equilibrium satisfyFrom the first equation of (10), we getSubstituting Eq. 11 into the second equation of (10) and canceling , we have thatApparently, if has a solution , then system (10) has an endemic equilibrium. Note that has the following properties:Consequently, we conclude the following theorem on the existence of system (2).
Theorem 3.3 If then has at least one positive solution In turn, system (2) has at least one positive endemic equilibrium
Proof If then and The Immediate Value Theorem ensures that has at least one positive point in
From Theorem 3.3, we assert that system (2) has the existence of the endemic equilibrium, but it does not guarantee the uniqueness of the positive solution when However, if α is too small or large enough, it is easy to find that the function F is a decreasing function associated with Then we can claim the uniqueness of the endemic equilibrium whereas However, for we do not catch the existence of such equilibria due to the complex structure of To overcome that difficulty, we select a sensitive analysis for model parameters to investigate the existence of a backward bifurcation of that model for
Now, we are concerned with the endemic curve which bifurcates backwards at If a backward bifurcation takes place, there always exist some additional conditions except To assess such a phenomena, we express from Substituting this expression into Eq. 12, one arrives atCalculating the derivative of Eq. 14 with respect to by the Implicit Function Theorem, we obtainwhereAfter a simple computation, we have that
L emma
. Let
be any feasible endemic equilibrium. The following statements are valid.
(1) If, thenis locally asymptotically stable;
(2) If, thenis unstable.
Proof To address the stability of case (1), we rewrite defined by Eq. 21 in the form ofwhere Suppose (6) has a solution with . Noting thatHence, if . Then contradicts with Eq. 6. This implies that Eq. 6 has no solution with nonnegative real parts. Consequently, the endemic equilibrium is locally stable if
On the contrary, if it follows from Eq. 21 thatThis, together with the Intermediate Value Theorem, ensures that Eq. 21 has at least one positive real solution. Hence, the endemic equilibrium is unstable. This completes the proof.
If a backward bifurcation takes place, Lemma 3.4 and Theorem 3.5 give the stability of two positive endemic equilibria.
. If , and , system (2) has two endemic equilibria. The one with a smaller quantitative of infected nodes is unstable; while the other, with a higher value of infected nodes, is locally asymptotically stable.
Theorem 3.7. SupposeIfor α is large enough, then the endemic equilibriumis globally asymptotically stable if
Proof Let us pick up a candidate Lyapunov function bywherewhere with
Taking the derivative of along the trajectory of Eq. 2, we obtainOn the contrary, differentiating along the solution of Eq. 2 leads toAdding Eqs 30, 31 together, one derives thatIf or it is large enough, then the termhere we have used the fact that and are both smaller than one. Therefore, the first and the last terms move towards to zero. Apparently, The equality holds if and only ifConsequently, the largest invariant set of contains only a singleton point . From LaSalle’s invariance principle, it follows that the endemic equilibrium is globally asymptotically stable.
4 Numerical Simulation
In this section, we will proceed with some numerical experiments to validate our theoretical results. We account for an epidemic spreading on a scale-free network. Hence, we assume that the degree distribution of that network is and maximum degree . Then Hence, the average degree, the second movement and the third movement of this network are , and
First, we fix If we take , we calculate and Theorem 3.2 ensures that the disease-free equilibrium is globally asymptotically stable. Let us pick up twenty different initial conditions for system (2). We find that the densities of infected edges associated with those initial data converge to zero (see Figure 1A). Enlarging , we estimate Since Theorem 3.7 ensures that the endemic equilibrium is globally asymptotically stable. Figure 1B displays that all the trajectories of system (2) move towards to a positive constant.
FIGURE 1
Second, we fix and hence We verify α and consider the existence and the uniqueness of the endemic equilibria. If we choose and it is large enough, Theorem 3.7 shows that the endemic equilibrium is globally asymptotically stable. Figure 2A depicts that the densities of infected edges from differential initial data converge to an endemic equilibrium . If we pick up which are mediated values, all the trajectories of system (2) converge to their corresponding positive equilibria. From Figure 2B we find that increasing values of α gradually advances the arrival times and increases the sizes of the associated endemic equilibria. This indicates that enriching the adequate medical resources effectively reduce the risk of an infection.
FIGURE 2
Third, if we fix the structure of the network, then a key value determines whether or not a backward bifurcation happens when If we select and , then From Lemma 3.4, system exhibits a backward bifurcation at From Figure 3, it follows that the trajectories of system (2) partially converge to a positive level or partially move towards to zero. However, when we take however, does not satisfy the conditions of Theorem 3.2. Figure 3B shows that the disease-free equilibrium is still asymptotically stable.
FIGURE 3
Finally, we want further insight into the existence of endemic equilibria when We take the parameters similar to case one except α. In this case, As we know, the endemic is a function with respect to α and . If we change values of α from 4.8 to 6.8 with step 0.5. Figure 4 suggests that there always exists a backward bifurcation; In addition, the lag effect α controls the depth of the backward bifurcation; the larger α, the bigger the depth of the bifurcation. α is bigger in case of strength of lag effects of medical resources. Figure 4B displays the distributions of the positive equilibria in a given region In the blue region, system (2) has only unique disease-free equilibrium , which is globally asymptotically stable; In the grey region, a backward bifurcation occurs, i.e, system (2) exhibits a bistable phenomena, one endemic equilibria with a large value is stable and the other one with a small value is unstable; Moreover, the disease-free equilibrium is locally stable. In the green region, system (2) has a unique endemic equilibrium and it is stable.
FIGURE 4
5 Discussion
In this paper, we considered a mean-field degree-based SIS epidemic model with a saturated treatment function. First, we adopted an edge-compartmental approach to simplify a pure degree-based model to a degree-edge-mixed model. Second, we proposed a novel method-a geometric approach to completely study the stability of each feasible equilibrium. The proposed model exhibits a backward bifurcation, i.e, does not sufficiently guarantee the eradication of an outbreak.
Compared with the results in [5, 25], a degree-based SIS epidemic model on complex networks shows a threshold dynamic, in the sense that, if then the disease-free equilibrium is globally asymptotically stable; Otherwise, the unique endemic equilibrium is globally asymptotically stable. However, in our system, a saturated treatment function radically altered such a threshold property. Lemma 3.4 suggests that the hysteretic effect leads to the occurrence of backward bifurcation.
The basic reproduction number has no relationship with the lag effect of α. However, it is a vital value which determines whether or not a backward bifurcation occurs. In view of most existing results about that phenomena, they usually showed the existence and did not point out the stability of each equilibrium except some simple models with lower dimensions [2, 3, 18, 19]. In this paper, we proposed a geometric approach to completely solve such an issue. Our results suggest that whereas and , one positive equilibrium with a larger value is stable; the other with small data is unstable. Their stabilities definitely depend on those signs of tangent slopes of the epidemic curve To establish it, we build a bridge between the local stability and the derivative of the epidemic curve, which is a universal result for any system even if it has a more complex structure [23].
Generally, the contact magnitude of an outbreak is characterized by the average degree and its heterogeneity is measured by the two movement degrees of a network. To assess how the heterogeneity of a network affects the transmission of an infection, we fixed all the parameters as in part one of Section 4. Then we take the degree distribution function with respectively. Then and 1.7924. As increases, the system undergoes three different phenomena: the extinction of an outbreak, the occurrence of a backward bifurcation, and the persistence of the disease (see Figure 5). This expounds that the system (2) undergoes a phase transition with the change of the network structure.
FIGURE 5
From an epidemiological viewpoint, the occurrence of a backward bifurcation implies that those control measures enabling do not efficiently ensures the eradication of an infection. In this case, the evolution of the disease heavily depends on the initial numbers of infected individuals. If it is small, the disease will die out; Otherwise, it will maintain a positive level and it cannot be radically eradicated from a region or a county. Hence, enhancing the management and allocations of medical resources plays an important role in slowing down a disease prevalence. Those control measures are in favor of effectively reducing the lag effect due to the limitations of medical resources as long as people face an emerging disease.
However, there are some limitations of this paper. First, we do not incorporate the population demography into the modeling process because introducing birth and death nodes essentially changes the topology of a network [27, 28]. This makes the model become too complex and then it has been become an unresolved issue to analyze its dynamical behaviors from mathematical view of points. Second, we do not couple individual contact data with some reported data for a realistic disease to study its evolutionary behaviors [29]. Third, we do not consider the convolution of information spread and epidemic transmission on multi-layered networks [30]. To carry out such a project, it may provide some valuable control suggestions for policymakers and public health government. We leave these for our future works.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
JY and XW conceived of the presented idea. JY developed the theory and performed the computations. XW verified the analytical methods. All authors discussed the results and contributed to the final manuscript. JY designed the modeling process and analyzed theoretical results. XW conceived of the study and helped to draft the manuscript. All the authors read and approved the final manuscript.
Funding
This work is partially supported by the National Natural Science Foundation of China (No.12001339, No.61573016), and the Shanxi Province Science Foundation for Youths (No. 201901D211413). Shanxi University of Finance and Economics Youth Research Fund Project (QN-2019017).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
KermackWMcKendrickA. A Contribution to the Mathematical Theory of Epidemics. Proc R Soc Lond A (1927) 115:700–21. 10.1098/rspa.1927.0118
2.
van den DriesschePWatmoughJ. A Simple Sis Epidemic Model with a Backward Bifurcation. J Math Biol (2000) 40:525–40. 10.1007/s002850000032
3.
Kribs-ZaletaCMVelasco-HernándezJX. A Simple Vaccination Model with Multiple Endemic States. Math Biosciences (2000) 164:183–201. 10.1016/s0025-5564(00)00003-1
4.
WangLDaiG-z.Global Stability of Virus Spreading in Complex Heterogeneous Networks. SIAM J Appl Math (2008) 68:1495–502. 10.1137/070694582
5.
d’OnofrioA. A Note on the Global Behaviour of the Network-Based Sis Epidemic Model. Nonlinear Anal Real World Appl (2008) 9:1567–72. 10.1016/j.nonrwa.2007.04.001
6.
ZhuLGuanGLiY. Nonlinear Dynamical Analysis and Control Strategies of a Network-Based Sis Epidemic Model with Time Delay. Appl Math Model (2019) 70:512–31. 10.1016/j.apm.2019.01.037
7.
XieYWangZLuJLiY. Stability Analysis and Control Strategies for a New Sis Epidemic Model in Heterogeneous Networks. Appl Math Comput (2020) 383:125381. 10.1016/j.amc.2020.125381
8.
LauMSYDalzielBDFunkSMcClellandATiffanyARileySet alSpatial and Temporal Dynamics of Superspreading Events in the 2014-2015 West Africa Ebola Epidemic. Proc Natl Acad Sci USA (2017) 114:2337–42. 10.1073/pnas.1614595114
9.
Pastor-SatorrasRVespignaniA. Epidemic Spreading in Scale-free Networks. Phys Rev Lett (2001) 86:3200–3. 10.1103/physrevlett.86.3200
10.
WangXYangJ. Dynamical Analysis of a Mean-Field Vector-Borne Diseases Model on Complex Networks: An Edge Based Compartmental Approach. Chaos (2008) 30:013103. 10.1063/1.5116209
11.
WangYJinZYangZZhangZ-KZhouTSunG-Q. Global Analysis of an Sis Model with an Infective Vector on Complex Networks. Nonlinear Anal Real World Appl (2012) 13:543–57. 10.1016/j.nonrwa.2011.07.033
12.
YangMChenGFuX. A Modified Sis Model with an Infective Medium on Complex Networks and its Global Stability. Physica A: Stat Mech its Appl (2011) 390:2408–13. 10.1016/j.physa.2011.02.007
13.
WuQFuXSmallMXuX-J. The Impact of Awareness on Epidemic Spreading in Networks. Chaos (2012) 22:013101. 10.1063/1.3673573
14.
YangJKuniyaTLuoX. Competitive Exclusion in a Multi-Strain Sis Epidemic Model on Complex Networks. Electron J Differ Equat (2019) 2019:1–30.
15.
van den DriesschePWatmoughJ. Reproduction Numbers and Sub-threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Math Biosciences (2002) 180:29–48. 10.1016/s0025-5564(02)00108-6
16.
KeelingMJGrenfellBT. Individual-based Perspectives on R0. J Theor Biol (2000) 203:51–61. 10.1006/jtbi.1999.1064
17.
DiekmannOHeesterbeekJAMetzJA. On the Definition and the Computation of the Basic Reproduction Ratio R0 in Models for Infectious Diseases in Heterogeneous Populations. J Math Biol (1990) 28:365–82. 10.1007/BF00178324
18.
WangW. Backward Bifurcation of an Epidemic Model with Treatment. Math Biosciences (2006) 201:58–71. 10.1016/j.mbs.2005.12.022
19.
CuiJMuXWanH. Saturation Recovery Leads to Multiple Endemic Equilibria and Backward Bifurcation. J Theor Biol (2008) 254:275–83. 10.1016/j.jtbi.2008.05.015
20.
BuonomoBLacitignolaD. On the Backward Bifurcation of a Vaccination Model with Nonlinear Incidence. Namc (2011) 16:30–46. 10.15388/na.16.1.14113
21.
GarbaSMGumelABAbu BakarMR. Backward Bifurcations in Dengue Transmission Dynamics. Math Biosciences (2008) 215:11–25. 10.1016/j.mbs.2008.05.002
22.
GumelAB. Causes of Backward Bifurcations in Some Epidemiological Models. J Math Anal Appl (2012) 395:355–65. 10.1016/j.jmaa.2012.04.077
23.
LiC-HYousefAM. Bifurcation Analysis of a Network-Based Sir Epidemic Model with Saturated Treatment Function. Chaos (2019) 29:033129. 10.1063/1.5079631
24.
HuangY-JLiC-H. Backward Bifurcation and Stability Analysis of a Network-Based Sis Epidemic Model with Saturated Treatment Function. Physica A: Stat Mech its Appl (2019) 527:121407. 10.1016/j.physa.2019.121407
25.
YangJXuF. The Computational Approach for the Basic Reproduction Number of Epidemic Models on Complex Networks. IEEE Access (2019) 7:26474–9. 10.1109/access.2019.2898639
26.
MartchevaM. Methods for Deriving Necessary and Sufficient Conditions for Backward Bifurcation. J Biol Dyn (2019) 13:538–66. 10.1080/17513758.2019.1647359
27.
JinZSunGSunGZhuH. Epidemic Models for Complex Networks with Demographics. Math Biosci Engin (2014) 11:1295–317. 10.3934/mbe.2014.11.1295
28.
YinQWangZXiaCDehmerMEmmert-StreibFJinZ. A Novel Epidemic Model Considering Demographics and Intercity Commuting on Complex Dynamical Networks. Appl Math Comput (2020) 386:125517. 10.1016/j.amc.2020.125517
29.
YangJWangGWangGZhangS. Impact of Household Quarantine on SARS-Cov-2 Infection in mainland China: A Mean-Field Modelling Approach. Math Biosci Engin (2020) 17:4500–12. 10.3934/mbe.2020248
30.
WangZXiaC. Co-evolution Spreading of Multiple Information and Epidemics on Two-Layered Networks under the Influence of Mass media. Nonlinear Dyn (2020) 102:3039–52. 10.1007/s11071-020-06021-7
Summary
Keywords
complex networks, an edge-compartmental approach, a geometric approach, backward bifurcation, global stability
Citation
Wang X and Yang J (2021) A Bistable Phenomena Induced by a Mean-Field SIS Epidemic Model on Complex Networks: A Geometric Approach. Front. Phys. 9:681268. doi: 10.3389/fphy.2021.681268
Received
17 March 2021
Accepted
03 May 2021
Published
17 June 2021
Volume
9 - 2021
Edited by
Mahdi Jalili, RMIT University, Australia
Reviewed by
Gui-Quan Sun, North University of China, China
Chengyi Xia, Tianjin University of Technology, China
Lin Wang, University of Cambridge, United Kingdom
Updates
Copyright
© 2021 Wang and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Junyuan Yang, yjyang66@sxu.edu.cn
This article was submitted to Social Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.