Abstract

Population demography can change the network structure, which further plays an important role in the spreading of infectious disease. In this paper, we study the epidemic dynamics in temporal clustered networks where the local-world structure and clustering are incorporated into the attachment mechanism of new nodes. It is found that increasing the local-world size of new nodes has little influence on the clustering coefficient but increases the degree heterogeneity of networks. Besides, when the network evolves faster, increasing the local-world size of new nodes leads to a faster initial growth rate and a larger steady density of infectious nodes, while it has small impacts on the steady density of infectious disease when the network evolves slowly. Furthermore, if the average degree is fixed, increasing the probability of triad formation enlarges the clustering coefficient of a network, which reduces the initial growth rate and steady density of infectious nodes in the network. This work could provide a theoretical foundation for the control of infectious disease.

1. Introduction

Infectious diseases, such as H1N1 influenza and the ongoing coronavirus disease 2019 (COVID-19), pose a significant health and economic burden on society. Mathematical models have been used as a powerful tool to mitigate the impacts of infectious diseases as they can provide valuable insights into the transmission mechanism of diseases and then guide policymakers implementing effective prevention and control measures. More recently, researchers have devoted a great deal of effort to mathematically model the spread of COVID-19 and seek for effective nonpharmaceutical interventions, such as social distancing, isolation, and intercity travel restrictions, to control the local outbreak of COVID-19 [15]. Besides, Li et al. used the signed networks to study the influence of structural balance on epidemic transmission [6]. Sun et al. investigated the effects of resource diffusion on epidemic propagation in multilayer networks with simplicial complexes, showing that increasing the resource diffusion on 2-simplexes can suppress the epidemic spreading [7]. It is also discovered that information transmission induced by the 2-simplex can suppress epidemic outbreak to a certain extent [8]. Also, a review of different phenomena on higher order networks was given by Majhi et al. [9]. All these studies have demonstrated that the population contact pattern which can be described as a network plays an important role in the epidemic spreading.

The population demography can change the network structure which further affects the disease spreading process on networks. Kamp introduced an SID model with HIV as a case study to investigate the interplay between epidemic spreading and the dynamics of network structure [10]. Jin et al. established an SIS heterogeneous mean-field model with a constant recruitment rate and derived the global dynamics of the model where the new members link to existing nodes randomly [11]. Piccardi et al. integrated the birth and death process into the SIRS model on heterogeneous networks and disentangled the underlying causes of the different outcomes obtained for SIR and SIRS processes by applying an age-degree model [12]. Based on conditional Markov chain and pairwise approximation, Luo et al. studied an SIS pairwise model on networks, finding that demographics of individuals can induce the extinction of an epidemic [13]. Jing et al. studied the effects of different moment closure methods [14], adaptive behaviors of new members [15], and the distribution of infection age [16] on epidemic spreading in networks with demographics. Leung et al. provided a new framework using digital proxies of population mobility and mixing to monitor viral transmissibility and effectiveness of social distancing interventions in the ongoing COVID-19 pandemic [17].

Currently, most research in network epidemic models with demographics has assumed that a new node links to existing nodes according to random or preferential attachment mechanism when it enters the network. In fact, there has been evidence to indicate that the preferential attachment mechanism works on a local world of each node [18, 19]. Using a local preferential attachment mechanism, Li and Chen constructed a local-world evolving network whose distribution is a transition between that of an exponential network and of a power-law scaling network [18]. Zhang et al. proposed a local-world evolving network model with changeable local-world size and tunable clustering by including a triad formation step in the local-world evolving network model [20]. Wen and Duan introduced extended links to mimic the weak relations between nodes in different local worlds, obtaining a local-world evolving model which exhibits scale-free behavior as well as the small-world property [21]. Wang et al. presented a local preferential attachment model where a local area consists of a node and all its neighbors [22]. They derived the stable degree distributions and clustering-degree correlations of the network. Besides, the epidemic spreading behavior [23], rumor spreading process [24], and cascading failures [25] in local-world evolving networks are studied. It shows that the attachment mechanism of new members has a significant impact on the network structure and the propagation process on networks.

Many real-world networks exhibit nonnegligible clustering which can have an important impact on epidemic dynamics. However, there is no clear consensus on the role of clustering. For example, it shows that the effect of clustering on epidemic spreading in Watts–Strogatz networks is stronger than the effect of clustering on epidemic spreading in Facebook network [26]. Considering the local-world structure and clustering characteristics, we will investigate the effects of attachment mechanism of new members on the network structure and epidemic dynamics in this paper.

The rest of this paper is organized as follows. In Section 2, we introduce some basic characteristics of networks. Then, the evolving network model and epidemic spreading process on networks are given in Section 3. In Section 4, numerical simulations are conducted to explore the impacts of mechanisms that new nodes connect existing nodes on the network structure and epidemic dynamics. Finally, some conclusions and discussions about this work are given in Section 5.

2. Basic Characteristics of Networks

Let be an unweighted, undirected network with representing the set of nodes and representing the set of links between nodes in the network. For a finite network with nodes, the set can be encoded in the adjacency matrix whose entry is 1 if there is a link between and , and 0 otherwise. A finite network is determined uniquely by its adjacency matrix.

Degree distribution is the most important characteristic of the network. The degree of one node is the number of links connecting it. Then, the degree distribution is given by , where is the number of nodes with degree . Some widely used networks include regular random, Poisson, bimodal, and scale-free networks.

Clustering, one frequently used statistical property of characterizing networks, measures the probability that two neighbors of a randomly chosen node share a link to form a triangle [27]. Specifically, the clustering coefficient of node is defined as the ratio of the number of triangles and triples:where and are the total number of existing links and the number of all possible links between the neighbors of node , respectively. Furthermore, the degree-dependent clustering coefficient of node with degree can be defined as the average clustering coefficient of nodes with degree , i.e.,where is the set of nodes with degree . Then, the clustering coefficient of the network is defined as the weighted average of degree-dependent clustering coefficient, which is also the average of clustering coefficient of nodes:where is the maximum degree of nodes.

3. The Model

3.1. The Evolving Network Model

The local-world evolving network was first proposed by Li and Chen [18]. It captures the localization of real-life networks. Then, Zhang et al. presented expanded local-world evolving network with tunable clustering by including a triad formation (TF) step [20]. Considering the local-world structure and clustering of networks, we propose the generating algorithm of the dynamic clustered network with demographics. It is summarized as follows:(i)Initial condition: the initial network consists of nodes and has a Poisson degree distribution with average degree .(ii)Network evolving: At each time step , we fist add a new node to the network and then delete a randomly selected node as well as its links from the network. The attachment mechanism between the new node and existing nodes is as follows:(1)Local-world establishment: randomly select nodes from the network, referred to as the “local world” of .(2)Local preferential attachment (LPA): The new node connects to different nodes in its local world according to the preferential attachment. Thus, the probability that node is connected to an existing node in the local world of node depends on the degree of node :(3)Triad formation (TF): When the new node connects an existing node , we then implement a TF step with probability . Specifically, if a link between and was added in the previous LPA step, then one more link between and a randomly chosen neighbor of is added. If all neighbors of have already been connected to , do one more LPA step instead.

A sketch of the network evolving process is shown in Figure 1. After steps (the designed total iteration steps which are large enough to reach a stationary network), the algorithm results in a connected network with nodes and average degree is . Obviously, the network size keeps a constant and it has and . There are two limiting cases for the evolving network: and . When , the stationary degree distribution is Poisson with mean , while when , the stationary degree distribution is “stretched exponential,” namely, [28]. Besides, the parameter allows us to include the clustering effect into network by building triads.

3.2. Epidemic Spreading on Evolving Networks

We consider an SIS epidemic spreading on evolving networks, where the initial network is the stationary network obtained by the generating algorithm in Section 3.1. During the epidemic, nodes can either be susceptible (S) or infectious (I). The susceptible nodes can be infected by their infectious neighbors with a per-contact transmission rate if they are connected to one or more infected nodes, while the infectious one can be cured at a recovery rate γ. When an infectious node recovers, it turns to susceptible and can be infected again. Meanwhile, existing nodes together with all their links can leave the network due to death, and new nodes enter the network due to birth. Here, the mechanism that nodes enter or leave the network is identical to that in Section 3.1. Thus, the population size is a constant. Given a short time interval , the algorithm of epidemic spreading on evolving networks can be summarized as follows:(i)Initialization: An undirected and unweighted network with nodes is created and set to the initial network by performing the network generating algorithm in Section 3.1. Then, randomly selected nodes in the network are initially infected and the other nodes are susceptible.(ii)Epidemic process: Given the state of each node at time , the states of nodes at time are updated as follows. Each infected node recovers and becomes susceptible again with probability . Each susceptible node, such as node , turns into infectious with probability , where is the number of infectious neighbors of .(iii)Network evolving: At time , (extracted from the binomial distribution ) randomly selected nodes as well as all their links are eliminated from the network. Then, nodes are added into the network sequentially where links are created according to the mechanism in Section 3.1.

To sum up, we can simulate the epidemic spreading on dynamic clustered networks by iterating the epidemic process and networking evolving process until there is no infectious node in the network or the number of infectious nodes reaches a steady state. During the iteration, record the number of infectious nodes in each time step.

3.3. The Heterogeneous Pairwise Model on Evolving Networks

In line with the above algorithms, we construct the heterogeneous pairwise model on evolving networks in this section. The notations of the model follow those established in [29], i.e., using square brackets to denote the expected numbers of nodes, pairs, or triples of a particular type. Specifically, let be the number of nodes in state and with degree , and be the number of nodes in state . represents the number of pairs with one node in state having degree and with the other node in state having degree . stands for the number of triples of type . If a disease state of a variable has no subscript, it implicitly contains the sum over all possible degrees, such as and . Let denote the probability that a link of new nodes is connected to an existing node in state with degree . It depends on the values of and . Furthermore, according to the attachment mechanism, the degree distribution of new nodes is

Using the notations above, the SIS pairwise model on evolving networks can be written as

In order to close the above equations, we need to give the expression of and an approximation of triples in terms of nodes and pairs. The expression of depends on the values of and , which are classified as the following three cases:(1), : In this special case, newcomers connect to all nodes in their local world. It means that the LPA is reduced to uniform attachment, that is, new nodes connect to existing nodes randomly. Therefore, it hasMoreover, the network clustering coefficient approached 0 such that the triples can be approximated by the triple closure in [30]:Substituting equations (7) and (8) into equation (6) yields a closed system.(2), : In this limiting case, the network clustering coefficient approached 0 and new nodes connect existing nodes exactly the same as global preferential attachment mechanism. Therefore, it hasSubstituting equations (8) and (9) into equation (6) gives a closed system.(3): In this general case, the explicit expressions of and are much more complicated and unable to be derived since they depend on the values of and . Thus, we do not give the exact formula of the model here. There are two remarks about the model.Remark 1. It is assumed in [18] that , finding that the local-world growing network also has a power-law degree distribution when . Similarly, the formulas of and are exactly the same as in equation (9) under this assumption. This reveals the main reason that as increases, the epidemic curves are closer to those when as shown in Figure 2.Remark 2. When , there is small but nonnegligible clustering in the network. Let be the network clustering coefficient; then, the triples can be approximated by the triple closure with clustering in [31]:

4. Numerical Simulations

4.1. Degree Distribution and Clustering Coefficient of Networks

We will investigate the degree distribution and clustering coefficient of networks by stochastic simulations. In these simulations, the initial network size is and each pair of nodes has a link independently of any others with probability . Thus, the initial network has a Poisson degree distribution with mean . The simulation lasts 20000 time steps according to the network generating algorithm described in Section 3.1.

Figure 3 gives the stationary degree distribution of networks under different values of and . It demonstrates that the local-world size and the probability of triad formation have a direct impact on the network structure. It can be seen from Figure 3(a) that when , the larger the value of , the flatter and wider the degree distribution. This means that increasing the local-world size reinforces the degree heterogeneity for networks. However, as both the values of and increase, this impact disappears as shown in Figure 3(c). If the values of and are large enough, new nodes link to existing nodes with approximately a global preferential attachment mechanism, that is, the local preferential attachment selection is not effective. Thus, there is relatively little change in the degree distribution of networks.

Figure 4 shows the evolution of clustering coefficients for different values of and . It suggests that the local-world size has little effect on the network clustering for a fixed value of . However, for fixed , the network clustering coefficient is approximately 0 if , whereas it converges rapidly to a stable value if . Furthermore, compared with the growing network without deletion of nodes [20], it can be found that node removal dramatically decreases the degree distribution and clustering in networks. But the network still has nonnegligible clustering where the clustering coefficient is increasing with in our model. Thus, the clustering coefficient can be tuned by changing the probability of triad formation.

4.2. Epidemic Dynamics

Figure 5 provides a comparison of stochastic simulations to numerical simulations of pairwise model for two limiting cases and , which correspond to uniform attachment and global preferential attachment, respectively. One can see that the agreement is very good in both cases. Besides, when new nodes connect existing nodes according to global preferential attachment, the epidemic has a faster initial growth rate and a larger steady density of infectious nodes compared to the case of uniform attachment. This phenomenon is more significant in Figure 5(a) where the average degree of the network is smaller. It indicates that the attachment of new nodes can significantly affect disease dynamics, and decreasing the degree heterogeneity of networks is helpful in mitigating epidemics.

Figure 2 investigates the effects of local-world size on epidemic dynamics. In Figure 2(a) where is large, higher values of correspond to faster initial growth rate and larger steady density of infectious nodes. However, when is small as shown in Figure 2(b), the value of has only a relatively small effect on the initial growth rate and does not affect the steady density of infectious nodes. Besides, it can be seen that a larger value of leads to a higher steady density of infectious nodes, and the epidemic reaches the steady state quicker by comparing Figures 2(a) and 2(b). This is due to the fact that the network evolves faster such that has a more significant effect on the network distribution as increases. Specifically, infectious nodes are removed faster from the network and susceptible nodes enter the network instead. Meanwhile, more infectious nodes leave the network before they infect others, reducing the average infection period. All these factors contribute to the decline in epidemic spreading. In contrast, when is small, the network evolves slowly and the network structure is relatively stable, such that the local-world size has little influence on the steady density of infectious diseases.

Figure 6 gives the simulation results based on the average of 200 stochastic simulations for different values of and . For comparison, and are determined such that the networks have the same average degree, i.e., . It shows that the initial growth rate and steady density of infectious nodes decrease with increasing values for fixed and average degree. On the other hand, combining with the result noted in Figure 4 that larger values of lead to a larger clustering coefficient of the network, it implies that clustering in networks is capable of suppressing infection of an epidemic. This is consistent with the previous work in static clustered networks [29, 32].

5. Conclusion and Discussion

The present study has explored the interplay between epidemic spreading and network evolution. In our model, new nodes fist connect existing nodes according the local preferential attachment mechanism and then connect to one randomly chosen neighbor of these nodes with a given probability . Such network evolution mechanism incorporates both the local-world structure and clustering into networks. Simulation results illustrate that the attachment mechanism of new nodes can influence the degree distribution, clustering, and epidemic dynamics in networks. It shows that increasing the local-world size leads to a stronger heterogeneity in the network degree distribution but has little impact on the clustering in networks. Furthermore, it corresponds to a faster initial growth rate as well as a larger steady density of infectious nodes when the network evolves faster. In contrast, the local-world size has a small impact on the steady density of infectious nodes if the network evolves slowly. On the other hand, for a given average degree and local-world size, increasing the probability of triad formation leads to a larger clustering coefficient of the network which further suppresses the infection of an epidemic. In conclusion, controlling the attachment mechanism of new nodes could be an effective measure to mitigate an epidemic.

A limitation of our study is that the probability that an edge emanating from new nodes connects a susceptible or an infectious node in networks cannot be calculated explicitly. Besides, the high dimensionality of the model poses a challenge for theoretical analysis of epidemic dynamics. How to derive an explicit formula for the probability ( and ) and investigate the steady states as well as their stability of the model needs more rigorous mathematical techniques. This remains an open question.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (nos. 12101373 and 11971278), the Natural Science Foundation of Shanxi Province (no. 20210302124261), and the Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi (no. 2020L0246).