# Cooperative and Competitive Dynamics Model for Information Propagation in Online Social Networks

Authors: Yaming Zhang; Chaosheng Tang; Li Weigang
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

Traditional empirical models of propagation consider individual contagion as an independent process, thus spreading in isolation manner. In this paper, we study how different contagions interact with each other as they spread through the network in order to propose an alternative dynamics model for information propagation. The proposed model is a novel combination of Lotka-Volterra cooperative model and competitive model. It is assumed that the interaction of one message on another is flexible instead of always negative. We prove that the impact of competition depends on the critical speed of the messages. By analyzing the differential equations, one or two stable equilibrium points can be found under certain conditions. Simulation results not only show the correctness of our theoretical analyses but also provide a more attractive conclusion. Different types of messages could coexist in the condition of high critical speed and intense competitive environment, or vice versa. The messages will benefit from the high critical speed when they are both competitive, and adopting a Tit-for-Tat strategy is necessary during the process of information propagation.


---

## Body

## 1. Introduction

With the rapid development of information technology and the widespread use of intelligent communication tools, we progress further into the age of information explosion and have been bombarded with massive amounts of data. For instance, there are more than 368 million users in Sina Weibo and 100 million messages are distributed every day [1], while over 500 million registered users and more than 200 million tweets are posted every signal day in Twitter [2]. As the Twitter and other online social networking (OSN) applications provide great facility by their convenient and efficient way of information dissemination, various studies focus on the issue of many complex phenomena during the process of the spreading of information over networks. Because of the complexity of the system, OSN is generally characterized as being complex networks, in which individuals are represented by nodes, and information interactions between individuals are represented by links between nodes [3]. Since the propagation of information and spreading of epidemic infection have a lot in common, many epidemic models have been used to describe and explore the evolution process of information. For instance, the process of rumor propagation can be viewed as a process of the individuals being infected. After being infected, the ignorants become spreader with a probability, and then the spreaders can become the stiflers, which is similar as the infected persons can recover over a period of time [4].Many researches have been devoted to the propagation dynamics. The significant model for information or rumor spreading was introduced by Daley and Kendall many years ago [5, 6]. Other variants based on Daley-Kendall (DK) model, such as Maki-Thompson (MK) model [7], have been widely used in the past for quantitative study of the mechanism of information dissemination. In DK model, a closed and homogeneously mixed population is classified into three main groups: ignorants, spreaders, and stiflers. Ignorants are susceptible to rumor or other information. Spreaders play an important role in accelerating information diffusion. Stiflers are the individuals who get tired of rumors and refuse to propagate them [8]. In DK model, the rumor is propagated through the population by pairwise contacts between spreaders and others, while it is spread by directed contacts of the spreaders with others in MK model.As the continuous development of research on complex networks, a number of recent studies confirm that the network topology has a marked influence on the whole diffusion process. In the implementation of the MK model, Zanette [9] showed that the simulation results exhibited a transition between regimes of localization and propagation at a finite value of the network randomness. Moreno et al. [10] studied the dynamics of epidemic spreading processes by Monte Carlo simulations and gave a numerical solution of rumor diffusion mean-field equations [11]. Zhou et al. [12] focused on the influence of network structure on rumor propagation and concluded that the total final infected nodes will decrease when the structure changes from random to scale-free network. Some studies described a formulation of the DK model on complex networks in terms of interacting Markov chains. Nekovee et al. [13] derived mean-field equations for the dynamics of rumor spreading on the arbitrary degree correlated networks and drew a conclusion that scale-free social networks are prone to the spreading of rumors. These studies revealed a complex interplay between the network topology and found that the final infected density of population with degree k is related to the network structure. When the structure changes from random networks to scale-free networks, the rumor is more easily spread.Besides the epidemic dynamics models mentioned above, other widely used models describing information propagation are diffusion models in social networks. For instance, the threshold models developed by Granovetter [14] described collective social behavior, and the simplest version of linear threshold model presented by Kempe et al. [15] has taken social network structure into account during the process of information propagation. Those models characterized the information propagation based on a directed graph which originated from social networks, and the inactive nodes could turn into active ones with a certain probability under the influence of a small number of initially activated nodes.Although extensive progress has been made in understanding the dynamics of the information spreading, most studies focused on one type of information or topic through a complex network in the context of dynamics of information spreading. Only a few of unsystematic studies concern the case that more than one type of information coexists and spreads among the population, such as two rumors with different probabilities of acceptance spreading among nodes [16], two types of information spreading among a population of individuals [17], competing ideas in complex social systems [18], and the propagation of competing memes on composite networks [19, 20].Nevertheless, all of these studies explored different contagions competing with each other as they spread over the network but neglected the phenomena of the cooperation as they mutually help each other in spreading. In fact, there are multiple pieces of information about an event spreading through the OSN simultaneously. Moreover, these pieces of information do not spread in isolation manner from all other information currently diffusing in the network. The built-in mechanism underlying these complex interactions between different types of contagions makes a noticeable impact on dynamic spreading. For instance, competing contagions decrease each other’s probability of spreading, while cooperating contagions help each other in being adopted throughout the network [21]. Thus, it is necessary for us to model the interaction between contagions and not just consider each contagion in isolation manner separated from others.As the main contribution of this paper, we propose a cooperative and competitive model for the different contagions of information diffusion in OSN. First of all, the proposed model is based on the population dynamics models to describe the development tendency and interaction between two different types of information. It provides a more realistic description of this process comparing with previous models by Lotka-Volterra and others. Secondly, we use differential equations of stability theory to analyze the proposed model. By means of the approximate analytical and exact numerical solutions of these equations, we examine the steady-state of information propagation in OSN. Finally, the analytic results with the real data collected from Sina Weibo indicate that the proposed model is effective for understanding and explaining some social phenomena presented in the process of online information spreading, especially the vital role of the cooperation and competition in the information propagation.The rest of this paper is organized as follows. In Section2 a formulation of the information spreading model within the framework of biological models is proposed. Section 3 is devoted to analyzing the cooperative and competitive model and explaining the results of the steady-state behavior. In Section 4, the numerical simulation with the data from Sina Weibo demonstrates the investigations of the steady state and dynamics of the model. The last section concludes the paper with the direction of further research.


## 2. Model

In this section, we describe the basic information spreading system and then establish a general co-competition model based on Lotka-Volterra model.


### 2.1. System Description

In OSNs websites, once a user publishes a message on Twitter or other microblog websites, the message will be transmitted to his/her followers. If the message is so interesting that some followers decide to repost it [3], this message can be read by more users and the information could thus reach beyond the network of the original author. We define the user behaviors, such as retweet, review, and mention, as a growth factor. The greater this value is, the faster the message spreads, and more users see the message. Meanwhile, the individuals may also cease spreading a message spontaneously at a certain rate influenced by time decaying effect [22]. This phenomenon can be considered as decay factor of the message.Considering the special time attributes of the information spreading in OSNs, we choose the user growth ratev
(
t
) as dynamical variables instead of the users themselves. In OSNs, the spread of messages is primarily determined by the repost number in time interval. For example, when we describe the spread of a tweet in Twitter, it is usually defined that this message spread more than N times per hour, or the message has been reposted a total number of nearly N times by time t. However, the latter is a cumulative value that cannot reflect a certain number of changes at some points. Furthermore, the message spread process has strong timeliness and variation in OSNs, and each message has a life cycle that spreads slowly at beginning and then breaks out and vanishes finally. The variables v
(
t
) can reflect the life circle curve of the message while the variables N cannot.Assume that there is a message published at timet in Twitter, then users repost this tweet and the message could spread at speed v
(
t
). As mentioned above, we set μ
(
v
) as a function of the speed v
(
t
). However, the speed could not immeasurably increase, and there should be a maximum value V. When the speed v
(
t
) achieves V, its value will decrease. We suppose that the spread of messages obeys logistic growth pattern in OSNs and set
(1)
μ
(
v
)
=
δ
·
(
1
-
v
(
t
)
V
)
.
Here δ is a constant, denoting the growth factor influence on speed. The speed v
(
t
) can be denoted as v
(
t
)
=
d
N
(
t
)
/
d
t, N
(
t
) denotes the total number of people who have heard the message up to time t, and v
(
t
) denote the variation number in time interval [
t
-
1
,
t
]. Meanwhile, the spread of message in OSNs will be influenced by growth factor θ
(
t
) and decay factor λ; we can establish an equation to describe this dynamic process
(2)
d
v
(
t
)
d
t
=
δ
·
(
1
-
v
(
t
)
V
)
·
θ
(
t
)
-
λ
v
(
t
)
.We discuss the main features of (2) as follows.(i) Whenθ
(
t
)
=
0 or the speed v
(
t
) reaches its maximum value V, (2) becomes
(3)
d
v
(
t
)
d
t
=
-
λ
v
(
t
)
.
In (3), v
(
t
) is a reduction function that the speed will decrease with time.(ii) Whenθ
(
t
) is a nonzero constant, we introduce variables γ and κ; set
(4)
γ
=
δ
·
θ
(
t
)
V
+
λ
κ
=
δ
·
θ
(
t
)
then (2) becomes
(5)
d
v
(
t
)
d
t
=
κ
-
γ
v
(
t
)
.
Set v
(
0
)
=
v
0,
(6)
v
(
t
)
=
κ
γ
(
1
-
e
-
γ
t
)
+
v
0
e
-
γ
t
.
Equation (8) denotes the function of v
(
t
) when θ
(
t
) is a nonzero constant.(iii) When the speed remains unchanged,d
v
(
t
)
/
d
t
=
0, (2) becomes
(7)
δ
·
(
1
-
v
(
t
)
V
)
·
θ
(
t
)
-
λ
v
(
t
)
=
0
then we obtain
(8)
v
(
t
)
=
δ
V
θ
(
t
)
δ
θ
(
t
)
+
λ
V
.
Based on the above analysis, we prove that the speed v
(
t
) is a function of the growth factor θ
(
t
). In next section, we further discuss the situation that different types of messages are spreading over networks.


### 2.2. Co-Competition Lotka-Volterra Model

The spreading of abundance of information to which we are exposed through online social networks is a complex sociopsychological process. In the real world, these contagions not only propagate at the same time but also interplay with each other as they spread over the networks. This phenomenon is similar to the mutualism-competition interaction among multiple species. During the spreading of pieces of information, the cooperation happens when the different types of information are at low speed, and the competition happens when they are at high speed. Hence, according to [23], the general cooperation-competition model of n types of information can be presented as follows: d
v
i
/
d
t
=
v
i
f
i
(
v
1
,
v
2
,
…
,
v
n
), where v
i denotes the speed of the ith messages and f
i
(
v
1
,
v
2
,
…
,
v
n
) represents its growth rate. The partial derivatives of the function f
i
(
v
1
,
v
2
,
…
,
v
n
) can be written as ∂
f
i
/
∂
v
j, i
,
j
=
1,2
,
…
,
n. The cooperation-competition interaction is set by (9) by satisfying the condition that for each i, 1
≤
i
≤
n, there is b
i
>
0. Consider
(9)
∂
f
j
∂
v
i
>
0
,
0
<
v
i
<
b
i
,
∂
f
j
∂
v
i
<
0
,
v
i
>
b
i
,
where i
≠
j, i
,
j
=
1,2
,
…
,
n. According to the ecological theory, we can infer that the ith message has a positive effect on other messages as it is at low density (
0
<
v
i
<
b
i
), while it has negative effect on them at high density (
v
i
>
b
i
).We study the common case that two different types of information interact with each other and then make a natural extension of the classical two-species Lotka-Volterra competitive model and Lotka-Volterra cooperative model. We assume that the spread of each message is affected by the internal attractiveness itself and the external pressure among messages simultaneously. Considering two different types of information on the same event spreading around the same time in Twitter or other OSNs, they will compete with each other for the reason that users only choose to believe one of them. In the meantime, if all the messages are so important or interesting, users would be more likely to adopt and share them. It is considered that they mutually help each other in spreading through the network. We term the two types of information as message 1 and message 2, and they will interact with each other when they spread at the same time. The propagation process and the interactions between the two types of information can be governed by the following set of rules.(i)
Each message has a relatively stable natural growth coefficientμ
i and a decay coefficient λ
i, which represent the growth factor and decay factor, respectively, and μ
i
>
λ
i.
(ii)
Since a signal message could not spread indefinitely, there exists a maximum speedV
i during the propagation process.
(iii)
Two messages will compete or cooperate with one another and take negative or positive effect, represented byσ
i
j.Supported by the rules (i) and (ii), we have the following model:(10)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
1
-
v
1
V
1
)
-
λ
1
v
1
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
1
-
v
2
V
2
)
-
λ
2
v
2
.Hereμ
i denotes the growth coefficient of the ith message and μ
i
∈
[
0,1
]; λ
i denotes the decay coefficient of the ith message and λ
i
∈
[
0,1
], i
=
1,2.According to the rule (iii), when the messages 1 and 2 coexist in the network, the interactive influence among the two messages produces positive effect or negative effect as they cooperate or compete with each other. Based on the analysis above, (10) can be written as
(11)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
1
+
σ
21
β
V
2
-
v
1
V
1
-
σ
21
|
v
2
-
β
|
*
V
2
)
-
λ
1
v
1
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
1
+
σ
12
α
V
1
-
v
2
V
2
-
σ
12
|
v
1
-
α
|
*
V
1
)
-
λ
2
v
2
.In order to facilitate the analysis, seta
=
1
-
λ
1
/
μ
1, b
=
1
-
λ
2
/
μ
2; (11) can be written as
(12)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
a
+
σ
21
β
V
2
-
v
1
V
1
-
σ
21
|
v
2
-
β
|
*
V
2
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
b
+
σ
12
α
V
1
-
v
2
V
2
-
σ
12
|
v
1
-
α
|
*
V
1
)
.Equation (12) can be considered as a competition and cooperation multisystem. Let v
2
=
0 in the first equation of the system (12); we have d
v
1
/
d
t
=
μ
1
v
1
(
a
V
1
-
v
1
)
/
V
1, then the parameter a
V
1 is the speed of the message 1 when it disseminates independently form message 2. Let C
1
(
v
1
,
v
2
)
=
a
+
σ
21
(
β
/
V
2
)
-
v
1
/
V
1
-
σ
21
(
|
v
2
-
β
|
*
/
V
2
); the function |
v
2
-
β
|
* can be defined as the absolute function |
v
2
-
β
| while the function |
v
2
-
β
| is smoothed in a very small neighborhood of its vertex (
0
,
β
). Thus, the isocline L
1
:
C
1
(
v
1
,
v
2
)
=
0 is smooth as shown in Figure 1. Then the function f
1
(
v
1
,
v
2
)
=
μ
1
v
1
C
1
(
v
1
,
v
2
) and satisfies
(13)
∂
f
1
∂
v
2
=
σ
21
V
2
μ
1
v
1
,
0
<
v
2
<
β
,
∂
f
1
∂
v
2
=
-
σ
21
V
2
μ
1
v
1
,
v
2
>
β
.Figure 1
The smoothed parts of the isoclinesL
1 and L
2.As shown in (13), the positive or negative effect acting on message 1 depends on the speed of message 2. When message 2 travels slowly (
v
2
<
β
), it takes positive effect on message 1; otherwise, it produces negative effect on message 1 at high speed (
v
2
>
β
). The parameter β signifies the critical speed of message 2 near which the effect of message 2 acting on message 1 would turn positive effect into negative effect. Furthermore, according to the Lotka-Volterra models, the parameter σ
21 represents both the cooperation (at low speed) and competition (at high speed) level of message 2 to message 1. Hence, the low speed of message 2, 0
<
v
2
<
β, represents the scope of cooperation of message 2 with message 1. A similar conclusion can be obtained that the low speed of message 1, 0
<
v
1
<
α, represents the region of net cooperation message 1 with message 2. The isocline L
2
:
C
2
(
v
1
,
v
2
)
=
b
+
σ
12
(
α
/
V
1
)
-
v
2
/
V
2
-
σ
12
(
|
v
1
-
α
|
*
/
V
1
)
=
0 is smooth as shown in Figure 1.Since the system (12) can exhibit the features of the Lotka-Volterra cooperative model and competitive model, it can be named co-competition Lotka-Volterra model. Obviously, when the messages are at low speed, for instance, v
1
<
α and v
2
<
β, the system (12) turns into the Lotka-Volterra cooperative system
(14)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
a
+
σ
21
v
2
V
2
-
v
1
V
1
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
b
+
σ
12
v
1
V
1
-
v
2
V
2
)
.Conversely, when the messages are at high speed, for instance,v
1
>
α and v
2
>
β, the system (12) turns into the Lotka-Volterra competitive model
(15)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
a
+
2
σ
21
β
V
2
-
v
1
V
1
-
σ
21
v
2
V
2
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
b
+
2
σ
12
α
V
1
-
v
2
V
2
-
σ
12
v
1
V
1
)
.By systematic comparison of the system (12) with the Lotka-Volterra competitive system (15), an important conclusion is reached that the cooperation at low speed promotes the coexistence level of the two messages in co-opetition Lotka-Volterra model. When message 2 is at a low speed (
v
2
<
β
), it follows the first equation of system (15) that the effect of message 2 acting on message 1 is negative (
-
σ
21
v
2
/
V
2
<
0
) while it follows the first equation of system (12) that the effect of message 2 acting on message 1 is positive (
σ
21
v
2
/
V
2
>
0
), which originates from the cooperation of message 2 at low speed. As a consequence, message 1 has survived for a long time when message 2 is at low speed. When message 2 is at high speed (
v
2
>
β
), the negative effect produced by message 2 acting on message 1 in system (15) is (
-
σ
21
v
2
/
V
2
) while that in system (12) is -
σ
21
(
v
2
-
β
)
/
V
2. Since -
σ
21
(
v
2
-
β
)
/
V
2
>
-
σ
21
v
2
/
V
2, the negative effect is lessened for the cooperation among messages in system (12). Similarly, the cooperation in system (12) promotes the existence of message 1. In a word, the cooperation in system (12) facilitates the coexistence of both the two messages.


## 2.1. System Description

In OSNs websites, once a user publishes a message on Twitter or other microblog websites, the message will be transmitted to his/her followers. If the message is so interesting that some followers decide to repost it [3], this message can be read by more users and the information could thus reach beyond the network of the original author. We define the user behaviors, such as retweet, review, and mention, as a growth factor. The greater this value is, the faster the message spreads, and more users see the message. Meanwhile, the individuals may also cease spreading a message spontaneously at a certain rate influenced by time decaying effect [22]. This phenomenon can be considered as decay factor of the message.Considering the special time attributes of the information spreading in OSNs, we choose the user growth ratev
(
t
) as dynamical variables instead of the users themselves. In OSNs, the spread of messages is primarily determined by the repost number in time interval. For example, when we describe the spread of a tweet in Twitter, it is usually defined that this message spread more than N times per hour, or the message has been reposted a total number of nearly N times by time t. However, the latter is a cumulative value that cannot reflect a certain number of changes at some points. Furthermore, the message spread process has strong timeliness and variation in OSNs, and each message has a life cycle that spreads slowly at beginning and then breaks out and vanishes finally. The variables v
(
t
) can reflect the life circle curve of the message while the variables N cannot.Assume that there is a message published at timet in Twitter, then users repost this tweet and the message could spread at speed v
(
t
). As mentioned above, we set μ
(
v
) as a function of the speed v
(
t
). However, the speed could not immeasurably increase, and there should be a maximum value V. When the speed v
(
t
) achieves V, its value will decrease. We suppose that the spread of messages obeys logistic growth pattern in OSNs and set
(1)
μ
(
v
)
=
δ
·
(
1
-
v
(
t
)
V
)
.
Here δ is a constant, denoting the growth factor influence on speed. The speed v
(
t
) can be denoted as v
(
t
)
=
d
N
(
t
)
/
d
t, N
(
t
) denotes the total number of people who have heard the message up to time t, and v
(
t
) denote the variation number in time interval [
t
-
1
,
t
]. Meanwhile, the spread of message in OSNs will be influenced by growth factor θ
(
t
) and decay factor λ; we can establish an equation to describe this dynamic process
(2)
d
v
(
t
)
d
t
=
δ
·
(
1
-
v
(
t
)
V
)
·
θ
(
t
)
-
λ
v
(
t
)
.We discuss the main features of (2) as follows.(i) Whenθ
(
t
)
=
0 or the speed v
(
t
) reaches its maximum value V, (2) becomes
(3)
d
v
(
t
)
d
t
=
-
λ
v
(
t
)
.
In (3), v
(
t
) is a reduction function that the speed will decrease with time.(ii) Whenθ
(
t
) is a nonzero constant, we introduce variables γ and κ; set
(4)
γ
=
δ
·
θ
(
t
)
V
+
λ
κ
=
δ
·
θ
(
t
)
then (2) becomes
(5)
d
v
(
t
)
d
t
=
κ
-
γ
v
(
t
)
.
Set v
(
0
)
=
v
0,
(6)
v
(
t
)
=
κ
γ
(
1
-
e
-
γ
t
)
+
v
0
e
-
γ
t
.
Equation (8) denotes the function of v
(
t
) when θ
(
t
) is a nonzero constant.(iii) When the speed remains unchanged,d
v
(
t
)
/
d
t
=
0, (2) becomes
(7)
δ
·
(
1
-
v
(
t
)
V
)
·
θ
(
t
)
-
λ
v
(
t
)
=
0
then we obtain
(8)
v
(
t
)
=
δ
V
θ
(
t
)
δ
θ
(
t
)
+
λ
V
.
Based on the above analysis, we prove that the speed v
(
t
) is a function of the growth factor θ
(
t
). In next section, we further discuss the situation that different types of messages are spreading over networks.


## 2.2. Co-Competition Lotka-Volterra Model

The spreading of abundance of information to which we are exposed through online social networks is a complex sociopsychological process. In the real world, these contagions not only propagate at the same time but also interplay with each other as they spread over the networks. This phenomenon is similar to the mutualism-competition interaction among multiple species. During the spreading of pieces of information, the cooperation happens when the different types of information are at low speed, and the competition happens when they are at high speed. Hence, according to [23], the general cooperation-competition model of n types of information can be presented as follows: d
v
i
/
d
t
=
v
i
f
i
(
v
1
,
v
2
,
…
,
v
n
), where v
i denotes the speed of the ith messages and f
i
(
v
1
,
v
2
,
…
,
v
n
) represents its growth rate. The partial derivatives of the function f
i
(
v
1
,
v
2
,
…
,
v
n
) can be written as ∂
f
i
/
∂
v
j, i
,
j
=
1,2
,
…
,
n. The cooperation-competition interaction is set by (9) by satisfying the condition that for each i, 1
≤
i
≤
n, there is b
i
>
0. Consider
(9)
∂
f
j
∂
v
i
>
0
,
0
<
v
i
<
b
i
,
∂
f
j
∂
v
i
<
0
,
v
i
>
b
i
,
where i
≠
j, i
,
j
=
1,2
,
…
,
n. According to the ecological theory, we can infer that the ith message has a positive effect on other messages as it is at low density (
0
<
v
i
<
b
i
), while it has negative effect on them at high density (
v
i
>
b
i
).We study the common case that two different types of information interact with each other and then make a natural extension of the classical two-species Lotka-Volterra competitive model and Lotka-Volterra cooperative model. We assume that the spread of each message is affected by the internal attractiveness itself and the external pressure among messages simultaneously. Considering two different types of information on the same event spreading around the same time in Twitter or other OSNs, they will compete with each other for the reason that users only choose to believe one of them. In the meantime, if all the messages are so important or interesting, users would be more likely to adopt and share them. It is considered that they mutually help each other in spreading through the network. We term the two types of information as message 1 and message 2, and they will interact with each other when they spread at the same time. The propagation process and the interactions between the two types of information can be governed by the following set of rules.(i)
Each message has a relatively stable natural growth coefficientμ
i and a decay coefficient λ
i, which represent the growth factor and decay factor, respectively, and μ
i
>
λ
i.
(ii)
Since a signal message could not spread indefinitely, there exists a maximum speedV
i during the propagation process.
(iii)
Two messages will compete or cooperate with one another and take negative or positive effect, represented byσ
i
j.Supported by the rules (i) and (ii), we have the following model:(10)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
1
-
v
1
V
1
)
-
λ
1
v
1
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
1
-
v
2
V
2
)
-
λ
2
v
2
.Hereμ
i denotes the growth coefficient of the ith message and μ
i
∈
[
0,1
]; λ
i denotes the decay coefficient of the ith message and λ
i
∈
[
0,1
], i
=
1,2.According to the rule (iii), when the messages 1 and 2 coexist in the network, the interactive influence among the two messages produces positive effect or negative effect as they cooperate or compete with each other. Based on the analysis above, (10) can be written as
(11)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
1
+
σ
21
β
V
2
-
v
1
V
1
-
σ
21
|
v
2
-
β
|
*
V
2
)
-
λ
1
v
1
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
1
+
σ
12
α
V
1
-
v
2
V
2
-
σ
12
|
v
1
-
α
|
*
V
1
)
-
λ
2
v
2
.In order to facilitate the analysis, seta
=
1
-
λ
1
/
μ
1, b
=
1
-
λ
2
/
μ
2; (11) can be written as
(12)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
a
+
σ
21
β
V
2
-
v
1
V
1
-
σ
21
|
v
2
-
β
|
*
V
2
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
b
+
σ
12
α
V
1
-
v
2
V
2
-
σ
12
|
v
1
-
α
|
*
V
1
)
.Equation (12) can be considered as a competition and cooperation multisystem. Let v
2
=
0 in the first equation of the system (12); we have d
v
1
/
d
t
=
μ
1
v
1
(
a
V
1
-
v
1
)
/
V
1, then the parameter a
V
1 is the speed of the message 1 when it disseminates independently form message 2. Let C
1
(
v
1
,
v
2
)
=
a
+
σ
21
(
β
/
V
2
)
-
v
1
/
V
1
-
σ
21
(
|
v
2
-
β
|
*
/
V
2
); the function |
v
2
-
β
|
* can be defined as the absolute function |
v
2
-
β
| while the function |
v
2
-
β
| is smoothed in a very small neighborhood of its vertex (
0
,
β
). Thus, the isocline L
1
:
C
1
(
v
1
,
v
2
)
=
0 is smooth as shown in Figure 1. Then the function f
1
(
v
1
,
v
2
)
=
μ
1
v
1
C
1
(
v
1
,
v
2
) and satisfies
(13)
∂
f
1
∂
v
2
=
σ
21
V
2
μ
1
v
1
,
0
<
v
2
<
β
,
∂
f
1
∂
v
2
=
-
σ
21
V
2
μ
1
v
1
,
v
2
>
β
.Figure 1
The smoothed parts of the isoclinesL
1 and L
2.As shown in (13), the positive or negative effect acting on message 1 depends on the speed of message 2. When message 2 travels slowly (
v
2
<
β
), it takes positive effect on message 1; otherwise, it produces negative effect on message 1 at high speed (
v
2
>
β
). The parameter β signifies the critical speed of message 2 near which the effect of message 2 acting on message 1 would turn positive effect into negative effect. Furthermore, according to the Lotka-Volterra models, the parameter σ
21 represents both the cooperation (at low speed) and competition (at high speed) level of message 2 to message 1. Hence, the low speed of message 2, 0
<
v
2
<
β, represents the scope of cooperation of message 2 with message 1. A similar conclusion can be obtained that the low speed of message 1, 0
<
v
1
<
α, represents the region of net cooperation message 1 with message 2. The isocline L
2
:
C
2
(
v
1
,
v
2
)
=
b
+
σ
12
(
α
/
V
1
)
-
v
2
/
V
2
-
σ
12
(
|
v
1
-
α
|
*
/
V
1
)
=
0 is smooth as shown in Figure 1.Since the system (12) can exhibit the features of the Lotka-Volterra cooperative model and competitive model, it can be named co-competition Lotka-Volterra model. Obviously, when the messages are at low speed, for instance, v
1
<
α and v
2
<
β, the system (12) turns into the Lotka-Volterra cooperative system
(14)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
a
+
σ
21
v
2
V
2
-
v
1
V
1
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
b
+
σ
12
v
1
V
1
-
v
2
V
2
)
.Conversely, when the messages are at high speed, for instance,v
1
>
α and v
2
>
β, the system (12) turns into the Lotka-Volterra competitive model
(15)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
a
+
2
σ
21
β
V
2
-
v
1
V
1
-
σ
21
v
2
V
2
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
b
+
2
σ
12
α
V
1
-
v
2
V
2
-
σ
12
v
1
V
1
)
.By systematic comparison of the system (12) with the Lotka-Volterra competitive system (15), an important conclusion is reached that the cooperation at low speed promotes the coexistence level of the two messages in co-opetition Lotka-Volterra model. When message 2 is at a low speed (
v
2
<
β
), it follows the first equation of system (15) that the effect of message 2 acting on message 1 is negative (
-
σ
21
v
2
/
V
2
<
0
) while it follows the first equation of system (12) that the effect of message 2 acting on message 1 is positive (
σ
21
v
2
/
V
2
>
0
), which originates from the cooperation of message 2 at low speed. As a consequence, message 1 has survived for a long time when message 2 is at low speed. When message 2 is at high speed (
v
2
>
β
), the negative effect produced by message 2 acting on message 1 in system (15) is (
-
σ
21
v
2
/
V
2
) while that in system (12) is -
σ
21
(
v
2
-
β
)
/
V
2. Since -
σ
21
(
v
2
-
β
)
/
V
2
>
-
σ
21
v
2
/
V
2, the negative effect is lessened for the cooperation among messages in system (12). Similarly, the cooperation in system (12) promotes the existence of message 1. In a word, the cooperation in system (12) facilitates the coexistence of both the two messages.


## 3. Model Analysis

In this section, we discuss the intersection of two isoclinesL
1 and L
2 by a symmetrical hypothesis and then analyze the dynamics of the model and the global stability of the system.


### 3.1. Basic Analysis

In order to show the existence of intersection of two isoclinesL
1 and L
2, we assume that the positive or negative effects produced by messages 1 and 2 are almost the same as each other. In the condition, we can set hypothesis as follows: a
=
b
=
z, σ
21
=
σ
12
=
σ, V
1
=
V
2
=
V, and α
=
β
=
ζ, then the values of μ
1 and μ
2 can reflect the difference of the growth factors among the two messages. For |
v
1
-
α
| and |
v
2
-
β
| are smoothed in a very small neighborhood of their vertexes (
0
,
α
) and (
0
,
β
), respectively, the system (12) can be written as
(16)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
z
+
σ
ζ
V
-
v
1
V
-
σ
ζ
V
|
v
2
ζ
-
1
|
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
z
+
σ
ζ
V
-
v
2
V
-
σ
ζ
V
|
v
1
ζ
-
1
|
)
;
set x
1
=
v
1
/
ζ and x
2
=
v
2
/
ζ; then
(17)
d
x
1
d
t
=
μ
1
v
1
(
z
+
σ
ζ
V
-
ζ
x
1
V
-
σ
ζ
V
|
x
2
-
1
|
)
,
d
x
2
d
t
=
μ
2
v
2
(
z
+
σ
ζ
V
-
ζ
x
2
V
-
σ
ζ
V
|
x
1
-
1
|
)
.
In int
⁡
R
+
2, the zero growth isoclines of x
1 and x
2 are as follows:
(18)
x
1
=
-
σ
|
x
2
-
1
|
+
C
x
2
=
-
σ
|
x
1
-
1
|
+
C
,
where C
=
V
z
/
ζ
+
σ.According to [24], the isoclines include three parts as shown in Figure 2: positive part (
+
), neutral part (0), and negative part (
-
). The positive part, neutral part, or negative part is the part of isoclines when the slope (for message 1, the slope is d
v
1
(
t
)
/
d
v
2
(
t
); for message 2, the slope is d
v
2
(
t
)
/
d
v
1
(
t
)) is positive part (
+
), neutral part (0), and negative part (
-
). The two isoclines are symmetry around the line x
1
=
x
2. Then, we only need to discuss
(19)
x
1
=
x
2
,
x
2
=
-
σ
|
x
1
-
1
|
+
C
.
Obviously, (19) always have solutions; we can also obverse the following aspects.(i)
There are two different roots(
σ
-
C
)
/
(
σ
-
1
) and (
σ
+
C
)
/
(
σ
+
1
), when σ
2
≠
0,1 and C
≠
1.
(ii)
There is double rootx
=
1 or x
=
C, when σ
2
≠
1 and C
=
1, or σ
=
0.
(iii)
There is root(
σ
-
C
)
/
(
σ
-
1
) or (
σ
+
C
)
/
(
σ
+
1
), when σ
2
=
1 and C
≠
1.The parabolic zero growth isoclines of message 1 (a) and message 2 (b).
(a)
(b)Based on the above analysis, we can conclude that trajectories in phase space always intersect isoclines either in the horizontal or in the vertical direction. In the next section, we make further analysis about the stability of intersection.


### 3.2. Equilibrium

In order to study the final results of interactive effect between two different types of messages, we need to conduct a comprehensive analysis for the stability of the equilibrium points of simultaneous differential equations.Theorem 1.
System (12) admits no periodic orbit.Proof.
Let(20)
B
(
v
1
,
v
2
)
=
1
v
1
v
2
(
v
1
,
v
2
)
∈
int
⁡
R
+
2
then we have
(21)
∂
∂
v
1
[
μ
1
(
a
+
σ
21
β
V
2
-
v
1
V
1
-
σ
21
|
v
2
-
β
|
*
V
2
)
]
+
∂
∂
v
2
[
μ
2
(
b
+
σ
12
α
V
1
-
v
2
V
2
-
σ
12
|
v
1
-
α
|
*
V
1
)
]
=
-
μ
1
v
2
V
1
-
μ
2
v
1
V
2
<
0
.According to the Bendixson-Dulac theorem [25], system (12) admits no periodic orbit.We study the evolution ofv
1
(
t
) and v
2
(
t
) as time t approaches infinity. As shown in Figure 1, there are three equilibrium points in system (12): O
1
(
0,0
), O
2
(
a
V
1
,
0
), and O
3
(
0
,
b
V
2
). To determine the stability of the equilibrium points, let
(22)
M
=
[
f
v
1
f
v
2
g
v
1
g
v
2
]
;
we have
(23)
p
=
-
(
f
v
1
+
g
v
2
)
|
p
i
,
q
=
det
⁡
M
|
p
i
.
Then, the stability conditions of equilibrium point O
i are p
>
0 and q
>
0. The Jacobin matrix of the system (12) at equilibrium state is as follows.(i) For the equilibrium pointO
1
(
0,0
), we get
(24)
M
(
O
1
)
=
[
a
μ
1
0
0
b
μ
2
]
p
1
=
-
(
a
μ
1
+
b
μ
2
)
,
q
1
=
a
b
μ
1
μ
2
.
Based on the assumption that a
=
1
-
λ
1
/
μ
1
>
0 and b
=
1
-
λ
2
/
μ
2
>
0, we have p
1
<
0 and q
1
>
0; the point O
1
(
0,0
) is an unstable node.Since our aim is to show the coexistence of competition and cooperation by the introduction ofα and β, we suppose that the condition (25) is always satisfied in the rest of the paper. Consider
(25)
b
+
σ
12
(
2
α
V
1
-
a
)
>
0
,
a
+
σ
21
(
2
β
V
2
-
b
)
>
0
.(ii) For the equilibrium pointO
2
(
a
V
1
,
0
), we get
(26)
M
(
O
2
)
=
[
-
a
μ
1
a
V
1
μ
1
σ
21
V
2
0
μ
2
[
b
+
(
2
α
V
1
-
a
)
σ
12
]
]
,
p
2
=
a
μ
1
-
μ
2
[
σ
12
(
2
α
V
1
-
a
)
+
b
]
,
q
2
=
a
μ
1
μ
2
(
-
2
α
σ
12
+
a
σ
12
V
1
-
b
V
1
)
V
1
.
The stability condition of equilibrium point O
2
(
a
V
1
,
0
) is σ
12
>
b
V
1
/
(
a
V
1
-
2
α
).Conditionσ
12
>
b
V
1
/
(
a
V
1
-
2
α
) means that message 1 is more competitive than message 2. It implies that the persons who should have adopted message 2 choose to believe message 1. Consequently, message 2 will eventually become extinct in the process of competition, while message 1 can continue to diffuse and reach the maximum speed V
1.(iii) For the equilibrium pointO
3
(
0
,
b
V
2
), we get
(27)
M
(
O
3
)
=
[
μ
1
[
a
+
(
2
β
V
2
-
b
)
σ
21
]
0
b
V
2
μ
2
σ
12
V
1
-
b
μ
2
]
,
p
3
=
b
μ
2
-
μ
1
[
a
+
σ
21
(
2
β
V
2
-
b
)
]
,
q
3
=
b
μ
1
μ
2
(
b
σ
21
V
2
-
a
V
2
-
2
β
σ
21
)
V
2
.
The stability condition of equilibrium point O
3
(
0
,
b
V
2
) is σ
21
>
a
V
2
/
(
b
V
2
-
2
β
).Conditionσ
21
>
a
V
2
/
(
b
V
2
-
2
β
) shows that message 2 is more competitive than message 1. This just means that people accepted message 2 even though they should have believed message 1. Therefore, during the process of competition, message 2 can survive and reach a maximum speed V
2, while message 1 will eventually become extinct finally.


### 3.3. Dynamics

Assume that the conditions in (25) are satisfied; the two messages in system (12) can coexist. The isocline L
1 of message 1 can be subdivided into L
11 and L
12 as follows:
(28)
L
11
:
a
+
σ
21
v
2
V
2
-
v
1
V
1
,
0
<
v
2
<
β
,
L
12
:
a
+
2
σ
21
β
V
2
-
v
1
V
1
-
σ
21
v
2
V
2
,
v
2
>
β
,
where L
1 is smoothed in a small neighborhood of the vertex D
1
(
a
V
1
+
β
σ
21
V
1
/
V
2
,
β
) as shown in Figure 1. Similarly, we divide the isocline L
2 of message 2 into L
21 and L
22 as follows:
(29)
L
21
:
b
+
σ
12
v
1
V
1
-
v
2
V
2
,
0
<
v
1
<
α
,
L
22
:
b
+
2
σ
12
α
V
1
-
v
2
V
2
-
σ
12
v
1
V
1
,
v
1
>
α
,
where L
2 is smoothed in a small neighborhood of the vertex D
2
(
α
,
b
V
2
+
α
σ
12
V
2
/
V
1
) as shown in Figure 1.While the line segmentsL
11 and L
21 are restricted to the intervals (
a
V
1
,
a
V
1
+
β
σ
21
V
1
/
V
2
) and (
0
,
α
), respectively, they are impossible to intersect in the interval α
≤
a
V
1. Then there is a maximum of three isolated positive equilibrium points in system (12). As shown in Figure 1, let K
1
(
k
11
,
k
12
) represent the intersection of L
11 and L
22, let K
2
(
k
21
,
k
22
) represent the intersection of L
21 and L
12, and let K
3
(
k
31
,
k
32
) represent the intersection of L
12 and L
22. Then we have
(30)
k
11
=
V
1
(
a
+
b
σ
21
)
+
2
α
σ
12
σ
21
σ
12
σ
21
+
1
,
k
12
=
V
2
(
V
1
(
b
-
a
σ
12
)
+
2
α
σ
12
)
(
σ
12
σ
21
+
1
)
V
1
,
k
21
=
V
1
(
V
2
(
a
-
b
σ
21
)
+
2
β
σ
21
)
(
σ
12
σ
21
+
1
)
V
2
,
k
22
=
V
2
(
a
σ
12
+
b
)
+
2
β
σ
12
σ
21
σ
12
σ
21
+
1
,
k
31
=
V
1
(
V
2
(
b
σ
21
-
a
)
-
2
β
σ
21
)
+
2
α
σ
12
σ
21
V
2
(
σ
12
σ
21
-
1
)
V
2
,
k
32
=
V
1
(
V
2
(
a
σ
12
-
b
)
+
2
β
σ
12
σ
21
)
-
2
α
σ
12
V
2
(
σ
12
σ
21
-
1
)
V
1
.Suppose thatG
1 and G
2 represent the slope of the lines L
12 and L
22, respectively. Suppose that G denotes the slope of the line D
1
D
2. Then we have
(31)
G
1
=
-
V
2
V
1
σ
21
<
0
,
G
2
=
-
V
2
σ
12
V
1
<
0
,
G
=
V
2
(
b
V
2
V
1
+
α
σ
12
V
2
-
β
V
1
)
V
1
(
a
V
1
V
2
-
α
V
2
+
β
σ
21
V
1
)
<
0
.Dynamics of system (12) can be displayed in five cases, while we discuss the situation of |
G
1
|
<
|
G
2
| in cases 1–4 and that of |
G
1
|
>
|
G
2
| in case 5.(i) Case
1
(
|
G
1
|
<
|
G
|
<
|
G
2
|
)
. In this case, there are three positive equilibria K
1, K
2, and K
3 as shown in Figure 3. The Jacobian matrix of system (12) at K
1 is
(32)
M
(
K
1
)
=
[
-
k
11
μ
1
V
1
k
11
μ
1
σ
21
V
2
-
k
12
μ
2
σ
12
V
1
-
k
12
μ
2
V
2
]
then p
(
M
(
K
1
)
)
=
k
11
μ
1
/
V
1
+
μ
2
k
12
/
V
2
>
0 and q
(
M
(
K
1
)
)
=
k
11
k
12
μ
1
μ
2
(
σ
12
σ
21
+
1
)
/
V
1
V
2
>
0. That is, K
1 is asymptotically stable. Then we obtain that the stability condition of equilibrium point K
1 is
(33)
0
<
σ
12
<
μ
1
V
1
(
λ
2
-
μ
2
)
μ
2
(
2
α
μ
1
+
λ
1
V
1
-
μ
1
V
1
)
,
k
k
k
k
k
k
k
k
k
σ
21
>
0
,
V
1
>
2
α
μ
1
μ
1
-
λ
1
.Figure 3
Three positive equilibriaK
1, K
2, and K
3.Condition0
<
σ
12
<
μ
1
V
1
(
λ
2
-
μ
2
)
/
μ
2
(
2
α
μ
1
+
λ
1
V
1
-
μ
1
V
1
) means that message 1 is less competitive than message 2. When (33) are satisfied, both message 1 and message 2 coexist, and the speed of message 1 and message 2 will tend to reach a nonzero finite value.Similarly, The Jacobian matrix of system (12) at K
2 is
(34)
M
(
K
2
)
=
[
-
k
21
μ
1
V
1
-
k
21
μ
1
σ
21
V
2
k
22
μ
2
σ
12
V
1
-
k
22
μ
2
V
2
]
then the p
(
M
(
K
2
)
)
=
k
21
μ
1
/
V
1
+
k
22
μ
2
/
V
2
>
0 and q
(
M
(
K
2
)
)
=
k
21
k
22
μ
1
μ
2
(
σ
12
σ
21
+
1
)
/
V
1
V
2
>
0. That is, K
2 is asymptotically stable. Then we obtain that the stability condition of equilibrium point K
2 is
(35)
0
<
σ
21
<
μ
2
V
2
(
λ
1
-
μ
1
)
μ
1
(
2
β
μ
2
+
λ
2
V
2
-
μ
2
V
2
)
,
k
k
k
k
k
k
k
k
k
σ
12
>
0
,
V
2
>
2
β
μ
2
μ
2
-
λ
2
.The condition0
<
σ
21
<
μ
2
V
2
(
λ
1
-
μ
1
)
/
μ
1
(
2
β
μ
2
+
λ
2
V
2
-
μ
2
V
2
) means that message 2 is less competitive than message 1. When (35) are satisfied, both message 1 and message 2 coexist, and their speed will change over time to reach a stable state.The Jacobian matrix of system (12) at K
3 is
(36)
M
(
K
3
)
=
[
-
k
31
μ
1
V
1
-
k
31
μ
1
σ
21
V
2
-
k
32
μ
2
σ
12
V
1
-
k
32
μ
2
V
2
]
then the p
(
M
(
K
3
)
)
=
k
31
μ
1
/
V
1
+
k
32
μ
2
/
V
2
>
0 and q
(
M
(
K
3
)
)
=
k
31
k
32
μ
1
μ
2
(
1
-
σ
12
σ
21
)
/
V
1
V
2. By |
G
1
|
<
|
G
2
|, we have σ
12
σ
21
>
1, and q
(
M
(
K
3
)
)
<
0. It means that K
3 is a saddle point. In int R
+
2, all orbits but K
3 and the separatrixes of K
3 converge to equilibria K
1 and K
2 as shown in Figure 3.(ii) Case
2
(
|
G
1
|
<
|
G
|
=
|
G
2
|
)
. As can be seen in Figure 4, there are two positive equilibrium points K
1 and K
2. It is important to note that, however, K
1 and K
3 in Figure 3 coincide. Similar to the analysis in case 1, K
2 is asymptotically stable, and the stability condition of equilibrium point K
2 satisfies (35).Figure 4
Two positive equilibriaK
1 and K
2.Sincep
(
M
(
K
1
)
)
=
k
11
μ
1
/
V
1
+
μ
2
k
12
/
V
2
>
0 and q
(
M
(
K
1
)
)
=
0, it can be determined that K
1 is a saddle node by the saddle node rule [26]. According to Theorem 1, we can draw the dynamics of system (12) as shown in Figure 4. The same approach also applies to the case |
G
1
|
=
|
G
|
<
|
G
2
|; the difference is that K
1 is asymptotically stable and K
2 is a saddle node equilibrium, and all orbits converge to equilibria K
1 in int R
+
2 except K
2 and the separatrixes of K
2.(iii) Case
3
(
|
G
1
|
≤
|
G
2
|
<
|
G
|
)
. In case 3, there is only one positive equilibrium point K
2. By Theorem 1, K
2 is asymptotically stable as shown in Figure 5, and the stability condition of equilibrium point K
2 satisfies (33). Similar discussions can be given for the case |
G
|
<
|
G
1
|
≤
|
G
2
| in which K
1 is asymptotically stable.Figure 5
The unique positive equilibriumK
2.(iv) Case
4
(
|
G
1
|
=
|
G
2
|
=
|
G
|
)
. In this case, the isoclines L
12 and L
22 coincide, and all the positive equilibrium points distribute on the line segment D
1
D
2. According to Theorem 1, in int R
+
2, all orbits of system (12) converge to the line segment as shown in Figure 6.Figure 6
Equilibrium points on the line segmentD
1
D
2.(v) Case
5
(
|
G
2
|
<
|
G
1
|
)
. In this situation, we have similar discussions with the case 3, and globally stable is K
2 as shown in Figure 5. The same principle applies to the condition of |
G
|
<
|
G
2
|
<
|
G
1
| that K
1 is globally stable. In the case of |
G
2
|
≤
|
G
|
<
|
G
1
| or |
G
2
|
<
|
G
|
≤
|
G
1
|, there is a unique positive globally stable point K
3 of system (12) in Figure 7.Figure 7
The unique positive equilibriumK
3.


## 3.1. Basic Analysis

In order to show the existence of intersection of two isoclinesL
1 and L
2, we assume that the positive or negative effects produced by messages 1 and 2 are almost the same as each other. In the condition, we can set hypothesis as follows: a
=
b
=
z, σ
21
=
σ
12
=
σ, V
1
=
V
2
=
V, and α
=
β
=
ζ, then the values of μ
1 and μ
2 can reflect the difference of the growth factors among the two messages. For |
v
1
-
α
| and |
v
2
-
β
| are smoothed in a very small neighborhood of their vertexes (
0
,
α
) and (
0
,
β
), respectively, the system (12) can be written as
(16)
d
v
1
(
t
)
d
t
=
μ
1
v
1
(
z
+
σ
ζ
V
-
v
1
V
-
σ
ζ
V
|
v
2
ζ
-
1
|
)
,
d
v
2
(
t
)
d
t
=
μ
2
v
2
(
z
+
σ
ζ
V
-
v
2
V
-
σ
ζ
V
|
v
1
ζ
-
1
|
)
;
set x
1
=
v
1
/
ζ and x
2
=
v
2
/
ζ; then
(17)
d
x
1
d
t
=
μ
1
v
1
(
z
+
σ
ζ
V
-
ζ
x
1
V
-
σ
ζ
V
|
x
2
-
1
|
)
,
d
x
2
d
t
=
μ
2
v
2
(
z
+
σ
ζ
V
-
ζ
x
2
V
-
σ
ζ
V
|
x
1
-
1
|
)
.
In int
⁡
R
+
2, the zero growth isoclines of x
1 and x
2 are as follows:
(18)
x
1
=
-
σ
|
x
2
-
1
|
+
C
x
2
=
-
σ
|
x
1
-
1
|
+
C
,
where C
=
V
z
/
ζ
+
σ.According to [24], the isoclines include three parts as shown in Figure 2: positive part (
+
), neutral part (0), and negative part (
-
). The positive part, neutral part, or negative part is the part of isoclines when the slope (for message 1, the slope is d
v
1
(
t
)
/
d
v
2
(
t
); for message 2, the slope is d
v
2
(
t
)
/
d
v
1
(
t
)) is positive part (
+
), neutral part (0), and negative part (
-
). The two isoclines are symmetry around the line x
1
=
x
2. Then, we only need to discuss
(19)
x
1
=
x
2
,
x
2
=
-
σ
|
x
1
-
1
|
+
C
.
Obviously, (19) always have solutions; we can also obverse the following aspects.(i)
There are two different roots(
σ
-
C
)
/
(
σ
-
1
) and (
σ
+
C
)
/
(
σ
+
1
), when σ
2
≠
0,1 and C
≠
1.
(ii)
There is double rootx
=
1 or x
=
C, when σ
2
≠
1 and C
=
1, or σ
=
0.
(iii)
There is root(
σ
-
C
)
/
(
σ
-
1
) or (
σ
+
C
)
/
(
σ
+
1
), when σ
2
=
1 and C
≠
1.The parabolic zero growth isoclines of message 1 (a) and message 2 (b).
(a)
(b)Based on the above analysis, we can conclude that trajectories in phase space always intersect isoclines either in the horizontal or in the vertical direction. In the next section, we make further analysis about the stability of intersection.


## 3.2. Equilibrium

In order to study the final results of interactive effect between two different types of messages, we need to conduct a comprehensive analysis for the stability of the equilibrium points of simultaneous differential equations.Theorem 1.
System (12) admits no periodic orbit.Proof.
Let(20)
B
(
v
1
,
v
2
)
=
1
v
1
v
2
(
v
1
,
v
2
)
∈
int
⁡
R
+
2
then we have
(21)
∂
∂
v
1
[
μ
1
(
a
+
σ
21
β
V
2
-
v
1
V
1
-
σ
21
|
v
2
-
β
|
*
V
2
)
]
+
∂
∂
v
2
[
μ
2
(
b
+
σ
12
α
V
1
-
v
2
V
2
-
σ
12
|
v
1
-
α
|
*
V
1
)
]
=
-
μ
1
v
2
V
1
-
μ
2
v
1
V
2
<
0
.According to the Bendixson-Dulac theorem [25], system (12) admits no periodic orbit.We study the evolution ofv
1
(
t
) and v
2
(
t
) as time t approaches infinity. As shown in Figure 1, there are three equilibrium points in system (12): O
1
(
0,0
), O
2
(
a
V
1
,
0
), and O
3
(
0
,
b
V
2
). To determine the stability of the equilibrium points, let
(22)
M
=
[
f
v
1
f
v
2
g
v
1
g
v
2
]
;
we have
(23)
p
=
-
(
f
v
1
+
g
v
2
)
|
p
i
,
q
=
det
⁡
M
|
p
i
.
Then, the stability conditions of equilibrium point O
i are p
>
0 and q
>
0. The Jacobin matrix of the system (12) at equilibrium state is as follows.(i) For the equilibrium pointO
1
(
0,0
), we get
(24)
M
(
O
1
)
=
[
a
μ
1
0
0
b
μ
2
]
p
1
=
-
(
a
μ
1
+
b
μ
2
)
,
q
1
=
a
b
μ
1
μ
2
.
Based on the assumption that a
=
1
-
λ
1
/
μ
1
>
0 and b
=
1
-
λ
2
/
μ
2
>
0, we have p
1
<
0 and q
1
>
0; the point O
1
(
0,0
) is an unstable node.Since our aim is to show the coexistence of competition and cooperation by the introduction ofα and β, we suppose that the condition (25) is always satisfied in the rest of the paper. Consider
(25)
b
+
σ
12
(
2
α
V
1
-
a
)
>
0
,
a
+
σ
21
(
2
β
V
2
-
b
)
>
0
.(ii) For the equilibrium pointO
2
(
a
V
1
,
0
), we get
(26)
M
(
O
2
)
=
[
-
a
μ
1
a
V
1
μ
1
σ
21
V
2
0
μ
2
[
b
+
(
2
α
V
1
-
a
)
σ
12
]
]
,
p
2
=
a
μ
1
-
μ
2
[
σ
12
(
2
α
V
1
-
a
)
+
b
]
,
q
2
=
a
μ
1
μ
2
(
-
2
α
σ
12
+
a
σ
12
V
1
-
b
V
1
)
V
1
.
The stability condition of equilibrium point O
2
(
a
V
1
,
0
) is σ
12
>
b
V
1
/
(
a
V
1
-
2
α
).Conditionσ
12
>
b
V
1
/
(
a
V
1
-
2
α
) means that message 1 is more competitive than message 2. It implies that the persons who should have adopted message 2 choose to believe message 1. Consequently, message 2 will eventually become extinct in the process of competition, while message 1 can continue to diffuse and reach the maximum speed V
1.(iii) For the equilibrium pointO
3
(
0
,
b
V
2
), we get
(27)
M
(
O
3
)
=
[
μ
1
[
a
+
(
2
β
V
2
-
b
)
σ
21
]
0
b
V
2
μ
2
σ
12
V
1
-
b
μ
2
]
,
p
3
=
b
μ
2
-
μ
1
[
a
+
σ
21
(
2
β
V
2
-
b
)
]
,
q
3
=
b
μ
1
μ
2
(
b
σ
21
V
2
-
a
V
2
-
2
β
σ
21
)
V
2
.
The stability condition of equilibrium point O
3
(
0
,
b
V
2
) is σ
21
>
a
V
2
/
(
b
V
2
-
2
β
).Conditionσ
21
>
a
V
2
/
(
b
V
2
-
2
β
) shows that message 2 is more competitive than message 1. This just means that people accepted message 2 even though they should have believed message 1. Therefore, during the process of competition, message 2 can survive and reach a maximum speed V
2, while message 1 will eventually become extinct finally.


## 3.3. Dynamics

Assume that the conditions in (25) are satisfied; the two messages in system (12) can coexist. The isocline L
1 of message 1 can be subdivided into L
11 and L
12 as follows:
(28)
L
11
:
a
+
σ
21
v
2
V
2
-
v
1
V
1
,
0
<
v
2
<
β
,
L
12
:
a
+
2
σ
21
β
V
2
-
v
1
V
1
-
σ
21
v
2
V
2
,
v
2
>
β
,
where L
1 is smoothed in a small neighborhood of the vertex D
1
(
a
V
1
+
β
σ
21
V
1
/
V
2
,
β
) as shown in Figure 1. Similarly, we divide the isocline L
2 of message 2 into L
21 and L
22 as follows:
(29)
L
21
:
b
+
σ
12
v
1
V
1
-
v
2
V
2
,
0
<
v
1
<
α
,
L
22
:
b
+
2
σ
12
α
V
1
-
v
2
V
2
-
σ
12
v
1
V
1
,
v
1
>
α
,
where L
2 is smoothed in a small neighborhood of the vertex D
2
(
α
,
b
V
2
+
α
σ
12
V
2
/
V
1
) as shown in Figure 1.While the line segmentsL
11 and L
21 are restricted to the intervals (
a
V
1
,
a
V
1
+
β
σ
21
V
1
/
V
2
) and (
0
,
α
), respectively, they are impossible to intersect in the interval α
≤
a
V
1. Then there is a maximum of three isolated positive equilibrium points in system (12). As shown in Figure 1, let K
1
(
k
11
,
k
12
) represent the intersection of L
11 and L
22, let K
2
(
k
21
,
k
22
) represent the intersection of L
21 and L
12, and let K
3
(
k
31
,
k
32
) represent the intersection of L
12 and L
22. Then we have
(30)
k
11
=
V
1
(
a
+
b
σ
21
)
+
2
α
σ
12
σ
21
σ
12
σ
21
+
1
,
k
12
=
V
2
(
V
1
(
b
-
a
σ
12
)
+
2
α
σ
12
)
(
σ
12
σ
21
+
1
)
V
1
,
k
21
=
V
1
(
V
2
(
a
-
b
σ
21
)
+
2
β
σ
21
)
(
σ
12
σ
21
+
1
)
V
2
,
k
22
=
V
2
(
a
σ
12
+
b
)
+
2
β
σ
12
σ
21
σ
12
σ
21
+
1
,
k
31
=
V
1
(
V
2
(
b
σ
21
-
a
)
-
2
β
σ
21
)
+
2
α
σ
12
σ
21
V
2
(
σ
12
σ
21
-
1
)
V
2
,
k
32
=
V
1
(
V
2
(
a
σ
12
-
b
)
+
2
β
σ
12
σ
21
)
-
2
α
σ
12
V
2
(
σ
12
σ
21
-
1
)
V
1
.Suppose thatG
1 and G
2 represent the slope of the lines L
12 and L
22, respectively. Suppose that G denotes the slope of the line D
1
D
2. Then we have
(31)
G
1
=
-
V
2
V
1
σ
21
<
0
,
G
2
=
-
V
2
σ
12
V
1
<
0
,
G
=
V
2
(
b
V
2
V
1
+
α
σ
12
V
2
-
β
V
1
)
V
1
(
a
V
1
V
2
-
α
V
2
+
β
σ
21
V
1
)
<
0
.Dynamics of system (12) can be displayed in five cases, while we discuss the situation of |
G
1
|
<
|
G
2
| in cases 1–4 and that of |
G
1
|
>
|
G
2
| in case 5.(i) Case
1
(
|
G
1
|
<
|
G
|
<
|
G
2
|
)
. In this case, there are three positive equilibria K
1, K
2, and K
3 as shown in Figure 3. The Jacobian matrix of system (12) at K
1 is
(32)
M
(
K
1
)
=
[
-
k
11
μ
1
V
1
k
11
μ
1
σ
21
V
2
-
k
12
μ
2
σ
12
V
1
-
k
12
μ
2
V
2
]
then p
(
M
(
K
1
)
)
=
k
11
μ
1
/
V
1
+
μ
2
k
12
/
V
2
>
0 and q
(
M
(
K
1
)
)
=
k
11
k
12
μ
1
μ
2
(
σ
12
σ
21
+
1
)
/
V
1
V
2
>
0. That is, K
1 is asymptotically stable. Then we obtain that the stability condition of equilibrium point K
1 is
(33)
0
<
σ
12
<
μ
1
V
1
(
λ
2
-
μ
2
)
μ
2
(
2
α
μ
1
+
λ
1
V
1
-
μ
1
V
1
)
,
k
k
k
k
k
k
k
k
k
σ
21
>
0
,
V
1
>
2
α
μ
1
μ
1
-
λ
1
.Figure 3
Three positive equilibriaK
1, K
2, and K
3.Condition0
<
σ
12
<
μ
1
V
1
(
λ
2
-
μ
2
)
/
μ
2
(
2
α
μ
1
+
λ
1
V
1
-
μ
1
V
1
) means that message 1 is less competitive than message 2. When (33) are satisfied, both message 1 and message 2 coexist, and the speed of message 1 and message 2 will tend to reach a nonzero finite value.Similarly, The Jacobian matrix of system (12) at K
2 is
(34)
M
(
K
2
)
=
[
-
k
21
μ
1
V
1
-
k
21
μ
1
σ
21
V
2
k
22
μ
2
σ
12
V
1
-
k
22
μ
2
V
2
]
then the p
(
M
(
K
2
)
)
=
k
21
μ
1
/
V
1
+
k
22
μ
2
/
V
2
>
0 and q
(
M
(
K
2
)
)
=
k
21
k
22
μ
1
μ
2
(
σ
12
σ
21
+
1
)
/
V
1
V
2
>
0. That is, K
2 is asymptotically stable. Then we obtain that the stability condition of equilibrium point K
2 is
(35)
0
<
σ
21
<
μ
2
V
2
(
λ
1
-
μ
1
)
μ
1
(
2
β
μ
2
+
λ
2
V
2
-
μ
2
V
2
)
,
k
k
k
k
k
k
k
k
k
σ
12
>
0
,
V
2
>
2
β
μ
2
μ
2
-
λ
2
.The condition0
<
σ
21
<
μ
2
V
2
(
λ
1
-
μ
1
)
/
μ
1
(
2
β
μ
2
+
λ
2
V
2
-
μ
2
V
2
) means that message 2 is less competitive than message 1. When (35) are satisfied, both message 1 and message 2 coexist, and their speed will change over time to reach a stable state.The Jacobian matrix of system (12) at K
3 is
(36)
M
(
K
3
)
=
[
-
k
31
μ
1
V
1
-
k
31
μ
1
σ
21
V
2
-
k
32
μ
2
σ
12
V
1
-
k
32
μ
2
V
2
]
then the p
(
M
(
K
3
)
)
=
k
31
μ
1
/
V
1
+
k
32
μ
2
/
V
2
>
0 and q
(
M
(
K
3
)
)
=
k
31
k
32
μ
1
μ
2
(
1
-
σ
12
σ
21
)
/
V
1
V
2. By |
G
1
|
<
|
G
2
|, we have σ
12
σ
21
>
1, and q
(
M
(
K
3
)
)
<
0. It means that K
3 is a saddle point. In int R
+
2, all orbits but K
3 and the separatrixes of K
3 converge to equilibria K
1 and K
2 as shown in Figure 3.(ii) Case
2
(
|
G
1
|
<
|
G
|
=
|
G
2
|
)
. As can be seen in Figure 4, there are two positive equilibrium points K
1 and K
2. It is important to note that, however, K
1 and K
3 in Figure 3 coincide. Similar to the analysis in case 1, K
2 is asymptotically stable, and the stability condition of equilibrium point K
2 satisfies (35).Figure 4
Two positive equilibriaK
1 and K
2.Sincep
(
M
(
K
1
)
)
=
k
11
μ
1
/
V
1
+
μ
2
k
12
/
V
2
>
0 and q
(
M
(
K
1
)
)
=
0, it can be determined that K
1 is a saddle node by the saddle node rule [26]. According to Theorem 1, we can draw the dynamics of system (12) as shown in Figure 4. The same approach also applies to the case |
G
1
|
=
|
G
|
<
|
G
2
|; the difference is that K
1 is asymptotically stable and K
2 is a saddle node equilibrium, and all orbits converge to equilibria K
1 in int R
+
2 except K
2 and the separatrixes of K
2.(iii) Case
3
(
|
G
1
|
≤
|
G
2
|
<
|
G
|
)
. In case 3, there is only one positive equilibrium point K
2. By Theorem 1, K
2 is asymptotically stable as shown in Figure 5, and the stability condition of equilibrium point K
2 satisfies (33). Similar discussions can be given for the case |
G
|
<
|
G
1
|
≤
|
G
2
| in which K
1 is asymptotically stable.Figure 5
The unique positive equilibriumK
2.(iv) Case
4
(
|
G
1
|
=
|
G
2
|
=
|
G
|
)
. In this case, the isoclines L
12 and L
22 coincide, and all the positive equilibrium points distribute on the line segment D
1
D
2. According to Theorem 1, in int R
+
2, all orbits of system (12) converge to the line segment as shown in Figure 6.Figure 6
Equilibrium points on the line segmentD
1
D
2.(v) Case
5
(
|
G
2
|
<
|
G
1
|
)
. In this situation, we have similar discussions with the case 3, and globally stable is K
2 as shown in Figure 5. The same principle applies to the condition of |
G
|
<
|
G
2
|
<
|
G
1
| that K
1 is globally stable. In the case of |
G
2
|
≤
|
G
|
<
|
G
1
| or |
G
2
|
<
|
G
|
≤
|
G
1
|, there is a unique positive globally stable point K
3 of system (12) in Figure 7.Figure 7
The unique positive equilibriumK
3.


## 4. Numerical Simulation

To better understand the proposed model, we carry out computer simulations with real data from Sina Weibo, a Chinese microblog. Our data collected by a crowd sourcing Weibo visual analytic system [27], and the interval time is from March 31th, 2014, to May 30th, 2014. We extract two messages published by famous Chinese actors Wen Zhang and Ma Yili about their marriage. In Sina Weibo, Wen Zhang has 53 706 926 fans, ranked in 7th place and Ma Yili has 34 759 267 fans, ranked in 35th place. At midnight of March 31th, Wen Zhang apologised to his families and devoted fans, and then his wife Ma Yili forgives him. For their action in saving marriage, amount of their fans reposted the two messages. As the repost actions have occurred mainly before 14:00, we specify a unit time in 15 minutes to statistical analysis among 11:00–14:00 in March 31th. In order to explore the dynamics process with varied parameters, we put the observations in our model and then further study the influence of the competition coefficient (σ
12 and σ
21) in case of messages propagation with low (or high) critical speed (α and β).We first consider the case that both messages have less competitiveσ
12
<
b
V
1
/
(
a
V
1
-
2
α
) and σ
21
<
a
V
2
/
(
b
V
2
-
2
β
) at low critical speed. Figure 8(a) shows the evolutionary trend of message 1 and message 2 over time. We find that the speed of the two messages can be viewed as a nonmonotonic function, and the velocity curve of message 1 can be divided into three distinctly different stages. In the initial stage, message 1 appears in OSNs and then attracts more and more users, leading to rapid spread of message 1. However, message 1 descends slowly followed by the end of initial stage, because of the natural dying out of its own as well as the competition coming from message 2. Specifically, since message 1 is less competitive than message 2, the growth of message 2 hinders the spread of message 1 and message 1 will continuously decay due to the competition. Finally, message 1 arrives at its equilibrium and then converges to a constant value. The evolution process of message 2 is similar to message 1. After a slow growth in the initial stage, message 2 spreads rapidly and then slows down till a steady state is reached.(a)v
1 and v
2 as a function of t, with v
1
(
0
)
=
2, v
2
(
0
)
=
1, V
1
=
3000, α
=
80, σ
12
8
(
a
)
=
0.6, μ
1
=
0.7, λ
1
=
0.2, V
2
=
1800, β
=
20, σ
21
8
(
a
)
=
0.8, μ
2
=
0.5, and λ
2
=
0.15. (b) The development trend of relationship between message 1 and message 2.
(a)
(b)The change trends of message 1 and message 2 can be observed from Figure8(b). At the beginning, the two messages grow quickly and then slow down and converge to a nonzero constant. This process reflects that, with less competitive and low critical speed, message 1 and message 2 will cooperate with each other in the propagation. Both of them can attract a certain amount of users, and they can coexist in a period of time.In the case ofσ
12
<
b
V
1
/
(
a
V
1
-
2
α
) and σ
21
>
a
V
2
/
(
b
V
2
-
2
β
), the message which is more competitive than the other will survive at the final state. As shown in Figure 9(a), along with the increasing number of people who have heard message 2, the speed of message 1 increases at the beginning and then decreases rapidly. The main reason for this phenomenon is that the two messages have cooperated with each other when their speed is lower than the critical speed at the beginning, and then competed against each other when their speed is higher than the critical speed after the initial phase. Finally, only the more competitive message (
σ
21
9
(
a
)
≫
σ
12
9
(
a
)
) survived in the system. When σ
12
>
b
V
1
/
(
a
V
1
-
2
α
) and σ
21
<
a
V
2
/
(
b
V
2
-
2
β
), we can get a similar conclusion that message 2 will vanish at stable state. That is, the more competitive message can draw the users attention and spread widely while the message with lesser competitive is difficult to attract users and vanished after a while.(a)v
1 and v
2 as a function of t, with v
1
(
0
)
=
2, v
2
(
0
)
=
1, V
1
=
3000, α
=
80, σ
12
9
(
a
)
=
0.6, μ
1
=
0.7, λ
1
=
0.2, V
2
=
1800, β
=
20, σ
21
9
(
a
)
=
2, μ
2
=
0.5, and λ
2
=
0.15. (b) The development trend of relationship between message 1 and message 2.
(a)
(b)The change trends of message 1 and message 2 in Figure9(b) is similar to that in Figure 8(b). The difference is that only message 2 can survive when a final state is reached. Under the constraint condition of parameters σ
21
9
(
a
)
≫
σ
21
8
(
a
), we can make a conclusion that message 2 attracts amount of users from message 1 at low critical speed and spreads widely in online social networks.As shown in Figure10(a), under the condition of σ
12
>
b
V
1
/
(
a
V
1
-
2
α
) and σ
21
>
a
V
2
/
(
b
V
2
-
2
β
), only message 1 survives and message 2 vanishes in steady state, while they coexist with each other at the final stage in Figure 10(b). It can be inferred that in strong competitive environment, the high critical speed promotes the coexistence of the two messages. At the same time, although restricted by σ
21
<
σ
12, the velocity values of both message 1 and message 2 in Figure 10(b) are greater than that in Figure 10(a); it can be concluded that when the two messages are both competitive, they will benefit from the high critical speed.Positive or negative effects for different critical speed withv
1
(
0
)
=
2, v
2
(
0
)
=
1, V
1
=
3000, σ
12
10
(
a
)
=
3, μ
1
=
0.7, λ
1
=
0.2, V
2
=
1800, σ
21
10
(
a
)
=
2, μ
2
=
0.5, and λ
2
=
0.15. (a) v
1 and v
2 as a function of t, with α
=
80, β
=
20. (b) v
1 and v
2 as a function of t, with α
=
2000, β
=
1200.
(a)
(b)Generally, more competitive messages disseminate more rapidly than less competitive messages over network. Unfortunately, the situation is often not the case. It can be observed that under the conditions with the sameσ
21, message 2 can survive and reach its maximum speed in Figure 9(a), but it vanished finally in Figure 10(a). The difference exists because the competition coefficient of message 1 (σ
12) is considerably larger in Figure 10(a) than that in Figure 9(a) (σ
12
10
(
a
)
≫
σ
12
9
(
a
)). It means that the advantage of the prior competitive message will quickly disappear as the competitiveness gap narrowed. Therefore, under the condition of low critical speed, it is necessary to adopt Tit-for-Tat-like strategy to improve their competitiveness when there already exists a competitive message in online social network.


## 5. Conclusions

The majority of studies on information propagation have focused on a single contagion spreading through the social network without considering the influence of other contagions. In this paper, we have explored the information spreading mechanism impacted by cooperative and competitive effectiveness during the process of the propagation of two types of messages over networks. We developed a new OSN information propagation model based on Lotka-Volterra models to demonstrate the dynamics of a specific cooperation-competition system of these two messages. The stability of the system was proved by the differential equations and the phase trajectory of the stability theory. The results of stability analysis further indicated that the differential equations admit no periodic orbit, and the system can reach steady state in certain condition.We proved that two types of messages usually cannot spread together peacefully, and only one type message survives till final state. Nevertheless, it produced stable points in the condition of high (or low) critical speed, which shows that it is possible for the coexistence of both messages in competitive propagation of different contagions. Using the real data collected from Sina Weibo, the results validated the theoretical analysis and presented another interesting conclusion that the messages will benefit from the high critical speed when they are both competitive. If there already is a competitive message in OSN, it is wise to adopt a Tit-for-Tat strategy.As the great massive data are generated every moment in OSN, the rapid progress of information technology provides us with an easy way to collect data and observe the phenomenon that various types of information flow in the form of cascades on Twitter or other OSNs. Naturally, the proposed dynamic information propagation model should be tested and applied in real OSN with large topology network and vast numbers of the nodes and connections. These issues will be addressed in future research.



---

*Source: 610382-2014-07-13.xml*