Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Scheduling Optimization of Compound Operations in Autonomous Vehicle Storage and Retrieval System
Next Article in Special Issue
Schwartz Symmetry Condition for Coherent Diffraction Imaging Patterns
Previous Article in Journal
Determining the Safest Anchoring Knot in a Fall Arrest System by Means of Static Tests
Previous Article in Special Issue
Distribution of Return Transition for Bohm-Vigier Stochastic Mechanics in Stock Market
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Propagation Dynamics of an Epidemic Model with Heterogeneous Control Strategies on Complex Networks

1
Department of Computer Science, Shanghai University of Engineering Science, Shanghai 201620, China
2
College of Media Engineering, Communication University of Zhejiang, Hangzhou 310018, China
3
Key Lab of Film and TV Media Technology of Zhejiang Province, Hangzhou 310018, China
4
Computational Communication Collaboratory, Nanjing University, Nanjing 210023, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2024, 16(2), 166; https://doi.org/10.3390/sym16020166
Submission received: 4 January 2024 / Revised: 26 January 2024 / Accepted: 29 January 2024 / Published: 31 January 2024
(This article belongs to the Special Issue Advances in Symmetry and Complex Systems)

Abstract

:
Complex network theory involves network structure and dynamics; dynamics on networks and interactions between networks; and dynamics developed over a network. As a typical application of complex networks, the dynamics of disease spreading and control strategies on networks have attracted widespread attention from researchers. We investigate the dynamics and optimal control for an epidemic model with demographics and heterogeneous asymmetric control strategies (immunization and quarantine) on complex networks. We derive the epidemic threshold and study the global stability of disease-free and endemic equilibria based on different methods. The results show that the disease-free equilibrium cannot undergo a Hopf bifurcation. We further study the optimal control strategy for the complex system and obtain its existence and uniqueness. Numerical simulations are conducted on scale-free networks to validate and supplement the theoretical results. The numerical results indicate that the asymmetric control strategies regarding time and degree of node for populations are superior to symmetric control strategies when considering control cost and the effectiveness of controlling infectious diseases. Meanwhile, the advantages of the optimal control strategy through comparisons with various baseline immunization and quarantine schemes are also shown.

1. Introduction

The prevalence and outbreak of infectious diseases are always important issues affecting the national economy and people’s livelihoods. Research on the transmission mechanism and dynamics of infectious diseases has attracted the attention of many scholars. Specifically, mathematical models have been a potent tool for forecasting the trajectory of infectious diseases and evaluating different prevention and control measures [1,2]. Through the application of optimal control theory and qualitative assessment of mathematical models, the effective management, reduction, and potential eradication of an infectious disease can be achieved [3,4,5,6]. Based on different mechanisms of disease transmission, Kermack and McKendrick introduced the SIS and SIR epidemic models [7,8,9]. Subsequent to their work, a variety of epidemic models have emerged to shed light on the progression of various diseases and to offer insightful control measures.
Any research progress on epidemic transmission dynamics and control measures may have a significant impact on the prevention and control of infectious diseases, and this topic has attracted widespread attention in scientific research [10,11,12]. During the spread of infectious diseases, the infected transmission mainly occurs through contact between susceptible and infected individuals. Therefore, quarantine, which prevents this contact, has become a commonly employed strategy for controlling disease. In the past two decades, some important methods [13,14,15,16] have been introduced to control outbreaks of infectious diseases on the assumption of a uniformly mixing population, which show that the quarantine rate might be a factor contributing to the observed sustained oscillations in some directly transmitted viral diseases. Since Barabási and Albert [17] proposed a scale-free network model, in which the degree distribution follows a power-law distribution, researchers have gradually focused on the study of transmission and control strategies of infectious diseases on complex networks and obtained a lot of useful and insightful results. The most influential results are the SIS and SIR epidemic models established by Pastor Satoras and Vespigani [18] through mean field theory, which indicate that the epidemic threshold of a disease will tend to 0 when the network scale is large enough. Furthermore, there are many results in theoretical research on disease control strategies on heterogeneous networks. Many researchers have introduced quarantine compartments Q into some epidemic models (e.g., SIS model [10,19,20], SIR model [10], SIRS model [21], SEIR model [22,23], etc.) on scale-free networks; these results all suggest that the epidemic threshold of a disease is closely linked to the network’s structural characteristics and quarantine rate, and the heterogeneity of networks and higher infectivity increase the risk of disease transmission, while quarantine measures help to prevent infection. Therefore, constructing a network-based transmission model can more realistically reflect the propagation laws of infectious diseases in real contact networks and explore the coupling effects between network structure, epidemic dynamics, and control strategies of diseases. However, those models either ignore the demographics or assume that infectious disease control strategies may maintain symmetry among individuals, that is, the control strategy for the nodes in the network is uniform.
Control measures for infectious diseases may maintain symmetry in time by taking the same measure with the same strength at different time points. This can include isolation, social distancing, or other control measures. Such strategies can effectively reduce the scale of infection. However, relying solely on the best uniform control measures feasible given economic, social, or other constraints may not be sufficient to achieve the disease control goals. An attractive alternative approach is the use of optimal control theory to employ time-varying control within certain bounds to strike a balance between control objectives and the associated costs. Optimal control applications in epidemic dynamics have primarily concentrated on homogeneous contact networks [24,25,26,27,28,29]. Recently, many intriguing works on the optimal control of epidemics on heterogeneous networks have emerged [21,30,31,32,33,34,35,36,37]. Li et al. presented a nonlinear SIQS epidemic model of networks and explored the issue of optimal quarantine control to minimize the cost of control measures [34]. Zhang et al. examined the optimal control of an SIQRS model that includes vaccination in a network and studied the effects of different control strategies [36]. Yang et al. investigated the stability and optimal control of SIS epidemic systems in directed networks [37]. However, these contributions were restricted to an optimal control strategy for isolation and did not compare results with other forms of heterogeneous isolation strategies.
The above works provide strong theoretical supports for the prevention of infectious diseases through isolation measures. However, appropriate vaccination for susceptible individuals is also crucial for preventing and controlling infection. Based on the above considerations, we study an SIS epidemic model with heterogeneous immunization for susceptible individuals and quarantine for infected individuals, while also considering the demographics of individuals on heterogeneous networks. The organization of this paper is as follows. In Section 2, we introduce an SIQS model with heterogeneous asymmetric control strategies. Traditional infectious disease control strategies maintain the symmetry principle that the quarantine rate of an infected node is constant across individuals. By heterogeneous asymmetric control strategies, we mean that the quarantine and immunization rates of an infected node are not constant but are related to the node degree. In other words, the control strategy on the network is not uniform, and instead, different control strategies are formulated based on degree, on complex networks, and in this section, we state our assumptions. In Section 3, we obtain the equilibria and basic reproduction number. The main results are shown in Section 4, where the uniform persistence and global stability of the disease-free and endemic equilibria of a system are analyzed by different mathematical methods. The optimal time-varying quarantine control for the model is considered in Section 5. In Section 6, some numerical simulations are conducted to illustrate and supplement the analysis results. Finally, a brief conclusion is given in Section 7.

2. Description and Formation of Epidemic Models

Guided by the symmetry of the contact mode on real contact networks when diseases spread, we consider a population situated on an undirected complex network N. Every node has three optional states: the susceptible state (S), quarantined state (Q), and infected state (I) [20]. We divided the population into n groups. Let S k ( t ) , I k ( t ) , and Q k ( t ) denote the densities of susceptible, infected, and quarantine nodes (individuals) with connectivity (degree) k at time t, respectively. The mechanism of our model mainly involves three factors. A schematic diagram is shown in Figure 1, and the meanings of the parameters are listed in Table 1. These factors are as follows:
(1) 
Birth and death: Each vacant node i in the network randomly selects a neighbor. If the neighbor is a vacant node, the state of i remains unchanged. If the neighbor is a nonvacant node, the vacant node i will be activated to generate a new susceptible node at the birth rate b. Each nonvacant node becomes a vacant node at a natural death rate d per unit time. We assume that each nonvacant node has the same birth contact ability A (where A = 1 ) due to physiological constraints.
(2) 
Immunization and quarantine ( I Q ): At each time step, susceptible individuals with degree k are immunized at the immune rate δ k . The infected nodes with degree k will be quarantined at rate β k . The quarantined individuals will recover to a susceptible node at rate η . Nodes with the same degree have identical quarantine and immunization strategies, while those with different degrees have different strategies.
(3) 
SIS epidemic frameworkInfection  S I : At the initial moment, some nodes are randomly selected as infected nodes. At each time step, the possibility that each infected node i will connect to its neighboring nodes is φ ( k i ) k i , where φ ( k ) represents the infectivity of infected nodes with degree k, and φ ( k )   = k [38,39,40], φ ( k )  = A [41], φ ( k )  =  k m [42], φ ( k )  =  a k m 1 + b k m [43]. If an infected node i interacts with a susceptible node j along a connecting edge, node j has a possibility of being infected by i at a transmission rate λ i j . For a node with degree k, the overall transmission rate is λ ( k ) .
Recovery  I S : Each infected node reverts to being susceptible at recovery rate γ .
According to the above assumptions, the model is as follows:
d S k ( t ) d t = b k [ 1 N k ( t ) ] Ψ k d S k ( t ) λ ( k ) ( 1 δ k ) S k ( t ) Θ ( t ) + γ I k ( t ) + η Q k ( t ) , d I k ( t ) d t = λ ( k ) ( 1 δ k ) S k ( t ) Θ ( t ) ( β k + γ + d ) I k ( t ) , d Q k ( t ) d t = β k I k ( t ) ( η + d ) Q k ( t ) .
with initial conditions 0 S k ( 0 ) , I k ( 0 ) , Q k ( 0 ) 1 , k = 1 , 2, …, n, where Θ ( t ) = i = 1 n p ( i | k ) φ ( i ) i I i ( t ) stands for the probability that a link which originates from a node with degree k points to an infected node. Ψ k ( t ) = i = 1 n p ( i | k ) A i N i ( t ) denotes the probability that a node of degree k has a reproductive contact with an individual neighbor. The 1 i denotes the possibility of selecting a neighbor with degree i of an empty node to activate it. On uncorrelated networks, the likelihood is independent of the connectedness of the originating node of the link. Thus, we note Θ ( t ) = 1 k i = 1 n φ ( i ) p ( i ) I i ( t ) , Ψ k ( t ) = 1 k i = 1 n p ( i ) N i ( t ) . The average density of susceptible, infected, and quarantined nodes is S ( t ) = k = 1 n S k ( t ) p ( k ) , I ( t ) = k = 1 n I k ( t ) p ( k ) , and Q ( t ) = k = 1 n Q k ( t ) p ( k ) .
Let N k ( t ) denote the density of individuals of degree k at time t, and N k = S k + I k + Q k . We develop the equation for N k ( t ) as
d N k ( t ) d t = b ( 1 N k ) Ψ k ( t ) d N k ( t ) .
According to Ref. [44], we obtain lim t N k ( t ) = 0 when b d , and no other dynamic patterns are present; when b > d , lim t + N k ( t ) = N k ( t ) * , where N k * = b k Ψ * d + b k Ψ * , Ψ * = 1 k i = 1 n i p ( i ) b Ψ * d + b i Ψ * .
Given that the original system and the limiting system exhibit identical long-term dynamic characteristics, in order to better analyze the stability of the model, we consider the limit system of (1) under the condition of b > d and N k * = S k ( t ) + I k ( t ) + Q k ( t ) based on the analysis of N k ( t ) . The corresponding limiting systems of model (1) are formulated as follows:
d S k ( t ) d t = b k [ 1 N k * ] Ψ * d S k ( t ) λ ( k ) ( 1 δ k ) S k ( t ) Θ ( t ) + γ I k ( t ) + η Q k ( t ) , d I k ( t ) d t = λ ( k ) ( 1 δ k ) S k ( t ) Θ ( t ) ( β k + γ + d ) I k ( t ) , d Q k ( t ) d t = β k I k ( t ) ( η + d ) Q k ( t ) .
with initial conditions 0 S k ( 0 ) , I k ( 0 ) , Q k ( 0 ) N k * < 1 , k = 1 , 2 , , n .
The following proposition shows that the state space of the solutions of system (3) is positively invariant.
Proposition 1. 
If ( S k ( t ) , I k ( t ) , Q k ( t ) ) k is a solution to system (3) that meets the initial conditions, then 0 S k ( t ) + I k ( t ) + Q k ( t ) N k * < 1 , and k N n for any t 0 , i.e., Γ is positively invariant.
Proof. 
x 0 Γ . Thus, t > 0 , and so we have N k 0 0 , Ψ ( 0 ) 0 .
d Ψ ( t ) d t = ( b d ) Ψ ( t ) b Ψ ( t ) k i i p ( i ) N i ( t ) .
We find that Ψ ( t ) = Ψ ( 0 ) e b d b k Σ i i p ( i ) N i ( t ) 0 . By N k ( t ) = S k ( t ) + I k ( t ) + Q k ( t ) and Equation (2), the evolution of N k ( t ) is derived by the following equation:
d N k ( t ) d t = b k [ 1 N k ( t ) ] Ψ d N k ( t ) [ b k Ψ + d ] N k ( t ) .
We have N k ( t ) N k ( 0 ) e ( b k Ψ + d ) t 0 . Then, from system (2), we also acquire
d ( 1 N k ( t ) ) d t = ( b k Ψ + d ) [ 1 N k ( t ) ] + d .
Obviously, N k ( t ) < 1 , and 0 S k ( t ) + I k ( t ) + Q k ( t ) = N k ( t ) < 1 .
Then, we verify I k 0 for any k = 1 , 2 , , n . If not, because I k ( 0 ) 0 , there exists a j { 1 , 2 , , n } and t 0 0 , satisfying
t 0 = inf { t | I j ( t ) = 0 , d I j ( t ) d t < 0 } .
It is evident that I j ( t 0 ) = 0 ,   I j ( t 0 ) < 0 for any t [ 0 , t 0 ] , I j ( t ) 0 ; therefore, we have
d I j ( t 0 ) d t = λ ( j ) ( 1 δ j ) S j ( t 0 ) Θ ( t 0 ) ( β j + γ + d ) I j ( t 0 ) = λ ( j ) ( 1 δ j ) S j ( t 0 ) Θ ( t 0 ) < 0 .
Because Θ ( t 0 ) = 1 k i = 1 n φ ( i ) p ( i ) I i ( t 0 ) 0 , we obtain S j ( t 0 ) < 0 . On the other hand, for any t [ 0 , t 0 ] , Q j ( 0 ) 0 , we attain
d Q j ( t ) d t = β j I j ( t ) ( η + d ) Q j ( t ) ( η + d ) Q j ( t ) .
That is to say, Q j ( t ) Q j ( 0 ) e ( η + d ) t 0 . When t [ 0 , t 0 ] ,
d S j ( t ) d t = b j [ 1 N j ( t ) ] Ψ j d S j ( t ) λ ( j ) ( 1 δ j ) S j ( t ) Θ ( t ) + γ I j ( t ) + η Q j ( t ) d S j ( t ) λ ( j ) ( 1 δ j ) S j ( t ) Θ ( t ) .
Then, we acquire S j ( t 0 ) S j ( 0 ) e 0 t 0 d + λ j ( 1 δ j ) Θ ( t ) d τ 0 . It is contradictory; hence, I k ( t ) 0 , for any k and t 0 . We also can obtain S k ( t ) , Q k ( t ) 0 , and S k ( t ) , I k ( t ) , Q k ( t ) < 1 . That completes the proof.  □

3. Equilibria and Basic Reproduction Number

In this section, the basic reproduction number R 0 is determined, which is the average number of new infections generated by a single newly infected individual throughout the entire infection period [19]. Considering the equivalent system,
d I k ( t ) d t = λ ( k ) ( 1 δ k ) ( N k * I k ( t ) Q k ( t ) ) Θ ( t ) ( β k + γ + d ) I k ( t ) , d Q k ( t ) d t = β k I k ( t ) ( η + d ) Q k ( t ) .
Let X denote the space of functions and x = ( I 1 ( t ) , Q 1 ( t ) , , I n ( t ) , Q n ( t ) ) X denote the state space for system (4).
Γ = { x X : 0 I k ( t ) , Q k ( t ) < 1 , k = 1 , 2 , , n } .
Proposition 1 can show that Γ is positively invariant for system (4). All viable steady states of system (4) satisfy the following equation:
λ ( k ) ( 1 δ k ) ( N k ( t ) I k ( t ) Q k ( t ) ) Θ ( t ) ( β k + γ + d ) I k ( t ) = 0 , β k I k ( t ) ( η + d ) Q k ( t ) = 0 .
In particular, there exists a disease-free equilibrium of (4) E 0 = ( 0 , 0 , , 0 ) ) , and the positive equilibrium E * = ( I k * , Q k * ) k , where E * satisfies
Q k * = β k η I k * , Θ * = 1 k i = 1 n φ ( i ) p ( i ) I i * , I k * = λ ( k ) ( 1 δ k ) N k * Θ * λ ( k ) ( 1 δ k ) ( 1 + β k η ) Θ * + β k + γ + d .
The basic reproduction number R 0 can be determined according to the following theorem:
Theorem 1. 
Define the basic reproduction number
R 0 = 1 k i = 1 n φ ( i ) p ( i ) λ ( i ) ( 1 δ i ) N i * β i + γ + d .
There always exists a disease-free equilibrium E 0 . If and only if R 0 > 1 , the system (4) has unique positive equilibrium point E * , which satisfies
Q k * = β k η I k * , Θ * = 1 k i = 1 n φ ( i ) p ( i ) I i * , I k * = λ ( k ) ( 1 δ k ) N k * Θ * λ ( k ) ( 1 δ k ) ( 1 + β k η ) Θ * + β k + γ + d .
Proof. 
Equation (6) implies that
I k = λ ( k ) ( 1 δ k ) N k Θ λ ( k ) ( 1 δ k ) ( 1 + β k η ) Θ + β k + γ + d ,
Substituting (8) into Θ , we obtain
Θ = 1 k i = 1 n φ ( i ) p ( i ) λ ( i ) ( 1 δ i ) N i Θ λ ( i ) ( 1 δ i ) ( 1 + β i η ) Θ + β i + γ + d Θ f ( Θ ) ,
where
f ( Θ ) = 1 k i = 1 n φ ( i ) p ( i ) λ ( i ) ( 1 δ i ) N i λ ( i ) ( 1 δ i ) ( 1 + β i η ) Θ + β i + γ + d .
Because f ( Θ ) < 0 and lim Θ + f ( Θ ) = 0 , it is clear that Equation (9) has a positive solution if and only if f ( 0 ) = 1 k i = 1 n φ ( i ) p ( i ) λ ( i ) ( 1 δ i ) N i β i + γ + d > 1 . That is to say, when R 0 > 1 , Equation (6) has a unique positive solution Θ * , namely, a unique positive equilibrium point E * of system (4) exists as described by Equation (7) if and only if R 0 > 1 . The proof is completed.  □
Remark 1. 
Without losing generality, we know that there always exists constant β m i n , β m a x , δ m i n , δ m a x , N m i n   ( 0 , 1 ) that satisfy
β m i n β k β m a x , δ m i n δ k δ m a x , N m i n N k * < 1 . Then, let φ ( k ) = k ,   λ ( k ) = λ k , we can obtain R 0 λ k 2 k ( 1 δ m a x ) N m i n β m a x + γ + d . Based on [18], we can find when the network size is large enough and the second-order moment of degree k 2 ; then, R 0 can always be guaranteed to be greater than 1 and the disease will persist. Therefore, the topology of the network, such as the average degree k and the second moment of degree k 2 , has an important impact on the epidemic threshold of infectious diseases.

4. Stability Analysis for SIQS Model

4.1. Stability Analysis of Disease-Free Equilibrium

Firstly, it is proved that the disease-free equilibrium (DFE) E 0 is locally asymptotically stable by the same method used in [45] by the following theorem.
Theorem 2. 
If R 0 < 1 , then the DFE of system (4) is locally asymptotically stable and unstable if R 0 > 1 .
Proof. 
The Jacobian matrix of system (4) at E 0 is
A = A 11 A 12 A 21 A 22 ,
where A 12 is the zero matrix of order n, A 21 = β k E n , A 22 = ( η + d ) E n , E n is the identity matrix, and
A 11 = ξ 1 φ 1 p 1 ( γ + β 1 ) ξ 1 φ 2 p 2 ξ 1 φ n p n ξ 2 φ 1 p 1 ξ 2 φ 2 p 2 ( γ + β 2 ) ξ 2 φ n p n ξ n φ 1 p 1 ξ n φ 2 p 2 ξ n φ n p n ( γ + β n ) ,
where ξ i = λ ( i ) ( 1 δ i ) N i * k ,   φ i p i = φ ( i ) p ( i ) . Let y = { y 1 , y 2 , , y n , y n + 1 , , y 2 n } Γ .
Subsequently, the system (4) can be reformulated in a vector format as
d y d t = A y + N ( y ) .
Let N ( y ) represent a column vector and N ( y ) 0 ,
N ( y ) = λ ( 1 ) ( 1 δ 1 ) ( I 1 + Q 1 ) Θ λ ( 2 ) ( 1 δ 2 ) ( I 2 + Q 2 ) Θ λ ( n ) ( 1 δ n ) ( I n + Q n ) Θ 0 0 0 .
The eigenvalues of matrix A are determined by | ρ E n A 11 | | ρ E n A 22 | = 0 . It is obvious that the characteristic equation has n multiple roots ( η + d ) from | ρ E n A 22 | = 0 . Other eigenvalues of matrix A are determined by | ρ E n A 11 | = 0 . Through some elementary transformations, | ρ E n A 11 | = 0 can be changed to
ρ ξ 1 φ 1 p 1 + ( β 1 + γ + d ) ξ 1 φ 2 p 2 ξ 1 φ n p n ξ 1 ξ 2 ( ρ + β 1 + γ + d ) ρ + ( β 2 + γ + d ) 0 ξ 1 ξ n ( ρ + β 1 + γ + d ) 0 ρ + ( β n + γ + d ) = 0 .
Case 1 
If ρ ( β i + γ + d ) , then we obtain
( ρ + β 1 + γ + d ) 1 1 k i = 1 n λ ( i ) ( 1 δ i ) N i * φ ( i ) p ( i ) ρ + β i + γ + d = 0 .
Note that G ( ρ ) = 1 k i = 1 n λ ( i ) ( 1 δ i ) N i * φ ( i ) p ( i ) ρ + β i + γ + d . Obviously, the solutions of Equation (12) are determined by G ( ρ ) = 1 . Since G ( ρ ) < 0 , G ( 0 ) = R 0 , lim ρ + G ( ρ ) = 0 ,   lim ρ ( β i + γ + d ) G ( ρ ) = ; hence, when R 0 < 1 , the solutions of equation G ( ρ ) = 1 are negative. Therefore, Equation (12) only has negative solutions.
Case 2 
If i = 2 n ( ρ + β i + γ + d ) = 0 , since the sum of all eigenvalues is equal to the trace of the matrix, when R 0 < 1 , ρ 1 = ξ 1 φ ( i ) p ( i ) ( β 1 + γ + d ) < 0 , and ρ i = ( β i + γ + d ) < 0 , i = 2 , , n .
Furthermore, because N ( y ) 0 , we have d y d t A y . Since all the eigenvalues of matrix A are negative, the equilibrium point E 0 of system (4) is locally asymptotically stable when R 0 < 1 .  □
Then, we get the global asymptotic stability of E 0 .
Theorem 3. 
When R 0 < 1 , the disease-free equilibrium E 0 of system (4) is global asymptotic stability.
Proof. 
From system (3), we obtain
d Θ ( t ) d t = 1 k k = 1 n φ ( k ) p ( k ) λ ( k ) ( 1 δ k ) ( N k * I k ( t ) Q k ( t ) ) Θ ( t ) ( β k + γ + d ) Θ ( t ) 1 k k = 1 n φ ( k ) p ( k ) λ ( k ) ( 1 δ k ) N k * Θ ( t ) ( β k + γ + d ) Θ ( t ) = ( β k + γ + d ) 1 k k = 1 n φ ( k ) p ( k ) λ ( k ) ( 1 δ k ) N k * β k + γ + d 1 Θ ( t ) = ( β k + γ + d ) ( R 0 1 ) Θ ( t ) .
From Proposition 1, we obtain Θ ( t ) Θ ( 0 ) e ( β k + γ + d ) ( R 0 1 ) t , while Θ ( t ) 0 . Therefore, if R 0 < 1 , we can acquire lim t + Θ ( t ) = 0 , and lim t + I k ( t ) = 0 . From the second equation of system (4), we obtain lim t + Q k ( t ) = 0 . This demonstrates that E 0 of system (4) is globally stable when R 0 < 1 . The proof is finished.  □

4.2. Global Stability of Endemic Equilibrium

In this section, we seek to study the conditions for the uniform persistence of disease, which is important in proving the global stability of endemic equilibrium E * . We obtain the uniform persistence based on Theorem 4.2 of Ref. [46]. The result of uniform persistence is shown in the following Lemma.
Lemma 1. 
When R 0 > 1 , the system (4) is permanent, i.e., there exists ε > 0 , such that
lim t + inf { I k ( t ) , Q k ( t ) } k = 1 n ε ,
where ( I k ( t ) , Q k ( t ) ) is any solution of system (4) satisfying (5).
Proof. 
Based on Theorem 4.2 of [46], we need to verify that all hypotheses for system (4) are satisfied. By setting C = Γ , Condition (1) is fulfilled. Condition (2) is evidently met. For Condition (3), notice that A 11 T is irreducible and a i j 0 whenever i j ; hence, there exists a positive eigenvector h ˜ = ( h 1 , h 2 , , h n ) of A 11 T , and its corresponding eigenvalue is s ( A 11 ) ( s ( A ) = R e i λ i ) when R 0 > 1 , s ( A 11 ) > 0 . Let h n + 1 = h n + 2 = = h 2 n = 0 ; then, we note h = ( h 1 , h 2 , , h 2 n ) , and it satisfies A T h = λ h . Let r = min 1 i n ( h i ) ; for any y Γ , one has h · y r i = 1 n y i r i = 1 n y i 2 . Thus, Condition (3) is verified. Since every part of N ( y ) is negative and h 0 , Condition (4) is satisfied. To confirm Condition (5), we set G = { y Γ | h · N ( y ) = 0 } . If y G , then
i = 1 n h i ( y i + y n + i ) λ ( i ) ( 1 δ i ) Θ = 0 .
Given that every term in the above sum is non-negative, one has y i = 0 for all i = 1 , 2 , n . Therefore, the only invariant set corresponding to (11) that is contained within G is y = 0 . Condition (5) is satisfied. All conditions are satisfied, and the proof of uniform persistence is completed.  □
From Lemma 1, it is clear that if R 0 > 1 , the infection will always exist. The following theorem shows the global stability of endemic equilibrium.
Theorem 4. 
Suppose that ( I k ( t ) , Q k ( t ) ) is a solution of system (4), satisfying the initial condition. If R 0 > 1 , then lim t + I k ( t ) = I k * , lim t + Q k ( t ) = Q k * , that is to say, the endemic equilibrium E * is globally attractive.
Proof. 
Firstly, from Lemma 1 and for any k { 1 , 2 , , n } , we find that there exists 0 < ε 1 , such that, for t that is sufficiently large, I k ε and Q k ε are satisfied. Hence, when t > 0 , we obtain ε φ ( k ) k Θ ( t ) φ ( k ) k . Then, from system (4), we obtain
d I k ( t ) d t λ ( k ) φ ( k ) k ( 1 δ k ) ( N k * I k ( t ) ) ( β k + γ + d ) I k ( t ) , t > 0 .
For any given constant 0 < ε 1 < min { ε , λ ( k ) ( 1 δ k ) φ ( k ) k ( 1 N k * ) + β k + γ + d λ ( k ) ( 1 δ k ) φ ( k ) k + β k + γ + d } , there exists a t 1 > 0 , when t > t 1 , such that
I k ( t ) A k 1 = λ ( k ) ( 1 δ k ) φ ( k ) k N k * λ ( k ) ( 1 δ k ) φ ( k ) k + β k + γ + d + ε 1 , t > t 1 .
Then, from the second equation of system (4), we have
d Q k ( t ) d t = β k ( N k * S k ( t ) Q k ( t ) ) ( η + d ) Q k ( t ) β k ( N k * Q k ( t ) ) ( η + d ) Q k ( t ) , t > 0 .
Hence, for any given constant 0 < ε 2 < min { 1 2 , ε 1 , η + d + β k β k N k * η + β k } , there exists a t 2 > t 1 , such that
Q k ( t ) B k 1 = β k N k * η + d + β k , t > t 2 .
On the other hand, from the first equation of system (4), it can be inferred that
d I k ( t ) d t λ ( k ) ε φ ( k ) k ( 1 δ k ) ( N k * I k ( t ) B k 1 ) ( β k + γ + d ) I k ( t ) .
Hence, for any given constant 0 < ε 3 < min { 1 3 , ε 2 , λ ( k ) ( 1 δ k ) ε φ ( k ) k ( N k * B k 1 ) λ ( k ) ( 1 δ k ) ε φ ( k ) k + β k + γ + d } , there exists a t 3 > t 2 , such that,
I k ( t ) a k 1 = λ ( k ) ( 1 δ k ) ε φ ( k ) k ( N k * B k 1 ) λ ( k ) ( 1 δ k ) ε φ ( k ) k + β k + γ + d ε 3 , t > t 3 .
Next, following from the second equation of system (4), we also obtain
d Q k ( t ) d t = β k I k ( t ) ( η + d ) Q k ( t ) β k a k 1 ( η + d ) Q k ( t ) .
Similarly, for any given constant 0 < ε 4 < min { 1 4 , ε 3 , β k a k 1 η + d } , there exists a t 4 > t 3 , such that
Q k ( t ) b k 1 = β k a k 1 η + d ε 4 , t > t 4 .
Since ε i is a small constant, it follows that
0 < a k 1 I k ( t ) A k 1 < 1 , 0 < b k 1 Q k ( t ) B k 1 < 1 , t > t 4 .
and we also obtain 0 < 1 k i = 1 n φ ( i ) p ( i ) a i 1 Θ ( t ) 1 k i = 1 n φ ( i ) p ( i ) A i 1 < φ ( k ) k , t > t 4 . Again, by system (4), one has
d I k ( t ) d t λ ( k ) ( 1 δ k ) ( N k * I k ( t ) b k 1 ) 1 k i = 1 n φ ( i ) p ( i ) A i 1 ( β k + γ + d ) I k ( t ) , t > t 4 .
Hence, 0 < ε 5 < min { 1 5 , ε 4 , ( 1 N k * + b k 1 ) λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) A i 1 + β k + γ + d λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) A i 1 + β k + γ + d } , there exists a t 5 > t 4 , such that,
I k ( t ) A k 2 = ( N k * b k 1 ) λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) A i 1 λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) A i 1 + β k + γ + d + ε 5 , t > t 5 .
Thus,
d Q k ( t ) d t = β k I k ( t ) ( η + d ) Q k ( t ) β k A k 2 ( η + d ) Q k ( t ) , t > t 5
and for any given constant 0 < ε 6 < min { 1 6 , ε 5 , η + d β k A k 2 η + d } , there exists a t 6 > t 5 , such that,
Q k ( t ) B k 1 = min { B k 1 , β k A k 2 η + d + ε 6 } , t > t 6 .
Consequently, one obtains that
d I k ( t ) d t λ ( k ) ( 1 δ k ) ( N k * I k ( t ) B k 2 ) 1 k i = 1 n φ ( i ) p ( i ) a i 1 ( β k + γ + d ) I k ( t ) , t > t 6 .
Hence, for any given constant 0 < ε 7 < min { 1 7 , ε 6 , ( N k * B k 2 ) λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) a i 1 λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) a i 1 + β k + γ + d } , there exists a t 7 > t 6 , such that
I k ( t ) a k 2 = ( N k * B k 2 ) λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) a i 1 λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) a i 1 + β k + γ + d ε 7 , t > t 7 .
Again, one has
d Q k ( t ) d t = β k I k ( t ) ( η + d ) Q k ( t ) β k a k 2 ( η + d ) Q k ( t ) , t > t 7 .
For any given constant 0 < ε 8 < min { 1 8 , ε 7 , β k a k 2 η + d } , there exists a t 8 > t 7 , such that
Q k ( t ) b k 2 = max { b k 1 , β k a k 2 η + d ε 8 } , t > t 8 .
It is clear that
0 < a k 1 < a k 2 I k ( t ) A k 2 < A k 1 < 1 , 0 < b k 1 < b k 2 Q k ( t ) B k 2 < B k 1 < 1 , t > t 8 .
Similarly, the calculation’s step m can be executed, resulting in four sequences. { A k m } , { B k m } , { a k m } , { b k m } are obtained as follows:
I k ( t ) A k m = ( N k * b k m 1 ) λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) A i m 1 λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) A i m 1 + ( β k + γ + d ) + ε 4 m 3 , t > t 4 m 3 , Q k ( t ) B k m = min { B k m 1 , β k A i m η + d + ε 4 m 2 } , t > t 4 m 2 , I k ( t ) a k m = ( N k * B k m ) λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) a i m 1 λ ( k ) ( 1 δ k ) 1 k i = 1 n φ ( i ) p ( i ) a i m 1 + ( β k + γ + d ) ε 4 m 1 , t > t 4 m 1 , Q k ( t ) b k m = max { b k m 1 , β k a k m η + d ε 4 m } , t > t 4 m .
It is clear that
0 < a k m I k ( t ) A k m < 1 , 0 < b k m Q k ( t ) B k m < 1 , t > t 4 m .
We find that { a k m } , { b k m } are monotone increasing sequences and { A k m } , { B k m } are strictly monotone decreasing sequences. It is clear that
B k m = β k A k m η + d + ε 4 m 2 , b k m = β k a k m η + d ε 4 m .
Since all are bounded monotonic sequences, the sequential limits exist. Let
lim t + A k ( m ) ( t ) = A k , lim t + B k ( m ) ( t ) = B k , lim t + a k ( m ) ( t ) = a k , lim t + b k ( m ) ( t ) = b k .
Since { ε m } satisfies that 0 < ε m < 1 m , one has ε m 0 as m . Therefore, from (14), as m , we obtain
A k = ( N k * b k ) λ ( k ) ( 1 δ k ) A λ ( k ) ( 1 δ k ) A + ( β k + γ + d ) , B k = β k A k η + d ,
a k = ( N k * B k ) λ ( k ) ( 1 δ k ) a λ ( k ) ( 1 δ k ) a + ( β k + γ + d ) , b k = β k a k η + d .
where A = 1 k i = 1 n φ ( i ) p ( i ) A i , a = 1 k i = 1 n φ ( i ) p ( i ) a i . Then, we obtain
A k = N k * ς k A ( η + d ) [ a ς k ( η + d β k ) + ( η + d ) ( β k + γ + d ) ] H k , a k = N k * ς k a ( η + d ) [ A ς k ( η + d β k ) + ( η + d ) ( β k + γ + d ) ] H k .
where ς k = λ ( k ) ( 1 δ k ) , H k = ( η + d ) 2 ( ς k A + ( β k + γ + d ) ) ( ς k a + ( β k + γ + d ) ) ς k 2 β k 2 A a . For any k, if η + d = β k , it is clear that A = a . If η + d β k , substituting (16) into A, we obtain
A = 1 k i = 1 n φ ( i ) p ( i ) N i * ς i A ( η + d ) [ a ς i ( η + d β i ) + ( η + d ) ( β i + γ + d ) ] H i , a = 1 k i = 1 n φ ( i ) p ( i ) N i * ς i a ( η + d ) [ A ς i ( η + d β i ) + ( η + d ) ( β i + γ + d ) ] H i .
We obtain 1 k i = 1 n φ ( i ) p ( i ) ς i N i * H i ( η + d ) 2 ( β i + γ + d ) ( A a ) = 0 . Therefore, A = a . Above all, A = a is equivalent to A i = a i for any i = 1 , 2 , , n . Then, from (13) and (14), it follows that
I k * = lim t + I k ( t ) = A k = a k = ( N k * b k ) λ ( k ) ( 1 δ k ) Θ * λ ( k ) ( 1 δ k ) Θ * + ( β k + γ + d ) ,
Q k * = lim t + Q k ( t ) = B k = b k = β k I k * η + d .
which is consistent with Theorem 2. This completes the proof.  □

5. The Optimal Control for the SIQS Model

In this section, we study the optimal control for system (1). To meet the control objective and minimize control expenses at the same time, a viable approach to targeting epidemic outbreaks is by implementing time-varying control using optimal control theory for various infectious diseases, for example, the results in Culshaw et al. [47], Chen et al. [48], and Abboubakar et al. [49]. The aim of this section is to identify the optimal quarantine measures for controlling the transmission of diseases.
First, we introduce a time-varying control rate u k ( t ) representing the percentage of infected individuals with degree k being quarantined. Hence, system (3) can be written as
d S k ( t ) d t = b k [ 1 N k * ] Ψ * d S k ( t ) λ ( k ) ( 1 δ k ) S k ( t ) Θ ( t ) + γ I k ( t ) + η Q k ( t ) , d I k ( t ) d t = λ ( k ) ( 1 δ k ) S k ( t ) Θ ( t ) ( u k ( t ) + γ + d ) I k ( t ) , d Q k ( t ) d t = u k ( t ) I k ( t ) ( η + d ) Q k ( t ) .
We note the quarantine functions u ( t ) = ( u 1 ( t ) , u 2 ( t ) , u n ( t ) ) are bounded, Lebesgue integrable functions. Since our aim is both to reduce the scale of infection and to control its strength, we consider the following objective functional J ( u k ) associated with the model in the control set Ω = { ( u 1 ( t ) , u 2 ( t ) , u n ( t ) ) 0 u i 0.9 } . B k denotes a positive weight parameter. The objective (cost) functional is given by
J ( u k ) = 0 T k = 1 n [ I k ( t ) + 1 2 B k u k 2 ( t ) ] d t .
The Lagrange function for optimal control problems (18) and (19) is defined as
L = i = 1 n [ I i + B i 2 u i 2 ( t ) ] .
In order to solve the optimal control problem, the existence of an optimal control must be assured. We applied Corollary 4.1 in Ref. [50] to prove the existence of the solution, which is shown in the proof of Theorem 5.
Theorem 5. 
For the objective functional J ( u k ) associated with model (18) defined in Ω, there exists an optimal quarantine strategy μ * = ( μ 1 , , μ n ) , minimizing J ( u ) .
Proof. 
(i) Clearly, for any 0 t T , the Ω is a nonempty set of Lebesgue integrable functions. (ii) It is obvious that the solutions are bounded, ensuring the boundedness and convexity of the admissible control set. (iii) Model (18) can be written as
d Φ ( t ) d t = F ( t , Φ , u k ( t ) )
where F ( t , Φ , u k ) = ( F 1 , F 2 , , F n ) is the right side of system (18), and Φ = ( Φ 1 , Φ 2 , , Φ n ) T , Φ k = ( S k ( t ) , I k ( t ) , Q k ( t ) ) T . It follows that
| F ( t , Φ 1 , u k ( t ) ) F ( t , Φ 2 , u k ( t ) ) | i | F i ( t , Φ 1 , u i ( t ) ) F i ( t , Φ 2 , u i ( t ) ) | .
For any k N n , since | Θ 1 ( t ) Θ 2 ( t ) | = 1 k i = 1 n φ ( i ) p ( i ) | I i 1 ( t ) I i 2 ( t ) | and the boundedness of the solution, there exists a constant M 1 such that | Θ 1 ( t ) Θ 2 ( t ) | M 1 | I k 1 ( t ) I k 2 ( t ) | . Then, there exists a positive constant M to guarantee the inequality | F k ( Φ 1 ) F k ( Φ 2 ) | M | Φ k 1 Φ k 2 | is satisfied. (iv) Letting 0 < ξ < 1 and m k ( t ) , v k ( t ) Ω , we obtain L ( t , Φ , ( 1 ξ ) m ( t ) + ξ v ( t ) ) ( 1 ξ ) L ( t , Φ , m ( t ) ) ξ L ( t , Φ , v ( t ) ) = i = 1 n B i 2 ξ ( ξ 1 ) ( v i m i ) 2 < 0 . (v) There exist c 1 > 0 , ρ > 1 , c 2 > 0 such that L ( t , Φ , u k ( t ) ) = i = 1 n [ I i + B i 2 u i 2 ( t ) ] i = 1 n B i 2 u i 2 ( t ) B m i n 2 i = 1 n u i 2 ( t ) C 1 | u | 2 C 2 .
By utilizing the results of Ref. [50], we find that there exists an optimal control μ * = ( μ 1 , , μ n ) , minimizing J ( μ ) . The proof is completed.  □
The solution to the optimal control problem is determined through Pontryagin’s minimum principle [50]. We note the Hamiltonian H as
H = k = 1 n [ I k ( t ) + 1 2 B k u k 2 ( t ) ] d t + i = 1 3 k = 1 n ω k i F k i ( t , Φ , u k ( t ) ) .
If ( Φ * , u k * ) is an optimal solution of the optimal control problem, then there exists a nontrivial vector function ω k = ( ω k 1 , ω k 2 , ω k 3 ) , k = 1 , 2 , , n , satisfying the following equalities:
d Φ ( t ) d t = H ( t , Φ * ( t ) , u k * , ω ( t ) ) ω , H ( t , Φ * ( t ) , u k * , ω ( t ) ) u k ( t ) = 0 , d ω ( t ) d t = H ( t , Φ * ( t ) , u k * , ω ( t ) ) Φ .
It follows from the derivation above that
ω * = 0 , H u k ( t ) < 0 , 0 < ω * < 0.9 , H u k ( t ) = 0 , ω * = 0.9 , H u k ( t ) > 0 .
Now the necessary conditions are implemented on the Hamiltonian H in (21) and we obtain the following results.
Theorem 6. 
For a given optimal control solution μ ( t ) and the corresponding system solution S k * ( t ) , Q k * ( t ) , I k * ( t ) , there is an adjacency function ω k 1 ( t ) , ω k 2 ( t ) , ω k 3 ( t ) that satisfies the following conditions:
d ω k 1 ( t ) d t = λ ( k ) ( 1 δ k ) Θ * ( t ) [ ω k 1 ( t ) ω k 2 ( t ) ] + d ω k 1 ( t ) 1 , d ω k 2 ( t ) d t = λ ( k ) ( 1 δ k ) S k * ( t ) φ ( k ) p ( k ) k ( ω k 1 ( t ) ω k 2 ( t ) ) γ ω k 1 ( t ) + ( u k + γ + d ) ω k 2 ( t ) u k ω k 3 ( t ) 1 , d ω k 3 ( t ) d t = η ω k 1 ( t ) + ( η + d ) ω k 3 ( t ) .
with conditions ω k 1 ( T )  =  ω k 2 ( T )  =  ω k 3 ( T )  = 0. Furthermore, the optimal control u k * = ( u 1 * , u 2 * , , u n * ) is given by
u k * = min { max [ 0 , I k * ( t ) ( ω k 2 ω k 3 ) B k ] , 0.9 }
Proof. 
We apply the necessary conditions in the Pontryagin maximum principle with Hamiltonian function (21) to derive the adjoint variables determined by the following equations.
d ω k 1 ( t ) d t = λ ( k ) ( 1 δ k ) Θ * ( t ) [ ω k 1 ( t ) ω k 2 ( t ) ] + d ω k 1 ( t ) 1 , d ω k 2 ( t ) d t = λ ( k ) ( 1 δ k ) S k * ( t ) φ ( k ) p ( k ) k ( ω k 1 ( t ) ω k 2 ( t ) ) γ ω k 1 ( t ) + ( u k + γ + d ) ω k 2 ( t ) u k ω k 3 ( t ) 1 , d ω k 3 ( t ) d t = η ω k 1 ( t ) + ( η + d ) ω k 3 ( t ) .
By the optimal conditions, we have
H u k ( t ) = B k u k * ω k 2 I k * + ω k 3 I k * = 0 .
It follows that μ k * = ( ω k 2 ω k 3 ) I k * B k . Using the property of the control space, we obtain
ω * = 0 , ( ω k 2 ω k 3 ) I k * B k < 0 , 0 < ω * < 0.9 , ( ω k 2 ω k 3 ) I k * B k = 0 , ω * = 0.9 , ( ω k 2 ω k 3 ) I k * B k > 0 .
Then, we have the optimal control u k * in compact notation,
u k * = min { max [ 0 , ( ω k 2 ω k 3 ) I k * B k ] , 0.9 }
That completes the proof.  □

6. Simulations

In this section, we conduct some numerical simulations to validate and supplement our theoretical results and to study the impact of the parameters on the dynamics of disease in order to find better control strategies. We first verify and supplement the stability results of system (1) by applying the Ode45 function in Matlab to numerically solve differential equations. It is based on the Runge Kutta method, which can efficiently solve ordinary differential equations and rigid differential equations. Next, we depict the numerical results of the optimal control problem by a combination of forward and backward difference approximation in Matlab to solve the optimal control problem. The numerical simulations include three parts. The first part is to verify the main theoretical results (Theorems 1–4) including the epidemic threshold and global stability, shown in Figure 2 and Figure 3. In the second part, we provide discussion and simulations to show the effectiveness of different immunity and quarantine strategies with degree, as shown in Figure 4 and Figure 5. The numerical simulation in the third part demonstrates the effectiveness of the optimal control method that was proven in Section 5, as shown in Figure 6 and Figure 7. We construct a BA (Barabási–Albert) scale-free network [17], which satisfies a power-law degree distribution p ( k ) C k 3 . The network size N is set as 1000, n = k m a x  = 80. This network evolved from the initial network with a size of m 0 = 3 , and during each time step, a new node with m = 3 edges was added to the network. We note that λ ( k ) = λ k .
(I) Firstly, we perform some simulations in Figure 2 to verify the important results of global stability of equilibria. We denote
φ ( k ) = k , k 30 ; 30 , k > 30 .
Figure 2a shows the relationship between I ( t ) and infectious rate λ on several different networks. Based on the basic reproduction number R 0 , we find that the epidemic threshold of the disease is λ c = k i = 1 n φ ( i ) p ( i ) i N i * ( 1 δ i ) β i + γ + d , where b = 0.2 , d = 0.16 , β k = 0.17 , η = 0.3 , δ k = 0.2 , γ = 0.36 . One can see that the epidemic threshold is consistent with our theoretical results. Figure 2b and Table 2 verify the results of Theorems 3 and 4, i.e., if R 0 < 1 , the disease-free equilibrium is globally stable, and if R 0 > 1 , then the endemic equilibrium persists and is globally stable, where b = 0.2 ; d = 0.1 ; β k = 0.2 ; η = 0.3 ; δ k = 0.2 ; λ  = 0.01, 0.038, 0.055, 0.21, 0.28, 0.35; and γ  = 0.36, 0.36, 0.36, 0.18, 0.18, 0.18. Figure 2b shows that the disease-free equilibrium cannot undergo a Hopf bifurcation.
Figure 3a shows a bifurcation diagram of model system (3). It clearly appears that model system (3) exhibits a forward bifurcation, that is, the disease-free equilibrium D F E is stable if R 0 1 , while if R 0 > 1 , the disease-free equilibrium D F E is unstable, and there exists a unique endemic equilibrium E E which is stable. This result is consistent with the result shown in Figure 2b.
We conduct a parameter sensitivity analysis on the basic reproduction number R 0 , where
R 0 = 1 k k = 1 n φ ( k ) p ( k ) λ ( k ) ( 1 δ k ) N k * β k + γ + d .
Then, let λ ( i ) = λ i , we can obtain
R 0 γ = 1 k i = 1 n φ ( i ) p ( i ) λ ( i ) ( 1 δ i ) N i * ( β i + γ + d ) 2 , R 0 λ = 1 k i = 1 n φ ( i ) p ( i ) i ( 1 δ i ) N i * β i + γ + d .
R 0 β i = 1 k i = 1 n φ ( i ) p ( i ) λ ( i ) ( 1 δ i ) N i * ( β i + γ + d ) 2 , R 0 δ i = 1 k i = 1 n φ ( i ) p ( i ) λ ( i ) N i * β i + γ + d , i = 1 , 2 , , n .
This indicates that the basic reproduction number R 0 increases with the decreases in quarantine rate β k , immunization rate δ k , and recovery rate γ , and R 0 increases with the increase in infection rate λ .
Without losing generality, let β k β , δ k δ . We depict the relationship between the R 0 and parameters in Figure 3b–d, which are consistent with the analysis results. Figure 3b shows that R 0 is positively correlated with the birth rate b and negatively correlated with the death rate d.
(II) Secondly, we depict the effectiveness of the immune and quarantining control strategies in Figure 4 and Figure 5. In this part, we also discuss the SIQS model (1) on the network with various immunization (proportional immunization strategy δ ˜ , target immunization strategy δ k t , acquaintance immunization strategy δ k a c ) schemes based on [51] and define quarantine schemes of infected individuals (proportional quarantine strategy β k , target quarantine strategy β k t , acquaintance quarantining strategy β k a c ) according to various immunization schemes. Then, we define the heterogeneous quarantine rate β k as follows:
Proportional quarantine strategy: Denote the average quarantine rate of proportional quarantine β ˜ , 0 < β ˜ < 1 . In randomly selecting one infected node for isolation, we find the quarantine rate is independent of the degree of the node, which is also a situation discussed in many papers.
Target quarantine strategy: We can devise a quarantine strategy for the infected nodes according to the definition of immunization [51]. Introduce an upper threshold κ , such that all infected nodes with connectivity k > κ are priority quarantined, i.e., we define the quarantine rate β k t by
β k t = 1 , k > κ ; c , k = κ ; 0 , k < κ .
where 0 < c 1 , and k β k t p ( k ) = β k t ˜ , where β k t ˜ is the average quarantine rate of target quarantine.
Acquaintance quarantining strategy: Select a random portion p from the N nodes. The likelihood of quarantining an infected node with degree k is given by k p ( k ) N k ; therefore, β k a c = p k p ( k ) / k . We note that β k a c ˜ denotes the average quarantine rate of acquaintance quarantine, where β k a c ˜ = k β k a c p ( k ) .
Figure 4a shows the I ( t ) under different quarantining control strategies; the target quarantine strategy is the most effective in controlling disease. Figure 4b compares the effectiveness of different quarantine strategies. From Figure 4a,b and Table 3, we find that the target quarantine strategy has better effectiveness than others. Figure 4c shows the effectiveness of different immunity strategies; it shows that all three immunization schemes are effective compared to the case without any immunization, and the targeted immunization scheme is more efficient than the proportional scheme discussed.
Figure 5a depicts the average infectious density with respect to quarantine rate under different immunity strategies. It also shows that the target strategy is better than the proportional and acquaintance strategies. We further show the average infectious density with respect to infected rate for the target and acquaintance immunization schemes and quarantine strategies in Figure 5b. We find that the targeted immunity and the target quarantine strategy have better effectiveness than others.
(III) Thirdly, we analyze and illustrate the optimal control and various quarantine schemes by numerical simulations, as shown in Figure 6 and Figure 7. In Figure 6, we depict I ( t ) under different control strategies. In order to provide a clearer representation of the findings in Figure 6a, we additionally compute the objective function J ( μ ) values for diverse control approaches, as shown in Table 4. In the absence of any control strategy (i.e., u k = 0 ), the infection ultimately breaks out and reaches a stable level of infection. We also set a fixed control strategy compared with the optimal control strategy, finding that if we want to obtain better results than the optimal control strategy, it will usually lead to a significant increase in control costs (i.e., u k m a x = 0.9 ). In Figure 6b, we compare, the effects of optimal control, without control, and acquaintance immunity, combining these with multiple quarantine strategies for infectious disease control. Compared with other strategies, the optimal control strategy has the best control effect.
Figure 7a illustrates the I ( t ) under optimal control under three cases: low cost at B k = 0.01 , moderate cost at B k = 1 , and high cost at B k = 100 . Even under high costs, as seen in the third scenario, the scale of infection can be notably curtailed through optimal control. Figure 7b displays the fluctuation of μ k across varying degrees under optimal regulation. It is evident that μ k tends to escalate as the degree rises. It indicates a higher level of infection, resulting in a higher value of μ k .

7. Conclusions

We proposed an SIQS epidemic model with heterogeneous immune and quarantine measures on heterogeneous networks. We performed a detailed mathematical analysis of our system, which showed that the basic reproduction number R 0 is collectively determined by the infection rate and the effectiveness of disease control measures. Specifically, when R 0 < 1 , the disease-free equilibrium E 0 exhibits global asymptotic stability. When R 0 > 1 , the endemic equilibrium E * demonstrates both local and global asymptotic stability. At the same time, we conducted an in-depth theoretical examination of system (1) to tackle the optimal quarantine control issue, confirming the presence of optimal solutions. We hope this result can provide some ideas for solving the optimal solution of asymmetric quarantine strategies on heterogeneous networks. We defined different forms of quarantine measures, which are asymmetric quarantine strategies related in degree, and numerically simulated the effectiveness of a combination of different quarantine strategies and immune schemes in controlling diseases. Our findings indicate that the implementation of any control strategy can have a favorable impact on the magnitude of an epidemic outbreak, with the optimal control being particularly effective. Among them, the effect of a target quarantine scheme can be comparable to the effect of optimal control. However, target isolation requires knowing the degree of all nodes in the network, which has a significant cost and workload. Therefore, in order to achieve the control goal and reduce the control cost at the same time, the optimal control strategy considering the asymmetry of control measures related in degree and time is shown to be superior to the symmetric control strategies.
A variety of control strategies, such as quarantine and immune schemes, can be employed for the effective management, mitigation, and potential eradication of infectious diseases. This paper can contribute to an understanding of the dynamics of infectious disease with control measures. However, our paper is not based on a specific infectious disease, so there are still certain limitations in providing a specific basis for public health prevention measures. Studying the control measures for a specific disease based on disease data will be meaningful work.

Author Contributions

Conceptualization, S.C. and K.-K.S.; methodology, S.C., Y.W. and L.L.; software, S.C., Y.W. and L.L.; validation, S.C. and Y.W.; formal analysis, S.C. and K.-K.S.; writing—original draft preparation, Y.W. and S.C.; writing—review and editing, Y.W., S.C., D.Y. and K.-K.S.; visualization, Y.W., S.C. and K.-K.S.; supervision, K.-K.S.; funding acquisition, D.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Social Science Funds of China (Grant No. 22BSH025), National Natural Science Foundation of China (Grant No. 62303303 and Grant No. 61803047), Major Project of The National Social Science Foundation of China (19ZDA149, 19ZDA324), and Fundamental Research Funds for the Central Universities (14370119, 14390110). Ke-Ke Shang is supported by the Jiangsu Qing Lan Project.

Data Availability Statement

No new data were created in this study. Data sharing is not applicable to this article.

Acknowledgments

We acknowledge the Theoretical Physics of Complex Systems Center (Nanjing University). We are deeply grateful to Jack Murdoch Moore, a native English speaker, for his significant contributions in revising and enhancing our paper. His Google Scholar profile: https://scholar.google.com/citations?user=AFDBPpYAAAAJ&hl=en (accessed on 26 January 2024). His personal website: https://jackmurdochmoore.github.io/JackMurdochMoore/ (accessed on 26 January 2024).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yang, B.; Shang, K.K.; Michael, S.; Chao, N.P. Information overload: How hot topics distract from news—COVID-19 spread in the US. Natl. Sci. Open 2023, 2, 20220051. [Google Scholar] [CrossRef]
  2. Jose, S.A.; Raja, R.; Omede, B.; Agarwal, R.P.; Alzabut, J.; Cao, J.; Balas, V. Mathematical modeling on co-infection: Transmission dynamics of Zika virus and Dengue fever. Nonlinear Dyn. 2023, 111, 4879–4914. [Google Scholar] [CrossRef]
  3. Peter, O.J.; Afolabi, O.A.; Victor, A.A.; Akpan, C.E.; Oguntolu, F.A. Mathematical model for the control of measles. J. Appl. Sci. Environ. Manag. 2018, 22, 571–576. [Google Scholar] [CrossRef]
  4. Fraser, C.; Donnelly, C.A.; Cauchemez, S.; Hanage, W.P.; Van Kerkhove, M.D.; Hollingsworth, T.D.; Griffin, J.; Baggaley, R.F.; Jenkins, H.E.; Lyons, E.J.; et al. Pandemic potential of a strain of influenza A (H1N1): Early findings. Science 2009, 324, 1557–1561. [Google Scholar] [CrossRef]
  5. Lipsitch, M.; Cohen, T.; Cooper, B.; Robins, J.M.; Ma, S.; James, L.; Gopalakrishna, G.; Chew, S.K.; Tan, C.C.; Samore, M.H.; et al. Transmission dynamics and control of severe acute respiratory syndrome. Science 2003, 300, 1966–1970. [Google Scholar] [CrossRef] [PubMed]
  6. Cai, L.M.; Li, Z.Q.; Song, X.Y. Global analysis of an epidemic model with vaccination. J. Appl. Math. Comput. 2018, 57, 605–628. [Google Scholar] [CrossRef] [PubMed]
  7. Kermack, W.O.; McKendrick, A.G. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. Ser. Contain. Pap. Math. Phys. Character 1927, 115, 700–721. [Google Scholar]
  8. Kermack, W.O.; McKendrick, A.G. Contributions to the mathematical theory of epidemics. II.The problem of endemicity. Proc. R. Soc. Lond. Ser. Contain. Pap. Math. Phys. Character 1932, 138, 55–83. [Google Scholar]
  9. Kermack, W.O.; McKendrick, A.G. Contributions to the mathematical theory of epidemics. III.Further studies of the problem of endemicity. Proc. R. Soc. Lond. Ser. Contain. Pap. Math. Phys. Character 1933, 141, 94–122. [Google Scholar]
  10. Chen, S.S.; Small, M.; Fu, X.C. Global stability of epidemic models with imperfect vaccination and quarantine on scale-free networks. IEEE Trans. Netw. Sci. Eng. 2019, 7, 1583–1596. [Google Scholar] [CrossRef]
  11. Peter, O.J.; Adebisi, A.F.; Ajisope, M.O.; Ajibade, F.O.; Abioye, A.I.; Oguntolu, F.A. Global stability analysis of typhoid fever model. Adv. Syst. Sci. Appl. 2020, 20, 20–31. [Google Scholar]
  12. Jose, S.A.; Raja, R.; Dianavinnarasi, J.; Baleanu, D.; Jirawattanapanit, A. Mathematical modeling of chickenpox in Phuket: Efficacy of precautionary measures and bifurcation analysis. Biomed. Signal Process. Control 2023, 84, 104714. [Google Scholar] [CrossRef]
  13. Hethcote, H.; Zhi En, M.; Sheng Bing, L. Effects of quarantine in six endemic models for infectious diseases. Math. Biosci. 2002, 180, 141–160. [Google Scholar] [CrossRef] [PubMed]
  14. Suo, Y.H. Asymptotical Stability of an SIQS Epidemic Model with Age Dependence and Generally Nonlinear Contact Rate. Appl. Mech. Mater. 2011, 58, 292–297. [Google Scholar] [CrossRef]
  15. Wei, F.Y.; Chen, F.X. Stochastic permanence of an SIQS epidemic model with saturated incidence and independent random perturbations. Phys. A Stat. Mech. Its Appl. 2016, 453, 99–107. [Google Scholar] [CrossRef]
  16. Zhang, X.B.; Liu, R.J. The stationary distribution of a stochastic SIQS epidemic model with varying total population size. Appl. Math. Lett. 2021, 116, 106974. [Google Scholar] [CrossRef]
  17. Barabási, A.L.; Albert, R. Emergence of scaling in random networks. Science 1999, 509–512, 286. [Google Scholar] [CrossRef]
  18. Pastor-Satorras, R.; Vespignani, A. Epidemic dynamics and endemic states in complex networks. Phys. Rev. E 2001, 63, 066117. [Google Scholar] [CrossRef]
  19. Cheng, X.X.; Wang, Y.; Huang, G. Global dynamics of a network-based SIQS epidemic model with nonmonotone incidence rate. Chaos Solitons Fractals 2021, 153, 111502. [Google Scholar] [CrossRef]
  20. Zhao, R.D.; Liu, Q.M.; Sun, M.C. Dynamical behavior of a stochastic SIQS epidemic model on scale-free networks. J. Appl. Math. Comput. 2022, 68, 813–838. [Google Scholar] [CrossRef]
  21. Li, T.; Wang, Y.M.; Guan, Z.H. Spreading dynamics of a SIQRS epidemic model on scale-free networks. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 686–692. [Google Scholar] [CrossRef]
  22. Kang, H.Y.; Liu, K.H.; Fu, X.C. Dynamics of an epidemic model with quarantine on scale-free networks. Phys. Lett. A 2017, 381, 3945–3951. [Google Scholar] [CrossRef]
  23. Wang, H.Y.; Moore, J.M.; Small, M.; Wang, J.; Yang, H.; Gu, C. Epidemic dynamics on higher-dimensional small world networks. Appl. Math. Comput. 2022, 421, 126911. [Google Scholar] [CrossRef]
  24. Iacoviello, D.; Stasio, N. Optimal control for SIRC epidemic outbreak. Comput. Methods Programs Biomed. 2013, 110, 333–342. [Google Scholar] [CrossRef]
  25. Buonomo, B.; Lacitignola, D.; Vargas-De-León, C. Qualitative analysis and optimal control of an epidemic model with vaccination and treatment. Math. Comput. Simul. 2014, 100, 88–102. [Google Scholar] [CrossRef]
  26. Kandhway, K.; Kuri, J. How to run a campaign: Optimal control of SIS and SIR information epidemics. Appl. Math. Comput. 2014, 231, 79–92. [Google Scholar] [CrossRef]
  27. Kandhway, K.; Kuri, J. Optimal control of information epidemics modeled as Maki Thompson rumors. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 4135–4147. [Google Scholar] [CrossRef]
  28. Jang, J.; Kwon, H.D.; Lee, J. Optimal control problem of an SIR reaction–diffusion model with inequality constraints. Math. Comput. Simul. 2020, 171, 136–151. [Google Scholar] [CrossRef]
  29. Wang, B.; Tian, X.H.; Xu, R.; Song, C.W. Threshold dynamics and optimal control of a dengue epidemic model with time delay and saturated incidence. J. Appl. Math. Comput. 2023, 69, 871–893. [Google Scholar] [CrossRef]
  30. Chen, L.J.; Sun, J.T. Global stability and optimal control of an SIRS epidemic model on heterogeneous networks. Phys. A Stat. Mech. Its Appl. 2014, 410, 196–204. [Google Scholar] [CrossRef]
  31. Chen, L.J.; Sun, J.T. Optimal vaccination and treatment of an epidemic network model. Phys. Lett. A 2014, 378, 3028–3036. [Google Scholar] [CrossRef]
  32. Xu, D.G.; Xu, X.Y.; Xie, Y.F.; Yang, C.H. Optimal control of an SIVRS epidemic spreading model with virus variation based on complex networks. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 200–210. [Google Scholar] [CrossRef]
  33. Jia, N.; Ding, L.; Liu, Y.J.; Hu, P. Global stability and optimal control of epidemic spreading on multiplex networks with nonlinear mutual interaction. Phys. A Stat. Mech. Its Appl. 2018, 502, 93–105. [Google Scholar] [CrossRef]
  34. Li, K.Z.; Zhu, G.H.; Ma, Z.J.; Chen, L.J. Dynamic stability of an SIQS epidemic network and its optimal control. Commun. Nonlinear Sci. Numer. Simul. 2019, 66, 84–95. [Google Scholar] [CrossRef]
  35. Wei, X.D.; Xu, G.C.; Zhou, W.S. Global stability of endemic equilibrium for a SIQRS epidemic model on complex networks. Phys. A Stat. Mech. Its Appl. 2018, 512, 203–214. [Google Scholar] [CrossRef]
  36. Zhang, L.; Liu, M.X.; Xie, B.L. Optimal control of an SIQRS epidemic model with three measures on networks. Nonlinear Dyn. 2021, 103, 2097–2107. [Google Scholar] [CrossRef]
  37. Yang, P.; Jia, J.B.; Shi, W.; Feng, J.W.; Fu, X.C. Stability analysis and optimal control in an epidemic model on directed complex networks with nonlinear incidence. Commun. Nonlinear Sci. Numer. Simul. 2023, 121, 107206. [Google Scholar] [CrossRef]
  38. Pastor-Satorras, R.; Vespignani, A. Epidemic dynamics in finite size scale-free networks. Phys. Rev. E 2002, 65, 035108. [Google Scholar] [CrossRef]
  39. Moreno, Y.; Pastor-Satorras, R.; Vespignani, A. Epidemic outbreaks in complex heterogeneous networks. Eur. Phys. J.-Condens. Matter Complex Syst. 2002, 26, 521–529. [Google Scholar] [CrossRef]
  40. Wang, L.; Dai, G.Z. Global stability of virus spreading in complex heterogeneous networks. SIAM J. Appl. Math. 2008, 68, 1495–1502. [Google Scholar] [CrossRef]
  41. Yang, R.; Wang, B.H.; Ren, J.; Bai, W.J.; Shi, Z.W.; Wang, W.X.; Zhou, T. Epidemic spreading on heterogeneous networks with identical infectivity. Phys. Lett. A 2007, 364, 189–193. [Google Scholar] [CrossRef]
  42. Chu, X.W.; Zhang, Z.Z.; Guan, J.H.; Zhou, S.G. Epidemic spreading with nonlinear infectivity in weighted scale-free networks. Phys. A Stat. Mech. Its Appl. 2011, 390, 471–481. [Google Scholar] [CrossRef]
  43. Zhang, H.F.; Fu, X.C. Spreading of epidemics on scale-free networks with nonlinear infectivity. Nonlinear Anal. Theory Methods Appl. 2009, 70, 3273–3278. [Google Scholar] [CrossRef]
  44. Zhu, G.H.; Chen, G.R.; Xu, X.J.; Fu, X.C. Epidemic spreading on contact networks with adaptive weights. J. Theor. Biol. 2013, 317, 133–139. [Google Scholar] [CrossRef] [PubMed]
  45. Liu, X.N.; Takeuchi, Y.; Iwami, S. SVIR epidemic models with vaccination strategies. J. Theor. Biol. 2008, 253, 1–11. [Google Scholar] [CrossRef] [PubMed]
  46. Lajmanovich, A.; Yorke, J.A. A deterministic model for gonorrhea in a nonhomogeneous population. Math. Biosci. 1976, 28, 221–236. [Google Scholar] [CrossRef]
  47. Culshaw, R.V.; Ruan, S.; Spiteri, R.J. Optimal HIV treatment by maximising immune response. J. Math. Biol. 2004, 48, 545–562. [Google Scholar] [CrossRef] [PubMed]
  48. Chen, P.Y.; Cheng, S.M.; Chen, K.C. Optimal control of epidemic information dissemination over networks. IEEE Trans. Cybern. 2014, 44, 2316–2328. [Google Scholar] [CrossRef]
  49. Abboubakar, H.; Guidzavai, A.K.; Yangla, J.; Damakoa, I.; Mouangue, R. Mathematical modeling and projections of a vector-borne disease with optimal control strategies: A case study of the Chikungunya in Chad. Chaos Solitons Fractals 2021, 150, 111197. [Google Scholar] [CrossRef]
  50. Fleming, W.; Rishel, R. Deterministic and Stochastic Optimal Control; Springer: Berlin/Heidelberg, Germany, 1975; p. 222. [Google Scholar]
  51. Fu, X.C.; Small, M.; Walker, D.M.; Zhang, H.F. Epidemic dynamics on scale-free networks with piecewise linear infectivity and immunization. Phys. Rev. E 2008, 77, 036113. [Google Scholar] [CrossRef]
Figure 1. The block diagram of the SIQS model. Here S, I, and Q represent susceptible, infected, and quarantined states.
Figure 1. The block diagram of the SIQS model. Here S, I, and Q represent susceptible, infected, and quarantined states.
Symmetry 16 00166 g001
Figure 2. (a) The I ( t ) with respect to λ and the structure of a network with different parameters m and m 0 . λ c denotes the theoretical epidemic threshold. Other parameter settings are as follows: b = 0.2 , d = 0.16 , β k = 0.17 , η = 0.3 , δ k = 0.2 , γ = 0.36 . (b) The average infected density I ( t ) under different parameters: λ = 0.01, 0.038, 0.055, 0.21, 0.28, 0.3; γ  = 0.36, 0.36, 0.36, 0.18, 0.18, 0.18; b = 0.2 ; d = 0.1 ; β k = 0.2 ; η = 0.3 ; and δ k = 0.2 , which correspond to R 0 = 0.1734, 0.6590, 0.9538 < 1; R 0 = 5.3671, 7.1561, 8.9451 > 1.
Figure 2. (a) The I ( t ) with respect to λ and the structure of a network with different parameters m and m 0 . λ c denotes the theoretical epidemic threshold. Other parameter settings are as follows: b = 0.2 , d = 0.16 , β k = 0.17 , η = 0.3 , δ k = 0.2 , γ = 0.36 . (b) The average infected density I ( t ) under different parameters: λ = 0.01, 0.038, 0.055, 0.21, 0.28, 0.3; γ  = 0.36, 0.36, 0.36, 0.18, 0.18, 0.18; b = 0.2 ; d = 0.1 ; β k = 0.2 ; η = 0.3 ; and δ k = 0.2 , which correspond to R 0 = 0.1734, 0.6590, 0.9538 < 1; R 0 = 5.3671, 7.1561, 8.9451 > 1.
Symmetry 16 00166 g002
Figure 3. (a) Bifurcation diagram of model system (1). EE is the endemic equilibrium E * , DFE is the disease-free equilibrium E 0 : b = 0.12 , d = 0.06 , γ = 0.28 , η = 0.3 , β k = 0.6 , δ k = 0.2 . (bd) The combined influence of parameters on R 0 . (b) The relationship between R 0 and b, d: λ = 0.15 , γ = 0.36 , η = 0.3 , β k = 0.2 , δ k = 0.2 . (c) The relationship between R 0 and β , δ : b = 0.12 , d = 0.06 , λ = 0.15 , γ = 0.3 , η = 0.3 . (d) The relationship between R 0 and λ , γ : b = 0.12 , d = 0.06 , η = 0.3 , β k = 0.75 , δ k = 0.8 .
Figure 3. (a) Bifurcation diagram of model system (1). EE is the endemic equilibrium E * , DFE is the disease-free equilibrium E 0 : b = 0.12 , d = 0.06 , γ = 0.28 , η = 0.3 , β k = 0.6 , δ k = 0.2 . (bd) The combined influence of parameters on R 0 . (b) The relationship between R 0 and b, d: λ = 0.15 , γ = 0.36 , η = 0.3 , β k = 0.2 , δ k = 0.2 . (c) The relationship between R 0 and β , δ : b = 0.12 , d = 0.06 , λ = 0.15 , γ = 0.3 , η = 0.3 . (d) The relationship between R 0 and λ , γ : b = 0.12 , d = 0.06 , η = 0.3 , β k = 0.75 , δ k = 0.8 .
Symmetry 16 00166 g003
Figure 4. (a) I ( t ) under different quarantining control strategies: no quarantine (cyan line), acquaintance quarantine (magenta line), proportional quarantine (blue line), target quarantine (green line): b = 0.2 , d = 0.08 , η = 0.3 , λ = 0.28 ,   δ k = 0.2 ,   γ = 0.18 ,   c = 0.8 ,   κ = 15 , p = 0.386 , β k = β k a c ˜ = β k t ˜ = β ˜ = 0.0499 . (b) The epidemic threshold λ c under different quarantine control strategies: no quarantine (green line), target quarantine (magenta line), acquaintance quarantine (blue line), proportional quarantine (red line): b = 0.2 , d = 0.18 , δ k = 0.2 , γ = 0.5 , η = 0.35 , c = 0.8 ,   κ = 10 , p = 0.8946 , β k = β k a c ˜ = β k t ˜ = β ˜ = 0.1150 . (c) The epidemic threshold λ c under different immunization strategies: no immunity (blue line), target immunity (green line), proportional immunity (cyan line), acquaintance immunity (red line): b = 0.2 , d = 0.18 , γ = 0.5 , η = 0.35 , β k = 0.21 , c = 1 , κ = 7 , p = 0.9836 , δ k = δ k a c ˜ = δ k t ˜ = δ ˜ = 0.218 .
Figure 4. (a) I ( t ) under different quarantining control strategies: no quarantine (cyan line), acquaintance quarantine (magenta line), proportional quarantine (blue line), target quarantine (green line): b = 0.2 , d = 0.08 , η = 0.3 , λ = 0.28 ,   δ k = 0.2 ,   γ = 0.18 ,   c = 0.8 ,   κ = 15 , p = 0.386 , β k = β k a c ˜ = β k t ˜ = β ˜ = 0.0499 . (b) The epidemic threshold λ c under different quarantine control strategies: no quarantine (green line), target quarantine (magenta line), acquaintance quarantine (blue line), proportional quarantine (red line): b = 0.2 , d = 0.18 , δ k = 0.2 , γ = 0.5 , η = 0.35 , c = 0.8 ,   κ = 10 , p = 0.8946 , β k = β k a c ˜ = β k t ˜ = β ˜ = 0.1150 . (c) The epidemic threshold λ c under different immunization strategies: no immunity (blue line), target immunity (green line), proportional immunity (cyan line), acquaintance immunity (red line): b = 0.2 , d = 0.18 , γ = 0.5 , η = 0.35 , β k = 0.21 , c = 1 , κ = 7 , p = 0.9836 , δ k = δ k a c ˜ = δ k t ˜ = δ ˜ = 0.218 .
Symmetry 16 00166 g004
Figure 5. The effectiveness of combination immunization and quarantine control strategies: b = 0.2 , d = 0.18 , c = 1 , κ = 7 , p = 0.9836 , δ k = δ k a c ˜ = δ k t ˜ = δ ˜ = 0.218 . (a) Different immunization strategies with respect to quarantine rate: λ = 0.2 , γ = 0.38 , η = 0.3 ; acquaintance immunization (magenta line), proportional immunization (blue line), target immunization (cyan line). (b) Comparison of different combination control strategies: γ = 0.5 , η = 0.35 , c = 0.8 , κ = 10 , p = 0.8946 , β k a c ˜ ( δ k a c ˜ ) = β k t ˜ ( δ k t ˜ ) = β ˜ ( δ ˜ ) = 0.1150 ; target immunization and target quarantine (red line); target immunization and acquaintance quarantine (magenta line); acquaintance immunization and target quarantine (cyan line); acquaintance immunization and acquaintance quarantine (blue line).
Figure 5. The effectiveness of combination immunization and quarantine control strategies: b = 0.2 , d = 0.18 , c = 1 , κ = 7 , p = 0.9836 , δ k = δ k a c ˜ = δ k t ˜ = δ ˜ = 0.218 . (a) Different immunization strategies with respect to quarantine rate: λ = 0.2 , γ = 0.38 , η = 0.3 ; acquaintance immunization (magenta line), proportional immunization (blue line), target immunization (cyan line). (b) Comparison of different combination control strategies: γ = 0.5 , η = 0.35 , c = 0.8 , κ = 10 , p = 0.8946 , β k a c ˜ ( δ k a c ˜ ) = β k t ˜ ( δ k t ˜ ) = β ˜ ( δ ˜ ) = 0.1150 ; target immunization and target quarantine (red line); target immunization and acquaintance quarantine (magenta line); acquaintance immunization and target quarantine (cyan line); acquaintance immunization and acquaintance quarantine (blue line).
Symmetry 16 00166 g005
Figure 6. (a) The average density I ( t ) under different control strategies. The lines with different colors correspond to different optimal control strategies: u k ( t ) =  optimal control (purple line); u k ( t ) = 0.9 (cyan line); u k ( t ) = 0 (green line). Other parameter settings are as follows: b = 0.09 , d = 0.012 , λ = 0.3 , γ = 0.01 , η = 0.03 , δ k = 0.1 . (b) Comparison of optimal control and combined heterogeneous control strategies, without control (green line); acquaintance immunization and acquaintance quarantine (blue line); acquaintance immunization and proportional quarantine (magenta line); acquaintance immunization and target quarantine (cyan line); optimal control (purple line): b = 0.09 , d = 0.012 , λ = 0.3 , γ = 0.01 , δ = 0.218 , c = 1 , κ = 7 , p = 0.9836 , β k = β k a c ˜ = β k t ˜ = β ˜ = 0.218 .
Figure 6. (a) The average density I ( t ) under different control strategies. The lines with different colors correspond to different optimal control strategies: u k ( t ) =  optimal control (purple line); u k ( t ) = 0.9 (cyan line); u k ( t ) = 0 (green line). Other parameter settings are as follows: b = 0.09 , d = 0.012 , λ = 0.3 , γ = 0.01 , η = 0.03 , δ k = 0.1 . (b) Comparison of optimal control and combined heterogeneous control strategies, without control (green line); acquaintance immunization and acquaintance quarantine (blue line); acquaintance immunization and proportional quarantine (magenta line); acquaintance immunization and target quarantine (cyan line); optimal control (purple line): b = 0.09 , d = 0.012 , λ = 0.3 , γ = 0.01 , δ = 0.218 , c = 1 , κ = 7 , p = 0.9836 , β k = β k a c ˜ = β k t ˜ = β ˜ = 0.218 .
Symmetry 16 00166 g006
Figure 7. (a) I ( t ) under optimal control with varying weights. Low costs: B k = 0.01 (red line); moderate costs: B k = 1 (blue line); high costs: B k = 100 (cyan line). Other parameters are fixed as b = 0.09 , d = 0.012 , λ = 0.3 , γ = 0.008 , η = 0.03 , δ k = 0.1 . (b) The quarantine rates u k ( t ) for any degree. Other parameters are fixed as b = 0.09 , d = 0.02 , λ = 0.3 , γ = 0.01 , η = 0.03 , δ k = 0.1 .
Figure 7. (a) I ( t ) under optimal control with varying weights. Low costs: B k = 0.01 (red line); moderate costs: B k = 1 (blue line); high costs: B k = 100 (cyan line). Other parameters are fixed as b = 0.09 , d = 0.012 , λ = 0.3 , γ = 0.008 , η = 0.03 , δ k = 0.1 . (b) The quarantine rates u k ( t ) for any degree. Other parameters are fixed as b = 0.09 , d = 0.02 , λ = 0.3 , γ = 0.01 , η = 0.03 , δ k = 0.1 .
Symmetry 16 00166 g007
Table 1. Symbols employed in models.
Table 1. Symbols employed in models.
SymbolsDescription
p ( k ) Proportion of nodes with degree k.
k Average degree   k = k k p ( k ) .
nMaximum degree.
bBirth rate.
dNatural death rate.
Ψ k ( t ) Fertile contact probability between a node with degree k and its neighbors.
λ ( k ) Transmission rate of infected nodes with degree k.
δ k Vaccination rate of susceptible nodes with degree k.
β k Quarantine rate of infected nodes with degree k.
γ Recovery rate of infected nodes.
η Recovery rate of quarantined nodes.
Table 2. Final infection scale I ( t ) under different R 0 in Figure 2b.
Table 2. Final infection scale I ( t ) under different R 0 in Figure 2b.
Basic reproduction number R 0 0.17340.65900.95385.36717.15618.9451
Final infection scale I ( t ) 0000.04490.05970.0724
Table 3. Final infection scale and λ c under different control strategies in Figure 4a,b.
Table 3. Final infection scale and λ c under different control strategies in Figure 4a,b.
Quarantining Control StrategiesFinal Infection Scale I ( t ) Epidemic Threshold λ c
No quarantine0.07130.0468
Acquaintance quarantine0.05680.0724
Proportional quarantine0.04240.0798
Target quarantine0.01820.1483
Table 4. The values of J ( μ ) under different control strategies.
Table 4. The values of J ( μ ) under different control strategies.
Optimal ControlMax ControlNo Control
J 1.7817 × 10 3 J 6.5164 × 10 4 J 5.1127 × 10 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, Y.; Chen, S.; Yu, D.; Liu, L.; Shang, K.-K. Propagation Dynamics of an Epidemic Model with Heterogeneous Control Strategies on Complex Networks. Symmetry 2024, 16, 166. https://doi.org/10.3390/sym16020166

AMA Style

Wang Y, Chen S, Yu D, Liu L, Shang K-K. Propagation Dynamics of an Epidemic Model with Heterogeneous Control Strategies on Complex Networks. Symmetry. 2024; 16(2):166. https://doi.org/10.3390/sym16020166

Chicago/Turabian Style

Wang, Yan, Shanshan Chen, Dingguo Yu, Lixiang Liu, and Ke-Ke Shang. 2024. "Propagation Dynamics of an Epidemic Model with Heterogeneous Control Strategies on Complex Networks" Symmetry 16, no. 2: 166. https://doi.org/10.3390/sym16020166

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop