Inhomogeneity of epidemic spreading with entropy-based infected clusters.

Considering the difference in the sizes of the infected clusters in the dynamic complex networks, the normalized entropy based on infected clusters (δ*) is proposed to characterize the inhomogeneity of epidemic spreading. δ* gives information on the variability of the infected clusters in the system. We investigate the variation in the inhomogeneity of the distribution of the epidemic with the absolute velocity v of moving agent, the infection density ρ, and the interaction radius r. By comparing δ* in the dynamic networks with δH* in homogeneous mode, the simulation experiments show that the inhomogeneity of epidemic spreading becomes smaller with the increase of v, ρ, r.

The epidemic spreading in human populations is a typical public health incident, and has also become one of the major public safety issues that humanity is facing until now. In order to conveniently study the epidemic spreading, researchers usually adopt the homogeneous mixing hypotheses. However, such a hypothesis is not strong. Recently, researchers have discussed the credibility of the hypothesis. In this paper, we further study this credibility and carry out simulation experiments in the dynamic network. The results indicate that the inhomogeneity of epidemic spreading is affected by the absolute velocity v, the infection density , and the interaction radius r.

INTRODUCTION

The studied models of the epidemic spreading have greatly improved the public's understanding of infection mechanisms. Among the most studied models for investigating the spread of epidemic, susceptible-infected-susceptible (SIS) and susceptible-infected-recovered (SIR) epidemic models are two famous epidemiological models. When studying the epidemic spreading with these models in complex networks, researchers often adopt homogeneous mixing hypotheses which each infected agent has same probability of contact with any susceptible agent. However, such a hypothesis is not strong and should be thoroughly discussed. Recently, researchers have discussed the credibility of the hypothesis. For instance, Liu et al. implement the distribution of an epidemic is always inhomogeneous after the transition process in Refs. 13 and 14, where they present the characteristic infected cluster size (CICS), which characterizes the inhomogeneity of the epidemic spreading. However, the CICS measurement studied on the typical size of the largest infected cluster, which cannot reflect the infected clusters.

Entropy is a fundamental concept in physics. It is related to the information content and order/disorder of a system. In statistical physics, cluster entropy has been used in the study of problems, such as percolation and complex systems. Cluster entropy was utilized within the Axelrod model to measures the number of cultural groups of different sizes. Cluster-size entropy was also used in the Axelrod's cultural adaptation model which gives information on the variability of the cultural cluster size present in the system. In the paper, the normalized entropy based on infected clusters () can reflect the variability of the infected clusters in dynamic model. In addition in order to characterize the inhomogeneity of the epidemic spreading, we denote the Homogeneous Mode where all the infected agents are randomly distributed and keep static in the system, is used to represent the case where the epidemic is homogeneously distributed. For the convenience of discussion, we denote the of homogeneous mode as . Therefore, for given infection density ρ, the larger is, the more inhomogeneous the epidemic spreading is.

THE DYNAMIC MODEL

In order to be closer to reality, we consider the case that agents may travel by transport to long distance. So, in this paper the agents not only may perform the local motion but also move into long position. It is assumed that agents M are distributed randomly in a square of linear size L with periodic boundary conditions initially. Each agent can jump into any position (i.e., long position) inside the square with the probability q or perform the local motion with the probability . Considering the fact that the jumping process is much shorter than the infection period, we assume that the status of the jumping agents will be kept when it move. At each time step , the local motion process is modeled as follows: where is the position of the agent i at time t. is characterized by v and . The absolute velocity v of all the agents is the same and remains constant in motion. is motion direction of the agent i at time t and follows the uniform distribution in the interval .

In the model, each agent can be in two distinct states: susceptible (S) or infected (I). Some agents are infected at the seed of the infection, and all the others are susceptible. And from then on, their positions and states are updated at each time step. Concretely, if the agent i 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 with probability at time where is the infection rate and is the infected neighbors set of the agent i. Here is defined as the infected agents in the spherical neighborhood of the interaction radius r centered on the agent i where denotes the Euclidean distance between i and j in two-dimensional space. If , , then i and j are neighbors and can infect each other. And we assume that each agent has the same interaction radius. Without lack of generality, let us set since it only affects the time scale of the infection evolution. Here, the above spread mechanism will run until the system reaches its stationary state or the epidemic dies out.

ENTROPY-BASED INFECTED CLUSTERS

In the dynamic model, each node represents an agent and each link represents a connection along which an epidemic can spread. 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 agents. The infected cluster size is the number of infected agents of the infected cluster. Let be the set of the infected clusters in the network. In this study, entropy-based infected clusters is defined as where is the probability that an infect agent belongs to of C. We can calculate by obviously, , which denotes the maximum value of , is equal to , when for each , where m is the number of the infected agents in the dynamic network. Similarly, the minimum value of for the dynamic network with m infected agents, denoted as , is equal to 0, which is obtained when the infected agents gather into one cluster.

The normalized entropy based on infected clusters is defined as From Eq. (5), when is approaching to 0, the infected agents will converge into a single cluster. When is approaching to 1, the infected agents will show a scattered distribution and the formed clusters appear with a similar size.

RESULTS AND DISCUSSIONS

In this section, we will discuss the inhomogeneity of the epidemic spreading while the agents walk randomly in the square (), and suppose . In all the simulation experiments, L is set to 10 and the number of agents M is fixed to 200. The initial proportions of the infected agents are 10%, while all the others start from the susceptible state. m is the number of infected agents. And then we update agents' state synchronously in the model for ( is set to 1) time steps. All data points shown in each figure are acquired by averaging over 50 runs.

In Fig. 1, one can find that the distributions of the infected clusters are different when the infection density is different and . When ρ is low, the infected agents can scatter small clusters with similar size as shown in Figs. 1(a)–1(c). When ρ is high, infected agents are prone to gather into large clusters in Figs. 1(d)–1(f). The CICS is cannot reflect the variability of the infected clusters in the system. However, can give information on the variability of the infected clusters. Next, we will discuss the inhomogeneity of the epidemic spreading by a comparable analysis on in the dynamic networks and in homogeneous mode.

FIG. 1.

The distribution of infected clusters in different infection density.

The moving of agents cannot affect the proportion of the infected agents in Ref. 13. Fig. 2 shows the combined effect of q and on the evolution of ρ for , . For a fixed v, one can find that ρ increases monotonously with . The reason for this is that when the infection rate is large, the probability that a susceptible agent is infected by its infected neighbors is also high. Thus, more susceptible agents get infected. In addition, the result of numerical simulation shows that ρ is close to 1 with the increase of and q.

FIG. 2.

ρ versus q and α.

When , , we can plot versus ρ with different q (as illustrated in Fig. 3). The comparable curve in homogeneous mode () is also given in Fig. 3 with black color. It shows that is close to when , which it means the infected agents are distributed homogeneously. Obviously, this conclusion is coincident with that of theoretical analysis. However, is apparently less than when and which it means the epidemic spreading is inhomogeneous. This conclusion by letting is coincident with that of previous results.

FIG. 3.

versus ρ, with the change of q.

To study the effect of the absolute velocity v on the inhomogeneity of the epidemic spreading when , and , Fig. 4 demonstrates the as a function of v for different values of q with 0.15, 0.25, and 0.35, respectively. To comparison, in Fig. 4(a) we plot six curves (0.3, 0.1, 0.5, 1, 2, 4, 10, respectively) of dynamic network and a straight line () of homogeneous mode when . Similarly, Figs. 4(b) and 4(c) show the situations when 0.5, 0.7, respectively. One can find that the variation of the inhomogeneity is related to v. When , is apparently less than for each ρ. That means, the epidemic spreads inhomogeneously. When , is close to for each ρ. Namely, the epidemic spreading is homogeneous. In fact, when the agents move with a high velocity, they have greater chance to jump out of the area that is covered by its infected neighbors than that in the case of small v. Therefore, the infected agents cannot gather into a large cluster. In addition, The infected density ρ affect the difference between the maximum and minimum of . When and (or ) when and (or )

FIG. 4.

distribution for different values of v.

Then we can say that the smaller is the infected density ρ, the stronger v affects the inhomogeneity of the epidemic spreading. And the finding is in accordance with the previous work.

Previous results have not discussed about the role of the interaction radius r on the inhomogeneity of the epidemic spreading. Fig. 5 presents as a function of for various v when , . One can find the epidemic spreads inhomogeneously when , . Next, suppose , we plot against r for different values of when , in Fig. 6. As can be seen in Fig. 6, is approaching to with the increase of r which it means that the epidemic spreading tends to a homogeneous state. Our results indicate the inhomogeneity decreases even disappears with the increase of r. And one can see that monotonously decreases with the increase of r. The reason for this is that the neighbors can infect each other in interaction neighborhoods which are determined by the radius r at each time step. As pairwise interactions increase with r, it is possible for a susceptible agent to be infected more by its infected neighbors. So the infected agents are prone to form large infected clusters.

FIG. 5.

versus with the change of v.

FIG. 6.

versus r for different values of α.

CONCLUSIONS

In this paper, we investigate inhomogeneity of the epidemic spreading by the normalized entropy based on infected clusters () in the dynamic network. can reflect the variability of the infected clusters in the system. The simulations show that decreases with the increase of the infection density ρ and the interaction radius r. That means, the infected agents prone to form large clusters as the infected agents always infect their neighbors. And the inhomogeneity of epidemic spreading decreases with the increase of v, r.