Intermittent social distancing strategy for epidemic control
We study the critical effect of an intermittent social distancing strategy on the propagation of epidemics in adaptive complex networks. We characterize the effect of our strategy in the framework of the susceptible-infected-recovered model. In our model, based on local information, a susceptible individual interrupts the contact with an infected individual with a probability and restores it after a fixed time . We find that, depending on the network topology, in our social distancing strategy there exists a cutoff threshold beyond which the epidemic phase disappears. Our results are supported by a theoretical framework and extensive simulations of the model. Furthermore we show that this strategy is very efficient because it leads to a “susceptible herd behavior” that protects a large fraction of susceptibles individuals. We explain our results using percolation arguments.
pacs:89.75.-k, 64.60.aq, 64.60.ah
The study of the topology of complex networks, and the dynamical processes that use these networks as substrate to spread, has recently generated great interest in the scientific community Boccaletti et al. (2006); Dorogovtsev and Mendes (2003); Romualdo Pastor Satorras (2004). In the past, studies of dynamic processes such as the spreading of rumors, opinions, and diseases on static networks, were concentrated on how their topology affects these processes Boccaletti et al. (2006); Dorogovtsev and Mendes (2003); Romualdo Pastor Satorras (2004). However, it is known that these processes evolve on top of networks where the topology changes with time Thilo Gross (2009). As a consequence, recently many researchers began to study dynamic networks, and the interplay between the dynamic process and the network dynamics. Those networks in which the topology changes regardless of the process taking place on top of them are called evolutive networks, while those networks that change their topology to mitigate or promote these processes are called adaptive networks Thilo Gross (2009); Lagorio et al. (2011). Adaptive networks have been investigated in many disciplines such as social sciences, epidemiology, biology, etc. Thilo Gross (2009); Schwartz and Shaw (2010); Vazquez et al. (2008). In these networks there is a coevolution between the link dynamics and the state of the nodes which leads to a collective phenomenon on adapting networks. As an example, in a network of routes, where the nodes are cities and the links are the routes connecting them, some overloaded paths become dysfunctional and new paths between cities are built to avoid the congestion. In the analysis of opinion formation on social networks, nodes usually tend to rewire or break their links with individuals with different opinions, leading in some models to a network fragmentation into components or clusters in which all members have the same opinion Holme and Newman (2006); Gil and Zanette (2006); Vazquez et al. (2008). In biological networks, such as the vascular system, after an arterial occlusion collateral vessels grows in order to increase the blood flow to neighboring tissues Schaper and Scholz (2003). Similarly, in the widely studied epidemic models on static networks, adaptive processes are used to model strategies that reduce the impact of the disease spreading Gross et al. (2006); Gross and Kevrekidis (2008); Lagorio et al. (2011).
One of the most popular models in epidemiology that reproduces seasonal diseases is the susceptible-infected-recovered model (SIR) Boccaletti et al. (2006); Anderson and May (1992), where individuals can be in one of three states, susceptible (), infected () or recovered (). In its classical formulation, an infected node infects a susceptible neighbor with probability and recovers with a certain fixed probability, which implies an exponential distribution of times for which individuals remain infected. However, this distribution is rarely realistic and for most seasonal diseases it has a sharp peak around an average value Karrer and Newman (2010). As a consequence, some studies Lagorio et al. (2009); Newman (2002); Parshani et al. (2010) have used a different version of the SIR model in which an infected individual recovers after a fixed time , called the recovery time.
It is well known that in the SIR model on static networks, the size of the infection is governed by the effective probability of infection or transmissibility of the disease, where . In turn, it was shown that this model can be mapped into a link percolation process Grassberger (1983); Newman (2002), where plays the role of the link occupancy probability in percolation. In a percolation process, there is a critical probability where the finite cluster size distribution behaves as in the thermodynamic limit. Above this threshold a “giant component” appears. As a consequence of the mapping between percolation and the SIR model; in the latter there is an epidemic threshold at below which the disease is an outbreak were the infection reaches a small fraction of the population, which is equivalent to having finite clusters in percolation, while above an epidemic develops corresponding to the emergence of a percolating giant component Miller (2007); Newman et al. (2001). This threshold, in uncorrelated static networks, depends only on the degree distribution , where is the degree or the number of links that a node can have. In particular, for Erdös-Rényi (ER) networks, where is the mean connectivity, the threshold is . However, in pure scale-free (SF) networks , where is the broadness of the distribution, in the thermodynamic limit for , which means that the epidemic spreads for any value of . This last result indicates that highly heterogeneous networks, such as many theoretical social networks, are very likely to develop an epidemic Newman (2002); Cohen et al. (2000).
While it is well known how the topology affects the SIR process in static networks, there is very little literature about this model on adaptive networks. Recently, Lagorio Lagorio et al. (2011) studied two different strategies to mitigate the spread of a disease in an adaptive SIR model. In that model a susceptible node disconnects the link with an infected individual with a rewiring probability , and creates a link with another susceptible node. The authors found that there is a phase transition at a critical rewiring threshold separating an epidemic from a non-epidemic phase, which can be related to static link percolation. A feature of this rewiring process is that links between susceptible individuals can be established independently of their previous relationship. In that strategy, the nodes have no memory because two individuals can be connected independently of their past. Even though this adaptive process could be representative of casual contacts between individuals such as those generated in public buildings like shopping, theaters, etc., that strategy will not work for other types of interaction such as friendship and working partners, where individuals preserve their closer contacts. Therefore, if an individual is separated from its closer neighbors, it will tend to reconnect with them, at some time, more often than with an unknown individual. On the other hand, recently Wang Wang et al. (2011) and Van Segbroeck Van Segbroeck et al. (2010) proposed strategies to stop the spread in the SIR model where the susceptible individuals driven by fear disconnect their links with their neighbors, infected or not, without creating a new link. This strategy reduces the number of contacts permanently, which is efficient but also very inconvenient from an economical point of view.
In this paper we propose a strategy based on “intermittent social distancing” in adaptive networks and study its efficiency in stopping the spread of diseases in theoretical and real networks. We found theoretically that in our model there exists a cutoff threshold that prevents an epidemic phase. Our results are supported by extensive simulations. Moreover, we found that the intermittent social distancing strategy is efficient to protect a large susceptible cluster. The paper is organized as following: in Sec. II, we derive the theoretical transmissibility for our model and show how our strategy diminishes the epidemic phase. In Sec. III, we show how the epidemic size is reduced with our strategy and the agreement between our theoretical approach and the simulations. In Sec. IV, we present a study of the “susceptible herd behavior” that we use as a criterion to evaluate the effectiveness of our strategy. In Sec. V, we present our conclusions.
Ii Analytical approach
We propose a SIR model in an adaptive network where an infected individual transmits the disease to a susceptible neighbor with probability and, if he fails, with probability the susceptible individual breaks the link with the infected one for a period . Thus, the effective probability of breaking a link is . After a time both nodes are reconnected and the process is repeated until the infected node recovers at a fixed time , , there is an intermittent connection between a susceptible node and its infected neighbor. This mimics a behavioral adaptation of the society to avoid contacts with infected individuals by imposing a social distancing during one or more periods of duration . Notice that in our model a susceptible node breaks its links using only local information and not global knowledge as provided by communication media. In our model, the time is increased by 1 after every infected node tries to infect its neighbors and the updates are done after each time step. In this process the dynamic transmissibility can be written as
where is given by
and denotes the integer part function. The first term of Eq. (1) is the probability that a node in the state transmits the disease to any neighbor node in the state (before it recovers), considering that the link - has never been broken. The second term represents the probability for a node in the state to transmit the disease to an node after that pair has been disconnected times for a period . The binomial coefficient of the second term takes into account the number of ways to arrange intermittent disconnected periods before the susceptible becomes infected individual at time . Notice that is the total times that the pair - is broken and represent a temporal social distancing. In Table 1, we illustrate the element of the second term of Eq. (1) with and .
It is known that the disease becomes an epidemic if the basic reproductive number , where is the number of secondary infections. For uncorrelated networks, is related to the connectivity distribution through the branching factor ,
where , and and are the first and the second moments of respectively. The branching factor is a measure of the heterogeneity of the network that diverges for SF networks with in the thermodynamic limit because Cohen et al. (2000). Then the critical cutoff , above which the disease dies, as a function of , , and , can be found through the condition , which yields
Cohen et al. (2000); Callaway et al. (2000) where is the critical transmissibility for the SIR model in static networks, and consequently our dynamic process in the steady state is related to a static topological property of the network. This is expected due to the fact that in our model the network topology does not change globally in the characteristic time scale of the disease spreading, and hence our strategy can be understood as an SIR model on a static network but with a transmissibility .
In Fig. 1, we plot the plane in order to show that with our strategy the epidemic phase is reduced compared to the static case. Notice that for the epidemic phase shrinks substantially compared to the case .
In the figure, the light-gray area, delimited between the curves corresponding to the blocking periods and , displays the region of parameters controlled by the intervention strategy. We can see that even the milder intervention expands the epidemic-free area compared to the static case. On the other hand, the stronger intervention drastically shrinks the epidemic phase as increases. From the figure, we can also see that for , as the social distancing increases the epidemic-free phase increases; meanwhile, for a theoretical SF network with , as , when there is no epidemic-free phase for any values of and . However, real networks are finite and they are not generally pure uncorrelated SF networks which implies that our intermittent social distancing strategy could be applied in these networks.
In order to determine how the heterogeneity of the network affects the strategy performance, in Fig. 2 we plot a phase diagram in the plane for fixed and different values of .
We can see that as the heterogeneity increases, the social distancing and have to increase in order to prevent the epidemic phase. Surprisingly, in high heterogeneous networks, we find that for , the critical cutoff probability is almost constant (see the Appendix A) with
, for very heterogeneous networks above , does not change the transmissibility; then the best and least expensive strategy is to reconnect at . This means that if we know the average duration of a disease , individuals can return to their activities with low risk just after half of the characteristic time. From Eq. (7), it is straightforward that for there exists an upper value of that depends only on , where is given by
where this strategy can be applied. From Eq. (6) the limit can also be expressed in terms of the minimum critical transmissibility,
As a consequence, for very low infection rates, our strategy predicts an epidemic-free phase even for highly heterogeneous networks with a finite epidemic threshold. Notice that since real networks have degree correlations and clustering, then the relation does not hold and the only magnitude that matters is , which has to be measured by the peak of the second giant component Wu et al. (2007). In those cases, for we have to use Eq. (9) instead of Eq. (8). We found that for the condensed matter coauthorship network Newman (2001) () and mathematics coauthorship network Palla et al. (2008) () the epidemic spread can be stopped in diseases with and , respectively. Similar values of were used for real networks by Kitsak Kitsak et al. (2010).
Iii Numerical results
In our model at the initial stage, all nodes are susceptible and we infect a randomly chosen node (patient zero) in the biggest component of our network. From the patient zero the disease spreads to its neighbors according to the rules of our model described above: an infected individual transmits the disease to a susceptible neighbor with probability and, if it fails, with probability the susceptible individual breaks the link with the infected one for a period . After a time both nodes are reconnected and the process is repeated until the infected node recovers at a fixed time . The time is increased by 1 after every infected node tries to infect its neighbors and the updates are done after each time step.
All our results are presented for but qualitatively all the results are the same for .
In Fig. 3, we plot the relative epidemic size as a function of for ER and SF networks for different values of and . From the plot we can see that decreases as and increase compared to the static case . We can also see that a critical probability exists, which can be obtained theoretically from Eq. (4), above which the disease dies. Then, depending on how virulent is the disease without intervention, we can control in order to stop the spread. In the insets of the figures, we collapse all the curves using the value of obtained from Eq. (4). The collapse close to the critical value shows the excellent agreement between the theoretical value of and the simulations.
Iv Susceptible herd behavior
With our intermittent social distancing strategy, the susceptible nodes dynamically reduce their contact with the nodes in the infected cluster, mitigating the spread of the disease, , our strategy produces a resistance to the disease which we call “susceptible herd behavior”. As a result of the coevolutive process, at the end of the spreading there is only one cluster composed of recovered individuals and, depending on , one or more susceptible clusters that give raise to a cluster size distribution of susceptible individuals or “voids” Exx (a). In our model, the cluster size distribution of voids or susceptible individuals is important since the formation of a susceptible herd behavior induced by the network dynamics also measures how effective our strategy is to preserve a whole part of the society safe from the disease. Next, we derive the value of the transmissibility for which a susceptible crowding develops, , the value below which our strategy is efficient.
Describing the growth of an epidemic cluster as a Leath process Leath (1976); Braunstein et al. (2007) for a value of the link occupancy probability and denoting by the probability that a cluster reaches the generation following a link, then the probability that a link leads to a giant component when is given by
where is the solution of
When the “epidemic” cluster grows, the size of the void clusters is reduced as in a node dilution process, since when a link is occupied a void cluster loses a node and all its edges. Then is the probability that a void cluster loses a node. If we denote by the fraction of void nodes removed the following relation holds Exx (b),
For node percolation, it is known that Stauffer and Aharony (1985) , where is the fraction of nodes in the giant void component and the number of finite void clusters of size . The fraction of remaining nodes, below which the giant void cluster is destroyed corresponds to the critical probability of node void percolation . At this value ; then . This means that at the void transition, only a fraction of the nodes belong to void clusters. As a consequence, the fraction of links needed to reach this point fulfills
where and is the solution of Eq. (14). Notice that at the void transition . A similar result was also found by Newman for a process of spreading of two pathogens Newman (2005), employing a different approach. Using the mapping between our adaptive SIR model and link percolation, the transmissibility needed to create a giant component of crowded susceptible individuals is . Thus, for a disease spreading in an ER network, the dynamical critical transmissibility for the giant susceptible cluster is given by .
In Fig. 4, we compare the cluster size distribution of susceptible individuals or voids for different values of the transmissibility where we distinguish outbreaks from epidemics Lagorio et al. (2009); Kenah and Robins (2007); Newman (2002) on an ER network with .
For we obtain only small clusters of susceptible nodes that decay faster than exponentially. On the other hand, for a giant susceptible cluster appears, meaning that our strategy is efficient producing large connected clusters of susceptible nodes that crowd in order to protect themselves. When with our strategy the transmissibility is reduced to the epidemic spreading slows down and the distribution decays as a power law with the same exponent as in a void node percolation transition, in contrast with the distribution of infected sizes of outbreaks which at criticality goes as as in a Leath process Leath (1976). These results confirms the importance of the interplay between void nodes and link percolation in our model. .
In this paper we present a novel adaptive strategy based on intermittent social distancing in an adaptive SIR model. In the intermittent social distancing strategy a susceptible individual breaks the link with the infected neighbor with a cutoff probability and then both individuals are reconnected after a time before the infected individual recovers at . Using the framework of percolation theory, we derive the dynamical transmissibility, and we find that there exists a critical cutoff where the epidemic spread is stopped for non-highly heterogeneous networks. We show that in some real networks our intermittent social distance strategy could stop the epidemic spreading for not very virulent diseases. We find an excellent agreement between the theory and the simulations. For heterogeneous networks, we find that a very high value of does not lead to a decrease of the epidemic, which implies that the less expensive strategy is to chose . Finally we verify that our strategy reduces the transmissibility below a value where a susceptible crowding is produced. Any clever strategy used to reduce the disease from spreading should protect the population at the least economic cost. We believe that our strategy, which allows us to control the disconnection of periods through , is a very convenient strategy because it creates a susceptible herd cluster. Our present findings could be used as a support and reference guidance for the development of further strategies to stop diseases in real networks.
Acknowledgements.This work was supported by UNMdP and FONCyT (Pict Grant No. 0293/2008). The authors thank Federico Vazquez and Camila Buono for useful discussions.
For and , in Eq. (2) the summation has only one term, because for we have only one period of disconnection, , . Intuitively, if then in order to compute the transmissibility, we have to consider at most only one break period since, otherwise the transmission of the disease is not possible. Denoting by , then Eq. (1) reduces to
Neglecting the higher powers of (which hold for , , or high values of ), we obtain
Notice that loses all the dependence on and as shown in Fig. 2. With this approximation and using the fact that , the critical cutoff is given by
If Eq. (4) holds, then
- Boccaletti et al. (2006) S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D. Hwang, Physics Reports 424, 175 (2006).
- Dorogovtsev and Mendes (2003) S. N. Dorogovtsev and J. F. F. Mendes, Evolution of Networks (Oxford, 2003).
- Romualdo Pastor Satorras (2004) A. V. Romualdo Pastor Satorras, Evolution and Structure of the Internet: a statistical approach (Cambridge, 2004).
- Thilo Gross (2009) H. S. Thilo Gross, Adaptive Networks: Theory, Models and Applications (Springer, 2009).
- Lagorio et al. (2011) C. Lagorio, M. Dickison, F. Vazquez, L. A. Braunstein, P. A. Macri, M. V. Migueles, S. Havlin, and H. E. Stanley, Phys. Rev. E 83, 026102 (2011).
- Schwartz and Shaw (2010) I. B. Schwartz and L. B. Shaw, Physics 3, 17 (2010).
- Vazquez et al. (2008) F. Vazquez, V. M. Eguíluz, and M. S. Miguel, Phys. Rev. Lett. 100, 108702 (2008).
- Holme and Newman (2006) P. Holme and M. E. J. Newman, Phys. Rev. E 74, 056108 (2006).
- Gil and Zanette (2006) S. Gil and D. H. Zanette, Phys. Lett. A 356, 89 (2006).
- Schaper and Scholz (2003) W. Schaper and D. Scholz, Thromb. Vasc. Biol. 23, 1143 (2003).
- Gross et al. (2006) T. Gross, C. J. D. D’Lima, and B. Blasius, Phys. Rev. Lett. 96, 208701 (2006).
- Gross and Kevrekidis (2008) T. Gross and I. G. Kevrekidis, EPL (Europhysics Letters) 82, 38004 (2008).
- Anderson and May (1992) R. M. Anderson and R. M. May, Infectious Diseases of Humans: Dynamics and Control (Oxford University Press, Oxford, 1992).
- Karrer and Newman (2010) B. Karrer and M. E. J. Newman, Phys. Rev. E 82, 016101 (2010).
- Lagorio et al. (2009) C. Lagorio, M. Migueles, L. Braunstein, E. López, and P. Macri, Physica A: Statistical Mechanics and its Applications 388, 755 (2009).
- Newman (2002) M. E. J. Newman, Physical Review E 66, 016128 (2002).
- Parshani et al. (2010) R. Parshani, S. Carmi, and S. Havlin, Phys. Rev. Lett. 104, 258701 (2010).
- Grassberger (1983) P. Grassberger, Math. Biosci. 63, 157 (1983).
- Miller (2007) J. C. Miller, Phys. Rev. E 76, 010101 (2007).
- Newman et al. (2001) M. E. J. Newman, S. H. Strogatz, and D. J. Watts, Phys. Rev. E 64, 026118 (2001).
- Cohen et al. (2000) R. Cohen, K. Erez, D. ben Avraham, and S. Havlin, Phys. Rev. Lett. 85, 4626 (2000).
- Wang et al. (2011) Y. Wang, G. Xiao, L. Wong, X. Fu, S. Ma, and T. H. Cheng, J. Phys. A: Math. Theor. 44, 355101 (2011).
- Van Segbroeck et al. (2010) S. Van Segbroeck, F. C. Santos, and J. M. Pacheco, PLoS Comput Biol 6, e1000895 (2010).
- Callaway et al. (2000) D. S. Callaway, M. E. J. Newman, S. H. Strogatz, and D. J. Watts, Phys. Rev. Lett. 85, 5468 (2000).
- Wu et al. (2007) Z. Wu, C. Lagorio, L. A. Braunstein, R. Cohen, S. Havlin, and H. E. Stanley, Physical Review E 75, 066110 (2007).
- Newman (2001) M. E. J. Newman, Proc. Natl. Acad. Sci. USA 98, 404 (2001).
- Palla et al. (2008) G. Palla, I. J. Farkas, P. Pollner, I. Derényi, and T. Vicsek, New J. Phys. 10, 123026 (2008).
- Kitsak et al. (2010) M. Kitsak, L. K. Gallos, S. Havlin, F. Liljeros, L. Muchnik, H. E. Stanley, and H. A. Makse, Nature Physics 6, 888 (2010).
- Exx (a) In Euclidean networks, the void or hole structure has been studied in many disciplines for example to characterize the morphology of the bones and to understand how habitat fragmentation affects animal movement processes.
- Leath (1976) P. Leath, Phys. Rev B 14, 5046 (1976).
- Braunstein et al. (2007) L. A. Braunstein, Z. Wu, Y. Chen, S. V. Buldyrev, T. Kalisky, S. Sreenivasan, R. Cohen, E. López, S. Havlin, and H. E. Stanley, I. J. Bifurcation and Chaos 17, 2215 (2007).
- Exx (b) Notice that is the fraction of links in a link percolation process, while is the fraction of void nodes for a void node percolation process.
- Stauffer and Aharony (1985) D. Stauffer and A. Aharony, Introduction to percolation theory (Taylor & Francis, 1985).
- Newman (2005) M. E. J. Newman, Phys. Rev. Lett. 95, 108701 (2005).
- Kenah and Robins (2007) E. Kenah and J. M. Robins, Phys. Rev. E 76, 036113 (2007).