Next Article in Journal
Protein Is an Intelligent Micelle
Previous Article in Journal
A Hybrid Parallel Balanced Phasmatodea Population Evolution Algorithm and Its Application in Workshop Material Scheduling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Community Connectivity on the Spreading Process of Epidemics

1
School of Mathematics and Statistics, Tianshui Normal University, Tianshui 741000, China
2
School of Mathematics and Data Science, Shaanxi University of Science & Technology, Xi’an 710021, China
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(6), 849; https://doi.org/10.3390/e25060849
Submission received: 12 April 2023 / Revised: 21 May 2023 / Accepted: 23 May 2023 / Published: 26 May 2023
(This article belongs to the Section Complexity)

Abstract

:
Community structure exists widely in real social networks. To investigate the effect of community structure on the spreading of infectious diseases, this paper proposes a community network model that considers both the connection rate and the number of connected edges. Based on the presented community network, a new SIRS transmission model is constructed via the mean-field theory. Furthermore, the basic reproduction number of the model is calculated via the next-generation matrix method. The results reveal that the connection rate and the number of connected edges of the community nodes play crucial roles in the spreading process of infectious diseases. Specifically, it is demonstrated that the basic reproduction number of the model decreases as the community strength increases. However, the density of infected individuals within the community increases as the community strength increases. For community networks with weak strength, infectious diseases are likely not to be eradicated and eventually will become endemic. Therefore, controlling the frequency and range of intercommunity contact will be an effective initiative to curb outbreaks of infectious diseases throughout the network. Our results can provide a theoretical basis for preventing and controlling the spreading of infectious diseases.

1. Introduction

Throughout history, viruses have been a major factor in the development of societies. Describing the mechanisms of epidemics spreading and predicting their epidemiological trends is the long-standing focus of research with the aim of effectively controlling the spread of diseases [1,2,3].
As we know, the study of complex networks has gradually become a hot issue in the field of complexity disciplines. Scholars have made significant contributions to the study of areas such as transportation, social, financial, and biological networks. As research into complex networks has continued, the spreading of computer viruses in computer networks, contagious diseases in social populations, and public opinion and rumors in social networks can have an egregious effect on the development of human society [4,5,6]. Therefore, the behavior of transmission dynamics on complex networks has become one of the research directions of great interest.
Following the introduction of small-world networks and scale-free networks [7,8,9,10], scientists have studied the many networks that exist in the real world and found that complex networks also have characteristics of community structure. Girvan and Newman have introduced community structure properties in many networks and proposed methods for detecting such structures [11]. The community structure means that the nodes in the network are divided into several groups; i.e., the individuals within the groups are relatively tightly connected, and the individuals between the groups are relatively sparsely connected. In real social networks, individuals form groups because of the same characteristics. In turn, scientists study the dynamics of disease transmission in networks by constructing community network models to mirror real networks [12,13,14,15,16,17,18,19,20,21,22,23,24,25].
Ref. [26] defined a parameter to represent the degree of communities to study the effect of community structure on propagation dynamics. Ultimately, it was found that community networks have wider degree distributions and smaller epidemic thresholds compared with random networks. Newman and Girvan defined the modularity coefficient to determine the strength of the community structure network [27]. Later, Salathe and Jones studied the impact of interventions on the spread of disease between communities and found that immunization interventions for individuals within communities in a strong community-structured network were a more effective way to control the spread of disease [28]. Authors investigated two models of community networks with different structures and concluded that the community structure can both prompt and inhibit the spread of viruses, while the intracommunity structure does not affect the spread behavior between communities [29]. Jean-Gabriel studied the emergence of community structures and the network model of structural evolution within communities [30]. Li and Jiang considered a disease model with community heterogeneity and found that the community heterogeneity affects the transmission threshold and disease prevalence rate [31].
Although network propagation models with community structures have been studied very extensively [32,33], most of the models ignore the effect of the number of connected edges within and among communities. To understand the influence of community connectivity on dynamic processes, this paper further studies the epidemic spread model on dynamic networks with community structures. Moreover, demographics play a crucial role in the disease-spreading process. Demographics change the number and internal relationships of individuals in different communities. Therefore, we investigate the effect of the connection rate and the number of connected edges among communities on epidemic propagation in the proposed model.
In fact, in real social networks, the physical contacts among individuals in some communities are closer than in other communities, reflecting the connections among communities. For example, in real social networks, the young-student community interacts with each other more frequently than the elderly-people community. Furthermore, studies have suggested that coupling strength in contact patterns among individuals in different communities has an important effect on epidemic spreading. We firstly construct a new community structure network model with the connection rate and the number of connected edges among communities. Then we establish a novel epidemic spreading model based on this network model and calculate the epidemic threshold.
This paper is organized as follows. In Section 2, a community growth network that considers the number of connected edges is presented. Moreover, we introduce a modularity coefficient to determine the strength of the community structure. Section 3 focuses on the SIR epidemic model via applying the mean-field approach and calculating the basic reproduction number. In Section 4, we give numerical simulations to verify the theoretical results. Finally, Section 5 concludes the paper.

2. Community Network Model

The network model generation algorithm with community structure can be constructed as follows [26,29].
Consider a network containing m communities, and each community contains n 1 , n 2 , , n m nodes, respectively. Here, we assume that the number of nodes in each community is n i .
At the initial moment, the nodes of the network are independent. Each node is connected to m 1 nodes with probability α of intracommunity links, and every two communities have m 2 links.
At each time step, each node in the i-th community is connected to m 3 nodes with probability p i n in the same community. Moreover, this node is connected to m 4 nodes with probability p o u t between different communities. Then, a community network is generated.
In what follows, we discuss the effect of the strength of the community structure on the spreading of the epidemic by introducing a modularity coefficient. The modularity coefficient is defined as [27]
Q = i [ e i i ( j e i j ) 2 ] ,
where e i j is the proportion of connected edges in community i and j to the total connected edges in the network.
According to the above network evolution, one can obtain
e i i = n α m 1 + n p i n m 3 n α m m 1 + 1 2 m ( m 1 ) m 2 + n p i n m m 3 + n p o u t m ( m 1 ) m 4 .
j e i j = e i i + j i e i j = n α m 1 + n p i n m 3 + ( m 1 ) m 2 + n p o u t ( m 1 ) m 4 n α m m 1 + 1 2 m ( m 1 ) m 2 + n p i n m m 3 + n p o u t m ( m 1 ) m 4 .
Substituting Equations (2) and (3) into (1), we obtain
Q = ( m 1 ) ( n α m 1 + n p i n m 3 ) 2 + ( n α m 1 + n p i n m 3 ) A ( ( m 1 ) m 2 + n p o u t ( m 1 ) m 4 ) 2 ( n α m m 1 + 1 2 m ( m 1 ) m 2 + n p i n m m 3 + n p o u t m ( m 1 ) m 4 ) 2 .
It can be seen from Equation (4) that the community strength is not only related to the connected rate p i n , p o u t but also to the number of intracommunity and intercommunity links. The following numerical simulations are performed. The parameters are taken as m = 4 , n = 500 , and α = 0.1 .
Figure 1 displays the relation between the community strength and the internal and external connection rates of the nodes. From Figure 1, we observe that the same modularity coefficient corresponds to different p i n and p o u t . Moreover, it can be seen that the internal connection rate p i n increases as the external connection rate p o u t increases.
Figure 2a shows that when the external connection rate p o u t is fixed, the modularity coefficient increases as the internal connection rate p i n increases. Figure 2b shows that when the internal connection rate p i n is fixed, the modularity coefficient decreases as the external connection rate p o u t increases. It implies that an increase in the frequency of connections between communities leads to a less distinct community structure in the network.
Figure 3a,c display the effect of the number of connected edges of each node within the community on the modularity coefficient. The result is that the higher the number of connected edges, the stronger the community strength. Moreover, Figure 3b,d reflect the increase in the number of connected edges between communities, thus making the community structure inconspicuous. However, Figure 3a,b show that the number of connected edges at the initial moment has a smaller effect on the modularity coefficient. This means that it will not play a decisive role in the strength of the community structure.

3. Infectious Disease Dynamics Model and Analysis

3.1. Model Description

In this subsection, we consider an epidemic network model with two communities. It is assumed that the individuals in each community will be in three different states: susceptible (S), infected (I), and recovered (R). Parameters S i , k , I i , k , and R i , k are defined as the density of susceptible, infected, and recovered individuals with degree k at time t, where i indicates the i - th community, i = 1 , 2 . Table 1 denotes the explanation of the key parameters in the transmission of infectious diseases. Here, we assume that all parameters are positive.
In order to further explain, the transformation between the states of the infectious disease in the two communities is described as follows. The Figure 4 shows the infectious disease transmission in complex network with two communities.
Based on the above community network model, the SIR epidemic transmission model is established as the following equations:
d S 1 , k ( t ) d t = β 1 p i n k S 1 , k Θ I 1 ( t ) β 1 α k S 1 , k Θ I 1 ( t ) γ 12 β 2 p o u t k S 1 , k Θ I 2 ( t ) + μ 1 R 1 , k ( t ) d I 1 , k ( t ) d t = β 1 p i n k S 1 , k Θ I 1 ( t ) + β 1 α k S 1 , k Θ I 1 ( t ) + γ 12 β 2 p o u t k S 1 , k Θ I 2 ( t ) δ 1 I 1 , k ( t ) d R 1 , k ( t ) d t = δ 1 I 1 , k ( t ) μ 1 R 1 , k ( t ) d S 2 , k ( t ) d t = β 2 p i n k S 2 , k Θ I 2 ( t ) β 2 α k S 2 , k Θ I 2 ( t ) γ 21 β 1 p o u t k S 2 , k Θ I 1 ( t ) + μ 2 R 2 , k ( t ) d I 2 , k ( t ) d t = β 2 p i n k S 2 , k Θ I 2 ( t ) + β 2 α k S 2 , k Θ I 2 ( t ) + γ 21 β 1 p o u t k S 2 , k Θ I 1 ( t ) δ 2 I 2 , k ( t ) d R 2 , k ( t ) d t = δ 2 I 2 , k ( t ) μ 2 R 2 , k ( t )
where n i , k represents the total number of individuals of degree k in the i - th community, k = 1 , 2 , , d , d is the maximum degree, n i denotes the whole number of individuals in the i - th community, and N is the total number of individuals in the network, satisfying N = n 1 + n 2 . We define the power-law degree distribution in the whole network as
p ( k ) = n 1 , k + n 2 , k n 1 + n 2 .
It is well known that the degree distribution in the i - th community satisfies
p i ( k ) = n i , k n i ,
and the average degree in the whole network is
k = k k p ( k ) = k k n 1 , k + n 2 , k n 1 + n 2 = 1 N k k ( n 1 , k + n 2 , k ) . = 1 N k k ( n 1 p 1 ( k ) + n 2 p 2 ( k ) )
In Equation (5), Θ I i ( t ) = 1 k k k p ( k ) I i , k ( t ) n i , k n 1 , k + n 2 , k denotes the probability that a randomly chosen link of a node is connected to an infected node at time t :
Θ I 1 ( t ) = 1 k k k p ( k ) I 1 , k ( t ) n 1 , k n 1 , k + n 2 , k = 1 k k k n 1 , k + n 2 , k n 1 + n 2 I 1 , k ( t ) n 1 , k n 1 , k + n 2 , k = 1 k k k I 1 , k ( t ) n 1 p 1 ( k ) N = n 1 N 1 1 N ( k 1 n 1 + k 2 n 2 ) k k p 1 ( k ) I 1 , k ( t ) = a k k p 1 ( k ) I 1 , k ( t )
where a = n 1 k 1 n 1 + k 2 n 2 . Similarly,
Θ I 2 ( t ) = b k k p 2 ( k ) I 2 , k ( t ) ,
where b = n 2 k 1 n 1 + k 2 n 2 .

3.2. The Basic Reproduction Number

This subsection focuses on the investigation of the basic reproduction number R 0 via the next-generation matrix scheme. In order to obtain the infection matrix F and the internal evolution matrix V , we translate Equation (5) into the following format:
d I 1 , k ( t ) d t = β 1 p i n k S 1 , k Θ I 1 ( t ) + β 1 α k S 1 , k Θ I 1 ( t ) + γ 12 β 2 p o u t k S 1 , k Θ I 2 ( t ) δ 1 I 1 , k ( t ) d I 2 , k ( t ) d t = β 2 p i n k S 2 , k Θ I 2 ( t ) + β 2 α k S 2 , k Θ I 2 ( t ) + γ 21 β 1 p o u t k S 2 , k Θ I 1 ( t ) δ 2 I 2 , k ( t ) d S 1 , k ( t ) d t = β 1 p i n k S 1 , k Θ I 1 ( t ) β 1 α k S 1 , k Θ I 1 ( t ) γ 12 β 2 p o u t k S 1 , k Θ I 2 ( t ) + μ 1 R 1 , k ( t ) d S 2 , k ( t ) d t = β 2 p i n k S 2 , k Θ I 2 ( t ) β 2 α k S 2 , k Θ I 2 ( t ) γ 21 β 1 p o u t k S 2 , k Θ I 1 ( t ) + μ 2 R 2 , k ( t ) d R 1 , k ( t ) d t = δ 1 I 1 , k ( t ) μ 1 R 1 , k ( t ) d R 2 , k ( t ) d t = δ 2 I 2 , k ( t ) μ 2 R 2 , k ( t )
where the infection matrix F and the internal evolution matrix V are given:
F = β 1 p i n k S 1 , k Θ I 1 ( t ) + β 1 α k S 1 , k Θ I 1 ( t ) + γ 12 β 2 p o u t k S 1 , k Θ I 2 ( t ) β 2 p i n k S 2 , k Θ I 2 ( t ) + β 2 α k S 2 , k Θ I 2 ( t ) + γ 21 β 1 p o u t k S 2 , k Θ I 1 ( t ) 0 0 0 0 V = δ 1 I 1 , k ( t ) δ 2 I 2 , k ( t ) β 1 p i n k S 1 , k Θ I 1 ( t ) + β 1 α k S 1 , k Θ I 1 ( t ) + γ 12 β 2 p o u t k S 1 , k Θ I 2 ( t ) μ 1 R 1 , k ( t ) β 2 p i n k S 2 , k Θ I 2 ( t ) + β 2 α k S 2 , k Θ I 2 ( t ) + γ 21 β 1 p o u t k S 2 , k Θ I 1 ( t ) μ 2 R 2 , k ( t ) δ 1 I 1 , k ( t ) + μ 1 R 1 , k ( t ) δ 2 I 2 , k ( t ) + μ 2 R 2 , k ( t ) .
Then, the infection matrix of the system at the disease-free equilibrium E 0 ( 1 , 0 , 0 , , 1 , 0 , 0 ) can be obtained:
F 11 = β 1 ( p i n + α ) a p 1 ( 1 ) β 1 ( p i n + α ) a 2 p 1 ( 2 ) β 1 ( p i n + α ) a d p 1 ( d ) β 1 ( p i n + α ) 2 a p 1 ( 1 ) β 1 ( p i n + α ) 2 a 2 p 1 ( 2 ) β 1 ( p i n + α ) 2 a d p 1 ( 1 ) β 1 ( p i n + α ) d a p 1 ( 1 ) β 1 ( p i n + α ) d a 2 p 1 ( 1 ) β 1 ( p i n + α ) d a d p 1 ( d ) F 12 = γ 12 β 2 p o u t b p 2 ( 1 ) γ 12 β 2 p o u t b d p 2 ( d ) γ 12 β 2 p o u t d b p 2 ( 1 ) γ 12 β 2 p o u t d b d p 2 ( d ) F 21 = γ 21 β 1 p o u t a p 1 ( 1 ) γ 21 β 1 p o u t a d p 1 ( d ) γ 21 β 1 p o u t d a p 1 ( 1 ) γ 21 β 1 p o u t d a d p 1 ( d ) F 22 = β 2 ( p i n + α ) b p 2 ( 1 ) β 2 ( p i n + α ) b d p 2 ( d ) β 2 ( p i n + α ) d b p 2 ( 1 ) β 2 ( p i n + α ) d b d p 2 ( d )
where F i j d × d , i = 1 , 2 , j = 3 , , 6 , and F i j d × d , i = 3 , , 6 , and j = 1 , 2 , , 6 are zero matrix.
The internal evolution matrix of the system at the disease-free equilibrium E 0 ( 1 , 0 , 0 , , 1 , 0 , 0 ) is given by
V 11 = δ 1 0 0 δ 1 V 22 = δ 2 0 0 δ 2 V 31 = β 1 ( p i n + α ) b p 1 ( 1 ) β 1 ( p i n + α ) b d p 1 ( d ) β 1 ( p i n + α ) d b p 1 ( 1 ) β 1 ( p i n + α ) d b d p 1 ( d ) V 32 = γ 12 β 2 p o u t b p 2 ( 1 ) γ 12 β 2 p o u t b d p 2 ( d ) γ 12 β 2 p o u t d b p 2 ( 1 ) γ 12 β 2 p o u t d b d p 2 ( d ) V 35 = μ 1 0 0 μ 1 V 41 = γ 21 β 1 p o u t a p 1 ( 1 ) γ 21 β 1 p o u t a d p 1 ( d ) γ 21 β 1 p o u t d a p 1 ( 1 ) γ 21 β 1 p o u t d a d p 1 ( d ) V 42 = β 2 ( p i n + α ) b p 2 ( 1 ) β 2 ( p i n + α ) b d p 2 ( d ) β 2 ( p i n + α ) d b p 2 ( 1 ) β 2 ( p i n + α ) d b d p 2 ( d ) V 46 = μ 2 0 0 μ 2   V 51 = δ 1 0 0 δ 1   V 55 = μ 1 0 0 μ 1 V 62 = δ 2 0 0 δ 2   V 66 = μ 2 0 0 μ 2
where V 1 j d × d , j = 2 , , 6 , V 2 j d × d , j = 1 , 3 , , 6 , V 3 j d × d , j = 3 , 4 , 6 , V 4 j d × d , j = 3 , 4 , 5 , V 5 j d × d , j = 2 , 3 , 4 , 6 , V 6 j d × d , and j = 1 , 3 , 4 , 5 are zero matrices according to the principle of the next-generation matrix, where R 0 = ρ ( F V 1 ) , which denotes the spectral radius of the matrix F V 1 .
Defining C = F V 1 , we have C 1 j d × d , ( j = 3 , , 6 ) , C 2 j d × d , ( j = 3 , , 6 ) , C i j d × d , and ( i = 3 , , 6 , j = 1 , , 6 ) are zero matrices,
C 11 d × d = β 1 ( p i n + α ) a δ 1 T 1 d × d , C 12 d × d = γ 12 β 2 p o u t b δ 2 T 2 d × d , C 21 d × d = γ 21 β 1 p o u t a δ 1 T 1 d × d , C 22 d × d = β 2 ( p i n + α ) b δ 2 T 2 d × d
where T i d × d = p i ( 1 ) d p i ( d ) d p i ( 1 ) d 2 p i ( d ) .
In the following, we determine the matrix C of the spectral radius:
C = C 11 d × d C 12 d × d C 21 d × d C 22 d × d .
It is found that the matrix C has 2 d 2 eigenvalues equal to 0 by applying similar transformation to the matrix C . The remaining two eigenvalues satisfy
β 1 ( p i n + α ) a k 2 1 δ 1 x γ 12 β 2 p o u t b k 2 2 δ 2 γ 21 β 1 p o u t a k 2 1 δ 1 β 2 ( p i n + α ) b k 2 2 δ 2 x = 0 .
From the determinant (12), the characteristic equation can be obtained as
β 1 ( p i n + α ) a k 2 1 δ 1 x β 2 ( p i n + α ) b k 2 2 δ 2 x γ 12 β 2 p o u t b k 2 2 δ 2 γ 21 β 1 p o u t a k 2 1 δ 1 = 0 .
Therefore, the discriminant of the roots of quadratic Equation (13) is
Δ = β 1 ( p i n + α ) a k 2 1 δ 1 + β 2 ( p i n + α ) b k 2 2 δ 2 2   4 β 1 ( p i n + α ) a k 2 1 δ 1 β 2 ( p i n + α ) b k 2 2 δ 2 γ 12 β 2 p o u t b k 2 2 δ 2 γ 21 β 1 p o u t a k 2 1 δ 1 . = β 1 ( p i n + α ) a k 2 1 δ 1 β 2 ( p i n + α ) b k 2 2 δ 2 2 + 4 γ 12 β 2 p o u t b k 2 2 δ 2 γ 21 β 1 p o u t a k 2 1 δ 1 > 0
The roots of the quadratic equation; that is, the eigenvalues of the matrix C , are
x 1 , 2 = β 1 ( p i n + α ) a k 2 1 δ 1 + β 2 ( p i n + α ) b k 2 2 δ 2 ± Δ 2 .
Thus, the basic reproduction number of the whole network is expressed as [34]
R 0 = 1 2 β 1 ( p i n + α ) a k 2 1 δ 1 + β 2 ( p i n + α ) b k 2 2 δ 2 + 1 2 β 1 ( p i n + α ) a k 2 1 δ 1 β 2 ( p i n + α ) b k 2 2 δ 2 2 + 4 γ 12 β 2 p o u t b k 2 2 δ 2 γ 21 β 1 p o u t a k 2 1 δ 1 .
Equation (16) shows that the probability of being connected in the inner and outer communities plays a positive role in the basic reproduction number. According to the above analysis, the following theorem is obtained.
Theorem 1:
For system (5), if the basic reproduction number   R 0 < 1 , the disease-free equilibrium is stable and the disease will die out in two communities. If the basic reproduction number R 0 > 1 , the disease-free equilibrium is unstable and the epidemic can occur in two communities.

4. The Impact of Community Structure on the Spreading of Infectious Diseases

In the following, we study the influence of community structure on the process of spreading infectious diseases. Here, we focus on the variation in the basic reproduction number and infection density.

4.1. The Influence of the Connection Rate on the Basic Reproduction Number

Firstly, we study the impact of the connection rate within and between the communities on the basic reproduction number. The parameters are taken as k 1 = 9.08 , k 2 = 8.83 , δ 1 = 0.05 , δ 2 = 0.1 , γ 12 = 0.5 , and γ 21 = 0.35 .
Fix parameter p o u t = 0.2 ; Figure 5 shows the spatial-temporal distribution of the basic reproduction number in the network with the infection and connection rates within the community. Let parameter p i n = 0.2 ; Figure 6 shows the spatial-temporal distribution of the basic reproduction number in the network with the infection rate within the community and the external connection rate. These figures represent that the internal connection rate and the external connection rate are positively correlated with the basic reproduction number. Comparing Figure 5 and Figure 6, it is observed that the external connection rate has a greater impact on the basic reproduction number for the network.

4.2. The Influence of the Connection Rate and Connected Edges on Infection Density

In this section, we discuss the impact of the connection rate and the number of connected edges within and outside the communities on the transmission of infectious diseases. The parameter takes the value n = 100 , and the initial infected node I ( 0 ) = 1 .
Figure 7 shows the variation in infection density with the internal connection rate in the community when p o u t = 0.2 . One can see that the greater the internal connection rate in the community, the more rapidly the density of infected individuals grows. The reason is that individual-to-individual contact becomes frequent within the community, and if an infected individual appears at this time, it is more likely for the infected individual to spread the virus to susceptible individuals, which leads to an increase in the size of the infection in the community as a whole. It is also found that p i n can affect the peak time in single communities.
Figure 8 describes the variation in the infection density with the external connection rate when p i n = 0.5 . It can be observed that if the external connection rate p o u t increases, it means that intercommunity contact increases, which then also increases the rate and scale of the spreading of infectious diseases throughout the network. In reality, the multiple cross-provincial movements of individuals lead to frequent external connections, with the result that infectious diseases spread on a large scale across the country. Thus, it seems that reducing the number of trips is an effective way to control the spreading of infectious diseases.
Observing the red solid line and the blue dashed part in Figure 7 and Figure 8, it is known that when p o u t increases by 0.05, the density of infection can increase by about 0.03. When p i n increases by 0.15, the density of infection only increases by less than 0.02. Based on the above data, one obtains that intercommunity connection plays a crucial role in the spread of infectious diseases.
In what follows, we investigate the effect of the number of connected edges within communities and between communities on the transmission of infectious diseases, respectively.
Figure 9 shows the variation in infection density with the number of connected edges within the community for the case of R 0 < 1 . It can be observed when the number of new connected edges m 3 of nodes within the community increases, the speed of infection and the size of the outbreak increase. Thus, it can be concluded that the more susceptible people an infected person in the community comes in contact with, the wider the spread of the infection within the network, and ultimately the greater the number of infected people.
Figure 10 indicates the variation in infection density with the number of connected edges between communities for the case of R 0 > 1 . From this figure, one can observe that the higher the number of connected edges between communities, the faster the speed of the infection. When infected individuals move between communities, they come into contact with many individuals from various other communities. The result shows that the speed of the spreading of infected individuals, as well as the final size of the infectious disease, also increases.
The above analysis implies that reducing the movement of individuals across provinces can control large-scale outbreaks of infectious diseases in certain conditions.
In addition, we consider the number of nodes within the communities is n 1 = n 2 = 100 . The number of new connected edges added by individuals at each time step within the two communities is m 3 1 = 8 , m 3 2 = 5 . Figure 11 gives the curve of the density of each state of the infectious disease with time for different community structures.
Comparing the infection density curves in Figure 11, the basic reproduction number R 0 > 1 , the infection density of the network with more connected nodes within the community, will continue to increase until it reaches a steady state. This implies that the virus is not spreading massively, so the curve of the infection density in the graph gradually decreases to a steady state.
As we know, intercommunity activities play an important role in the spreading of infectious diseases. Here we combine Equations (4) and (16) to investigate the effect of community strength on the spread of infectious diseases. Figure 12 shows the variation curve of community strength with the basic reproduction number, where α = 0.1 and p i n = 0.2 .
From Figure 12, we can find that the basic reproduction number decreases as the modularity coefficient Q increases. This implies that the stronger the community structure, the less widespread the spread of infectious diseases will be. In fact, this conclusion holds, but only if the internal connection rate p-in and the number of inner edges are constant and the external connection rate p-out and the number of intracommunity edges are decreasing. Therefore, the possibility of a virus outbreak at the social level as a whole is small. When the community structure is weak enough, the basic regeneration number is greater than 1, and eventually the epidemic evolves into an endemic disease.
From the above analysis, we can summarize that the strength of the community can affect the final form of the infectious disease, i.e., whether it is eradicated or forms an endemic disease.

5. Conclusions

This paper presents a modified community network model considering the connection rate of nodes and the number of connected edges between communities. Based on this model, an SIR infectious disease transmission model with two community structures is constructed, which is closer to real networks. Furthermore, the transmission dynamics of epidemics are analyzed via theoretical analysis and numerical simulation. Our results show that the higher the frequency of susceptible nodes being infected due to the frequent contact of individuals within the community, the easier the infectious disease will spread within the community. When the external connection rate of the community increases, the movement of people between communities becomes frequent, which leads to the transmission of the infection throughout the network. In addition, when the external connected edges of communities increase, the infectious disease spreads rapidly between communities, and the number of infected individuals in the entire network increases massively.
In conclusion, controlling unnecessary movement of people and minimizing visits to crowded spaces are effective measures to curb the transmission of infectious diseases. This is also in line with our current epidemic prevention policy of not going out of town unless necessary and reducing the number of trips. In the future, we intend to investigate mathematical properties of the epidemic model and further explore the impact of data-driven models for multiplex networks.

Author Contributions

Z.G. (Zhongshe Gao): Revise and writing. Z.G. (Ziyu Gu): Software, Writing—original draft. L.Y.: Visualization, Revising. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (No. 11702195).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Gong, Y.W.; Song, Y.R.; Jiang, G.P. Global dynamics of a novel multi-group model for computer worms. Chin. Phys. B 2013, 22, 040204. [Google Scholar] [CrossRef]
  2. Pastor-Satorras, R.; Castellano, C.; Van Mieghem, P.; Vespignani, A. Epidemic processes in complex networks. Rev. Mod. Phys. 2015, 87, 925. [Google Scholar] [CrossRef]
  3. Song, H.T.; Wang, R.F.; Liu, S.Q.; Jin, Z.; He, D. Global stability and optimal control for a COVID-19 model with vaccination and isolation delays. Results Phys. 2022, 42, 106011. [Google Scholar] [CrossRef] [PubMed]
  4. Borge-Hohhoefer, J.; Meloni, S.; Goncalves, B.; Moreno, Y. Emergence of influential spreaders in modified rumor models. J. Stat. Phys. 2013, 151, 383–393. [Google Scholar] [CrossRef]
  5. Zhao, L.J.; Wang, Q.; Chen, J.J.; Chen, Y.; Wang, J.; Huang, W. Rumor spreading model with consideration of forgetting mechanism: A case of online blogging Live Journal. Physica A 2011, 390, 2619–2625. [Google Scholar] [CrossRef]
  6. Carvalho, S.A. Modeling and analysis of dengue epidemic: Control methods and vaccination strategies. Theory Biosci. 2019, 138, 223–239. [Google Scholar] [CrossRef]
  7. Pastor-Satorras, R.; Vespignani, A. Epidemic spreading in scale-free networks. Phys. Rev. Lett. 2000, 86, 3200–3203. [Google Scholar] [CrossRef]
  8. Pastor-Satorras, R.; Vespignani, A. Epidemic spreading in finite size scale-free networks. Phys. Rev. E 2002, 65, 035108. [Google Scholar] [CrossRef]
  9. Dickison, M.; Havlin, S.; Stanley, H.E. Epidemics on interconnected networks. Phys. Rev. E 2012, 85, 066109. [Google Scholar] [CrossRef]
  10. Zhang, M.Y.; Huang, T.; Guo, Z.X.; He, Z.G. Complex-network-based traffic network analysis and dynamics: A comprehensive review. Phys. A Stat. Mech. Its Appl. 2022, 607, 128063. [Google Scholar] [CrossRef]
  11. Girvan, M.; Newman, M.E.J. Community structure in social and biological networks. Proc. Natl. Acad. Sci. USA 2002, 99, 7821–7826. [Google Scholar] [CrossRef]
  12. Huang, W.; Li, C.G. Epidemic spreading in scale-free networks with community structure. J. Stat. Mech. Theory Exp. 2007, 2007, P01014. [Google Scholar] [CrossRef]
  13. Zhang, J.P.; Jin, Z. Epidemic spreading on complex networks with community structure. Appl. Math. Comput. 2012, 219, 2829–2838. [Google Scholar] [CrossRef]
  14. Chen, J.; Hu, M.B.; Li, M. Traffic-driven epidemic spreading in multiplex networks. Phys. Rev. E 2020, 101, 012301. [Google Scholar] [CrossRef]
  15. Li, J.X.; Wang, J.; Jin, Z. SIR dynamics in random networks with communities. J. Math. Biol. 2018, 77, 1117–1151. [Google Scholar] [CrossRef]
  16. Wang, S.F.; Gong, M.G.; Liu, W.F.; Wu, Y. Preventing epidemic spreading in networks by community detection and memetic algorithm. Appl. Soft Comput. J. 2020, 89, 106–118. [Google Scholar] [CrossRef]
  17. Arenas, A.; Cota, W.; Gomez-Gardenes, J.; Gomez, S.; Granell, C.; Matamalas, J.T.; Soriano-Panos, D.; Steinegger, B. Derivation of the effective reproduction number R for COVID-19 in relation to mobility restrictions and confinement. medRxiv 2020, 11, 2020-04. [Google Scholar]
  18. Matamalas, J.T.; Arenas, A.; Gomez, S. Effective approach to epidemic containment using link equations in complex networks. Sci. Adv. 2018, 4, eaau4212. [Google Scholar] [CrossRef]
  19. Zhua, P.C.; Wang, X.; Zhi, Q.; Ma, J.Z.; Guo, Y.M. Analysis of epidemic spreading process in multi-communities. Chaos Solitons Fractals 2018, 109, 231–237. [Google Scholar] [CrossRef]
  20. Fu, X.H.; Song, Q.L.; Li, S.Q.; Shen, Y.; Yue, S. Dynamic changes in bacterial community structure are associated with distinct priming effect patterns. Soil Biol. Biochem. 2022, 169, 108671. [Google Scholar] [CrossRef]
  21. Shao, F.; Jiang, G.P. Traffic driven epidemic spreading in homogeneous networks with community structure. J. Netw. 2012, 7, 850–855. [Google Scholar] [CrossRef]
  22. Gong, Y.W.; Song, Y.R.; Jiang, G.P. Time-varying human mobility patterns with metapopulation epidemic dynamics. Physica 2013, 392, 4242–4251. [Google Scholar] [CrossRef] [PubMed]
  23. Saramfiki, J.; Kaski, K. Modelling development of epidemics with dynamic small world networks. J. Theor. Biol. 2005, 234, 413–421. [Google Scholar] [CrossRef] [PubMed]
  24. Wei, L.X.; Zhang, J.G.; An, X.L.; Nan, M.; Qiao, S. Stability and Hopf bifurcation analysis of flux neuron model with double time delays. J. Appl. Math. Comput. 2022, 68, 4017–4050. [Google Scholar] [CrossRef]
  25. Torri, G.; Giacometti, R. Financial contagion in banking networks with community structure. Commun. Nonlinear Sci. Numer. Simul. 2022, 117, 106924. [Google Scholar] [CrossRef]
  26. Liu, Z.H.; Hu, B.B. Epidemic spreading in community networks. Europhys. Lett. 2005, 72, 315–321. [Google Scholar] [CrossRef]
  27. Newman, M.E.J.; Girvan, M. Finding and evaluating community structure in networks. Phys. Rev. E 2004, 69, 026113. [Google Scholar] [CrossRef]
  28. Salathe, M.; Jones, J.H. Dynamics and control of diseases in networks with community structure. PLoS Comput. Biol. 2010, 6, 1000736. [Google Scholar] [CrossRef]
  29. Stegehuis, C. Epidemic spreading on complex networks with community structures. Sci. Rep. 2016, 6, 29748. [Google Scholar] [CrossRef]
  30. Young, J.-G.; Hébert-Dufresne, L.; Allard, A.; Dubé, L.J. Growing networks of overlapping communities with internal structure. Phys. Rev. E 2016, 94, 022317. [Google Scholar] [CrossRef]
  31. Li, C.C.; Jiang, G.P.; Song, Y.R.; Xia, L.L.; Li, Y.W.; Song, B. Modeling and analysis of epidemic spreading on community networks with heterogeneity. J. Parallel Distrib. Comput 2018, 119, 136–145. [Google Scholar] [CrossRef]
  32. Xu, Z.P.; Li, K.Z.; Sun, M.F.; Fu, X.C. Interaction between epidemic spread and collective behavior in scale-free networks with community structure. J. Theor. Biol. 2019, 462, 122–133. [Google Scholar] [CrossRef]
  33. Dansu, E.J.; Seno, H. A model for epidemic dynamic in a community with visitor subpopulation. J. Theor. Biol. 2019, 478, 115–127. [Google Scholar] [CrossRef]
  34. Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef]
Figure 1. The relationships between Q and p i n , p o u t .
Figure 1. The relationships between Q and p i n , p o u t .
Entropy 25 00849 g001
Figure 2. (a,b) The relationship between Q and the internal connection rates.
Figure 2. (a,b) The relationship between Q and the internal connection rates.
Entropy 25 00849 g002
Figure 3. The relationship between the modularity coefficient Q and the connected edges.
Figure 3. The relationship between the modularity coefficient Q and the connected edges.
Entropy 25 00849 g003
Figure 4. The flow chart of the infectious disease transmission in two community structure networks.
Figure 4. The flow chart of the infectious disease transmission in two community structure networks.
Entropy 25 00849 g004
Figure 5. (a,b) The spatial-temporal distribution of R 0 with β and p i n .
Figure 5. (a,b) The spatial-temporal distribution of R 0 with β and p i n .
Entropy 25 00849 g005
Figure 6. (a,b) The spatial-temporal distribution of R 0 with β and p o u t .
Figure 6. (a,b) The spatial-temporal distribution of R 0 with β and p o u t .
Entropy 25 00849 g006
Figure 7. The variation curves of infection density with time for different p i n , p o u t = 0.2 .
Figure 7. The variation curves of infection density with time for different p i n , p o u t = 0.2 .
Entropy 25 00849 g007
Figure 8. The evolutions of infection density with time for different p o u t , p i n = 0.5 .
Figure 8. The evolutions of infection density with time for different p o u t , p i n = 0.5 .
Entropy 25 00849 g008
Figure 9. The variation curves of infection density with time for different m 3 and R 0 < 1 .
Figure 9. The variation curves of infection density with time for different m 3 and R 0 < 1 .
Entropy 25 00849 g009
Figure 10. The variation curves of infection density with time for different m 4 and R 0 > 1 .
Figure 10. The variation curves of infection density with time for different m 4 and R 0 > 1 .
Entropy 25 00849 g010
Figure 11. The variation in S ( t ) , I ( t ) , and R ( t ) with time.
Figure 11. The variation in S ( t ) , I ( t ) , and R ( t ) with time.
Entropy 25 00849 g011
Figure 12. The variation in R 0 with Q .
Figure 12. The variation in R 0 with Q .
Entropy 25 00849 g012
Table 1. The explanation of all the system parameters.
Table 1. The explanation of all the system parameters.
ParameterExplanation
β i Infection rate in the i - th community
α Probability of nodes being connected in the inner community at the initial moment
p i n Probability of nodes being connected in the inner community at each time step
p o u t Probability of nodes being connected in the outer community at each time step
δ i Recovery rate of infected in the i - th community
μ i Probability of recovered reverting to susceptible
γ i j Coefficient factor affecting the infection rate in the outer community
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

Gao, Z.; Gu, Z.; Yang, L. Effects of Community Connectivity on the Spreading Process of Epidemics. Entropy 2023, 25, 849. https://doi.org/10.3390/e25060849

AMA Style

Gao Z, Gu Z, Yang L. Effects of Community Connectivity on the Spreading Process of Epidemics. Entropy. 2023; 25(6):849. https://doi.org/10.3390/e25060849

Chicago/Turabian Style

Gao, Zhongshe, Ziyu Gu, and Lixin Yang. 2023. "Effects of Community Connectivity on the Spreading Process of Epidemics" Entropy 25, no. 6: 849. https://doi.org/10.3390/e25060849

APA Style

Gao, Z., Gu, Z., & Yang, L. (2023). Effects of Community Connectivity on the Spreading Process of Epidemics. Entropy, 25(6), 849. https://doi.org/10.3390/e25060849

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