Loading [MathJax]/jax/element/mml/optable/MathOperators.js
Editorial Special Issues

Effects of ionizing radiation in biomolecules, cells and tissue/organs: basic mechanisms and applications for cancer therapy, medical imaging and radiation protection

  • Citation: Francesca Ballarini, Mario P. Carante, Alessia Embriaco, Ricardo L. Ramos. Effects of ionizing radiation in biomolecules, cells and tissue/organs: basic mechanisms and applications for cancer therapy, medical imaging and radiation protection[J]. AIMS Biophysics, 2022, 9(2): 108-112. doi: 10.3934/biophy.2022010

    Related Papers:

    [1] Morten Risstad, Mathias Holand . On the relevance of realized quarticity for exchange rate volatility forecasts. Data Science in Finance and Economics, 2024, 4(4): 514-530. doi: 10.3934/DSFE.2024021
    [2] Manveer Kaur Mangat, Erhard Reschenhofer, Thomas Stark, Christian Zwatz . High-Frequency Trading with Machine Learning Algorithms and Limit Order Book Data. Data Science in Finance and Economics, 2022, 2(4): 437-463. doi: 10.3934/DSFE.2022022
    [3] Kasra Pourkermani . VaR calculation by binary response models. Data Science in Finance and Economics, 2024, 4(3): 350-361. doi: 10.3934/DSFE.2024015
    [4] Obinna D. Adubisi, Ahmed Abdulkadir, Chidi. E. Adubisi . A new hybrid form of the skew-t distribution: estimation methods comparison via Monte Carlo simulation and GARCH model application. Data Science in Finance and Economics, 2022, 2(2): 54-79. doi: 10.3934/DSFE.2022003
    [5] Jin Zeng, Yijia Zhang, Yun Yin, Peter G Moffatt . The effect of the Covid pandemic on stock market volatility: Separating initial impact from time-to-recovery. Data Science in Finance and Economics, 2024, 4(4): 531-547. doi: 10.3934/DSFE.2024022
    [6] Maria Francesca Carfora, Albina Orlando . Application of the Gordon Loeb model to security investment metrics: a proposal. Data Science in Finance and Economics, 2024, 4(4): 601-614. doi: 10.3934/DSFE.2024025
    [7] Yue Yuin Lim, Sie Long Kek, Kok Lay Teo . Efficient state estimation strategies for stochastic optimal control of financial risk problems. Data Science in Finance and Economics, 2022, 2(4): 356-370. doi: 10.3934/DSFE.2022018
    [8] Dominic Joseph . Estimating credit default probabilities using stochastic optimisation. Data Science in Finance and Economics, 2021, 1(3): 253-271. doi: 10.3934/DSFE.2021014
    [9] Chong Li, Guoqiong Long, Shuai Li . Research on measurement and disequilibrium of manufacturing digital transformation: Based on the text mining data of A-share listed companies. Data Science in Finance and Economics, 2023, 3(1): 30-54. doi: 10.3934/DSFE.2023003
    [10] Nitesha Dwarika . Asset pricing models in South Africa: A comparative of regression analysis and the Bayesian approach. Data Science in Finance and Economics, 2023, 3(1): 55-75. doi: 10.3934/DSFE.2023004


  • With the development of society, the number of deaths caused by influenza infection has had gradually increased. For example, the global pandemic in 1918 caused more deaths in the world than in the First World War[1,2]. According to statistical analysis, the global influenza pandemic is caused by H2N2 virus, which is called "Asian influenza" because it first occurred in Asia. The incidence rate was about 15–30% [3,4]. Recently, COVID-19 also has a great impact on human life. According to the latest real-time statistics of the World Health Organization, as of October 10, 2022, the cumulative number of confirmed cases worldwide exceeded 600 million, and the cumulative number of deaths was about 6.5 million[5]. Therefore, in order to better control and reduce the number of deaths, how to effectively and quickly control the spread of disease is worthy of further study.

    In the face of the global pandemic, many countries have responded to and mitigated the impact of influenza on human life through media reports, vaccination, disinfection, the use of protective equipment (such as masks), quarantine and other measures. In addition, some biological mathematicians have also established many mathematical models in order to better understand influenza infection and quantify the effectiveness of various control measures [6,7,8,9,10,11,12,13]. For example, considering social factors such as vaccination, media reports and protective measures, [14,15,16,17,18,19,20,21] established some different forms of disease models to analyze the impact of disease transmission. Practice shows that frequent hand washing, disinfection or wearing protective masks, as well as reducing or avoiding close contact with infected persons, can effectively reduce transmission [22,23,24,25,26,27,28,29]. However, media reports on influenza cases and deaths may have a great impact on the public, because the public will take some protective measures after media reports. In addition, quarantine and other measures will greatly reduce the effective contact rate between vulnerable people and infected people, thus reducing the spread of disease.

    Recently, many different control strategies have been proposed to apply to the control of some disease models. For example, by using multiple optimal control, Ndii and Adi [12] have studied the effects of individual awareness and vector controls on Malaria transmission dynamics. By using a non-smooth control strategy, Li et al.[13] have considered the bifurcations and dynamics of a plant disease system. By using the two thresholds control strategy, Li et al.[30] have considered the global dynamics of a Filippov predator-prey model. Chen et al. [31] have considered a two-thresholds policy for a Filippov model in combating influenza. Zhou et al. [32] have discussed a two-thresholds policy to interrupt transmission of West Nile Virus to birds. Meanwhile, based on the above mentioned transmission factors, some scholars have also done a lot of meaningful work. For example, Dong et al. [33] studied a kind of nonlinear incidence Filippov epidemic model to describe the impact of media in the process of epidemic transmission. They found that choosing an appropriate threshold value and control intensity can prevent the outbreak of infectious diseases, and media coverage can reduce the burden of disease outbreaks and shorten the duration of disease outbreaks. Can et al.[25] have studied a Filippov model describing the effects of media coverage and quarantine on the spread of human influenza. Xiao et al. [28] have discussed a media impact switching surface during an infectious disease outbreak. In real life, when the influenza started to spread, that is, when the epidemic was not serious, the public paid little attention to the influenza. Generally, only when the number of people who have already felt it reaches and exceeds a certain threshold level will the mass media begin to coverage in large numbers and the public begin take some countermeasures [11,17,19,20]. Currently, many control managements have been developed, including threshold [30], harvesting control management[34,35,36], impulsive control management[37,38] and other controls [39,40,41,42,43,44,45]. However, the implementation of this measure may cause some social impacts, such as socioeconomic decline and psychological impact on the quarantined people. It is important to consider when to take these control measures. Therefore, an appropriate threshold policy is needed to deal with the influenza outbreak, or at least reduce the number of infected people to an acceptable level. As far as we know, under logistic source and broken line control strategy, the sliding dynamics and bifurcations of a human influenza system have been seldom reported in existing work.

    Motivation and inspiration come from the discussion above. In this paper, we propose a human influenza system under logistic source and broken line control strategy. The main contributions of this paper include three points: First, choose two thresholds IT and ST as a broken line control strategy: Once the number of infected people exceeds IT, the media influence comes into play, and when the number of susceptible individuals is greater than ST, the control by quarantine of susceptible individuals control is open. Furthermore, by choosing different thresholds IT and ST, the existence and stability of all possible equilibria considered, and then the Filippov system tends to the pseudo-equilibrium on sliding mode domain or one endemic equilibrium or bistability endemic equilibria under some conditions. It is worth pointing out that in our paper there exists a new phenomenon of two real equilibrium points co-existing.

    In fact, from the perspective of control strategy, our paper puts forward a broken line control strategy. Different from the previous single-stage control strategy [30,34,35,37,38], our control strategy is divided into two stages and can control and simulate the real life situation well.

    This paper is structured as follows. In Section 2, we propose a non-smooth model under logistic source and broken line control strategy. By changing the infection threshold values and susceptibility threshold values, in Section 3, we consider the global dynamics of the system under Case 1: ST<S1<S2=S3. In Section 4, we study the global dynamics of the system under Case 2: S1<ST<S2=S3. Section 5 considers the global dynamics of the system under Case 3: S2=S3<ST.The regular/virtulal equilibrium bifurcations are given in Section 6. Finally, we summarize the main results of this paper and discuss some biological conclusions in Section 7.

    In [25], a Filippov human influenza model with effects of media coverage and quarantine was considered. Because a logistic source can better depict the real situation, motivated by the above discussion and the existing multiple threshold conditions [25,30,31,32,45], we consider the logistic source factor with the susceptible population, and then the new human influenza model can be described by

    {St=γS(1SK)βSI,It=βSI(d+δ+r)I,Rt=rIdI, (2.1)

    where S is susceptible individuals, I is infected individuals, and R is recovered individuals. γ is an intrinsic growth rate of susceptible individuals, and K is carrying capacity. β is the transmission rate. d is the natural death rate. δ is the death rate caused by disease, and r is the recovered rate.

    To better control the spread of the disease, we give the following broken line control strategy. If the infected populations is lower than IT, no control is required. When the infected populations is greater than IT, the media coverage control is open. That is, the mass media being coverage of information about the disease, including the route of transmission, the number of infected cases and the number of deaths.

    Therefore, the public realizes the harm of the disease, and they changed their behavior, leading to a decline in the contact rate β. In this paper, we consider that the positive number v1(0v11) is the reduction amount of contact rate. In addition, we consider the quarantine control as follows. If the number of susceptible individuals is less than ST, we do not quarantine susceptible individuals. If S>ST holds, the quarantine control is open, and this time we assume that the positive number v2(0v21) is the quarantine rate of susceptible individuals. Figure 1 shows the broken line control strategy. Based on system (2.1) and the broken line control strategy, in this paper, we propose a Filippov influenza system with media control and quarantine of susceptible populations as follows:

    {St=γS(1SK)β(1v1)SIv2S,It=β(1v1)SI(d+δ+r)I,Rt=rIdI (2.2)
    Figure 1.  Schematic diagram of the broken line control strategy.

    with

    (v1,v2)={(0,0), for I<IT,(p,0), for I>IT and S<ST,(p,q), for I>IT and S>ST. (2.3)

    Notice that the third equation R does not contain the variables S and I of (2.2), so we can not consider the third equation R of the system. Then, in this paper, we investigate a Filippov system

    {St=γS(1SK)β(1v1)SIv2S,It=β(1v1)SI(d+δ+r)I. (2.4)

    Then, the first quadrant is divided by the following five regions:

    Γ1={(S,I)R2+:I<IT},Γ2={(S,I)R2+:I>IT and S<ST},Γ3={(S,I)R2+:I>IT and S>ST},Π1={(S,I)R2+:I=IT},Π2={(S,I)R2+:I>IT and S=ST}.

    The non-smooth system in region Γi for i=1,2,3 is described by

    (StIt)=(γS(1SK)βSIβSI(d+δ+r)I)=F1(S,I),(S,I)Γ1;(StIt)=(γS(1SK)β(1p)SIβ(1p)SI(d+δ+r)I)=F2(S,I),(S,I)Γ2;(StIt)=(γS(1SK)β(1p)SIqSβ(1p)SI(d+δ+r)I)=F3(S,I), (S,I)Γ3. (2.5)

    The normal vector of the system is defined as n1=(0,1)T in Π1, while the normal vector of the system is n2=(1,0)T in Π2. Denote the right-hand side of system (2.4) by f. The following definitions are necessary in this paper[14,46,47].

    Definition 1. [14] A point E is called a real equilibrium of system (2.4) if there exists i{1,2,3} such that Fi(E)=0 and EΓi, denoted by ER.

    Definition 2. [14] A point E is called a virtual equilibrium of system (2.4) if there exists i{1,2,3} such that Fi(E)=0 and E¯Γi, where ¯Γi is the closure of Γi, denoted by EV.

    Denote the equations that describe the sliding mode dynamics on the sliding mode domain iΠi by Fsi(S,I),i{1,2}.

    Definition 3. [14] A point E is called a pseudo-equilibrium of system (2.4) if point E is an equilibrium of Fsi(S,I) on sliding mode domain . That is, Fsi(E)=0, and point EΠi,i{1,2}.

    Lemma 1. The solution (S(t),I(t)) of system (2.4) with the positive initial value S(0) and I(0) satisfy S(t)>0 and I(t)>0 for t[0,+).

    Proof. First, we prove that S(t)>0. If not, we assume that there is a time t1 satisfying the solution S(t1)0. Then we know that there exists the other time t>0 such that the solution S(t)=0 and S(t)>0 for t[0,t). From the first equation of system (2.4), one has

    dSdt=γS(1SK)β(1v1)SIv2S=S[r(1SK)β(1v1)Iv2].

    Further, we obtain

    S(t)=S(0)et0[r(1SK)β(1v1)Iv2]ds>0,

    which is a contradiction with S(t)=0. Thus, S(t)>0, for all t>0.

    Next, we prove that I(t)>0. If not, we assume there is a time t2 satisfying the solution I(t2)0. Then, there exists the other time ˜t>0 such that the solution I(˜t)=0, and I(t)>0 for time t[0,˜t). From the second equation of system (2.4), one gets

    dIdt=I[(β(1v1)S(d+δ+r)](d+δ+r)I,t[0,˜t].

    Then, if t=˜t, we have

    I(˜t)I0e(d+δ+r)˜t>0

    which is a contradiction with I(˜t)=0. Thus, I(t)>0, for all t>0. In a word, we have that (S(t),I(t)) of (2.4) with the positive initial value S(0) and I(0) satisfies S(t)>0 and I(t)>0 for t[0,+). The proof is finished.

    Lemma 2. The solution (S,I) of system (2.4) is bounded.

    Proof. Notice the first equation of system (2.4). Then,

    dSdtx=K=β(1v1)KIv2S<0,anddSdtx>K<0.

    Hence, there is a positive number T such that the solution S(t)<K for tT.

    Let N=S+I, and then for tT, we obtain

    dNdt=γS(1SK)v2S(d+δ+r)IγK(1SK)v2S(d+δ+r)IγKmin{γ+v2,d+δ+r}N.

    Moreover,

    lim suptNγKmin{γ+v2,d+δ+r}.

    Therefore, under Lemma 1, we know that the solution of (2.4) is bounded. The proof is finished.

    Next, we give the basic reproduction number R0i,i=1,2,3, of the system. For example, in region Γ1, the basic reproduction number is defined as R01=Kβd+δ+r. In region Γ2, the basic reproduction number is defined as R02=Kβ(1p)d+δ+r. In region Γ3, the basic reproduction number is defined as R03=(γp)Kβ(1p)γ(d+δ+r).

    For system (2.4), clearly, in region Γ1, (2.4) has three equilibria, i.e., E10=(0,0), E11=(K,0) and E1=(S1,I1)=(d+δ+rβ,γβ(1d+δ+rKβ)), which is a stable node if R01>1.

    In region Γ2, system (2.4) has three equilibria, i.e., E20=(0,0), E21=(K,0) and E2=(S2,I2)=(d+δ+rβ(1p),γβ(1p)(1d+δ+rKβ(1p))), which is a stable node (focus) if R02>1.

    In region Γ3, system (2.4) has three equilibria, i.e., E30=(0,0), E31=(Kγqγ,0) and E3=(S3,I3)=(d+δ+rβ(1p),γβ(1p)(1d+δ+rKβ(1p)qr)), which is a stable node (focus) if R03>1.

    Theorem 1. Suppose that R0i<1, and the disease free equilibrium Ei1(i=1,2,3) of (2.4) is globally asymptotically stable. In addition, the endemic equilibrium Ei(i=1,2,3) of (2.4) is globally asymptotically stable if R0i>1.

    Proof. Since R0i<1,i=1,2,3, a Lyapunov function is considered as

    L(t)=SSi1Si1lnSSi1+I.

    Applying LaSalle's invariance principle[13,14], we know that Ei1(i=1,2,3) of system (2.4) is globally asymptotically stable.

    When R0i>1,i=1,2,3, the following Lyapunov function is considered:

    L(t)=SSiSilnSSi+IIiIilnIIi.

    Using LaSalle's invariance principle [13,14], we know that points Ei(i=1,2,3) of system (2.4) areglobally asymptotically stable. The proof is finished.

    This paper only considers the global dynamics of (2.4) under case R0i>1 (i=1,2,3). Next, we aim to address the richness of the possible equilibria and sliding modes on Π1 and Π2 that the system with can exhibit.

    From (2.4), when S1<S2=S3 and I3<I2, we consider the following three cases: ST<S1,S1<ST<S2=S3, and S2=S3<ST with varied IT. Further, according to the dynamics in each case, the biological phenomena of (2.4) are described in this section.

    Throughout the paper, the S-nullclines and I-nullclines of (2.4) are represented by the dashed curves and dash-dot lines, respectively. Thus, S=S1,S2 and S3 are the I-nullclines of F1,F2 and F3, denoted by L12,L22 and L32, respectively. That is, the curves

    L12:={(S,I)Γ1:γS(1SK)βSI=0},
    L22:={(S,I)Γ2:γS(1SK)β(1p)SI=0}

    and

    L32:={(S,I)Γ3:γS(1SK)β(1p)SIqS=0}

    are the S-nullclines of systems F1,F2 and F3, denoted by L11,L21 and L31, respectively.

    In this part, we first consider sliding mode dynamics of (2.4) on Π1 under Case 1: ST<S1<S2=S3. Second, the sliding mode dynamics on Π2 are also given under Case 1: ST<S1<S2=S3. In addition, we investigate the bifurcations of (2.4) under Case 1: ST<S1<S2=S3. Finally, some numerical simulations are displayed to confirm the results.

    This part investigates the existence of the sliding mode region on Π1. Based on Definition 3, if n1,F1>0 and n1,F3<0 hold, then we know that there is the sliding mode region 1, which is expressed as

    1={(S,I)Π1:S1<S<S3}.

    Using the Filippov convex method [13,14], we obtain

    (StIt)=λ1F1(S,I)+(1λ1)F3(S,I),where λ1=n1,F3n1,F3F1.

    So, we have the differential equations describing the sliding mode dynamics along the manifold 1 for system (2.4):

    (StIt)=(γS(1SK)qpS(d+δ+r)IT+(d+δ+r)qβp0). (3.1)

    Next, we analyze the existence of the positive equilibriums on 1 of (3.1). Let

    Δ1=(γqp)24γK[(d+δ+r)IT(d+δ+r)qβp],IT=qβp+(γqp)2K4γ(d+δ+r).

    Proposition 1. For varied IT, we have the following results.

    If IT>IT, then system (3.1) has no equilibrium.

    If IT>IT>qβp, system (3.1) has two positive equilibria E±s1=(S±s1,IT), where S±s1=(γqp)±Δ12γK.

    If IT<qβp, then system (3.1) has a unique positive equilibrium Es2=(Ss2,IT), where Ss2=(γqp)+Δ12γK.

    In addition, when the sliding mode 1 has a pseudo-equilibrium, we have

    ST=γqp2γK. (3.2)

    Proposition 2. Under the condition ST<S1, Es11, and the following results are given.

    If IT<I3, we have E+s11.

    If I3<IT<I1, we have E+s11.

    If IT>I1, we have E+s11.

    Proposition 3. Under the condition S1<ST<S3, the following results hold.

    (1) Assume that I1<I3, and then

    if I1<IT<I3, we have Es11, E+s11;

    if I3<IT<IT, we have Es11, E+s11.

    (2) Assume that I1>I3, and then

    if I3<IT<I1, we have Es11, E+s11;

    if I1<IT<IT, we have Es11, E+s11.

    Proposition 4. Under the condition ST>S3, E+s11, and the following conclusions hold.

    If IT<I1, we have Es11.

    If I1<IT<I3, we have Es11.

    If IT>I3, we have Es11.

    Theorem 2. If IT>IT>qβp, then a stable pseudo-equilibriums E+s1 is located on the sliding mode 1, and the unstable pseudo-equilibriums Es1 is located on the sliding mode 1.

    Proof. Notice that

    S(rS(1SK)qpS(d+δ+r)IT+(d+δ+r)qβp)|S±s1=Δ1.

    Therefore, the point E+s1 is attracting, and the point Es1 is repelling.

    Proposition 5. Assume that γqpγK>d+δ+rβ(1p), and then the pseudo-equilibrium Es2 is not located on 1.

    Proof. From Proposition 1, Eq (3.2) and the condition of Proposition 5, one has

    Ss2=(γqp)+Δ12γK>2ST=2γqpγK>d+δ+rβ(1p)=S3,

    which implies that Ss22ST>S3. Based on the definition of sliding mode region 1, we know that Ss21. Then, 1 does not have a pseudo-equilibrium Es2. The proof is completed.

    Let

    J1=γβ(1p)(1STKqγ),J2=γβ(1p)(1STK).

    Clearly J1<J2. A subset of Π2 is a sliding mode domain if F2,n2F3,n2<0. If IT>J2, there is no sliding mode domain on Π2. If IT<J2, then the sliding mode domain of (2.4) on Π2 is given as

    2={(S,I)Π2:max{IT,J1}<I<J2}. (3.3)

    Using the Filippov convex method[13], we have the differential equations describing the sliding mode dynamics along the manifold 2 for system (2.4) with (2.3):

    (StIt)=(0β(1p)STI(d+δ+r)I). (3.4)

    Clearly, there is not a positive equilibrium. Thus, if there exists a sliding domain 2 on Π2, we know the system does not have a pseudo-equilibrium.

    Due to ST<S1<S2=S3, we know that E2Γ2 is a virtual equilibrium, denoted by EV2. However, point E1 and point E3 may exist depending on IT, and we have the following three Propositions:

    Proposition 6. Assume that ST<S1, and the following assertions hold.

    If IT<I3, we have Es11,E+s11,E1Γ1,E3Γ3.

    If I3<IT<I1, we have Es11,E+s11,E1Γ1,E3Γ3.

    If I1<IT<IT, we have Es11,E+s11,E1Γ1,E3Γ3.

    If IT>IT, we have Es1andE+s1don't exist,E1Γ1,E3Γ3.

    Proposition 7. Under the condition S1<ST<S3, the following assertions hold.

    (1) Assume that I1<I3, and further,

    if IT<I1, we have Es11,E+s11,E1Γ1,E3Γ3;

    if I1<IT<I3, we have Es11,E+s11,E1Γ1,E3Γ3;

    if I3<IT<IT, we have Es11, E+s11,E1Γ1,E3Γ3;

    if IT>IT, we have Es1andE+s1don't exist,E1Γ1,E3Γ3.

    (2) Assume that I1>I3, and further,

    if IT<I3, we have Es11,E+s11,E1Γ1,E3Γ3;

    if I3<IT<I1, we have Es11,E+s11,E1Γ1,E3Γ3;

    if I1<IT<IT, we have Es11, E+s11,E1Γ1,E3Γ3;

    if IT>IT, we have Es1andE+s1don't exist,E1Γ1,E3Γ3.

    Proposition 8. Under the condition ST>S3, the following assertions hold.

    If IT<I1, we have Es11,E+s11,E1Γ1,E3Γ3.

    If I1<IT<I3, we have Es11,E+s11,E1Γ1,E3Γ3.

    If I3<IT<IT, we have Es11,E+s11,E1Γ1,E3Γ3.

    If IT>IT, we have Es1andE+s1don't exist,E1Γ1,E3Γ3.

    Based on Propositions 6–8, we have the following summary:

    B1. Let Es11,E+s11,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set B11. Then, we conclude that there does not exist a pseudo-equilibrium, and all trajectories of the system (2.4) will converge to ER3, as shown in Figure 2(a), where

    B11={ΥR2+:ST<S1,IT<min{I1,I3}}.
    Figure 2.  The trajectory of (2.4) under Case 1: ST<S1<S2=S3, d=0.01,δ=0.01,r=0.01,γ=0.2,K=3, (a) where β=0.7,p=0.2,q=0.2,ST=1.5,IT=1; (b) where β=0.7,p=0.2,q=0.2,ST=1.5,IT=2.8; (c) where β=0.7,p=0.2,q=0.2,ST=1.5,IT=3.5; (d) where β=1.5,p=0.5,q=0.2,ST=1,IT=2.8; (e) where β=1.2,p=0.5,q=0.2,ST=1,IT=3.

    B2. Let Es11,E+s11,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set B21. Then, we know that E+s11Π1 is a stable pseudo-equilibrium, and all solutions of the system will approach E+s1, as shown in Figure 2(b), where

    B21={ΥR2+:ST<S1,I3<IT<I1}.

    B3. Let Es11,E+s11,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set B31B32B33. Then, we conclude that the system (2.4) does not have a pseudo-equilibrium, and all trajectories of the system (2.4) will converge to the equilibrium point ER1, as shown in Figure 2(c), where

    B31={ΥR2+:ST<S1,I1<IT<IT, if ST<S1},B32={ΥR2+:ST<S1,IT>IT, if S1<ST<S3},B33={ΥR2+:ST<S1,I3<IT<IT, if ST>S3}.

    B4. Let Es11,E+s11,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set B41. Then, we know that Es11Π1 is an unstable pseudo-equilibrium, and all solutions will approach ER1 or ER3. The result of this numerical simulation is shown in Figure 2(d), where

    B41={ΥR2+:ST<S1,I1<IT<I3}.

    B5. Let Es11,E+s11,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set B51. Then, we conclude that Es11Π1 is an unstable pseudo-equilibrium, and all solutions of the system (2.4) will tend to ER1 or E+s1. The result of this numerical simulation is shown in Figure 2(e), where

    B51={ΥR2+:ST<S1,max{I1,I3}<IT<IT, if S1<ST<S3}.

    For the B1 of case 1, when system (2.4) has a unique equilibrium ER3, we have the following.

    Theorem 3. If ST<S1<S3=S2 and IT<min{I1,I3}, then the point ER3 of system (2.4) is globally asymptotically stable.

    Proof. Suppose that the system (2.4) has a closed orbit U (shown in Figure 3(a)) that surrounds the real equilibrium ER32 and the sliding mode 1. Define U=U1+U2+U3, where Ui=UΓi,i=1,2,3. Let Ω be the bounded region delimited by U and Ωi=ΩΓi for i=1,2,3. Considering the Dulac function D=1SI. Three Steps are given as follows:

    Figure 3.  The schematic diagram of limit cycle.

    Step 1: System (2.4) does not have a closed orbit in region Γi,i=1,2,3.

    U((Df1)S+(Df2)I)dSdI=3i=1Ui((DFi1)S+(DFi2)I)dSdI=3γK3i=1Ui1IdSdI<0,

    where f1 is the first component of f, and f2 is the second component of f.Fi1 is the first component of Fi and Fi2 is the second component of Fi,i=1,2,3. Let ˜Ωi be the region bounded by ˜Ui,Pi and Qi, where ˜Ωi and ˜Ui depend on ϵ and converge to Ωi and Ui as ϵ approaches 0.

    Step 2: System (2.4) does not have a closed trajectory in region Γ.

    We can get

    Ωi((DFi1)S+(DFi2)I)dSdI=limϵ0˜Ωi((DFi1)S+(DFi2)I)dSdI.

    Since dS=F11dt and dI=F12dt along ˜U1 and dI=0 along P1, by using Green's theorem, for region ˜Ω1, we have

    ˜Ω1((DF11)S+(DF12)I)dSdI=˜Ω1DF11dIDF12dS=˜U1DF11dIDF12dS+P1DF11dIDF12dS=P1DF12dS. (3.5)

    Similarly, we obtain

    ˜Ω2((DF21)S+(DF22)I)dSdI=P2DF22dS+Q2DF21dI (3.6)

    and

    ˜Ω3((DF31)S+(DF32)I)dSdI=P3DF32dS+Q3DF31dI. (3.7)

    From Eqs (3.5)–(3.7), we have

    0>3i=1Ωi((DFi1)S+(DFi2)I)dSdI=limϵ03i=1˜Ωi(DFi1S+DFi2I)dSdI=limϵ0(P1DF12dSP2DF22dS+Q2DF21dIP3DF32dS+Q3DF31dI). (3.8)

    Denote the intersection points of the closed trajectory U and the line I=IT by A1 and A2 and the intersection point of U and line S=ST if I>IT by A3. In addition, denote the intersection point of the line I=IT and the line S=ST by ET. Note that A1S<ST<A2S and A3I>IT. Then, Eq (3.6) becomes

    0>A1SA2S(βd+δ+rS)dSSTA1S(β(1p)d+δ+rS)dS+A3IIT(γ(1SK)Iβ(1p))dIA2SST(β(1p)d+δ+rS)dS+ITA3I(γ(1SK)qIβ(1p))dI=A2SA1SpdS+ITA3I(qI)dI=p(A2SA1S)+qln(A3IIT)>0, (3.9)

    which is a contradiction. Therefore, we know that there does not have a closed orbit U surrounding the sliding mode 1 and the real equilibrium ER3.

    Step 3: System (2.4) does not have a closed trajectory in regions Γi and Γj(ij). With a similar proof procedure to Step 2, it is easy to that there is no closed trajectory in Γ1 and Γ2 (see Figure 3(d)), Γ1 and Γ3 (see Figure 3(b)), Γ2 and Γ3 (see Figure 3(c)), respectively.

    Therefore, based on Steps 1–3, we know that the point ER3 of (2.4) is globally asymptotically stable if ST<S1<S3=S2 and IT<min{I1,I3}. This completes this theorem.

    With a similar proof procedure to Theorem 3, we have the following.

    Theorem 4. If ST<S1<S2=S3 and the conditions of (B3) hold, then the point ER1 of system (2.4) is globally asymptotically stable.

    In this part, we first consider sliding mode dynamics of (2.4) on Π1 under Case 2: S1<ST<S2=S3. Second, we study the sliding mode dynamics on Π2 under Case 2: S1<ST<S2=S3. In addition, the bifurcations of (2.4) are investigated under Case 2: S1<ST<S2=S3. Finally, some numerical simulations are displayed to confirm the results.

    For Case 2, S1<ST<S2=S3, two sliding domains on Π1 are given as

    3={(S,I)Π1:S1<S<ST},4={(S,I)Π1:ST<S<S3}.

    The dynamics on 3 are governed by

    (StIt)=(γS(1SK)(d+δ+r)IT0). (4.1)

    Now, we investigate the existence of a positive equilibrium on 3 of system (4.1). Let

    Δ2=γ24γ(d+δ+r)ITK,IT=γK4(d+δ+r).

    Proposition 9. For varied IT, we have the following results.

    (1) If IT>IT, then system (4.1) does not have an equilibrium;

    (2) If 0<IT<IT, system (4.1) has two positive equilibria E±s3=(S±s3,IT), where S±s3=γ±Δ22γK.

    Next, we find the conditions of the pseudo-equilibrium on the sliding mode 3. Let

    H1=γKd+δ+rST2+γd+δ+rST,ST=K2.

    Proposition 10. Under the condition ST<S1<ST, Es33, and the following assertions hold.

    (1) If IT>I1, we have E+s33;

    (2) If H1<IT<I1, we have E+s33;

    (3) If IT<H1, we have E+s33.

    Proposition 11. Under the condition S1<ST<ST, we have the following results.

    (1) If IT<min{I1,H1}, we have E±s33.

    (2) If min{I1,H1}<IT<max{I1,H1}, we have

    Es33,E+s33ifI1<H1;

    Es33,E+s33ifI1>H1.

    (3) If max{I1,H1}<IT<IT, we have E±s33.

    Proposition 12. Under the condition S1<ST<ST, E+s33, and the following assertions hold.

    (1) If IT<I1, we have Es33;

    (2) If I1<IT<H1, we have Es33;

    (3) If IT>H1, we have Es33.

    Theorem 5. If 0<IT<IT, then a stable pseudo-equilibrium E+s3 of (2.4) is located on the sliding mode 3, and the unstable pseudo-equilibrium Es3 of (2.4) is located on the sliding mode 3.

    Proof. Notice that

    S(γS(1SK)(d+δ+r)IT)|E±s3=Δ2.

    Therefore, the point E+s3 is attracting, and the point Es3 is repelling.

    The dynamics on region 4 are described by (3.1). Let

    H2=qβpγK(d+δ+r)ST2+(γd+δ+rqp(d+δ+r))ST.

    From Proposition 1, we have the following.

    Proposition 13. Under the condition ST<ST<S3, Es14, and the following assertions hold.

    (1) If IT<I3, we have E+s14;

    (2) If I3<IT<H2, we have E+s14;

    (3) If IT>H2, we have E+s14.

    Proposition 14. Under the condition ST<ST<S3, we have the following results.

    (1) If IT<min{I3,H2}, we have E±s14.

    (2) If min{I3,H2}<IT<max{I3,H2}, we have

    Es14,E+s14ifI3>H2;

    Es14,E+s14ifI3<H2.

    (3) If max{I3,H2}<IT<IT, we have E±s14.

    Proposition 15. Under the condition ST<S3<ST, E+s14, and the following assertions hold.

    (1) If IT<H2, we have Es14;

    (2) If H2<IT<I3, we have Es14;

    (3) If IT>I3, we have Es14.

    Theorem 6. If 0<IT<IT, then a stable pseudo-equilibrium E+s1 is located on 4, and the unstable pseudo-equilibrium Es1 is located on 4.

    If S1<ST<S2=S3, the sliding mode dynamics on region Π2 are the same as Section 3.2. We omit it here.

    For this case, we conclude that the point E2 is a virtual equilibrium. E1 and E3 are changeable depending on IT, and we have the following five situation.

    Proposition 16. Under the conditions ST<S1<ST and ST<ST<S3, the following assertions hold.

    (1) If IT<I3, we have E+s33,E+s14,E1Γ1,E3Γ3;

    (2) If I3<IT<H2, we have E+s33,E+s14,E1Γ1,E3Γ3;

    (3) If H2<IT<H1, we have E+s33,E+s14,E1Γ1,E3Γ3;

    (4) If H1<IT<I1, we have E+s33,E+s14,E1Γ1,E3Γ3;

    (5) If I1<IT<IT, we have E+s33,E+s14,E1Γ1,E3Γ3.

    Proposition 17. Under the conditions S1<ST<ST and ST<ST<S3, the following assertions hold.

    (1) Assume that I1>H1, and further,

    if IT<I3, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I3<IT<H2, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H2<IT<H1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H1<IT<I1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I1<IT<IT, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3.

    (2) Assume that H2<I1<H1, and further,

    if IT<I3, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I3<IT<H2, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H2<IT<I1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I1<IT<H1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H1<IT<IT, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3.

    (3) Assume that I3<I1<H2, and

    if IT<I3, then Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I3<IT<I1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I1<IT<H2, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H2<IT<H1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H1<IT<IT, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3.

    (4) Assume that I1<I3, and further,

    if IT<I1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I1<IT<I3, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if I3<IT<H2, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H2<IT<H1, we have Es33,E+s33,E+s14,E1Γ1,E3Γ3;

    if H1<IT<IT, we have Es33E+s33,E+s14,E1Γ1,E3Γ3.

    Proposition 18. Under the conditions S1<ST<ST and ST<ST<S3, we have the following results.

    (1) Assume that I1>H2, and further,

    if IT<I3, we have Es33,E+s14,E1Γ1,E3Γ3;

    if I3<IT<H2, we have Es33,E+s14,E1Γ1,E3Γ3;

    if H2<IT<I1, we have Es33,E+s14,E1Γ1,E3Γ3;

    if I1<IT<H1, we have Es33,E+s14,E1Γ1,E3Γ3;

    if H1<IT<IT, we have Es33,E+s14,E1Γ1,E3Γ3.

    (2) Assume that I3<I1<H2, and further,

    if IT<I3, we have Es33,E+s14,E1Γ1,E3Γ3;

    if I3<IT<I1, we have Es33,E+s14,E1Γ1,E3Γ3;

    if I1<IT<H2, we have Es33,E+s14,E1Γ1,E3Γ3;

    if H2<IT<H1, we have Es33,E+s14,E1Γ1,E3Γ3;

    if H1<IT<IT, we have Es33,E+s14,E1Γ1,E3Γ3.

    (3) Assume that I1<I3, and further,

    if IT<I1, we have Es33,E+s14,E1Γ1,E3Γ3;

    if I1<IT<I3, we have Es33,E+s14,E1Γ1,E3Γ3;

    if I3<IT<H2, we have Es33,E+s14,E1Γ1,E3Γ3;

    if H2<IT<H1, we have Es33,E+s14,E1Γ1,E3Γ3;

    if H1<IT<IT, we have Es33,E+s14,E1Γ1,E3Γ3.

    Proposition 19. Under the conditions S1<ST<ST and ST<ST<S3, we have the following results.

    (1) Assume that I3>H1, I1<H2<H1<I3, and

    if IT<I1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if I1<IT<H2, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if H2<IT<H1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if H1<IT<I3, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if I3<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3.

    (2) Assume that H2<I3<H1, I1<H2<I3<H1, and further,

    if IT<I1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if I1<IT<H2, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if H2<IT<I3, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if IT>I3, then

      ◇ if IT>H1, and further,

        ◇ if I3<IT<H1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if H1<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3.

      ◇ if IT<H1, and

        ◇if I3<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if IT<IT<H1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if H1<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3.

    (3) Assume that I1<I3<H2, I1<I3<H2<H1, and

    if IT<I1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if I1<IT<I3, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if I3<IT<H2, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if IT>H2,

      ◇ if IT>H1, and further,

        ◇ if H2<IT<H1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if H1<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3.

      ◇ if IT<H1, and further,

        ◇ if H2<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if IT<IT<H1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if H1<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3.

    (4) Assume that I3<I1, further, I3<I1<H2<H1, and

    if IT<I3, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if I3<IT<I1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if I1<IT<H2, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

    if IT>H2, and further,

      ◇ if IT>H1, and

        ◇if H2<IT<H1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇if H1<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3.

      ◇ if IT<H1, and further

        ◇ if H2<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if IT<IT<H1, we have Es33,Es14,E+s14,E1Γ1,E3Γ3;

        ◇ if H1<IT<IT, we have Es33,Es14,E+s14,E1Γ1,E3Γ3.

    Proposition 20. Suppose S1<ST<ST and ST<S3<ST, E+s33 and E+s14Π1, we have the following results.

    (1) Assume that I3>H1, and further,

    if IT<I1, we have Es33,Es14,E1Γ1,E3Γ3;

    if I1<IT<H2, we have Es33,Es14,E1Γ1,E3Γ3;

    if H2<IT<H1, we have Es33,Es14,E1Γ1,E3Γ3;

    if H1<IT<I3, we have Es33,Es14,E1Γ1,E3Γ3;

    if I3<IT<IT, we have Es33Es14,E1Γ1,E3Γ3.

    (2) Assume that I3<H1, and further,

    if IT<I1, we have Es33,Es14,E1Γ1,E3Γ3;

    if I1<IT<H2, we have Es33,Es14,E1Γ1,E3Γ3;

    if H2<IT<I3, we have Es33,Es14,E1Γ1,E3Γ3;

    if I3<IT<H1, we have Es33,Es14,E1Γ1,E3Γ3;

    if H1<IT<IT, we have Es33,Es14,E1Γ1,E3Γ3.

    Based on Propositions 16–20, the following summary is given.

    C1. Let Es33,E+s33,Es14,E+s14,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set C11. Then, we conclude that the system (2.4) does not have a pseudo-equilibrium, and all trajectories of the system (2.4) will converge to the equilibrium point ER3. The result of this numerical simulation is shown in Figure 4(a), where

    C11={ΥR2+:S1<ST<S3,IT<min{I1,I3}}.
    Figure 4.  The trajectory of system (2.4) under Case 2: S1<ST<S2=S3, d=0.1,δ=0.1,r=0.1,γ=0.2,K=1, (a) where β=0.7,p=0.2,q=0.2,ST=3,IT=1; (b) where β=0.7,p=0.2,q=0.2,ST=3,IT=2.8; (c) where β=0.7,p=0.2,q=0.2,ST=3,IT=3.1; (d) where β=0.7,p=0.2,q=0.2,ST=3,IT=3.21; (e) where β=0.7,p=0.2,q=0.2,ST=3,IT=3.4; (f) where β=1.8,p=0.6,q=0.4,ST=2.3,IT=2.5.

    C2. Let Es33,E+s33,Es14,E+s14,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set C21. Then, we know that E+s14Π1 is a stable pseudo-equilibrium, and all solutions of the system (2.4) will approach the point E+s1, as shown in Figure 4(b), where

    C21={ΥR2+:S1<ST<S3,I3<IT<min{H2,I1}}.

    C3. Let Es33,E+s33,Es14,E+s14,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set C31. Then, we conclude that the system (2.4) does not have a pseudo-equilibrium, and all trajectories of the system (2.4) will converge to the point ET=(ST,IT). The result of this numerical simulation is shown in Figure 4(c), where

    C31={ΥR2+:S1<ST<S3,H2<IT<min{H1,I1}}.

    C4. Let Es33,E+s33,Es14,E+s14Π1,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set C41. Then, we know that the point E+s33 is a stable pseudo-equilibrium, and all solutions of the system (2.4) will approach the point E+s3, as shown in Figure 4(d), where

    C41={ΥR2+:S1<ST<S3,H1<IT<I1}.

    C5. Let Es33,E+s33,Es14,E+s14,E1Γ1,E3Γ3, and the value Υ=(ST,IT) belongs to the set C51C52C53C54. Then, we know that the system (2.4) does not have a pseudo-equilibrium, and all trajectories of the system (2.4) will converge to the point ER1. The result of this numerical simulation is shown in Figure 4(e), where

    C51={ΥR2+:S1<ST<S3,I1<IT<IT, if ST<S1<ST and ST<ST<S3},C52={ΥR2+:S1<ST<S3,H1<IT<IT, if S1<ST<ST and ST<ST<S3},C53={ΥR2+:S1<ST<S3,H1<IT<IT, if S1<ST<ST and ST<ST<S3,I3<H1},C54={ΥR2+:S1<ST<S3,max{H1,I3}<IT<IT, if S1<ST<ST and ST<S3<ST}.

    C6. Let , and the value belongs to the set . Then, we conclude that is an unstable pseudo-equilibrium, and the solution of the system (2.4) will approach the point or or . The result of this numerical simulation is shown in Figure 4(f), where

    C7. Let , and the value belongs to the set . Then, we know that is an unstable pseudo-equilibrium, where is stable. The solution of the system (2.4) will tend to or or . The result of this numerical simulation is shown in Figure 5(a), where

    Figure 5.  The trajectory of system (2.4) under Case 2: , , (a) where ; (b) where ; (c) where ; (d) where ; (e) where ; (f) where .

    C8. Let , and the value belongs to the set . Then, we conclude that is an unstable pseudo-equilibrium. The solution of the system (2.4) will approach the equilibrium point or , as shown in Figure 5(b), where

    C9. Let , and the value belongs to . Then, we conclude that is a unstable pseudo-equilibrium, where is stable. The solution of the system (2.4) will tend to the equilibrium point or or . The results of this numerical simulation are shown in Figure 5(c), where

    C10. Let , and the value belongs to the set . Then, we conclude that is an unstable pseudo-equilibrium. The solution of the system (2.4) will converge to or or , as shown in Figure 5(d), where

    C11. Let , and the value belongs to the set . Then, we show that is an unstable pseudo-equilibrium. The solution of the system (2.4) will tend to or or , as shown in Figure 5(e), where

    C12. Let , , and the value belongs to the set . Then, we know that the point and the value are unstable pseudo-equilibriums, where is stable. The solution of the system (2.4) will approach or or . The result of this numerical simulation is shown in Figure 5(f), where

    C13. Let , and the value belongs to the set . Then, we conclude that is an unstable pseudo-equilibrium, where is stable. The solution of the system (2.4) will converge to or or , as shown in Figure 6, where

    Figure 6.  Dynamics of system (2.4) under Case 2: . Here, the parameter values are and .

    In this part, we first consider sliding mode dynamics of (2.4) on under Case 3: . Second, the sliding mode dynamics on are given under Case 3: . In addition, we investigate the bifurcations of (2.4) under Case 3: . Finally, some numerical simulations are displayed to confirm the results.

    If and on , then is described as

    Next, the conditions of a pseudo-equilibrium on the sliding mode are given as follows.

    Proposition 21. Under the condition , and the following assertions hold.

    (1) If , we have ;

    (2) If , we have ;

    (3) If , we have .

    Proposition 22. Under the condition , the following assertions hold.

    (1) Assume that , and further,

    if , we have , ;

    if , we have , .

    (2) Assume that , and further,

    if , we have , ;

    if , we have , .

    Proposition 23. Under the condition , and the following assertions hold.

    (1) If , we have ;

    (2) If , we have ;

    (3) If , we have .

    Theorem 7. If , the sliding mode has a stable pseudo-equilibrium , and the sliding mode has an unstable pseudo-equilibrium.

    When , the sliding mode dynamics on are the same as Section 3.2. We omit it here.

    For this Case 3, the point is a virtual equilibrium, denoted by . Points and are changeable depending on , and then we have the following.

    Proposition 24. Under the condition , the following assertions hold.

    (1) If , we have ;

    (2) If , we have ;

    (3) If , we have ;

    (4) If , we have that .

    Proposition 25. Under the condition , the following assertions hold.

    (1) Assume that , and further,

    if , we have ;

    if , we have ;

    if , we have , ;

    if , we have that .

    (2) Assume that , and further,

    if , we have ;

    if , we have ;

    if , we have , ;

    if , we have that .

    Proposition 26. Under the condition , the following assertions hold.

    (1) If , we have ;

    (2) If , we have ;

    (3) If , we have ;

    (4) If , we have that .

    Based on Propositions 24–26, the following summary is given.

    D1. Let , and the value belongs to the set . Then, we conclude that the system (2.4) does not have a pseudo-equilibrium, and all trajectories of the system (2.4) will converge to . The result of this numerical simulation is shown in Figure 7(a), where

    Figure 7.  The trajectory of the system (2.4) under Case 3: , , (a) where ; (b) where ; (c) where ; (d) where ; (e) where .

    D2. Let , and the value belongs to the set . Then, we show that is an unstable pseudo-equilibrium. The solution of the system (2.4) will approach or , as shown in Figure 7(b), where

    D3. Let , and the value belongs to the set . Then, we know that the system (2.4) does not have a pseudo-equilibrium, and all trajectories of the system (2.4) will converge to , as shown in Figure 7(c), where

    D4. Let , and the value belongs to the set . Then, we show that is an unstable pseudo-equilibrium, and the solution of system (2.4) will approach or , as shown in Figure 7(d), where

    D5. Let , and the value belongs to the set . Then, we conclude that is a stable pseudo-equilibrium. All solutions of the system (2.4) will tend to . The result of this numerical simulation is shown in Figure 7(e), where

    Remark 1. For smooth system (2.1), we have discussed the three equilibrium points, that is, , the disease free equilibrium point and the endemic equilibrium point of (2.4). By constructing a Lyapunov function, we obtain the global stability of system (2.1) in Theorem 1. For the non-smooth system (2.4), we investigate the non-smooth system (2.4) with two threshold control strategies. Using a Filippov analysis method, Green's formula, the comparison theorem and numerical simulation method, the rich dynamics of the system are given, such as the bistability phenomenon, the globally stable pseudo-equilibrium and the regular/virtulal equilibrium bifurcations. Through the two control strategies, we can control the disease individuals to the appropriate balance. In particular, Theorems 3, 6 and 7 in this paper cannot appear in the smooth system (2.1); please see Figure 4(b), (c). There is bistability in the system (2.4).

    Remark 2. In [25], a Filippov model describing the effects of media coverage and quarantine on the spread of human influenza was considered, and the threshold conditions for stability switches were obtained analytically. The discontinuous system (2.4) considered in our paper is a logistic source, and [25] considered a linear source. Second, the dynamics are different. Our paper employs the Green's theorem and a Dulac function. Then, we show that two real equilibria occur simultaneously in our paper. Using numerical simulation methods, the sliding dynamics and bifurcations of a human influenza system under logistic source and broken line control strategy are given. The results of this paper are new with respect to[25].

    Based on the previous discussion, it is shown that system (2.4) will exhibit multiple equilibriums and sliding modes. In order to better construct the bifurcation diagram, we choose and as bifurcation parameters, and the other parameters are fixed as shown in Figure 8. With the expressions of equilibria found in Section 2.3, the lines to divide the relevant parameter plane are given as follows:

    Figure 8.  The regular/virtulal equilibrium bifurcations where .

    The three solid lines , and divide the two-dimensional plane space into four regions in the first quadrant. Suppose that the control value satisfies and (that is, regions and region ; see Figure 8), the points and are virtual equilibria points (denoted by , and , respectively), and exists with the sliding mode domain. If the control value satisfies (that is, , as shown Figure 8), is a regular equilibrium, while the point is a virtual equilibrium (denoted by the equilibrium and the equilibrium , respectively), and point does not exist with the sliding mode domain. If (that is, ; see Figure 8), is a regular equilibrium, while point is a virtual equilibrium (denoted by and ), and point does not exist with the sliding mode domain.

    Next, we choose and as bifurcation parameters, and the other parameters are fixed. From Proposition 1 in this paper, the lines to divide the relevant parameter plane are given as follows:

    The two solid lines and divide the two-dimensional plane space into three regions in the first quadrant . Suppose that the control value satisfies (that is, region ; see Figure 9(a)), then system (3.1) has no equilibrium. If the control value satisfies (that is, region , as shown in Figure 9(a)), system (3.1) has two positive equilibria and , where . When the control value satisfies (that is, region ; see Figure 9(a)), system (3.1) has a unique positive equilibrium , where .

    Figure 9.  The pseudo-equilibrium bifurcations, (a) where ; (b) where .

    We choose and as bifurcation parameters, and the other parameters are fixed. With Proposition 9 in this paper, the line to divide the relevant parameter plane is given as

    The solid line divides the two-dimensional plane space into two regions in the first quadrant . Suppose that the control value satisfies (that is, region , as shown in Figure 9(b)), and then system (4.1) does not have an equilibrium. If the control value satisfies (that is, region , as shown in Figure 9(b)), system (4.1) has two positive equilibria and , where .

    Remark 3. Notice that[30] considered the global dynamics of a Filippov predator-prey model with two thresholds for integrated pest management. By using Filippov theory, the sliding mode dynamics and global dynamics were established. Different from [30], our paper shows the dynamic behavior of the Filippov model with respect to all possible equilibria. It is shown that the Filippov system tends to the pseudo-equilibrium on sliding mode domain or one endemic equilibrium or two endemic equilibria under some conditions. Second, although both this paper and [30] discuss the global dynamics of a Filippov model with two thresholds, this paper first gives different control strategies. In particular, the two real equilibria occur simultaneously using methods such as Green's theorem and a Dulac function.

    In this paper, we have established a non-smooth system to determine whether it is necessary to adopt the control strategy of media coverage and quarantine of susceptible individuals according to the number of infected and susceptible individuals. Media coverage changes the transmission mode of influenza. Further, in order to reduce the spread of influenza, when the number of cases exceeds the larger infection threshold , and the number of susceptible individuals is greater than , we will quarantine the susceptible individuals. It is worth noting that there are two difficulties in this paper. First, the traditional continuity theory cannot be applied due to the non-smooth system with the broken line control strategy. For example, when proving the global stability of discontinuous systems, the traditional Lyapunov function cannot be similarly constructed. Second, Green's formula of continuous systems cannot be used to prove the existence of global stability of the pseudo equilibriums in discontinuous systems. In this paper, by choosing different thresholds and and using Filippov theory, we study the dynamic behavior of the Filippov model with respect to all possible equilibria. The regular/virtulal equilibrium bifurcations are given. It is shown that the Filippov system tends to the pseudo-equilibrium on sliding mode domain or one endemic equilibrium or two endemic equilibria under some conditions.

    Next, we summarize the corresponding biological results of Tables 13. These results show that the choice of values of and is very important, and it determines whether to adopt control strategies.

    Table 1.  Dynamics of Filippov model (2.4).
    Condition 1 Condition 2 Result
    Figure 2(a)
    Figure 4(a)
    Figure 7(a)

     | Show Table
    DownLoad: CSV
    Table 2.  Dynamics of Filippov model (2.4).
    Condition 1 Condition 2 Result
    Figure 2(b)
    Figure 2(c)
    Figure 4(b)
    Figure 4(c)
    Figure 4(d)
    Figure 4(e)
    Figure 7(c)
    Figure 7(e)

     | Show Table
    DownLoad: CSV
    Table 3.  Dynamics of Filippov model (2.4).
    Condition 1 Condition 2 Result
    Figure 2(d)
    Figure 2(e)
    Figure 4(f)
    Figure 5(a)
    Figure 5(b)
    Figure 5(c)
    Figure 5(e)
    Figure 5(e)
    Figure 5(f)
    Figure 6
    Figure 7(b)
    Figure 7(d)

     | Show Table
    DownLoad: CSV

    ● In Table 1, we know that the infection threshold value is chosen to be small enough, i.e., , and then the number of infected individuals will reach the equilibrium of system (2.4).

    ● From Table 2, system (2.4) has a unique globally asymptotically stable pseudo-equilibrium or if or admits a unique globally asymptotically stable equilibrium when . Our control goal can be achieved finally, and there is no need to adjust the threshold strategy.

    ● In Table 3, the solution of system (2.4) will converge to a locally asymptotically stable equilibrium if or tends to a locally asymptotically stable equilibrium when or pseudo-equilibrium if . We show that it may be necessary to adjust the threshold policy according to the initial number of susceptible individuals and infected individuals. The results obtained have certain guiding significance for choosing thresholds and designing a corresponding threshold strategy.

    Next, we consider the effect of key parameters in the subsystem on the basic regeneration number as follows.

    The three-dimensional diagram of the parameter space is shown in Figure 10(a) under the parameter values of , , . It observe when the parameter increase from to , the basic reproduction number decreases correspondingly and is less than unit 1. The trajectory of the subsystem (2.4) will converge globally to the free equilibrium (see Theorem 1), implying that the infected individuals extinct and then a stable free steady state occurs.

    Figure 10.  (a) Variation of the basic reproduction number with the effect of parameters and . (b) Variation of with the effect of parameters and . (c) Variation of with the effect of parameters and . (d) Variation of with the effect of parameters and .

    The three-dimensional diagram of the parameter space is shown in Figure 10(b) under the parameter values of , , , . It is easy to observe when fixing the , increase from to , decreases correspondingly and is less than unit 1. By using Theorem 1, the infected individuals persist and the trajectory of the subsystem (2.4) will converge globally to the endemic steady state.

    The three-dimensional diagram of the parameter space is shown in Figure 10(c) under the parameter values of , , . We observe when the parameters increase from to , also increases correspondingly and is greater than the unit 1. By using Theorem 1, the infected individuals persist and the trajectory of the subsystem (2.4) will converge globally to the endemic steady state.

    The three-dimensional diagram of the parameter space is shown in Figure 10(d) under the parameter values of , , , , . It is easy to observe when fixing the parameter with a transmission rate increase from to , increases correspondingly and is greater than the unit 1. By using Theorem 1, the infected individuals persist, and the trajectory of the subsystem (2.4) will converge globally to the endemic steady state. When fixing the parameters transmission rate increase from to , decreases correspondingly and is less than the unit 1. The trajectory of the subsystem (2.4) will converge globally to the free equilibrium (see Theorem 1 of our paper), implying that the infected individuals become extinct, and then a stable free steady state occurs.

    In addition, this model has not been validated by actual influenza data, and we only have analyzed theoretically. In the next stage, we will verify and simulate the validity of the conclusions in actual time from some websites and statistics of health departments. However, the paper has studied the non-smooth system of two threshold control strategies and has validated the correctness of the theory through numerical simulation. Due to the serious lack of current influenza data from SARS-CoV-2 infections, the verification of the work is extremely difficult. However, the theory of this paper can provide appropriate guidance for the current influenza by SARS-CoV-2 infection. In this paper, we only consider the dynamics of the system (2.4) if the basic reproduction number . However, under the saturation rate , and the conditions and , the dynamical behaviors of system (2.4) and the method of proving global stability are not yet fully clear and would be our further topic.

    We sincerely thank the anonymous referees for their very detailed and helpful comments on which improved the quality of this paper. This work is supported in part by the Yunnan Fundamental Research Projects (No: 202101BE070001-051).

    The authors have no conflict of interest to declare in carrying out this research work.



    [1] Hall E, Giaccia A (2018) Radiobiology for the Radiologist.Lippincott Williams & Wilkins.
    [2] Ottolenghi A, Ballarini F, Merzagora M (1999) Modelling radiation induced biological lesions: from initial energy depositions to chromosome aberrations. Radiat Environ Bioph 38: 1-13. https://doi.org/10.1007/s004110050132
    [3] Cornforth MN, Bedford JS (1987) A quantitative comparison of potentially lethal damage repair and the rejoining of interphase chromosome breaks in low passage normal human fibroblasts. Radiat Res 111: 385-405. https://doi.org/10.2307/3576926
    [4] Cotran RS, Kumar V, Robbins SL (1989) Pathological basis of disease. Philadelphia: WB Saunders.
    [5] Ballarini F, Ottolenghi A (2004) A model of chromosome aberration induction and chronic myeloid leukaemia incidence at low doses. Radiat Environ Bioph 43: 165-171. https://doi.org/10.1007/s00411-004-0246-7
    [6] Ballarini F, Alloni D, Facoetti A, et al. (2006) Modelling radiation-induced bystander effect and cellular communication. Radiat Prot Dosim 122: 244-251. https://doi.org/10.1093/rpd/ncl446
    [7] Particle Therapy Co-Operative Group. Available from: http://www.ptcog.ch
    [8] Scholz M, Kellerer AM, Kraft-Weyrather W, et al. (1997) Computation of cell survival in heavy ion beams for therapy. Radiat Environ Bioph 36: 59-66. https://doi.org/10.1007/s004110050055
    [9] Inaniwa T, Kanematsu N, Matsufuji N, et al. (2015) Reformulation of a clinical-dose system for carbon-ion radiotherapy treatment planning at the National Institute of Radiological Sciences, Japan. Phys Med Biol 60: 3271. https://doi.org/10.1088/0031-9155/60/8/3271
    [10] Ballarini F, Carante MP (2016) Chromosome aberrations and cell death by ionizing radiation: Evolution of a biophysical model. Radiat Phys Chem 128: 18-25. https://doi.org/10.1016/j.radphyschem.2016.06.009
    [11] Ballarini F, Battistoni G, Campanella M, et al. (2006) The FLUKA code: an overview. J Phys: Conf Series 41: 151-160. https://doi:10.1088/1742-6596/41/1/014
    [12] Carante MP, Ballarini F (2016) Calculating variations in biological effectiveness for a 62 MeV proton beam. Front Oncol 6: 76. https://doi.org/10.3389/fonc.2016.00076
    [13] Carante MP, Aimè C, Cajiao JJT, et al. (2018) BIANCA, a biophysical model of cell survival and chromosome damage by protons, C-ions and He-ions at energies and doses used in hadrontherapy. Phys Med Biol 63: 075007. https://orcid.org/0000-0002-6629-3382
    [14] Carante MP, Aricò G, Ferrari A, et al. (2019) First benchmarking of the BIANCA model for cell survival prediction in a clinical hadron therapy scenario. Phys Med Biol 64: 215008. https://doi.org/10.1088/1361-6560/ab490f
    [15] Carante MP, Aricò G, Ferrari A, et al. (2020) In vivo validation of the BIANCA biophysical model: Benchmarking against rat spinal cord RBE data. Int J Mol Sci 21: 3973. https://doi.org/10.3390/ijms21113973
    [16] Carante MP, Embriaco A, Aricò G, et al. (2021) Biological effectiveness of He-3 and He-4 ion beams for cancer hadrontherapy: a study based on the BIANCA biophysical model. Phys Med Biol 66: 195009. https://doi.org/10.1088/1361-6560/ac25d4
    [17] Kozlowska W, Carante M, Aricò G, et al. (2022) First application of the BIANCA model to carbon-ion patient cases. Phys Med Biol . In press
    [18] Demaria S, Coleman CN, Formenti SC (2016) Radiotherapy: changing the game in immunotherapy. Trends Cancer 2: 286-294. https://doi.org/10.1016/j.trecan.2016.05.002
    [19] Colangelo NW, Azzam EI (2020) The importance and clinical implications of FLASH ultra-high dose-rate studies for proton and heavy ion radiotherapy. Radiat Res 193: 1-4. https://doi.org/10.1667/RR15537.1
    [20] Durante M, Cucinotta FA (2011) Physical basis of radiation protection in space travel. Rev Mod Phys 83: 1245. https://doi.org/10.1103/RevModPhys.83.1245
    [21] Ballarini F, Battistoni G, Cerutti F, et al. (2006) GCR and SPE organ doses in deep space with different shielding: Monte Carlo simulations based on the FLUKA code coupled to anthropomorphic phantoms. Adv Space Res 37: 1791-1797. https://doi.org/10.1016/j.asr.2006.03.007
    [22] Campa A, Alloni D, Antonelli F, et al. (2009) DNA fragmentation induced in human fibroblasts by 56Fe ions: experimental data and Monte Carlo simulations. Radiat Res 171: 438-445. https://doi.org/10.1667/RR1442.1
    [23] Ottolenghi A, Ballarini F, Biaggi M (2001) Modelling chromosomal aberration induction by ionising radiation: the influence of interphase chromosome architecture. Adv Space Res 27: 369-382. https://doi.org/10.1016/S0273-1177(01)00004-7
    [24] Ballarini F, Ottolenghi A (2003) Chromosome aberrations as biomarkers of radiation exposure: modelling basic mechanisms. Adv Space Res 31: 1557-1568. https://doi.org/10.1016/S0273-1177(03)00091-7
    [25] Embriaco A, Ramos R, Carante M, et al. (2021) Healthy tissue damage following cancer ion therapy: a radiobiological database predicting lymphocyte chromosome aberrations based on the BIANCA biophysical model. Int J Mol Sci 22: 10877. https://doi.org/10.3390/ijms221910877
    [26] Stram DO, Sposto R, Preston D, et al. (1993) Stable chromosome aberrations among A-bomb survivors: An update. Radiat Res 136: 29-36. https://doi.org/10.2307/3578636
    [27] Nakano M, Kodama Y, Ohtaki K, et al. (2001) Detection of stable chromosome aberrations by FISH in A-bomb survivors: comparison with previous solid Giemsa staining data on the same 230 individuals. Int J Radiat Biol 77: 971-977. https://doi.org/10.1080/09553000110050065
    [28] Bauchinger M, Schmid E, Braselmann H (2001) Time-course of translocation and dicentric frequencies in a radiation accident case. Int J Radiat Biol 77: 553-557. https://doi.org/10.1080/09553000010022382
    [29] George K, Willingham V, Wu H, et al. (2002) Chromosome aberrations in human lymphocytes induced by 250 MeV protons: Effects of dose, dose rate and shielding. Adv Space Res 30: 891-899. https://doi.org/10.1016/S0273-1177(02)00406-4
    [30] Durante M, Yamada S, Ando K, et al. (2000) X-rays vs. carbon-ion tumor therapy: cytogenetic damage in lymphocytes. Int J Radiat Oncol Biol Phys 47: 793-798. https://doi.org/10.1016/S0360-3016(00)00455-7
  • This article has been cited by:

    1. Wenjie Li, Yajuan Guan, Jinde Cao, Fei Xu, Global dynamics and threshold control of a discontinuous fishery ecological system, 2024, 182, 09600779, 114817, 10.1016/j.chaos.2024.114817
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1970) PDF downloads(204) Cited by(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog