Inhomogeneity of epidemic spreading.

In this study, we use the characteristic infected cluster size to investigate the inhomogeneity of the epidemic spreading in static and dynamic complex networks. The simulation results show that the epidemic spreads inhomogeneously in both cases. Also, the inhomogeneity of the epidemic spreading becomes smaller with increasing speed of moving individuals and almost disappears when the speed is high enough. (c) 2010 American Institute of Physics.

The epidemic has always been a threat to the human health. The rule of the epidemic spreading as an important issue has been quite extensively studied for many years, but has not been mastered well. The epidemic spreads by the interaction between people, and complex networks, whose nodes represent individuals or organizations and link the interactions among them, can characterize this interpersonal interaction perfectly. In order to conveniently study the epidemic spreading, researchers usually adopt the homogeneous network structure and the homogeneous mixing hypotheses. Recently, researchers began to discuss the credibility of the homogeneous network structure hypothesis and consider the inhomogeneity of the network structure. However, the credibility of the homogeneous mixing hypothesis has not been well discussed until now. In this paper, we try to discuss this credibility and carry out numerical simulations on an epidemic spreading model. The results show that the epidemic spreading is inhomogeneous whether the individuals are stationary or moving. However, the inhomogeneity decreases with the increase of individuals’ speed, and almost disappears when the speed is high enough.

INTRODUCTION

In recent years, physicists and biologists have been attracted to the spreading of the epidemic. Undoubtedly, the successive large-scale outbreaks of SARS (Severe Acute Respiratory Syndrome) and H1N1 make them more enthusiastic about this study. Complex networks as a new branch of statistical physics provide a reliable model for the intensive study of the epidemic spreading. When studying the epidemic spreading in complex networks, researchers usually adopt the homogeneous network structure and homogeneous mixing hypotheses. As research progressed, researchers began to discuss the credibility of the homogeneous network structure hypothesis and take into account the inhomogeneity of the network structure. In this study, we try to check the credibility of the homogeneous mixing hypothesis, meaning that each infected individual has the same probability of contacting with any susceptible (healthy) individual, which has not been well discussed until now.

We will focus on the susceptible-infected-susceptible (SIS) model on complex networks. Here each node in the network represents an individual, and each link represents a connection along which the epidemic can spread. A susceptible individual can be infected by its infected neighbors with some probability, and the infected individual can become susceptible with other probability. As the infected individuals always infect their neighbors, we predict that the spreading of the epidemic may actually be inhomogeneous. That is, the probabilities of infected individuals connecting with susceptible ones differ significantly. To investigate the inhomogeneity, we perform large scale numerical simulations on static and dynamic networks. In these simulations, we study two spreading modes of the SIS model. In one mode, the probability that a susceptible individual is infected is unrelated to the number of its infected neighbors. In another mode, the probability increases with the number.

In this work, the inhomogeneity of the epidemic spreading is characterized by the characteristic infected cluster size (hereinafter referred to as “CICS”), where the cluster is the subnet whose nodes are connected (i.e., from any node, one can reach any other node along links in the subnet), the infected cluster only includes infected individuals, and the CICS is the typical size of the largest infected cluster, namely, the number of infected individuals of the largest infected cluster. The simulation results show that the infected individuals are always distributed inhomogeneously and prone to gather into large clusters, even if they walk randomly. More interestingly, the inhomogeneity of the epidemic spreading decreases with increasing speed of the individuals, and the epidemic nearly spreads homogeneously when the moving speed is high enough.

MODEL

In our model, individuals walk randomly in a square of linear size with periodic boundary conditions, and they are distributed randomly in the square initially. The and are the position and motion direction of the individual at time . Then, at time , they are updated according towhere is the velocity of the individual at time , is the speed which is the modulus of the velocity. The speed is same for all individuals and remains constant in motion. follows the uniform distribution in . is the update interval and is set to 1. If , the individual is called a neighbor of the at time , where is the interaction radius. This means that there is a link between the and nodes in the corresponding network.

We assume that the update of the individuals’ positions and states is simultaneous with the fixed time interval . Each individual has two states: susceptible and infected . The next state of each individual depends on its current state and its neighbors’ states. Concretely, if the individual is infected currently, then it is cured and becomes susceptible at the next time step with probability ; if it is susceptible currently, then it can be infected by its infected neighbors in either of two modes. In mode 1, the individual is infected with probability if it has infected neighbors. In mode 2, it is infected with probability . Apparently, each infected neighbor of the individual infects it independently.

Researchers have analyzed the epidemic thresholds of mode 1 in various complex networks and have shown that the degree distribution of the network has a significant impact on the epidemic spreading. Liu et al. investigated the epidemic spreading of mode 2 in community networks. The epidemic spreading of mode 2 in dynamic networks has been explored by Frasca et al. recently. In this work, the mode 1 and mode 2 are both studied. Without lack of generality, we set , as it only affects the time scale of the infection evolution. Considering the reality, we will not set to a large value. As a result, the infected individuals can gather into several clusters rather than one. Then, the homogeneity of the epidemic spreading can be explored by analyzing and comparing the CICSs of the different cases.

In the following, the number of individuals is fixed to 100. is denoted as the density of the network. The size of the square is measured in units of the interaction radius . As we can see, is critical for determining the network structure. If approaches to , most individuals are connected together. Then, too many individuals will be infected persistently. By contrast, if is much larger than , individuals have few opportunities to interact with each other. The network is divided into many small clusters, then not many individuals will be infected. As a result, the epidemic will disappear soon.

In this work, we only discuss the cases when is proper. When (static networks), the percentage of the infected individuals is plotted as a function of with different in Fig. 1. As shown in Fig. 1(a) [Fig. 1(b)], when in mode 1 (mode 2), only when is above 0.6 (0.45), the epidemic can spread in the network. When in mode 1 (mode 2), the epidemic persists as long as is above 0.15 (0.07); about one third (half) of the individuals are infected when equals to 0.5. Thus, we choose to ensure that neither the epidemic will disappear, nor too many individuals will be infected, with a large range of value. Moreover, Fig. 6(b) in Sec. IV indicates that the moving of individuals cannot affect the percentage of the infected individuals in mode 2 and has little influence on the percentage in mode 1. So, we also choose when (dynamic networks).

FIG. 1.

The percentage of infected individuals as a function of with (circles), (squares), and (diamonds) in mode 1 (a) and mode 2 (b). Results are averaged over 15 runs.

FIG. 6.

(a) When , vs , with the change of : mode 1 (circles) and mode 2 (squares). Also, vs (diamonds). (b) vs , for : mode 1 (the bottom curves) and mode 2 (the top curves). Results are averaged over 15 runs.

In static networks, the individuals are distributed randomly and keep still over time, so each individual connects to another with the same probability. Thus the degree distribution , as the probability that a node is linked to other nodes, of static networks is Poissonian. In dynamic networks, the individuals are randomly distributed initially and walk randomly in the square. Therefore, the probability that an individual connects to another is a constant, which equals to the probability of the static case. Then the degree distribution of dynamic networks is also Poissonian. In Sec. III, we investigate the epidemic spreading in static networks . Section IV is devoted to the study of the epidemic spreading in dynamic networks .

EPIDEMIC SPREADING IN STATIC NETWORKS

In order to characterize the homogeneity of the epidemic spreading, we introduce the “homogeneous mode.” In this mode, the infected individuals, whose number is , are randomly scattered in a square of linear size . Then their distribution is homogeneous. When the infection density is not quite large, the distribution of the infected cluster size in homogeneous mode is similar to the distribution of the cluster size of the site percolation. That is, there exists a . When , the number of infected clusters decays approximately as a power law in , while decaying much faster than a power law for (see Fig. 2). In reality, the infected clusters whose sizes are larger than barely exist, so the infected cluster with the critical infected cluster size is regarded as the largest infected cluster that can be observed, and is denoted as the CICS.

FIG. 2.

Log-log plot of the infected cluster number as a function of the infected cluster size in a square of size in homogeneous mode. For (a), 40 (b), 50 (c), and 60 (d), the infected cluster number decays as a power law for and rapidly for . Results are averaged over 500 runs.

With the increase of , a critical infection density emerges, above which it is possible that all the individuals gather into one cluster. Our experiments show that is about 1.4. As shown in Fig. 1, the maximum number of infected individuals is about 50, that is, the maximum infection density is about 1 . In such cases, infected individuals always gather into more than one cluster. Then we can characterize the homogeneity of their distribution with the CICS. For the convenience of discussion, we denote the CICS of homogeneous mode as .

Correspondingly, the number and density of the infected individuals and the CICS in mode 1 (mode 2) are denoted by , , and , respectively. We will see later that when , , , and differ significantly, that is, the epidemic spreads inhomogeneously in mode 1 and mode 2. At , a given number of individuals are taken as the seeds of the infection (the proportion is 10%), while all the others start from the susceptible state. Also, all individuals are scattered randomly in the square, namely, the infected and susceptible individuals mix well. After an initial transient process, the systems stabilize in a steady state with a constant average infection density.

For each value of , we can get the values of . With changing , is plotted as a function of [the middle curve in Fig. 3(a)]. In the same way, is plotted as a function of [the top curve in Fig. 3(a)]. In homogeneous mode, each (or each ) corresponds to one , then we can plot as a function of [the bottom curve in Fig. 3(a)]. For the sake of getting value corresponding to a given , we plot versus , in Fig. 3(b).

FIG. 3.

(a) With the change of , vs , : mode 1 (circles) and mode 2 (squares); vs (diamonds) for homogeneous mode. The dashed line is for . (b) vs , : mode 1 (circles) and mode 2 (squares). Results are averaged over 15 runs.

The initial infection density is 10%, so the epidemic will be persistent with time, only when the infection density is larger than 0.1 after a long spreading. Therefore, we make the epidemic spread in the square for a long time and then discuss the situation when below. As is shown in Fig. 3(a), we draw the following conclusions by contrasting , with . When , and are larger than significantly, which means that the inhomogeneity exists. Concretely, when is near 0.2, and are up to 2–3 times , which is a significant difference. With increasing , the multiple becomes less although the difference between and has little change. Then we can say that the smaller the infection density , the larger the inhomogeneity of the epidemic spreading.

In these cases, when , can reach two or more times and [see the line of Fig. 3(a)], that is, the number of infected individuals of homogeneous mode is two or more times that of mode 1 and mode 2. It illustrates that, from another point of view, the epidemic spreading in mode 1 and mode 2 is inhomogeneous. When is about 0.8, and are similar, that is, infected individuals are distributed homogeneously in mode 1 and mode 2. However, the values of mode 1 and mode 2 reach 0.5 and 0.8, respectively [see Fig. 3(b)]. They are so large that we do not discuss this case. In fact, the epidemic always spreads near the infected individuals, that is, the infected individuals only infect their susceptible neighbors, and thus it is a natural thing for epidemic to spread inhomogeneously when the infection density is not so large.

Figure 3(a) also shows that when and , . It means that the epidemic spreading in mode 2 is more inhomogeneous than in mode 1, namely, the infected individuals are much easier to gather into large clusters in mode 2 than in mode 1. In addition, when the values of are the same, , as can be seen from Fig. 3(b). It also demonstrates the same conclusion. The reason is that the probability that a susceptible individual is infected increases with the number of its infected neighbors in mode 2, i.e., the susceptible individuals are more easily infected in mode 2 than in mode 1.

EPIDEMIC SPREADING IN DYNAMIC NETWORKS

In this section, we discuss the inhomogeneity of the epidemic spreading while the individuals walk randomly in the square , and let . As has been argued, is set to 7, the initial proportions of the infected individuals in mode 1 and mode 2 are 10%. In Fig. 4, we plot the evolution of the infection density and in mode 1 and mode 2 in the case that , respectively. As shown in Fig. 4, after the initial transient process, fluctuates narrowly around a value, which means that the SIS model reaches the steady state. Figure 5 shows the evolution of and in the case that and .

FIG. 4.

(a) The evolution of the infection density in mode 1. (from top to bottom). Each straight line corresponds to the average infection density after the initial transient process. (b) Same with (a) but for mode 2. Results are averaged over 15 runs.

FIG. 5.

(a) The evolution of the CICS when in mode 1. (from top to bottom). The straight line corresponds to with in homogeneous mode. (b) Same with (a) but for mode 2. Results are averaged over 15 runs.

For each value of , we can plot as a function of by changing the value. Curves of versus , at and the curve of versus [same with that in Fig. 3(a)] are given in Fig. 6(a). It shows that when , the epidemic spreading is inhomogeneous; when is about 0.8, the infected individuals are distributed homogeneously. Clearly, this conclusion is consistent with that of static networks, which are drawn in Fig. 3(a). The correspondences between and for different are given in Fig. 6(b). Interestingly, the different (including ) corresponds to the similar for a given , as shown in Fig. 6(b). This means that the moving of individuals cannot affect the proportion of the infected individuals in mode 2. This conclusion can be obtained by letting the delay of the model in Ref. 18.

In order to further study the impact of speed on the inhomogeneity of the epidemic spreading, we will investigate the variation of with increasing of for different . As the epidemic spreading emerges the inhomogeneity when , are set to 0.2, 0.4, and 0.6 in the following discussion. At first, we plot as a function of by changing the value of for each . Then a set of functions of and is obtained by changing the value. Corresponding to , a set of can be gotten. Then is plotted as a function of when in the bottom subplot of Fig. 7(a). For the sake of comparison, we plot two straight lines in the same subplot, which respectively correspond to (the bottom dashed line) and (the top one) for in Fig. 3(a). Similarly, the curves and straight lines corresponding to are also plotted in Fig. 7(a) (the middle subplot for , the top subplot for ). In the same way, we plot as a function of for , respectively, and the corresponding straight lines in Fig. 7(b).

FIG. 7.

(a) In mode 1, the CICS as a function of speed for the infection density (from bottom to top). Each bottom (top) dashed line corresponds to for the same in Fig. 3(a). (b) Same as (a) but for mode 2. Results are averaged over 15 runs.

As can be seen in Fig. 7, when , is apparently larger than for each , which means that the epidemic spreading is inhomogeneous. Concretely, when , is close to that of the static network . With increasing of , decreases evidently and moves toward (homogeneous mode) progressively. This means that the inhomogeneity of the epidemic spreading becomes smaller with increasing speed . Generally speaking, it is easier to maintain a structure in the static environment. Interestingly, our simulation results indicate that when the individuals walk randomly and the speed is not very high, their distribution is inhomogeneous. That is, the inhomogeneity of the epidemic spreading is kept in the dynamic environment. When , is always near , which means that the distribution of the infected individuals is approaching to that in homogeneous mode. This is because when is large enough, any infected individual can easily jump out the area that is covered by its infected neighbors. Therefore, the large infected clusters cannot form.

As also shown in Fig. 7, the infection density affects the difference between the maximum and minimum of . When , the maximum of is 1.54 (2.26) times the minimum, while when , the maximum of is 1.16 (1.34) times the minimum. Then we can say that the smaller the infection density , the stronger affects the inhomogeneity of the epidemic spreading. Besides, comparing Fig. 7(a) and Fig. 7(b) represents that when and are the same, are always larger than . This means that the inhomogeneity is more obvious in mode 2 than in mode 1 at the same speed , which is in accord with the case of static networks.

CONCLUSIONS

In this paper, the inhomogeneity of the epidemic spreading in two spreading modes of the SIS model is investigated. The simulations in the static and dynamic networks show that the infected individuals are usually prone to gather into large clusters as the infected individuals always infect their neighbors. For such a reason, the epidemic usually spreads inhomogeneously. Even in dynamic networks, the inhomogeneity can be kept well. And, the smaller the infection density, the more inhomogeneously the epidemic spreads. However, the inhomogeneity decreases with the increase of the individuals’ speed in the dynamic networks, and the epidemic spreading becomes almost homogeneous when the speed is large enough.