#### Abstract

In a recent paper by J. Ren et al. (2012), a novel computer virus propagation model under the effect of the antivirus ability in a real network is established. The analysis there only partially uncovers the dynamics behaviors of virus spread over the network in the case where around bifurcation is local. In the present paper, by mathematical analysis, it is further shown that, under appropriate parameter values, the model may undergo a global B-T bifurcation, and the curves of saddle-node bifurcation, Hopf bifurcation, and homoclinic bifurcation are obtained to illustrate the qualitative behaviors of virus propagation. On this basis, a collection of policies is recommended to prohibit the virus prevalence. To our knowledge, this is the first time the global bifurcation has been explored for the computer virus propagation. Theoretical results and corresponding suggestions may help us suppress or eliminate virus propagation in the network.

#### 1. Introduction

More and more applications based on computer networks go into our daily life. While bringing convenience to us, computer networks are exposed to various threats. Computer viruses, which are programs developed to attempt to attach themselves to a host and spread to other computers mainly through Internet, can damage network resources. Consequently, understanding the law governing the spread of computer virus is of considerable interest.

Due to the high similarity between computer viruses and biological viruses [1], the classical SIR (Susceptible-Infected-Recovered) computer virus propagation model borrowed from its epidemic counterpart was proposed [24], which is formulated as the following system of differential equations: Here it is assumed that all the computers connected to the network in concern are classified into three categories: susceptible, infected, and recovered computers. Let , , and denote their corresponding numbers at time , respectively. This model involves four positive parameters: denotes the rate at which external computers are connected to the network, denotes the recovery rate of infected computers due to the antivirus ability of the network, denotes the rate at which one computer is removed from the network, and denotes the rate at which, when having connection to one infected computer, one susceptible computer can become infected. For some variants of this model or its extensions, see [519].

The use of antivirus software is regarded as one of the most effective approaches to recovering infected computers [20]. In reality, the ability of an antivirus software is usually proportional to its cost. Due to the limited software cost, the antivirus ability of a network is limited. So, it is natural to consider the following recovery function: where is the recovery rate when the antivirus ability is not fully utilized; . The susceptible individuals are assumed to have the logistic growth with carrying capacity as well as intrinsic growth rate .

Based on the above mentioned work, we establish a novel computer virus propagation model in [19]: Because the first two equations in (3) are independent of , we can consider the following reduced model:

In [19], we carry out the global qualitative analysis of the model and obtain rich dynamical properties. First, we give the threshold value determining whether the virus dies out completely. Second, we study the existence of equilibria and investigate their local and global asymptotic stability. Next, we find that, depending on the antivirus ability, the system may undergo a local backward bifurcation or a Hopf bifurcation. Finally, we prove that, under appropriate conditions, the system may admit bistable states: a stable virus-free equilibrium and a stable virus equilibrium, or two stable virus equilibria. In this case, the initial condition is critical for the eventual steady state of the system.

In consideration of the fact that our former work on the dynamics behaviors of virus spread in model (4) is incomplete. In this paper, we continue the previous work and give detailed theoretical analysis of the presented model so as to supplement its rich dynamics. In fact, by our mathematical analysis, we have not only theoretically proved the existence of global Bogdanov-Takens bifurcation in the model, but also obtained the curves of saddle-node bifurcation, Hopf bifurcation, and homoclinic bifurcation, even for unstable limit cycle and stable virus equilibrium cases and especially found stability switches for the virus equilibria near the Bogdanov-Takens bifurcation point so that the qualitative behaviors of virus spread have been illustrated, all of which are instructive for us to choose appropriate virus-controlling strategy. On this basis, a set of policies is recommended for eradicating viruses spreading across the Internet effectively.

The remaining materials of this paper are organized this way: Section 2 studies the existence of equilibria of model (4); Section 3 examines the global Bogdanov-Takens bifurcation; Section 4 illustrates the qualitative behaviors of the bifurcation. We end the paper with a brief concluding remark in Section 5.

#### 2. Equilibria

To obtain equilibria, model (4) can be written as Equation (5) has a trivial equilibrium and a virus-free equilibrium . Let and if , then (5) admits a unique positive solution , where Obviously, is a virus equilibrium of (5) if and only if From (6) we get a quadratic equation in : where . When and (which is equivalent to ), (6) has the following positive solutions: Accordingly, , are virus equilibria if , .

Let us consider the inequality , which is equivalent to . By solving it, we get Likewise, we can derive if or Notice that is equivalent to . To summarize, we get the following.

Theorem 1 1 (see [19]). Consider model (4). Let If , then and are the only possible equilibria. Moreover, (i)if and , then both and are equilibria;(ii)if and , then is an equilibrium but is not;(iii)if , then does not exist. Furthermore, exists if and only if .

#### 3. Global Bifurcation

In this section, we will study the global bifurcation of the above model. In fact, if the condition (which is equivalent to   +  ), (6) has the following unique positive solution: Accordingly, the virus equilibrium is Solving the inequality , we can obtain To summarize, we get the following.

Theorem 2. Consider model (4). Let If holds, then is a unique virus equilibrium.

Theorem 3. If the virus equilibrium exists, then is a cusp of codimension 2; that is, it is a Bogdanov-Takens singularity.

Proof. To translate to the origin, let , ; then (4) becomes That the Jacobin matrices at the virus equilibrium have the double-zero eigenvalue implies that Let ; then (21) can become the following: where , , , and .
Again by letting , we can obtain the following: where are smooth functions in at least of third order and . It follows from [21, 22] that model (4) admits a Bogdanov-Takens bifurcation.

Next, we will explore the approximating BT bifurcation curves by choosing and as bifurcation parameters. First we fix the parameters , and and let , where and are parameters which vary in a small neighborhood of the origin. Let ; then is translated to and model (4) becomes where , , , and .

Let ; then (24) becomes where and is a smooth function of at least of third order.

Let ; then (26) becomes where , and is a smooth function of at least of third order.

A new time variable is introduced which satisfies . Still write as ; then (27) can be written as follows: where is a smooth function of at least of third order.

Let , , rewrite as , and then (28) becomes where , , , and is a smooth function of at least of third order.

Let , , and , rewrite as , and then (29) can become the following: where , , and is a smooth function of at least of third order.

According to the theorems of [23], we can obtain the following theorem:

Theorem 4. If the virus equilibrium exists, then model (4) admits the following bifurcation behavior:(1)there is a saddle-node bifurcation curve ;(2)there is a Hopf bifurcation curve ;(3)there is a homoclinic curve .

#### 4. Further Discussions and Suggestions

From the above theorem, we have found that, depending on the parameters that vary in a small neighborhood of the original, model (4) undergoes a B-T bifurcation consisting of the curves of saddle-node bifurcation, Hopf bifurcation, and homoclinic bifurcation where the virus exhibits the different spread behaviors. In the following, we will illustrate the qualitative behaviors and the bifurcation diagram as sketched in Figure 1.

Now, Let us make a roundtrip near the Bogdanov-Takens point . When the point lies in region 1 where there are no virus equilibria (and thus no limit cycles are possible) and the corresponding flow is parallelizable. When the point lies on the curve , model (4) has exactly one virus equilibrium , which is saddle-node, and has three separatrices: one of them tends to as to (unstable separatrix) and two are stable separatrices (tend to as to ). Entering from region 2 into region 3 through the component of the fold curve yields two virus equilibria: a saddle and a stable node and has no periodic solution. It is worth noting that when the points pass from one side of the cure to the other side, the number of equilibria of model (4) changes from zero to two. This hints that model (4) undergoes a saddle-node bifurcation of codimension 1, and the surface is called a saddle-node bifurcation surface. Then, the node turns into a focus and loses stability when crossing the Hopf bifurcation boundary . While entering into region 3, there exist exactly two virus equilibria: one saddle point and another stable focus surrounded by one limit cycle, which is unstable. As we move clockwise, it “grows” and approaches the saddle, turning into a homoclinic orbit at HL. For our remainder journey clockwise, no limit cycles must remain. Therefore, there must be global bifurcations “destroying” the cycle somewhere between and . When reaching curve HL, two virus equilibria appear, one saddle point and one stable focus. Two separatrices of the saddle point form a closed separatrix loop surrounding the stable focus. There is no periodic solution and if we trace the homoclinic orbit along the curve HL toward the Bogdanov-Takens point, the loop-like orbit shrinks and disappears. To complete our roundtrip, note that there are no cycles in region 4 located below the curve HL, and a stable node and a saddle, existing for the parameter values in these regions, collide and disappear, consequently, resulting in only one virus equilibrium which has three separatrices, one unstable and two stable at the fold curve . Let us also point out that at the critical equilibrium with a double-zero eigen value has exactly two asymptotic orbits (one tending to the equilibrium for and one approaching it as ).

We know that the virus propagation behaviors mainly depend on the dynamics, specifically, the existence and stability of virus-free equilibrium or virus equilibrium and the limit cycles of model. Based on the above analysis, we can see that the dynamics of B-T bifurcation can broadly fall into three types between different regions and curves: unstable virus equilibrium, no virus equilibrium, and no limit cycle; stable virus equilibrium; stable virus equilibrium and unstable limit cycle. Next, we provide some suggestions aiming at each type to suppress or eliminate the virus spread in the real network.

Case 1. Existence of an unstable virus equilibrium, no equilibrium, and no limit cycle: this implies that any positive orbit of model (4) meets the positive -axis in finite time; that is, there is no infected computer in the network, and therefore, the virus will become extinct. Thus, the virus should be left alone in a network.

Case 2. Existence of a stable virus equilibrium: this implies that the virus cannot be eradicated, but can be suppressed. Thus, virus prevalence, which is defined as the infected computers, independent of any initial state, would attain a certain level. So, it is necessary for us to carry out certain steps to inhibit the virus prevalence to an acceptable level. A proper suggestion is to increase , which corresponds to the increase of the antivirus ability . Say, it is sufficient to increase to the extent where a virus equilibrium becomes unstable so as to eliminate the virus in a real network. This leads us to perceive that, on the one hand, we should understand the structure of currently prevalent viruses, because the new antivirus software developed based on this understanding would possess a powerful capability to detect and clear existing viruses; toward this direction, we give the structure of the virus in Figure 2. On the other hand, we should invest more on developing antivirus software or updating it timely so as to enhance its antivirus ability.

Case 3. Coexistence of a stable virus equilibrium and an unstable limit cycle, that is, a unique unstable limit cycle inside which the positive orbits of model (4) tend to a virus equilibrium as tends to infinity: this implies that virus is persistent inside the cycle, whereas the virus becomes disappeared outside this cycle, which leads us to perceive that, on the one hand, we should enhance the antivirus ability of the network; on the other hand, we should take proper measures to make the initial position lie outside the cycle so as to eliminate the virus.

#### 5. Concluding Remarks

In our previous work, we presented a novel computer virus model describing the virus spread behaviors under the effect of the antivirus ability in a real network and analyzed its dynamics for the case where around bifurcation is local. In order to supplement these partial dynamical behaviors, in the present paper, it is further shown that, under appropriate parameter values, the model may undergo a global B-T bifurcation, and the curves of saddle-node bifurcation, Hopf bifurcation, and homoclinic bifurcation are obtained at the virus equilibrium. To our knowledge, this is the first time the global bifurcation has been explored for the computer virus propagation. Based on the theoretical results, corresponding suggestions may help us take practical measures to suppress or eliminate virus in the network.

Towards this direction, there are some further issues that are yet to be studied. First, the inherent structure of a computer virus represents the resource of its breakout. Therefore, it is crucial to dissect the logical structure of virus and analyze the function of each of the modules and their interrelations so as to uncover the resource of virus burst. Second, it is well known that a computer virus may spread through the Internet, which can be reflected in two-layer structure: router layer and autonomous systems layer. The former processes different extents of the antivirus ability due to the limit cost of antivirus and the latter emerges power-law distribution topology features [24], both of which will be combined to model and analyze the virus propagation behaviors on the Internet in our next work.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The work is supported by the National Natural Science Foundation of China under Grant no. 61304117, the Natural Science Foundation of the Jiangsu Higher Education Institutions of China under Grant no. 13KJB520008, and the doctorate teacher support project of Jiangsu Normal University under Grant no. 12XLR021.