Abstract

Virus spreading on the Internet will negatively affect cybersecurity. An intermittent quarantine immunization strategy to control virus spreading when containing information diffusion is proposed herein. In this model, information and virus spread on different subnetworks and interact with each other. We further develop a heterogeneous mean-field approach with time delays to investigate this model and use Monte Carlo simulations to systematically investigate the spreading dynamics. For a relatively short intermittent period, the optimal information transmission probability of the virus will be significantly suppressed. However, when the intermittent period is extremely long; increasing the probability of information transmission can control the virus spreading as well as suppress the increase in the intermittent period. Finally, it is shown that the average degree of the two subnetworks does not qualitatively affect the spreading dynamics.

1. Introduction

Computer viruses are spreading widely on the Internet, thereby significantly affecting the cyberspace security [14]. For instance, the WannaCry virus, one of the most influential viruses since Panda Burned Incense, has affected more than 100 countries and regions. Researchers from different fields, including computer science, network science, physics, and mathematics, have attempted to eliminate the spread of the virus [511]. In this study, we focused on network science to address this problem, in which the ultimate aim is to develop effective measures to control virus spreading using mathematical modeling and network theory.

Network science researchers have developed many effective models to investigate the spread of viruses on the Internet. Data analyses revealed that the topology of the Internet exhibits the heavy-tail degree distribution [1214]; hence, PastorCSatorras and Vespignani proposed a mathematical susceptible-infected-susceptible (SIS) model on scale-free complex networks [5, 15]. They used a heterogeneous mean-field theory to describe the dynamics and revealed that a few hubs resulted in vanishingly low values of infection transmission probability that triggered virus spreading on the Internet. Based on this finding, we can design an effective approach to control virus spread. Target immunization [1519] is one of the most popular strategies, which involve the immunization of hubs to suppress viruses significantly. To implement the target immunization strategy, the topology of the Internet must be determined, which limits its applications. Therefore, researchers proposed the acquaintance immunization strategy [2023], which necessitates only information regarding local network structures. The progress in state-of-the-art network immunization has been reported in a recent review [24].

Computer users typically communicate with their friends through social networks [2528]. Hence, the information spread by a virus will diffuse in the social network when the virus is spreading on the Internet. Therefore, researchers have developed models to describe the dynamics of virus-information spread. Scholars typically use epidemic spreading models to model the spread of viruses, such as SIS and susceptible-infected-removed (SIR) models. Therefore, we herein present some of the advancements in the field of epidemic-information spreading dynamics that can be used to investigate the spread of virus information. Funk et al. [29] assume that an epidemic is spread on homogeneous networks, whereas information regarding the epidemic is also spread within the same network. Using a mean-field approach, they discovered that information spreading can significantly suppress the spread of an epidemic but cannot alter the outbreak threshold. Granell et al. [30, 31] used the unaware-aware-unaware (UAU)-SIS model, which assumes UAU information spreading dynamics on social networks and SIS epidemic spreading dynamics on face-to-face contact networks. They used a discrete Markovian approach to describe the interacting dynamics and revealed that the epidemic outbreak threshold was determined by the topologies of social networks, contact networks, and dynamic information spreading parameters. Subsequently, researchers investigated the effects of network topologies and interacting mechanisms on spreading dynamics [3237]. Wang et al. [3840] used an asymmetrical interacting SIR-SIRV model and discovered that the interlayer degree correlation facilitated the control of epidemic spread.

To immunize virus spreading on the Internet, different strategies can be adopted, e.g., barring Internet access, implementing patches, and turning on the firewall. Additionally, advanced information-epidemic spreading dynamics can be adopted to describe more realistic scenarios. Zhao et al. [41] developed a virus-information spreading model to design an optimal allocation strategy for patches. Once a virus spreads on the Internet, the simplest approach to stop the spreading is to bar Internet access until the virus is eliminated [4244]. Because mathematical studies pertaining to information spreading dynamics are scarce, we herein propose a mathematical model to describe an interaction model based on multiplex networks. We assume that information regarding the virus spreads in the social network, i.e., the virus spreads on the Internet and then performs extensive numerical simulations to determine the optimal intervention strength.

2. Model Descriptions

In this section, an interacting spreading dynamics model based on the intermittent quarantine immunization strategy on multiplex networks with nodes is described. The multiplex networks comprise two subnetworks, denoted as and , and each subnetwork represents a communication platform. Each node in two subnetworks is matched one-to-one, which implies that a node exists in different communication platforms. The two subnetworks were constructed based on an uncorrelated configuration model with degree distributions and , separately. When the network size is , no inter- and interlayer degree correlations exist in the multiplex subnetwork (as shown in Figure 1(a)).

Assume that subnetworks and represent the social and virus networks, respectively. Information is typically spread on social network , whereas virus spreads on subnetwork . To describe information spreading, we used the classic SIR epidemiological model [4547], in which each node can be susceptible, informed, or recovered at a specified time. At each time step, each informed node first transmits the information to every susceptible neighbor with probability and then recovers with probability . It is noteworthy that the susceptible node in subnetwork is informed once its counterpart in subnetwork becomes infected. Hence, virus spread promotes information spreading (as shown in Figure 1(b)).

For virus spreading on subnetwork , we used a generalized SIR model. The virus spreading dynamics is the same as the information spreading dynamics on subnetwork but with different infection and recovery probabilities, i.e., and , respectively. We assume that the susceptible node in subnetwork adopts an intermittent quarantine strategy to control virus spread. Specifically, the susceptible node in subnetwork will be quarantined for a period if its corresponding node in subnetwork is in the informed state. After period , the susceptible node is reconnected its neighbors, i.e., the intermittent quarantine strategy is induced (as shown in Figure 1(c)). From the above descriptions, we know the state of the node depends on the information spreading and thus further determines the adoption intermittent quarantine strategy for the node . Hence, virus spreading on subnetwork is suppressed. According to the above descriptions, the differences between the information diffusion and virus spreading are listed as follows: (i) different dynamical parameters, i.e., transmission and recovery probabilities, and (ii) an addition intermittent quarantine state is induced in the virus spreading.

We randomly selected one node as the virus seed in subnetwork ; subsequently, we set its counterpart in subnetwork as the information seed. The remaining nodes were set in a susceptible state. The effective information transmission and virus transmission probabilities are expressed as and , respectively. The interacting spreading dynamics terminate once no node is available in the informed or infected state. In Table 1, we illustrate the definitions of the parameters in our model.

3. Results

The results of this study are presented in this section.

3.1. Theoretical Analysis

To analyze the results of this study quantitatively, we used a heterogeneous mean-field approach with time delays. We denote , , and as the fractions of the susceptible, informed, and recovered nodes with degree in subnetwork at time , respectively. Similarly, we denote , , , and to represent the fraction of nodes in the susceptible, infected, recovered, and vaccinated states with degree in subnetwork at time , respectively. Once the parameters above are obtained, the fraction of nodes in each state can be determined using the degree distributions of the two networks, denoted as and , separately. For instance, the fraction of nodes in the susceptible state in subnetwork is .

Inspired by [38, 39], we developed a heterogeneous mean-field approach to investigate virus-information spreading dynamics. For the information spreading dynamics on subnetwork , a susceptible node (e.g., node ) becomes infected via one of two approaches. (i) Node can receive the information from its infected neighbors with probability , where is the probability that node is connected to an infected neighbor. For uncorrelated networks, is expressed aswhere is the average degree of subnetwork . (ii) Node receives the infection once its corresponding node becomes infected, and the infection probability is , where is the probability that a susceptible node is connected to an infected neighbor in subnetwork . is expressed aswhere is the average degree of subnetwork . Combining the two approaches above, the evolution of can be expressed as follows:

Once the susceptible node is infected by the information, it is regarded as infected. Hence, the evolution of is expressed aswhere is the fraction of nodes recovered from the infected state. Therefore, the evolution of is expressed as

For virus spreading on subnetwork , the situation becomes more complex. A susceptible node may be quarantined for a period if its corresponding node becomes the infected state. Mathematically, we use a new state to represent the quarantined nodes. We use to denote the fraction of nodes with degree in the quarantined state. The fraction of nodes becomes a quarantined state with probability . Hence, for a randomly selected susceptible node, it becomes an infected state with probability . The evolution of is expressed as

The evolution of is

The evolution of occurs in one of two processes: (i) the susceptible nodes become infected by their infected neighbors with probability and (ii) the quarantined node experiences more time steps than and is infected by neighbors with probability . Therefore, the evolution of is expressed as

The evolution of is

3.2. Numerical Simulations

We performed extensive Monte Carlo simulations on uncorrelated configuration networks using the method proposed in [48]. Specifically, we generated an uncorrelated configuration network as follows. (i) The network size was determined based on the average degree and degree distributions of the two subnetworks. We set the network sizes of subnetworks and as . The degree distributions of the two networks were set based on a power law as and , respectively, where and are the degree exponents of the two subnetworks, separately: and . In the numerical simulations, we set , and the recovery probability . The multiplex networks are denoted as SFCSF multiplex networks. (ii) Degree sequences were generated for each subnetwork based on the degree distributions, and stubs were assigned for each node. (iii) For each, two stubs were selected randomly, and an edge was built until no stubs remained in the subnetwork. (iv) For each node in the two subnetworks, we performed one-to-one matching randomly.

We first investigated the final information spreading size and the virus spreading size for a specified quarantine value for period , as shown in Figure 2. For a specified value of , increased monotonically with (see Figure 2(a)) because the nodes in subnetwork had a higher probability of accessing the information. In addition, increased with because the spread of the virus promoted information spreading, as shown in Figures 2(a) and 2(b). In general, the virus spreading size increased with for a specified . However, we discovered an optimal information transmission probability for which the virus spreading suppressed significantly, as illustrated in Figures 2(c) and 2(d). Specifically, first decreased with because information was obtained by more nodes, and the quarantine strategy was adopted. However, as increased, more nodes adopted the quarantine strategy for period . If (i.e., the virus spreading terminal time) is sufficiently large, those nodes will be protected from being infected by the virus. Otherwise, the quarantined nodes will participate in virus spreading. Therefore, increases with . It is noteworthy that when is sufficiently large, many nodes are quarantined, and only a fraction of quarantined nodes will participate in virus spreading when their quarantine period is completed. Hence, decreases with .

We further investigated the effects of the quarantine period on the virus-information spreading dynamics, as shown in Figure 3. For the case of (see Figures 3(a)3(d)), i.e., when the quarantine period was relatively small, a similar phenomenon to that shown in Figure 2 was observed. For a large value of , i.e., when (as shown in Figures 3(e)3(h)), a different phenomenon was observed, and the optimal information disappeared. Specifically, increased monotonically with or , whereas decreases monotonically. When was extremely large, the quarantined nodes did not participate in the virus spreading dynamics when they completed their quarantine period. Hence, we conclude that increasing the information transmission probability can contain virus spreading for a long quarantine period.

Next, we investigate whether an optimal quarantine period can significantly suppress virus spreading. Hence, we systematically investigated the spread of virus information for different values of and , as shown in Figure 4. We discovered that decreased with , regardless of the values of and . In terms of virus spreading, decreased monotonically with because the larger value of resulted in fewer nodes participating in virus spreading when they complete their quarantine period.

Finally, we investigated the effects of the average degree of the multiplex networks on virus-information spreading, as shown in Figure 5). It is noteworthy that varying the values of and do not qualitatively affect the phenomena presented in Figure 2. In other words, an optimal information transmission probability exists at which the spread of virus will be significantly suppressed.

4. Conclusions

Controlling the spread of virus on the Internet is vital to cyberspace security. The spread of virus on the Internet triggers information spreading on social networks. Hence, we proposed an interacting virus-information spreading dynamics model for multiplex networks, in which a node receiving information is intermittently quarantined for a specified period. The spreading dynamics were described using a time-delay heterogeneous mean-field approach. By performing Monte Carlo simulations, we investigated the spreading dynamics of SFCSF networks and identified two situations. For relatively small intermittent periods, virus spreading suppressed significantly at the optimal information transmission probability. For large intermittent periods, i.e., the intermittent period was larger than the virus spreading time. Therefore, increasing the information diffusion probability will likely facilitate the control of viruses. Furthermore, we discovered that increasing the intermittent period suppressed virus spreading. Finally, we changed the average degrees of the two subnetworks and discovered that the above mentioned phenomena were not qualitatively affected.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.