It is well observed that some infectious diseases show a feature of periodicity; that is, the disease may prevail during a certain season and vanish afterwards. In the effort of understanding this specific phenomenon, in this paper, we propose a spreading model with a time-variant infectivity function considering both the features of periodicity and agent variability. We apply the modified model to a scale-free network frame with tunable power-law coefficient to find out the characteristics of periodic spreading dynamics. Our work consists of both theoretical derivation and a series of corresponding numerical simulations in order to find out the influence of the parameters of both network topology and the modified model on spreading dynamics in the underlined networks. The experiments prove the success of our model in producing a periodic behavior of spreading, and the results agree well with theoretical calculations.
National Science Foundation of China114471471. Introduction
The spreading process on top of complex systems is one of the most focused subjects in the realm of network science and intensively studied [1–6]. A large range of realistic systems such as communication network, transportation web, power grid, or large-scale supply chains [7–9] could undergo some sort of spreading processes and cause potential economic losses. The last fifteen years have seen the successful application of complex network theory on the analysis of dynamical behaviors on top of complex systems [10–14]. From the perspective of this theory, the individuals of a real-world complex system are modeled by agents in a network and the contacts among them are considered as the links connecting agents. A healthy agent may become infectious through the already infected agents that are linked to it. On this common basis, researchers have integrated the network view into a series of compartmental models like SI (Susceptible-Infected), SIS (Susceptible-Infected-Susceptible), and SIR (Susceptible-Infected-Recovered) among many others and have unveiled some exciting features of networked spreading dynamics [15–17]. Reference [18] considered the effect of medium by means of a modified SIS model and recovered that spreading medium could accelerate epidemic spreading without changing the final spreading scale. Reference [19] established R-SIR-R (Risk-SIR-Random) and R-SIR-A (Risk-SIR-Acquaintance) on the basis of cellular automata to study the safety risk spreading of complex information system on top of nearest coupled network, random network, WS small-world network, and BA scale-free network and concluded that the outbreak threshold decreases with the strengthening of network’s heterogeneity. Reference [4] built a SVIR (Susceptible-Vaccinated-Infected-Recovered) model and compared various vaccination strategies; the results showed that both the vaccination level and infection level are lower when adopting continuous strategy other than pure strategy.
Nevertheless, the compartmental models are too simplified to describe various spreading processes in real-world systems partly because the infectivity is treated as linear. To this end, nonlinear infectivity is proposed and many delicate models are put forward. Chu et al. in [20] proposed a epidemic spreading model in weighted SF network; they investigated the nonlinear infectivity and weighted transmission rate and found the strong effect of the infectivity parameter on the epidemic threshold and the prevalence. Reference [21] investigated a new SIS model with nonlinear infectivity as well as birth and death of agents and links, and the reproductive number is calculated mathematically.
In spite of the intensive research, many valuable problems still remain open. Why do some popular diseases (such as flu) show a notable feature of periodicity? What is the mechanism of this phenomenon? Note that a large number of spreading processes in practical systems display a feature of periodicity [22, 23]. Reference [22] studied the spreading in community networks, and Zhou et al. found periodic behaviors in both SIS and SIRS models when the infection rate is larger than the threshold which may be helpful in understanding the spatial-temporal pattern of DHF incidence in Thailand. Reference [23] proposed a modified SIRS (Susceptible-Infected-Recovered-Susceptible) model and found that the time-delay of the assessment of epidemic risk can cause periodic outbreak of diseases. In the above papers, researchers tried to demonstrate the mechanism of periodicity from their own perspectives; however, they all adopted a constant effective infectivity rate in the spreading model regardless of the fact that the contagion rate of epidemics could be time-variant. Moreover, agent variability is also advised to be considered in the confirmation of effective infectivity rate.
In this paper, we take both the factors into account. We apply the new model to topologically alterable scale-free networks and discuss the influence of network’s structural parameters, agent variability, and periodicity on the epidemic spreading behaviors. The results show that our model is successful in describing the periodic feature, and some influential parameters are proposed and investigated.
2. Modified SIS Model
There are two sorts of agents in SIS model (S and I) which are susceptible and infected. In this case, a susceptible agent can be infected when it is in contact with an already infected one with probability η. Once infected, agents recover with some probability β per unit time. The dynamics of SIS model can be described by the following system of ordinary differential equations:(1)dstdt=-ηitst+βit,ditdt=ηitst-βit,where s(t) (i(t)) indicates the number of susceptible (infected) agents at moment t. Without loss of generality, let β=1; hence the effective infectivity rate λ=η/β=η. Based on mean field theory, we can get the main equation:(2)dρktdt=-ρkt+λk1-ρktΘρt,where the first term on the right side is the extinction term and the second one is the production term. The meaning of each parameter is as follows: (a) ρkt is the probability that the agent of degree k gets affected; (b) λ is the probability that a healthy agent becomes infected when being in contact with an infected agent; (c) ρt=∑kPkρkt is the average proportion of all the infected agents in the network and can be used to measure the spreading prevalence. Here, we call it spreading scale in the remaining part of this paper and take it as an important meter of spreading dynamics; (d) Θρt is the probability that a given edge of a certain agent is pointed to an infected neighbor.
When the degree-degree interrelationship is weak, we have(3)Θ=∑kkPkρktk,where 〈k〉 is the mean degree of the network. Apparently, Θ∈[0,1]; the sufficient condition that the spreading process sustains in the system is dρkt/dt≥0. We take the critical value dρkt/dt=0; then we have(4)ρk=λkΘ1+λkΘ.
Taking agent variability into consideration, it requires(5)λk=λkkαkα,where α<0, and λk is the average of λk. We name α the agent variability weight for it defines the extent of agent difference. To make sure that λ∈0,1, we need two conditions to be satisfied:
Condition 1: λ<γ-1/γ-α-1.
Condition 2: kmin>γ-α-3/γ-2,
where kmin is the smallest degree of all agents and γ the power-law coefficient. The above inequations are easy to be satisfied in practice.
After agent variability, we take periodicity into consideration by modifying (5) as(6)λk,t=λk,tkαkα1+εft.
ft is a continuous function with period T, -1≤ft≤1; ε is the amplitude, 0≤ε≤1. Also, we have ∫aa+Tftdt=0,∀a∈R.
Revising (4), we can get(7)ρk=λk,tkΘ1+λk,tkΘ.
Integrating (7) into (3), we will have an autonomous equation of Θ:(8)Θ=∑kPkk·λk,tk2Θ1+λk,tk2Θ.
There exists a trivial solution Θ=0. To allow for nontrivial solutions, there must be(9)ddΘ∑kPkk·λk,tk2Θ1+λk,tk2ΘΘ=0≥1.
Clearly, there is(10)λkα/kα1+εftk2k≥1.
With λkα/〈kα〉1+εftk2/k≥λkα/kα1-εk2/k, we obtain the sufficient condition of (10):(11)λkα/kα1-εk2k≥1.
After solving it, we obtain(12)λ≥kαk1-εkα+2.
The outbreak threshold that sustains the spreading process on top of the network comes as(13)λ0c=kαk1-εkα+2.
λ0c is related with α, ε, and k. Comparing it with the threshold of traditional SIS model on top of scale-free network which is λc=k/k2, it is quite clear that λ0c>λc. In other words, epidemics are more difficult to sustain on the network when adopting our modified model. Specially, ε must to be set to a small value to avoid error enlargement.
3. Experiments and Discussion
In the above sections, we proposed the modified SIS spreading model which involves an abstract periodic function ft (-1≤ft≤1). Without loss of generality, let ft=sinσπt for simulative convenience’s sake. We conduct the simulations by distributing initiative infected agents randomly on the network with probability 20%.
3.1. Network Model
There was a prevailing consumption for a long period of time that all complex systems could be considered mathematically as random networks until Barabási and Albert [11] falsified that consumption. Their work on a WWW project indicates that the degree distribution of many practical networks submits power-law distribution and they proposed the BA scale-free network model with power-law coefficient fixed to 3 by preference attachment and growth mechanism. Although BA model is able to picture some realistic networks such as the Internet, most real-world complex systems which can be abstracted as scale-free networks do not submit a fixed coefficient power-law distribution. As a result, [24, 25] built a network construction mechanism which allows for a controllable power-law coefficient. The details are as follows.
By distributing a normalized weight to each of the N isolated agents and adding links among the agents with some probability, we get a tunable power-law coefficient and mean degree by adjusting the total edge number. The algorithm is as follows.
Step 1.
Distribute weight qi=i-ω to agent i, i=1,2,…,N, where N is the total number of agents; the tunable parameter ω is within interval 0≤ω<1. Normalize the weight as(14)qi∗=qi∑kqk.
Step 2.
Add edges among agents. If agent i and j satisfy(15)qi∗≥r1>qi+1∗,qj∗≥r2>qj+1∗,where r1 and r2 are random numbers in 0,1, and if i and j are not connected, add an edge between them. By this rule, isolated nodes with high weights are more likely to be linked.
Step 3.
Repeat Step 2 till there are mN edges.
The mean degree of this network is 2m, and the degree distribution satisfies power-law Pk∝k-γ, where k is the degree of agent and γ satisfies(16)γ=1+1ω.
Apparently, we are able to adjust the network topology by tuning parameters m and ω.
3.2. Time Evolution
To reveal the new characteristics of the modified spreading model built in Section 2, we compare it with three comparative models in the following part and analyze their simulative performance. The models we are going to discuss are listed as follows:
Model A: SIS model with no periodic infectivity and no agent variability.
Model B: SIS model with periodic infectivity and no agent variability.
Model C: SIS model with agent variability and no periodic infectivity.
Model D: SIS model with both periodicity and agent variability.
The interactive scale-free network has the size N=500 and mean degree k=10. Parameters in comparative spreading models are set as follows: σ=0.01, ε=0.2, α=-0.5, and λ=0.5 for periodic infectivity models and λ=0.5 for constant infectivity models.
The three curves shown in Figure 1 correspond, respectively, to the time evolution of Model A, Model B, and Model B’s infectivity rate λ. Similarly, Figure 2 shows the time evolution of Models C and D and λ. In particular, the vertical location of λ is adjusted to get a better view of comparison.
Evolution of Model A, Model B, and its infectivity rate λ. (a) For Model A, the spreading scale oscillates violently within the interval (0.46, 0.5) on a cycle of two time steps; nevertheless, Model B shows a more interesting feature; the evolution curve of Model B displays a dual-cycle characteristic where the spreading scale evolutes on a minor cycle and a major cycle simultaneously. (b) is a magnified view of the revolution of the two models from which we can clearly see the minor cycles.
Evolution of Model C, Model D, and infectivity rate λ. In (a), it is shown that the evolution of both Models C and D is asymptotically stable where Model D curve is synchronized with the curve of infection rate. (b) is a magnified view of Models C and D.
As shown in Figure 1, the evolution curve of Model A oscillates substantially within interval (0.44, 0.51) on a minor cycle which is about 2 time steps. Meanwhile, when considering periodicity in Model B, a dual-cycle effect is observed on the curve: on one hand, the curve oscillates on a short cycle as in Model A; on the other hand, Model B curve synchronizes with the infectivity rate on a major cycle which is 200 time steps according to our initial settings. By a simple modification of the SIS spreading model such as involving periodicity in infectivity, we can observe a realistic phenomenon in spreading dynamics that the spreading scale shows a periodic feature. In our modification, the infectivity rate λ is a sine function instead of a constant number as in most researches, which causes the spreading scale to display a pseudo-sine function as well. Moreover, we have tested cosine function and several other functions for λ; we found the tendency of spreading scale during evolution in accordance with that of the infectivity rate λ in all the results. To this end, we have reason to believe that our modified model on the infectivity rate provides an explanation for the periodic phenomenon in real-world spreading processes.
As to Models C and D, the simulation results are shown in Figure 2. After entering the steady stage, Model C curve is surprisingly steady and keeps in an almost horizontal state. Besides, Model D curve quite assembles that of Model A only with a smaller amplitude of the minor cycle. Moreover, if we take a close look at the curves in the figure, we can find that the amplitudes of both the models reduce as time goes on, but with a relatively low speed. For a better view, we compared Model A with C in Figure 3 where Model C involves agent variability and Model A does not and Model B with D in Figure 4 where Model D involves agent variability and Model B does not; by the clear difference shown in the figures, we can conclude with evidence that agent variability minishes the periodic amplitude.
Comparison of Model A with Model C. In Model A, the spreading model considers no periodicity of infectivity or agent variability, while in Model C agent variability is taken into consideration.
Comparison of Model B with Model D. In Model B, only infectivity periodicity is taken into consideration while Model D considers both infectivity periodicity and agent variability.
3.3. The Influence of Network Topology
In the art of complex network, researchers’ main focus lies in the interplay of network structure and the dynamics which occur on the network. In this subsection, we adopt the network model proposed in Section 2 and simulate the influence of power-law coefficient, mean degree, and so forth on the spreading dynamics. The initial setting of the modified SIS model is as follows: σ=0.01, ε=0.2, α=-0.5.
3.3.1. Power-Law Coefficient
The power-law coefficient is an indicator of homogeneity of a scale-free network. A big value of the coefficient corresponds to a high level of uniformity of degree distribution, and a more negative value of the coefficient corresponds to a low level of uniformity of degree distribution.
Figure 5 displays the plots of spreading scale ρ with respect to λ under several values of γ. The simulation results indicate that they basically agree with the theoretical calculation by comparison of Figure 5 and (13), but the simulation thresholds are a bit larger than the theoretical ones; this is because the simulation curve is an average of 200 independent simulations; in some simulations the epidemic does not survive. Let us take the situation of γ=2 as an example, λ0c=0.0609 according to (13); meanwhile, in the figure, the black line starts to ascend in the interval [0.06, 0.07]. It is also observed from the same figure that ρ is negatively related with γ, which means that the spreading scale descends with the homogeneity of the underlined network. This is because, in a more homogeneous network, the number of hub nodes (those nodes with rather high degrees) is reduced; as a result, the influence of “super spreaders” on their neighbors is restrained which leads to a lower level of spreading.
ρ-λ curve with variable γt→∞. As the power-law coefficient increases, the spreading outbreaks more difficultly, indicating the enhancement of network robustness against epidemic propagation.
3.3.2. Mean Degree of the Network
The mean degree indicates the average thickness of a network; this metric reflects the average connection number that a node has. The parameter setting of the spreading model is as follows: ε=0.2, σ=0.01, and α=-0.5.
We can see that the influence of mean degree on the outbreak threshold is monotonous from Figure 6. The theoretical outbreak thresholds corresponding to k=14, k=18 and k=22 are, respectively, 0.0797, 0.0683, and 0.0562 which are a bit smaller than simulation results due to the same averaging effect as in Section 3.3.1. The outbreak threshold is negatively related with mean degree of the network; in other words, the more connected the agents are to each other, the easier the diseases break out and stay active in the system. The reason for this tendency is that a more highly connected network provides more routes for the spreading process to expand, which leads to a low outbreak threshold.
ρ-λ curve with variable kt→∞. The more inner-connected the network is, the more easily the epidemic breaks out.
3.4. The Influence of Modified SIS Model
On the basis of the intensively studied linear infectivity, we propose a nonlinear periodic infectivity rate model with the consideration of agent variability and apply it to our network model. The modified spreading model we put forward includes two adjustable parameters: the cycle T and agent variability weight α.
3.4.1. Infectivity Rate Cycle
From the above discussion, we know the cycle of infectivity rate affects that of the spreading scale. In this subsection, we discuss its influence on the outbreak threshold and the spreading scale through numerical simulation. The parameter settings of the network and modified SIS model are as follows: N=500, k=10, ε=0.2, and α=-0.5.
The curves of ρ with respect to λ when T=50, T=20, and T=10 are displayed in Figure 7. Clearly, the cycle of infectivity rate has little effect on the outbreak threshold and the final spreading scale. To illustrate this with more evidence, we simulated the time evolution of ρ when T=200 and T=50 shown in Figure 8. It can be seen that after t≥400, the evolution enters a steady state. Both the curves are sine-wave shaped and oscillate within the same interval; the only difference lies in the oscillation cycle.
ρ-〈λ〉 curve with variable T (t→∞).
Evolution of ρ with two kinds of periods. There is no difference in the oscillation amplitude between the cases when T=200 and when T=50.
3.4.2. Agent Variability Weight
In a realistic system, each agent’s ability to avoid contagions varies because different agents have different local properties such as the degree property. Based on that idea, it is all natural to consider the infectivity rate as degree related. Obviously, it is reasonable to assume that agents with high degrees are more likely to contact with infected neighboring agents.
The interactive network we use has the minimum degree kmin=2 which satisfies Condition 2 in Section 2. We plot ρ versus λ when α=0, α=-0.5, α=-1, and α=-1.5. The results are demonstrated in Figure 9. We compare the simulated thresholds to theoretical results which are 0.0677, 0.0887, 0.1152, and 0.1511 and find that they fit well. As shown in (5), as α goes further from 0, high degree agents are less likely to infect neighbors in accordance with (5) which causes the enhancement of network robustness against epidemic spreading outbreak.
ρ-〈λ〉 curve with variable α (t→∞).
4. Conclusion
The variation of infectivity rate of epidemic is a realistic consideration. Based on that, we modified the conventional model by adding periodic feature and agent variability to the infectivity function. We discussed the spreading dynamics of our model on top of topologically alterable scale-free networks by both theoretical and simulative measures. It is found that the periodic infectivity rate causes the final spreading scale to oscillate with the infectivity. Also, a dual-cycle phenomenon is observed; meanwhile, agent variability decreases the amplitude of the minor cycle. The network topology influences the spreading dynamics in several aspects. The power-law coefficient which indicates the heterogeneity of network has positive relationship with the spreading threshold but is positively related with the final spreading scale and the mean degree of network is negatively related with the threshold. Furthermore, the cycle of infectivity rate affects only the cycle of final spreading scale. Agent variability weight and infectivity periodic amplitude are both positively related with the spreading threshold.
Competing Interests
The authors declare that they have no competing interests.
Authors’ Contributions
Z. Wang and H. Han wrote the paper and did the theoretical analysis, C. Ding did the simulation, and H. Yao and J. Du reviewed the paper.
Acknowledgments
This research is funded by the National Science Foundation of China (Grant no. 11447147).
ChenS.HuangW.CattaniC.AltieriG.Traffic dynamics on complex networks: a survey201220122373269810.1155/2012/732698MR28374542-s2.0-81555199726DickisonM.HavlinS.StanleyH. E.Epidemics on interconnected networks201285606610910.1103/physreve.85.0661092-s2.0-84862226520LiM.DuX.VillaruzA. E.DiepB. A.WangD.SongY.TianY.HuJ.YuF.LuY.OttoM.MRSA epidemic linked to a quickly spreading colonization and virulence determinant201218581681910.1038/nm.26922-s2.0-84860690538CaiC.-R.WuZ.-X.GuanJ.-Y.Effect of vaccination strategies on the dynamic behavior of epidemic spreading and vaccine coverage201462-63364310.1016/j.chaos.2014.04.0052-s2.0-84899636751GuoQ.LiL.ChenY.YangY.NiuX.Modeling dynamics of disaster spreading in community networks2011641-215716510.1007/s11071-010-9854-8MR27829662-s2.0-79953029186WangZ.YaoH.YangG.DuJ.Failure risk propagation and protection schemes in coupled systems201580627510.1016/j.chaos.2015.06.005MR3407003LüL.MedoM.YeungC. H.ZhangY.-C.ZhangZ.-K.ZhouT.Recommender systems2012519114910.1016/j.physrep.2012.02.0062-s2.0-84866757975HolmeP.SaramäkiJ.Temporal networks201251939712510.1016/j.physrep.2012.03.0012-s2.0-84866939709TayurS.GaneshanR.MagazineM.201217Berlin, GermanySpringer Science & Business MediaWattsD. J.StrogatzS. H.Collective dynamics of ‘small-world’ networks1998393668444044210.1038/309182-s2.0-0032482432BarabásiA.-L.AlbertR.Emergence of scaling in random networks1999286543950951210.1126/science.286.5439.509MR20916342-s2.0-0038483826LiuQ.YiJ.The research of the social network evolution based on the evolutionary game theory2013622323890210.7498/aps.62.2389022-s2.0-84889775654WilkinsonR. R.SharkeyK. J.Message passing and moment closure for susceptible-infected-recovered epidemics on finite networks2014892602280810.1103/PhysRevE.89.022808PengX.-Z.YaoH.DuJ.DingC.ZhangZ.-H.Study on cascading invulnerability of multi-coupling-links coupled networks based on time-delay coupled map lattices model201463707890110.7498/aps.63.0789012-s2.0-84897407939SongY.-R.JiangG.-P.GongY.-W.Epidemic propagation on adaptive coevolutionary networks with preferential local-world reconnecting strategy201322404020510.1088/1674-1056/22/4/0402052-s2.0-84876146972XiaC. Y.LiuZ. X.ChenZ. Q.YuanZ. Z.SIRS epidemic model with direct immunization on complex networks2008234468472MR2406683ZanetteD. H.Dynamics of rumor propagation on small-world networks200265404190810.1103/physreve.65.0419082-s2.0-41349109719ShiH. J.DuanZ. S.ChenG. R.An SIS model with infective medium on complex networks20083878-92133214410.1016/j.physa.2007.11.0482-s2.0-38349171523LiZ.XuG.-A.BanX.-F.ZhangY.HuZ.-M.Complex information system security risk propagation research based on cellular automata2013622020020310.7498/aps.62.2002032-s2.0-84886636830ChuX.ZhangZ.GuanJ.ZhouS.Epidemic spreading with nonlinear infectivity in weighted scale-free networks2011390347148110.1016/j.physa.2010.09.0382-s2.0-78649951759ZhuG.FuX.ChenG.Global attractivity of a network-based epidemic SIS model with nonlinear infectivity20121762588259410.1016/j.cnsns.2011.08.039MR28777032-s2.0-84855247936ZhouY.-Z.LiuZ.-H.ZhouJ.Periodic wave of epidemic spreading in community networks200724258110.1088/0256-307x/24/2/0782-s2.0-34247199991ZhangH.ZhangW. Y.Time-delayed information can induce the periodic outbreaks of infectious diseases201242631638HolmeP.KimB. J.Growing scale-free networks with tunable clustering2002652402610710.1103/physreve.65.0261072-s2.0-41349104828GohK.-I.KahngB.KimD.Universal behavior of load distribution in scale-free networks20018727427870110.1103/PhysRevLett.87.2787012-s2.0-0035981014