# The Dynamical Modeling Analysis of the Spreading of Passive Worms in P2P Networks.

1. IntroductionPeer-to-Peer (P2P) networks are composed of connected computers, which can function as both clients and servers. Each computer in the P2P network can request information and also offer information to other computers. P2P networks provide a very convenient way for people to share the information and gain widespread popularity. The rapid development of P2P networks has attracted the attention of the creators of viruses, worms, and other security threats. A number of viruses and worms are prone to spreading through P2P networks. A worm in P2P networks is a program that replicates itself and propagates through the network from one computer to another to infect healthy computers. According to spreading approach [1, 2], P2P worms can be divided into three categories: passive worms, reactive worms and proactive worms. Passive worms can hide themselves in files and spread as the files are downloaded and executed on healthy computers. Reactive worms propagate by using security vulnerabilities, and this propagation can occur through legitimate network behaviors. Proactive worms can automatically connect and infect neighbor computers by exploiting topological information acquired from infected ones. In this paper, we only focus on passive worms.

The propagation behavior of passive worms is similar to that of biological viruses. There have been many studies on modeling the propagation of passive worms in recent years. In 2006, the epidemiological models in P2P networks were proposed by Thommes and Coates in [3]. The models characterized the P2P virus propagation and pollution dissemination, respectively. In 2007, Zhou et al. [4] proposed a mathematical model for the propagation of passive worms based on the two-factor model [5] and analyzed the dynamics of passive worms. In 2008, Feng et al. [6] proposed three models of passive worm propagation: the SI model, SIS model, and SIR model. The key difference among the three models is the state transitions of peers. In 2009, Wang et al. [7] proposed a passive worm propagation model in unstructured P2P networks. And defense method was also studied based on healthy file dissemination. In 2010, Fan et al. [8] presented a logic matrix approach for modeling the spreading of P2P worms. This proposed model was essentially discretetime deterministic spreading model. In 2012, Rasheed [9] defined an SEI model based on the study of epidemiology and used a P2P simulator to implement this model in P2P networks. In 2014, Chen et al. [10] proposed a four-factor propagation model for passive worms. In this model, four factors were considered: address hiding, configuration diversity, online/offline behaviors, and download duration. In the same year, Yang et al. [11] proposed two propagation models of passive worms: the static model and the dynamic model. In 2015, Feng et al. [12] proposed an analytical model for modeling the propagation of passive worms by adopting epidemiological approaches. In the proposed model, they considered the dynamic characteristics of the P2P network. In 2018, Rguibi et al. [13] presented a propagation model of passive worms in P2P networks. In this proposed model, the hesitation to open a new downloaded file was considered.

However, the above models ignore the influences of network heterogeneity [14,15] on the propagation of passive worms. In a P2P network, each node represents a peer (i.e., a computer), and each edge between two nodes stands for a connection. Different peers have different connections with others, which exhibits an obvious heterogeneity. Additionally, the node degree distribution of unstructured P2P networks follows a power-law distribution, so it is of practical significance to study the propagation of passive worms in heterogeneous networks. On the other hand, a healthy peer that can be infected within a distance of TTL (Time to Live) hops from it [7,16]. Here TTL is the threshold representing the number of hops a search can reach. For convenience, we use d to stand for the number of hops in the following sections. In this paper, we take the above two aspects into account and propose a novel dynamical model to study the dynamics of passive worm propagation in P2P networks. For the proposed model, we first derive the basic reproduction number [R.sub.0]. Additionally, by using epidemiological approaches, it is proved that if [R.sub.0] < 1, the worm-free equilibrium is globally asymptotically stable; when [R.sub.0] >1, there exists a positive equilibrium, which is globally asymptotically stable. And then we verify the rationality of the proposed model by comparing the stochastic simulation with the numerical simulation and investigate the effects of network topology and hop count on passive worm propagation. Finally, some immunization strategies are discussed.

The rest of this paper is organized as follows. In Section 2, we propose a novel network-based SIS passive worm propagation model. In Section 3, the basic reproduction number [R.sub.0] is derived by employing the existence of the positive equilibrium. In Section 4, we investigate the global stabilities of the worm-free equilibrium and positive equilibrium. In Section 5, the numerical simulations are given to interpret the theoretical results, and the stochastic simulations are carried out to verify the rationality of the model. In Section 6, we discuss some immunization strategies and make a conclusion.

2. The Model Formulation

2.1. Modeling Background. To further study the propagation of passive worms, we first need to know the search mechanisms employed by P2P networks. In unstructured P2P networks, such as Gnutella, the most common mechanism is flooding. According to this mechanism, a peer searching for a file sends a query to each of its neighbors. And then each neighbor checks that it has the file matching the query. If it does, it responds to this request and then checks the number of hops the query can reach. The query is forwarded once, and the number of hops is reduced by 1. When hop count is greater than zero, it forwards the query to its neighbors; otherwise, the query is stopped. In particular, when we search for the files in the infected peers, they always respond for all the queries received. Figure 1 illustrates this process visually when hop count d = 2. From Figure 1, it can be seen that there are infected peers in the responding ones. Thus worms can spread easily from infected peers to uninfected ones. In order to model the spreading of passive worms, we need to quantify the search neighborhood.

Now, we adopt the generating function to calculate the number of peers available. The generating function for the degree distribution is denoted by [G.sub.0](x) and defined as

[G.sub.0](x) = [[infinity].summation over (k=0)] p(k) [x.sub.k], (1)

where p(k) is the probability that a randomly chosen node in the network has degree k. When studying the propagation problem on the network, we are also concerned about excess degree distribution. Since the probability that we reach a node of degree k by following an edge is kp(k)/<k>, the probability q(k) of having excess degree k is q(k) = (k + 1) p(k + 1)/<k> [17]. Therefore, the generating function for the excess degree distribution is defined by

[G.sub.1](x) = [[infinity].summation over (k=0)] q(k) [x.sub.k], (2)

In this paper, we are particularly concerned about the number of d-hop neighbors. Firstly, the probability [p.sup.(2).sub.k] that a node has k second neighbors (i.e., 2-hop neighbors) is given by

[p.sup.(2).sub.k] = [[infinity].summation over (m=0)] q(m) [P.sup.(2)] (k | m), (3)

where [P.sup.(2)](k | m) is the probability of having k second neighbors given that there are m first neighbors (i.e., 1-hop neighbors) and p(m) is the degree distribution of the network. From Figure 1, we can see that the number of 2-hop neighbors of a node is equal to the sum of the excess degrees of each of the 1-hop neighbors. Suppose the excess degrees of m first neighbors are [j.sub.1], [j.sub.2], [j.sub.3], ..., [j.sub.m], respectively. They are independent random numbers drawn from the distribution q(k). Then the probability distribution of the sum of m integers has generating function [[G.sub.1](x)].sup.m]. That is,

[[infinity].summation over (m=0)] [P.sup.(2)] (k | m) [x.sup.k] = [[[G.sub.1](x)].sup.m]. (4)

Therefore, substituting the formulas (3) and (4) into the generating function for [P.sup.(2).sub.k], we get

[G.sup.(2)] (x) = [[infinity].summation over (k=0)] [p.sup.(2).sub.k] [x.sup.k] = [G.sub.0] ([G.sub.1] (x)). (5)

Similarly, if there are m neighbors at distance d -1, then the generating function for the number of d-hop neighbors can be represented by

[mathematical expression not reproducible]. (6)

Differentiating (5) and letting x = 1, we get the average number [C.sub.2] of 2-hop neighbors of a node, that is [C.sub.2] = [G'.sub.0](1)[G'.sub.1](1), where [G'.sub.0](1) = (k) and [G'.sub.1](1) = [[summation].sup.[infinity].usb.k=0] kq(k) = (1/{k))(([k.sup.2]) - (k)). Hence the average number of 2-hop neighbors is rewritten as [C.sub.2] = ([k.sup.2]) - (k). Similarly, the average number [C.sub.d] of d-hop neighbors is given by

[C.sub.d] = [([C.sub.2]/[C.sub.1]).sup.d-1] [C.sub.1], (7)

where [C.sub.1] = <k>.

2.2. The Model. Based on modeling idea of epidemic diseases on networks, we propose a novel network-based SIS model to study the dynamics of passive worm propagation. We divide all peers in the network into n groups of sizes [N.sub.k](k = 1,2, ..., n), and each peer in group k has exactly degree k, which is the number of edges emanating from a peer. If the total number of the peers is N, i.e., N = [N.sub.1] + [N.sub.2] + ... + [N.sub.n], the probability that a randomly chosen peer has degree k is p(k) = [N.sub.k]/N, which is referred to as the degree distribution of the network. Additionally, according to the SIS transmission process, each peer can be in one of only two states: susceptible and infected. Let [S.sub.k](t) and [I.sub.k](t) denote the number of susceptible and infected peers of degree k at time t, respectively. In order to reflect the spreading behaviors of passive worms and simplify analysis, we make the following assumptions. (1) The total number of the peers is invariable. (2) The total number of files in all computers is the same. (3) The passive worm always generates an equal number of infected files on each peer. Therefore, the network-based SIS passive worm propagation model can be described by the following differential equations:

[mathematical expression not reproducible], (8)

where

[mathematical expression not reproducible]. (9)

The interpretations of the variables and parameters are shown in Table 1.

Remark 1. For uncorrelated networks, the conditional probability p(l | k) that a node with degree k is connected to a node with degree I is independent of k, and it is proportional to lp(l) [18], i.e., p(l | k) = lp(l)/(k).

3. Basic Reproduction Number

Let [[rho].sub.k](t) = [I.sub.k](t)/[N.sub.k] be the relative density of infected nodes with degree k at time i, and let [theta](t) = [[summation].sup.n.sub.m=1] (mp(m)[[rho].sub.m](t)/(k)) be the probability that any connection points to an infected node. So (8) becomes

d[[rho].sub.k](t)/dt = -[gamma][[rho].sub.k] (t) + [beta](k + H) (1 - [[rho].sub.k] (t)) [theta] (t), (10)

where

[mathematical expression not reproducible]. (11)

The equality (11) can be reduced to

[mathematical expression not reproducible]. (12)

The equilibrium of the system (10) satisfies

[[rho].sub.k] = [beta] (k + H)[theta]/ [gamma] + [beta] (k + H) [theta]. (13)

Substituting (13) into [theta], we get

[mathematical expression not reproducible]. (14)

Obviously, [theta] = 0 is a zero solution of (14), and 0 [less than or equal to] [theta] [less than or equal to] 1. Let

[mathematical expression not reproducible] (15)

By computing, we get

[mathematical expression not reproducible] (16)

Therefore, F([theta]) is a concave function for 0 [less than or equal to] [theta] [less than or equal to] 1 and F(0) = 0,

F(1) = 1 - 1/(k) [n.summation over (m=1)] mp (m) [beta] (m + H)/[gamma] + [beta] (m + H) > 0. (17)

So the sufficient and necessary condition for the existence of a unique positive solution of F([theta]) for 0 < 9 < 1 is obtained as follows:

dF([theta])/d[theta] | [sub.[theta]=0] = 1 - [beta]/[gamma] (([k.sup.2])/(k) + H) < 0. (18)

Then, we derive the basic reproduction number

[R.sub.0] = [beta]/[gamma] (([k.sup.2])/(k) + H). (19)

Remark 2. If d = 1, then H = 0. In this scenario, [R.sub.0] is exactly the basic reproduction number of Pastor-Satorras and Vespignani model in [19].

4. Global Stability of Equilibria

4.1. Global Stability of Worm-Free Equilibrium. It is easy to check that the set

[OMEGA] = {[[rho].sub.1], [[rho].sub.2],..., [[rho].sub.k], 1 [less than or equal to] k [less than or equal to] n, 0 [less than or equal to] [[rho].sub.k] [less than or equal to] 1} (20)

is an invariant set of the system (10). Multiplying the equation of the system (10) by kp(k) and then summing over k, we obtain

[mathematical expression not reproducible]. (21)

Let us consider the Lyapunov function

[mathematical expression not reproducible]. (22)

If [R.sub.0] < 1, then dV/dt [less than or equal to] 0 for all [[rho].sub.k] [greater than or equal to] 0. And V'(t) = 0 if and only if [I.sub.k] = 0. Therefore, the worm-free equilibrium is globally asymptotically stable.

4.2. Global Stability of Positive Equilibrium

Lemma3. Assume that [mathematical expression not reproducible].

Proof. Multiplying the equation of the system (10) by kp(k)/(k) and then summing over k, [theta](t) satisfies

[mathematical expression not reproducible]. (23)

Since [theta](0) = (1/(k)) [[summation].sup.n.sub.k=1] kp(k) [[rho].sub.k] (0) > 0, according to the existence and uniqueness of solutions of differential equations, we know that [theta](t) > 0, for all t > 0.

The equation of the system (10) can be rewritten as

[mathematical expression not reproducible]. (24)

Since [theta](t) > 0, we get

d[[rho].sub.k](t)/dt + [[gamma] + [beta](k + H) [theta] (t)] [[rho].sub.k] (t) > 0. (25)

From the inequality (25), we have

[mathematical expression not reproducible], (26)

for all t > 0.

Additionally, the function 1 - [[rho].sub.k](t) satisfies

[mathematical expression not reproducible]. (27)

Through the similar proof, we have 1 - [[rho].sub.k](t) > 0. Therefore, 0 < [[rho].sub.k](t) < 1 for all t > 0.

Lemma 4. Assume that 0 [less than or equal to] [[rho].sub.k](0) [less than or equal to] 1 satisfies [[summation].sup.n.sub.k=1] kp(k) [[rho].sub.k](0) > 0. If [R.sub.0] > 1, then the solution [[rho].sub.k](t) of the system (10) satisfies [mathematical expression not reproducible].

Proof. By Lemma 3, we know that [theta](t) > 0 for t [greater than or equal to] 0. According to (23), we get

[mathematical expression not reproducible]. (28)

From (28), we have

[mathematical expression not reproducible]. (29)

We will prove that there exists [t.sub.1] > 0. If t [greater than or equal to] [t.sub.1], then

[mathematical expression not reproducible]. (30)

Using reduction to absurdity, we assume that the following holds for all t > 0:

[mathematical expression not reproducible]. (31)

By integrating, for all t > 0, we have

[mathematical expression not reproducible]. (32)

This contradicts [theta](t) < 1. Therefore, there exists [t.sub.1] > 0, and the inequality (30) holds when t = [t.sub.1]. Next, we further prove that the inequality (30) holds when t [greater than or equal to] [t.sub.1]. Supposing that this is not true, then there exists [t.sub.2] > [t.sub.1] and it follows that

[mathematical expression not reproducible]. (33)

Therefore, the inequality (30) holds for 0 < t < [t.sub.2]. If t = [t.sub.2], then

[mathematical expression not reproducible]. (34)

By (29) and (34), we get

[mathematical expression not reproducible]. (35)

According to the definition of derivatives, we know that [theta](t) < [theta]([t.sub.2]) and the following holds:

[mathematical expression not reproducible]. (36)

This contradicts the definition of [t.sub.2]. Thus, the inequality (30) holds when t [greater than or equal to] [t.sub.1], i.e., the following inequality holds for t [greater than or equal to] [t.sub.1]:

[theta] (t) > 1/2[beta] (n + H) ([beta] ([k.sup.2])/(k) + [beta]H - 1). (37)

According to Lemma 3 and continuity of [theta](t), we know that

[mathematical expression not reproducible]. (38)

Let [sigma] = [inf.sub.t[greater than or equal to] 0] [theta](t), and the following holds for t > 0:

[mathematical expression not reproducible]. (39)

For the inequality (39), we get

[mathematical expression not reproducible]. (40)

By the inequality (40), we easily know that [inf.sub.t[greater than or equal to][tau]] [[rho].sub.k](t) > 0 for any [tau] > 0.

Lemma 5. If the solution [[rho].sub.k](t) of the system (10) satisfies [mathematical expression not reproducible], then

[mathematical expression not reproducible], (41)

where [u.sub.k] [greater than or equal to] 0 and [l.sub.k] [greater than or equal to] 0.

Proof. The way to prove it comes from [20]. Since lim [sup.sub.t[right arrow]+[infinity]] [[rho].sub.k](t) [less than or equal to] [u.sub.k], we know that for all [epsilon] > 0, there exists [tau] > 0 and [[rho].sub.k](t) [less than or equal to] [u.sub.k] + [epsilon] for t [greater than or equal to] [tau].

By Lemma 3, we have 1 - [[rho].sub.k](t) > 0. For t [greater than or equal to] r, we obtain

[mathematical expression not reproducible]. (42)

From the inequality (42), we have

[mathematical expression not reproducible], (43)

where [mathematical expression not reproducible]. Therefore, we get

[mathematical expression not reproducible]. (44)

Let [epsilon] [right arrow] 0, and the first inequality in (41) holds.

Additionally, since lim [inf.sub.t[right arrow]+[infinity]] [[rho].sub.k](t) [greater than or equal to] [l.sub.k], for all [eta] > 0, there exists [phi] > 0 such that as t [greater than or equal to] [[rho].sub.k] (t) [greater than or equal to] [l.sub.k] + [eta].

If t [greater than or equal to] [phi], then

[mathematical expression not reproducible]. (45)

By the inequality (45), we obtain

[mathematical expression not reproducible], (46)

where [mathematical expression not reproducible]. Therefore, the following inequality holds:

[mathematical expression not reproducible]. (47)

Let [eta] [right arrow] 0, and we get the second inequality in (41).

Theorem 6. Assume that [mathematical expression not reproducible] is the coordinate component of positive equilibrium of the system (10).

Proof. It is first proven that the following relation holds:

[mathematical expression not reproducible]. (48)

We let [u.sup.(1).sub.k] = 1 and define the sequence

[mathematical expression not reproducible]. (49)

By Lemma 3, we know that lim [sup.sub.t[right arrow]+[infinity]] [[rho].sub.k](t) [less than or equal to] 1 = [u.sup.(1).sub.k] for 1 [less than or equal to] k [less than or equal to] n. According to Lemma 5, we have

[mathematical expression not reproducible]. (50)

It can be easily proved that the sequence [u.sup.(j).sub.k] is decreasing by induction. Therefore, its limit exists and is expressed as [lim.sub.j[right arrow][infinity]] [u.sup.(j).sub.k] = [u.sub.k]. Setting j [right arrow] [infinity] in (49) and (50), we have

[mathematical expression not reproducible]. (51)

Additionally, the following function is considered:

[mathematical expression not reproducible]. (52)

By computing, we get

f(0) = 0, f' (0) = [beta]/[gamma] (([k.sup.2])/(k) + H) - 1 > 0. (53)

If x > 0 and x is sufficiently small, then f(x) > f(0) = 0 according to the definition of derivatives. On the basis of Lemma 5, let [l.sup.(1).sub.k] satisfy the following relation.

[mathematical expression not reproducible]. (54)

And the sequence is defined as

[mathematical expression not reproducible]. (55)

By (54) and Lemma 5, we obtain

[mathematical expression not reproducible]. (56)

On the basis of (52), (54) and (55), the following inequality holds:

[mathematical expression not reproducible]. (57)

According to (55) and (57), we get

[mathematical expression not reproducible]. (58)

If [l.sup.(j+1).sub.k] > [l.sup.(j).sujb.k] for 1 [less than or equal to] k [less than or equal to] n, then

[mathematical expression not reproducible]. (59)

Therefore, the sequence [l.sup.(j).sub.k] is increasing for j [greater than or equal to] 2 and its limit exists. Let [lim.sub.j[right arrow][infinity]] = [l.sub.k]. Setting j [right arrow] [infinity] in (55) and (56), we have

[mathematical expression not reproducible]. (60)

By (51) and (60), we know that [lim.sub.j[right arrow][infinty]] [u.sup.(j).sub.k] = [u.sub.k] and [lim.sub.j[right arrow][infinity]] [l.sup.(j).sub.k] = [l.sub.k] are the coordinate components of positive equilibrium of the system (10). According to the uniqueness of the positive equilibrium of the system (10), we deduce that [u.sub.k] = [l.sub.k] = [[rho].sup.*.sub.k] and

[mathematical expression not reproducible]. (61)

Therefore, [lim.sub.t[right arrow][infinity]](t) = [[rho].sup.*.sub.k]. Through the above proofs, we know that the positive equilibrium is globally asymptotically stable.

5. Simulation Results

In this section, we make a series of simulations to validate the proposed model in this paper. Our simulations are based on a network with power-law distribution p(k) = C[k.sup.-[mu]], where C = 1.92 and [mu] = 2.14. And let N = 5000.

In Figure 2, we examine the validity of the model through the numerical simulations and the stochastic simulations in the case of different d. In Figure 2(a), the parameters are chosen as [beta] = 0.002 and [gamma] = 0.0015. We give the results of the numerical simulation and the stochastic simulation about the model for d = 1. In Figure 2(b), when [beta] = 0.004 and [gamma] = 0.0065, the numerical simulation about the model for d = 2 is carried out to visually show the change of the proportion of infected peers with time. At the same time, the stochastic simulation for d = 2 is carried out to verify the rationality of the model for d = 2. Similarly, when [beta] = 0.003 and [gamma] = 0.0055, Figure 2(c) shows the results of the numerical simulation and the stochastic simulation about the model for d = 3. Ina word, it can be seen from Figure 2 that the fitting results of stochastic simulations and numerical simulations are more consistent, which shows that the model established in this paper is reasonable.

In Figure 3, when the basic reproduction number [R.sub.0] > 1, we examine the impact of d on the evolution of infection prevalence, which is defined as the proportion of infected peers with respect to the total number of the peers. The parameters are chosen as [beta] = 0.0021 and [gamma] = 0.0019. Obviously, the number of hops has a significant effect on the spreading of passive worms. From Figure 3, it can be seen that the propagation speed of passive worms is accelerated with the increase of the value of d. It is easy to understand that the larger the value of d is, the larger the search neighborhood becomes. Thus, the probability that there exist infected peers becomes larger, which accelerates the propagation of passive worms. Hence, the infection becomes more prevalent. From Figure 3, we can also see that the proportion of infected peers can converge to a positive value, which shows that passive worms exist in the network and reach a stable state eventually. This is consistent with the theoretical analysis.

In Figure 4, we show the tendency of passive worm propagation when the basic reproduction number [R.sub.0] = 0.8914 < 1. From Figure 4, it can be seen that the tendency of the spread of passive worms is in decline and reaches a worm-free equilibrium state eventually, which is consistent with the analysis in previous sections. Additionally, we find that the larger the value of d, the faster the descent.

In Figure 5, when the basic reproduction number [R.sub.0] = 3.9685 > 1, we examine the change of the relative densities of infected peers with different degrees in the case of d = 1,2 and 3, respectively. Here, the relative density of infected peers with degree k is defined as the proportion of infected peers with degree k with respect to the total number of the peers with degree k, i.e., the solutions of the system (10). In particular, Figure 5 gives the relative densities of the infected nodes with degree k = 4,20 and 30 when d = 1,2 and 3, respectively. From Figure 5, we find that the nodes with large degrees are easily to be infected, and the relative densities of the infected nodes with small degrees increase obviously with the increase of d. It is easy to understand that with the increase of d, the number of the nodes with small degrees increases more significantly than that of the nodes with large degrees. Therefore, the probability that the nodes with small degrees are infected becomes larger.

6. Discussions and Conclusions

In this paper, after understanding the mechanism of passive worm propagation, we propose a novel network-based SIS model to study the dynamics of its propagation in P2P networks and predict the trend of the propagation. For the presented model, the basic reproduction number [R.sub.0] is derived by employing the existence of the positive equilibrium. And then we prove that if [R.sub.0] < 1, the worm-free equilibrium is globally asymptotically stable. Additionally, if [R.sub.0] > 1, there exists a positive equilibrium, which is globally asymptotically stable. In order to interpret the theoretical results, we carry out the numerical simulations. We find that the simulation results are consistent with the theoretical analysis. On the other hand, we verify the rationality of the model established by comparing the stochastic simulation with the numerical simulation. We find that their fitting results are more consistent. Moreover, we examine the effect of hop count on the spread of passive worms. From simulation results, we can see that if [R.sub.0] > 1, the propagation speed of passive worms is accelerated with the increase of the value of d, which leads to an increase in the number of infected peers; if [R.sub.0] < 1, the number of infected peers decreases rapidly with the increase of the value of d and drops to zero eventually. To control the spread of passive worms, we need to restrict the number of hops a search can reach in addition to changing the network topology.

The appropriate immunization strategy is also very important for controlling the spread of the virus. Usually, there are two kinds of immunization strategies: uniform immunization and targeted immunization [21, 22]. Now, let us discuss their effects on the system (10).

Uniform immunization is a very simple immunization strategy. This strategy is to immune a fraction of the nodes in the network randomly. Let [psi] be the immunization rate and 0 < [psi] < 1. By substituting [beta](1 - [psi]) for [beta], the system (10) becomes

[mathematical expression not reproducible]. (62)

For system (62), we get its basic reproduction number

[[??].sub.0] = (1 - [psi]) [beta]/[gamma] (<[k.sup.2]>/(k) + H) = (1 - [psi]) [R.sub.0]. (63)

In the formula (63), if [psi] = 0, that is, no immunization is performed, [[??].sub.0] = [R.sub.0]. While 0 < [psi] < 1, [[??].sub.0] < [R.sub.0], which means that the immunization strategy is beneficial to control the spread of virus.

Targeted immunization is another effective strategy to control the spread of the virus, and it is to immune the nodes with large degrees. To this end, lower and upper thresholds [[kappa].sub.1] and [[kappa].sub.2] are introduced. Therefore, the immunization rate [[psi].sub.k] is a piecewise function about k and is defined as

[mathematical expression not reproducible], (64)

where 0 < [p.sub.k] < 1 and [bar.[psi]] = [[summation].sub.k] [[psi].sub.k]p(k) is the average immunization rate. And then the system (10) becomes

[mathematical expression not reproducible]. (65)

For system (65), we get its basic reproduction number

[mathematical expression not reproducible]. (66)

Note that

[mathematical expression not reproducible]. (67)

For appropriately small [[kappa].sub.2], if [[phi].sub.1] > 0 and [[phi].sub.2] > 0, then

[[mathematical expression not reproducible]. (68)

By observing (68), we know that the targeted immunization is more effective than the uniform immunization in the case of the same average immunization rate.

On the other hand, we compare the study in this work with those in [11, 12]. All three papers use dynamics to study the spread of passive worms, but the research focus is different. First, the modeling mechanisms of passive worm spread considered in model are different. In [11], the model mainly describes the effects of downloading files, executing files and deleting files on passive worm propagation, while [12] examines the impacts of downloading files, executing files, online and offline on passive worm spread. In comparison, this paper focuses on the effects of network heterogeneity and the number of hops on the spread of passive worms. Secondly, the basic reproduction number is different due to different models. In [11], the factors considered in the model include downloading, execution and deletion. Besides these factors, [12] also considers online and offline. The basic reproduction number in [11] is determined by the rate at which files are downloaded (i.e., [[lambda].sub.d]), the probability of successfully downloading files (i.e., [[lambda].sub.d]), the rate at which files are executed (i.e., [[lambda].sub.e]), the probability of successfully executing files (i.e., [p.sub.e]), the rate at which infected files are deleted (i.e., [[lambda].sub.s]), the probability of successfully deleting infected files (i.e., [p.sub.s]), and the number of uninfected files. If online and offline are not considered, the basic reproduction number in [12] is determined by [[lambda].sub.d], [p.sub.d], [[lambda].sub.e], [p.sub.e], and [[lambda].sub.r] (i.e., recovery rate) and the number of uninfected files. The parameters [[lambda].sub.d], [p.sub.d], [[lambda].sub.e], and [p.sub.e] are constants related to the occurrence of infection, which can cause the susceptible nodes to become the infected ones, while the parameters [[lambda].sub.s], [p.sub.s] and [[lambda].sub.r] are constants related to the occurrence of recovery, which can cause the infected nodes to become the susceptible ones. In this paper, transmission rate (i.e., [beta]) and recovery rate (i.e., [gamma]) are constants. In essence, we use [beta] to reflect the quantity associated with the occurrence of infection, similar to [[lambda].sub.d], [p.sub.d], [[lambda].sub.e], and [p.sub.e] in [11, 12]. The parameter y is used to reflect the quantity associated with the occurrence of recovery, similar to [[lambda].sub.s] and [p.sub.s] in [11] and [[lambda].sub.r] in [12]. Besides, the basic reproduction number in (19) is related to the values of <[k.sup.2]>, <k> and d; that is to say, the network topology and the number of hops can affect the value of [R.sub.0]. However, this information is not reflected in the basic reproduction number proposed in [11, 12]. Therefore, the basic reproduction number in (19) contains those proposed in [11, 12].

In this paper, we study the spread of passive worms in static networks. However, the propagation of passive worms can be affected by other factors, such as users' behaviors, which exhibit the dynamic characteristics. Therefore, we will focus on the spread of passive worms in dynamic networks in future work.

https://doi.org/10.1155/2018/1656907

Data Availability

No real data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant nos. 11331009 and 11601292).

References

[1] E. Mojahedi and M. A. Azgomi, "Modeling the propagation of topology-aware P2P worms considering temporal parameters," Peer-to-Peer Networking and Applications, vol. 8, no. 1, pp. 171-180, 2013.

[2] H. Tang, J. Huang, and W. Wang, "A game based passive worm defense model for P2P networks," ACM SIGAPP Applied Computing Review, vol. 14, no. 1, pp. 20-29, 2014.

[3] R. Thommes and M. Coates, "Epidemiological modelling of peer-to-peer viruses and pollution," in Proceedings of the INFOCOM 2006: 25th IEEE International Conference on Computer Communications, Spain, April 2006.

[4] H. Zhou, Y. Wen, and H. Zhao, "Passive Worm Propagation Modeling and Analysis," in Proceedings of the 2007 International Multi-Conference on Computing in the Global Information Technology (ICCGI'07), pp. 32-32, Guadeloupe, French Caribbean, March 2007.

[5] C. C. Zou, W. Gong, and D. Towsley, "Code red worm propagation modeling and analysis," in Proceedings of the 9th ACM Conference on Computer and Communications Security, pp. 138-147, Washington, DC, USA, November 2002.

[6] C. Feng, Z. Qin, L. Cuthbet, and L. Tokarchuk, "Propagation modeling of passive worms in P2P networks," in Proceedings of the 2008 IEEE International Conference on Cybernetics and Intelligent Systems, CIS 2008, pp. 1027-1031, China, September 2008.

[7] F. Wang, Y. Zhang, and J. Ma, "Defending passive worms in unstructured P2P networks based on healthy file dissemination," Computers & Security, vol. 28, no. 7, pp. 628-636, 2009.

[8] X. Fan and Y. Xiang, "Modelingthe propagation of Peer-to-Peer worms," Future Generation Computer Systems, vol. 26, no. 8, pp. 1433-1443, 2010.

[9] M. F. Rasheed, "Modelling virus propagation in P2P networks," International Journal of Computer Science, no. 9, pp. 580-587, 2012.

[10] T. Chen, X.-S. Zhang, and Y. Wu, "FPM: Four-factors Propagation Model for passive P2P worms," Future Generation Computer Systems, vol. 36, pp. 133-141, 2014.

[11] W. Yang, Y.-P. Gao, Z.-L. Zhu, G.-R. Chang, and Y. Yao, "Modelling, analysis and containment of passive worms in P2P networks," International Journal of Internet Protocol Technology, vol. 8, no. 2-3, pp. 130-142, 2014.

[12] C.-S. Feng, J. Yang, Z.-G. Qin, D. Yuan, and H.-R. Cheng, "Modeling and analysis of passive worm propagation in the P2P file-sharing network," Simulation Modelling Practice and Theory, vol. 51, pp. 87-99, 2015.

[13] M. A. Rguibi and N. Moussa, "The impact of hesitation on passive worm spreading in P2P networks," Peer-to-Peer Networking and Applications, vol. 11, no. 3, pp. 397-408, 2018.

[14] H. Kang, Y. Lou, G. Chen, S. Chu, and X. Fu, "Epidemic spreading and global stability of a new SIS model with delay on heterogeneous networks," Journal of Biological Systems, vol. 23, no. 4, pp. 577-595, 2015.

[15] C.-H. Li, C.-C. Tsai, and S.-Y. Yang, "Analysis of epidemic spreading of an SIRS model in complex heterogeneous networks," Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 4, pp. 1042-1054,2014.

[16] K. Ramachandran and B. Sikdar, "Modeling malware propagation in Gnutella type peer-to-peer networks," in Proceedings of the 20th IEEE International Parallel and Distributed Processing Symposium, IPDPS 2006, Greece, April 2006.

[17] M. E. J. Newman, Networks: An Introduction, Oxford University Press, Oxford, UK, 2010.

[18] M. Catanzaro, M. Boguna, and R. Pastor-Satorras, "Generation of uncorrelated random scale-free networks 027103," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 71, no. 2, 2005.

[19] R. Pastor-Satorras and A. Vespignani, "Epidemic spreading in scale-free networks," Physical Review Letters, vol. 86, no. 14, pp. 3200-3203, 2001.

[20] L. Wang and G.-Z. Dai, "Global stability of virus spreading in complex heterogeneous networks," SIAM Journal on Applied Mathematics, vol. 68, no. 5, pp. 1495-1502, 2008.

[21] X. Fu, M. Small, D. M. Walker, and H. Zhang, "Epidemic dynamics on scale-free networks with piecewise linear infectivity and immunization," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics,vol. 77,no. 3, ArticleID 036113, pp. 420-427, 2008.

[22] Y. Wang and J. Cao, "A note on global stability of the virose equilibrium for network-based computer viruses epidemics," Applied Mathematics and Computation, vol. 244, pp. 726-740, 2014.

Shujuan Zhang, (1,2) Zhen Jin (iD), (1,2,3,4) and Juan Zhang (3)

(1) School of Mechatronic Engineering, North University of China, Taiyuan 030051, Shanxi, China

(2) School of Data Science and Technology, North University of China, Taiyuan 030051, Shanxi, China

(3) Complex Systems Research Center, Shanxi University, Taiyuan 030006, Shanxi, China

(4) Shanxi Key Laboratory of Mathematical Techniques and Big Data Analysis on Disease Control and Prevention, Shanxi University, Taiyuan 030006, Shanxi, China

Correspondence should be addressed to Zhen Jin; jinzhn@263.net

Received 22 May 2018; Accepted 3 September 2018; Published 20 September 2018

Academic Editor: Zhengqiu Zhang

Caption: Figure 1: A schematic diagram of the working process using the flooding mechanism.

Caption: Figure 2: The numerical simulations and the stochastic simulations with different d. (a) d=1; (b) d=2; (c) d=3.

Caption: Figure 3: The impact of d on the infected peers with [beta] = 0.0021 and [gamma] = 0.0019.

Caption: Figure 4: The impact of d on the infected peers with [R.sub.0] = 0.8914.

Caption: Figure 5: The change of the relative densities of infected peers with degree k = 4,20 and 30 for different d when [R.sub.0] = 3.9685. (a) d=1; (b) d=2; (c) d=3.

Table 1: Descriptions of variables and parameters in model (8). Variables and Comments parameters [beta] Transmission rate. [gamma] Recovery rate. d The number of hops. (k) The mean degree. ([k.sup.2]) The second order moment. The probability that a susceptible peer with [[theta].sub.1](t) degree k contacts with infected files which are from infected peers in its neighbors. The probability that a susceptible peer with [[theta].sub.2])(t) degree k contacts with infected files which are from infected peers in the 2-hop neighbors. The probability that a susceptible peer with [[theta].sub.3])(t) degree k contacts with infected files which are from infected peers in the 3-hop neighbors. The probability that a susceptible peer with [[theta].sub.4]) (t) degree k contacts with infected files which are from infected peers in the 4-hop neighbors. The probability that a susceptible peer with [[theta].sub.d]) (t) degree k contacts with infected files which are from infected peers in the tl-hop neighbors.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Zhang, Shujuan; Jin, Zhen; Zhang, Juan |

Publication: | Discrete Dynamics in Nature and Society |

Geographic Code: | 9CHIN |

Date: | Jan 1, 2018 |

Words: | 6475 |

Previous Article: | Distributed Tracking Control for Discrete-Time Multiagent Systems with Novel Markovian Switching Topologies. |

Next Article: | On the Solutions of a System of Third-Order Rational Difference Equations. |

Topics: |