Abstract
COVID-19 has infected millions of people, with deaths in more than 200 countries. It is therefore essential to understand the dynamic characteristics of the outbreak and to design effective strategies to restrain the large-scale spread of the epidemic. In this paper, we present a novel framework to depress the epidemic spreading, by leveraging the decentralized dissemination of information. The framework is equivalent to finding a special minimum dominating set for a duplex network which is a general dominating set for one layer and a connected dominating set for another layer. Using the spin glass and message passing theory, we present a belief-propagation-guided decimation (BPD) algorithm to construct the special minimum dominating set. As a consequence, we could immediately recognize the epidemic as soon as it appeared, and rapidly immunize the whole network at minimum cost.
Introduction
In recent years, many different epidemic types have frequently erupted and spread worldwide, leading to not only great economic losses, but also widespread social disruption [7, 11, 14, 19]. The ongoing COVID-19 pandemic in particular has infected more than 10 million people causing deaths in more than 200 countries [5, 6, 10, 11, 24]. It is therefore essential to understand the dynamic characteristics of epidemic outbreaks and to design effective strategies to restrain the large-scale spread of the epidemic [4, 9, 13, 16–18].
The most traditional method of controlling an epidemic outbreak is network immunization, which cuts the path of the spread, by immunizing parts of the network nodes [2, 20, 23]. The common solution of network immunization uses the centrality index of the network, in which the nodes with large degree, k-shell, articulation points, betweenness, or spectral coefficients are selected as the immunized nodes. Recently, the control strategy of an epidemic, based on the competing spread on top of a duplex network, has attracted a lot of attention [1, 3, 15–17]. It is known that the information propagation in a virtual social network plays an important role in controlling epidemic spreading. A representative scenario is that the infectious diseases spread in a real physical contact network of individuals, will naturally lead to the propagation of crisis awareness information on virtual social networks. The individuals who obtained the awareness information will take measures to protect themselves, resulting in the suppression of the spread of the infectious disease. This is the framework of the competing spreading process of disease and awareness information on the duplex network coupled with the physical network and virtual network. Similarly, [21] introduced the competing spreading framework of the computer virus and patch, in which the patch dissemination restrains the virus propagation while they spread through different channels [21]. However, both of these competing spreading frameworks suffer the problem that it is surely improper to assume all network nodes have the ability to monitor their neighbors and recognize the epidemic and would like to distribute the awareness information or patch.
To compensate for the shortcomings of the competing spreading framework, Takaguchi et al. proposed a novel control method which takes advantage of both traditional network immunization and competing spreading [12]. They introduced the observer node into the networks, which could monitor the state of its neighbors and could therefore recognize the epidemic if its neighbors were infected. In addition, the observer node could send awareness information to its uninfected neighbors (i.e. immunize its uninfected neighbors) if it recognized the epidemic. However, such a method focusses on a single network which ignores the fact that physical monitoring and information propagation usually takes place in different channels. As an extension, [8] introduced an improved model of deploying an observer node on multiplex networks, consisting of an epidemic spreading layer and an information spreading layer. Though the observer node-based method could monitor all the network nodes and efficiently recognize the infection, it cannot immunize all the network nodes when infection spreading is confirmed.
In this sense, Zhao et al. introduced a strategy of controlling an epidemic through the minimum dominating set (MDS) of the multiplex network [22]. They proposed a framework coupled by a duplex network (one layer is the immunization distribution network, another is the epidemic spreading network layer) and a central node and let the nodes of the MDS of the duplex network be the observer nodes. When one neighbor is infected in the epidemic spreading network layer, the observer could immediately recognize the infection and send the information to the central node. The central node would then distribute the crisis information or patches to all the observers and the observers then spread them to their neighbors in the immunization distribution network. Under the definition of MDS of the multiplex network, such a method could recognize the infection immediately as it happened and could also immediately immunize all the network nodes using the smallest number of observer nodes.
However, Zhao’s method belongs to the centralization strategy. As an improvement, in this paper we propose an MDS-based decentralized strategy. We introduce a special MDS of the duplex network which is a general dominating set for one layer and a connected dominating set for another layer. Using the spin glass and message passing theory, we present a belief-propagation-guided decimation (BPD) algorithm to construct the special MDS. As a consequence, we could monitor all the network nodes and immediately recognize the infection, and rapidly immunize all the network nodes when infection spreading is confirmed, while the cost is minimized.
Deployment of Observer Node Based on Special MDS of Duplex Network
Different from the Zhao’s method, our method only contains one duplex network where one layer is the immunization distribution network (IL), and the other is the epidemic spreading network layer (EL). Here, we introduce a novel definition of special MDS of the duplex network (see Figure 1). It is a minimum node set which is a general dominating set for the EL and a connected dominating set for the IL; at the same time, the size of the node set is minimum.
FIGURE 1
To control the epidemic spreading, we let the nodes of the special MDS of the duplex network be the observer nodes. The observer nodes could monitor its neighbors’ states in EL, thus could immediately recognize the infection after the neighbor was infected. In addition, the observer nodes could diffuse the information to its neighbors in IL after they recognized the infection or after receiving the information from its neighbors. Therefore, under the definition of the special MDS of the duplex network, the epidemic could be immediately recognized as soon as it appeared since all network nodes were monitored, and all network nodes could receive the related information due to the observer nodes being fully connected in IL.
Construction of Special MDS of Duplex Network
When considering the discussion in the above section, the problem of controlling epidemic spreading is converted to construct a special MDS of the duplex network. Obviously, the construction of the special MDS is NP-hard, and so we can only try to obtain a near optimal solution. Therefore, all the MDS discussed in the rest of the paper are actually a near-optimal special dominating set. In [22]; Zhao et al. proposed a framework of the construction method of the MDS of the multiplex network, but no detailed implementations are given. In this section, we first show the detailed procedure of the construction of general MDS of the duplex network, then a simple greedy strategy is given to construct the special MDS of the duplex network.
Construction of the General MDS of the Duplex Network
Let be a duplex network, where denotes the network layer i. In our framework, and represent the IL and EL, respectively. The two layers have the same set V of N nodes and a distinct set of edges. If D is a dominating set of G, it must be a dominating set of each network layer of G, which means every node in must be connected to at least one node of D in every layer. The state of node is called occupied, denoted by . Otherwise is called empty but observed and denoted as . For constructing the MDS of the duplex network, we introduce the following partition functionwhere represents a configuration of node states, and denotes the set of neighbors of i in . x is regarded as the positive re-weighting parameter. Obviously, in Eq. 1 implies i is occupied or observed in both layers; otherwise the term is equal to 0. Therefore, Z is only contributed by the configuration of the dominating set. In particular, only the MDS of the duplex network contributes to Z in the case of .
Solving Eq. 1 using the message passing theory, we have the expression of marginal probability of i being in state as followswhere represents the union of set and , i.e. the set of all neighbors of nodes i in both layers. denotes the joint probability of j being in state and i in under the assumption that node i is deleted from the network. is given bywhere
Eq. 9 is called the belief propagation equation which has been widely used in different kinds of optimization problems.
Based on Eqs.1–9, we can use the classical Belief-Propagation-Guided Decimation Algorithm (BPD) to identify the MDS of a multiplex network. As discussed in [22]; we first simplify the calculation of .
In case j has not been occupied but observed in at least one layer, the states of its neighbors in the corresponding layer could not be considered. Therefore, and can be further simplified aswhere represents the set of nodes which have been observed in , and denotes the intersection of and . In addition, can also be simplified if i has been observed in at least one layer
The procedure of BPD is now given by:
The states of all nodes are initialized to be unoccupied and unobserved. The initial values of all joint probabilities could set to be random numbers.
Perform iterative calculations of the joint probabilities based on Eq. 5 for rounds. Then calculate the occupation probability for all nodes using the obtained according to Eq. 2.
Select the first r unoccupied nodes with the highest occupation probabilities occupied. Then update the states of all remaining nodes.
Remove the edges between the nodes that have been observed and remove all isolated observed nodes in every layer.
If the remaining network still has unobserved nodes, calculate Eq. 5 for rounds and then calculate Eq. 2. Whether using Eq. 3 and Eq. 4 or Eq. 3 and Eq. 12 in Eq. 2, depends on whether i is unobserved or observed.
Repeat operations (3)–(5) until all nodes are observed in all layers.
Greedy Construction of Special MDS of Duplex Network
In this section, we construct the special MDS of the duplex network based on the general MDS constructed via the BPD algorithm. We delete all of the nodes of layer IL which do not belong to the general MDS and denote the remaining network as IL’. If IL’ is not fully connected, it must have more than one cluster. We then remove nodes form
into IL’ to reduce the number of clusters based on the following operations.
If the insertion of one node into the IL’ can reduce the number of clusters, we select the node whose insertion can reduce the number of clusters to be inserted into the IL’ the most.
If the insertion of one node into the IL’ cannot reduce the number of clusters, we select two linked nodes whose insertion can reduce the number of clusters to be inserted into the IL’ the most.
The operations are performed repeatedly until the IL’ is fully connected. In this way, a special MDS of the duplex network is constructed.
Result
In this section, we verify the effectiveness of our method by comparing it with the random algorithm and the simulated annealing algorithm.
The random algorithm is very simple. It randomly selects an unoccupied node and changes the node’s state to be occupied, until a special dominating set is constructed. The simulated annealing algorithm is widely used to solve many different kinds of optimization problems. To construct a special MDS of the duplex network, we define the following energy function
In the energy function,
indicates that node
iis occupied but has no occupied neighbor in layer IL, otherwise
.
indicates that node
iis observed in layer α, that is it has at least one occupied neighbor in layer α, otherwise
. In addition,
u,
v, and
ware tunable parameters. At the beginning of the SA algorithm, we construct a special dominating set for the duplex network based on the random algorithm. Then we calculate the corresponding energy
Ebased on
Eq.13. At each time step of SA, we perform the following operations:
Randomly select one node from the network. If the selected node is occupied, change its state to unoccupied. We then calculate the new energy .
Randomly select one node from the network. If the selected node is unoccupied, change its state to occupied. We then calculate the new energy .
The operation is accepted with probability , otherwise, the new occupation should be omitted. In this paper, we set . β represents the inverse of the temperature. Its initial value is set to be 0.5 and then is increased by at each time step. The SA algorithm terminates when β equals . It should be noted that the SA algorithm could not ensure that the constructed dominating set is a special dominating set, therefore, a greedy process is still needed after the termination of the SA algorithm.
We performed experiments on four different kinds of duplex networks which are coupled by the SF network and ER random network. The results are given in Table 1. We can see that our method constructs the smallest special dominating set for any duplex network. Therefore, we could use the smallest observer nodes to control the epidemic spreading based on our deployment scheme.
TABLE 1
| Network | Layers | N | Random | SA | Our | |
|---|---|---|---|---|---|---|
| 1 | ER (IL) | 10,000 | 40,085 | 9,784 | 4,462 | 3,407 |
| ER (EL) | 40,130 | |||||
| 2 | ER (IL) | 10,000 | 40,085 | 9,723 | 4,564 | 3,402 |
| SF (EL) | 39,987 | |||||
| 3 | SF (IL) | 10,000 | 39,986 | 8,810 | 4,669 | 2,384 |
| SF (EL) | 39,987 | |||||
| 4 | SF (IL) | 10,000 | 39,986 | 9,211 | 4,532 | 3,052 |
| ER (EL) | 40,130 |
The size of constructed dominating set of duplex network coupled by different network layers
Conclusion
In this paper, we designed a deployment scheme of observer nodes to control epidemic spreading based on the special MDS of the duplex network. Through the nodes of the special MDS of the duplex network, we could immediately recognize the epidemic as soon as it appeared and could immunize all network nodes—at the same time, the cost remains minimal. The deployment problem of the observer nodes is then converted to the construction problem of the special MDS of the duplex network. We used the BPD algorithm and a greedy algorithm to construct the special MDS of the duplex network. Our solution is verified to be very efficient on different kinds of duplex networks.
Funding
This work was supported in part by the National Natural Science Foundation of China (61702309), the Innovation and entrepreneurship training program of Civil Aviation University of China (201910059044), the National Key Research and Development Project (2018YFE0119700) and the Key Research and Development Project of Shandong Province (2019JZZY010132).
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
LZ and DZ proposed the framework. JW, WZ, YZ, WW, and XX developed the mathematical theory. JW and WZ performed the experiment. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
GaoBDengZZhaoD. (2016) Competing spreading processes and immunization in multiplex networks. Chaos, Solit Fractals93, 175–181. 10.1016/j.chaos.2016.10.013
2.
GaoCLiuJZhongN. (2010) Network immunization with distributed autonomy-oriented entities. IEEE Trans Parallel Distr Syst22, 1222–1229. 10.1109/TPDS.2010.197
3.
GranellCGómezSArenasA. (2014) Competing spreading processes on multiplex networks: awareness and epidemics. Phys Rev E - Stat Nonlinear Soft Matter Phys90, 012808. 10.1103/PhysRevE.90.012808
4.
GuoZ-GSunG-QWangZJinZLiLLiC. (2020) Spatial dynamics of an epidemic model with nonlocal infection. Appl Math Comput377, 125158. 10.1016/j.amc.2020.125158
5.
KrögerMSchlickeiserR. (2020) Gaussian doubling times and reproduction factors of the covid-19 pandemic disease. Frontiers in Physics8, 276. 10.3389/fphy.2020.00276
6.
LiH-JBuZWangZCaoJ. (2019a) Dynamical clustering in electronic commerce systems via optimization and leadership expansion. IEEE Transactions on Industrial Informatics16, 5327–5334. 10.1109/TII.2019.2960835
7.
LiH-JWangLZhangYPercM. (2020) Optimization of identifiability for efficient community detection. New J Phys22, 063035. 10.1088/1367-2630/ab8e5e
8.
LiZZhuPZhaoDDengZWangZ. (2019b) Suppression of epidemic spreading process on multiplex networks via active immunization. Chaos: An Interdisciplinary Journal of Nonlinear Science29, 073111. 10.1063/1.5093047
9.
LiuCZhouNZhanX-XSunG-QZhangZ-K. (2020) Markov-based solution for information diffusion on adaptive social networks. Appl Math Comput380, 125286. 10.1016/j.amc.2020.125286
10.
SaberiMHamedmoghadamHMadaniKDolkHDMorganAMorrisJKet al (2020) Accounting for underreporting in mathematical modelling of transmission and control of covid-19 in Iran. Frontiers in Physics8, 289. 10.3389/fphy.2020.00289
11.
SunGQWangSFLiMTLiLZhangJZhangWet al (2020) Transmission dynamics of covid-19 in wuhan, China: effects of lockdown and medical resources. Nonlinear Dyn101, 1981–1993. 10.1007/s11071-020-05770-9
12.
TakaguchiTHasegawaTYoshidaY. (2014) Suppressing epidemics on networks by exploiting observer nodes. Phys Rev E - Stat Nonlinear Soft Matter Phys90, 012807. 10.1103/PhysRevE.90.012807
13.
WangXJiaDGaoSXiaCLiXWangZ. (2020a) Vaccination behavior by coupling the epidemic spreading with the human decision under the game theory. Appl Math Comput380, 125232. 10.1016/j.amc.2020.125232
14.
WangZBauchCTBhattacharyyaSd’OnofrioAManfrediPPercMet al (2016) Statistical physics of vaccination. Phys Rep664, 1–113. 10.1016/j.physrep.2016.10.006
15.
WangZGuoQSunSXiaC. (2019) The impact of awareness diffusion on sir-like epidemics in multiplex networks. Appl Math Comput349, 134–147. 1016/j.amc.2018.12.045
16.
WangZXiaCChenZChenG. (2020b) Epidemic propagation with positive and negative preventive information in multiplex networks. IEEE Transactions on Cybernetics. Available online at: https://ieeexplore.ieee.org/abstract/document/8957067.
17.
XiaCWangZZhengCGuoQShiYDehmerMet al (2019) A new coupled disease-awareness spreading model with mass media on multiplex networks. Inf Sci471, 185–200. 10.1016/j.ins.2018.08.050
18.
YinQWangZXiaCDehmerMEmmert-StreibFJinZ. (2020) A novel epidemic model considering demographics and intercity commuting on complex dynamical networks. Appl Math Comput386, 125517. 10.1016/j.amc.2020.125517
19.
ZhaoDLiLPengHLuoQYangY. (2014a) Multiple routes transmitted epidemics on multiplex networks. Phys Lett378, 770–776. 10.1016/j.physleta.2014.01.014
20.
ZhaoDWangLLiSWangZWangLGaoB. (2014b) Immunization of epidemics in multiplex networks. PLoS One9, e112018. 10.1371/journal.pone.0112018
21.
ZhaoDWangLWangZXiaoG. (2018) Virus propagation and patch distribution in multiplex networks: modeling, analysis, and optimal allocation. IEEE Trans Inf Forensics Secur14, 1755–1767. 10.1109/TIFS.2018.2885254
22.
ZhaoDXiaoGWangZWangLXuL. (2020a) Minimum dominating set of multiplex networks: definition, application, and identification. IEEE Transactions on Systems, Man, and Cybernetics: Systems99, 1–15. 10.1109/TSMC.2020.2987163
23.
ZhaoDYangSHanXZhangSWangZ. (2020b) Dismantling and vertex cover of network through message passing. IEEE Transactions on Circuits and Systems II: Express Briefs99, 1–15. 10.1109/TCSII.2020.2973414
24.
ZhuangZCaoPZhaoSLouYWangWYangSet al (2020) Estimation of local novel coronavirus (covid-19) cases in wuhan, China from off-site reported cases and population flow data from different sources. Frontiers in Physics8, 336. 10.3389/fphy.2020.00336
Summary
Keywords
COVID-19, epidemic spreading, minimum dominating set, multiplex network, control
Citation
Wang J, Zhang L, Zhu W, Jiang Y, Wu W, Xu X and Zhao D (2021) The Suppression of Epidemic Spreading Through Minimum Dominating Set. Front. Phys. 8:588513. doi: 10.3389/fphy.2020.588513
Received
29 July 2020
Accepted
19 October 2020
Published
14 January 2021
Volume
8 - 2020
Edited by
Hui-Jia Li, Beijing University of Posts and Telecommunications (BUPT), China
Reviewed by
Chengyi Xia, Tianjin University of Technology, China
Gui-Quan Sun, North University of China, China
Updates
Copyright
© 2021 Wang, Zhang, Zhu, Jiang, Wu, Xu and Zhao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lei Zhang, zhanglei@sdas.org
This article was submitted to Social Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.