Abstract

Coevolution spreading dynamics on complex networks is a hot topic, which attracts much attention in network science. This paper proposes a mathematical model to describe the two competing complex information spreading dynamics on multiplex networks. An individual can only accept one of the two pieces of information. A heterogeneous mean-field theory is developed to describe the spreading dynamics. We reveal different regions through Monte Carlo simulations of the competing complex information spreading dynamics: no global information, one information dominant, and two information coexistence. We finally find that the heterogeneity of the multiplex networks’ degree distributions does not qualitatively affect the results.

1. Introduction

Social networks, e.g., WeChat, Twitter, and Facebook, are the convenient way to express and share information with friends [15]. Students may retransmit an article to support or oppose the opinions of the authors. To describe the interaction of two types of information spreading dynamics on social networks, researchers used the framework of coevolution spreading dynamics on complex networks [6]. Generally speaking, scholars used complex networks to describe the social network, in which the nodes represent the individuals and edges stand for the relationships among individuals. Extensive empirical analyses revealed that social network exhibits degree heterogeneity, high clustering, multilayer, and temporality [718]. Including some of those characteristics, many research studies have been designed (see some important reviews in this field, e.g., reference [1921]).

Based on the shared information’s content, we can classify the information spreading dynamics into simple and complex information [2224]. The simple information means that a single contact between the informed and susceptible individuals can trigger the infection transmission and widely used to describe epidemic spreading and simple news diffusion. In contrast, single contact cannot eliminate some complex information’s legality and reliability (e.g., political networks). Therefore, social reinforcement is included in the dynamics of complex information spreading. In the following content, we present the progress for simple and complex information spreading detailedly.

The study about the simple information spreading on complex networks independently is the most widely investigated and revealed some critical results. For instance, simple information can always spread in the social network regardless of the values of information transmission probability when the heterogeneity of degree distribution is extremely strong [2527]. Besides, the high clustering promotes the breakout of the information while suppresses the spreading size for large transmission probability [28]. Once several information are transmitted simultaneously, scholars revealed more novel phenomena. When two pieces of information are spreading on networks consecutively, the second information’s outbreak threshold is larger than the threshold of the first information [29]. Karrer and Newman [30] further considered two competing spreading dynamics on the single network simultaneously. By using the competing percolation theory, they found that two epidemic cannot coexist in the thermodynamic limit (i.e., when the network size is extremely large). Sahneh and Scoglio [31] investigated two competing information on multiplex networks and computed the absolute-dominance and coexistence regions. When two information collaborate, i.e., if the acceptance probability for another information is enlarged when an individual accepts one information, the system exhibits a discontinuous phase transition [3234]. Scholars considered another situation, asymmetric interacting, i.e., one spreading dynamic promotes the other. In contrast, the other one suppresses the spreading and reveals that the global epidemic will be greatly suppressed once the information spreading is included [3537].

To include the social reinforcement in the dynamics of complex information spreading, the linear threshold model is the most widely used, in which an individual accepts the information only when the received information is larger than a threshold [22, 3840]. Setting the final spreading size as the order parameter, researchers revealed the discontinuous phase transition in the system. More detailedly, the final spreading size versus the transmission probability exhibits a discontinuous growth. Liu et al. [41] investigated two complex information spreading consecutively on single complex network and revealed a transition for continuous and discontinuous phase transitions. Furthermore, Liu et al. [42] studied two complex synergetic information spreading dynamics on multiplex networks and revealed that the synergetic induced the discontinuous phase transition. To the best of our knowledge, there is still lacking a mathematical model to investigate two competing complex information spreading dynamics on multiplex networks. We will propose a mathematical model in section 2 and then develop a heterogeneous mean-field theory in section 3. We further study the complex competing information spreading model by using Monte Carlo simulation in section 4. Finally, we conclude in section 5.

2. Competing Complex Information Spreading Model

This section proposes the competing information spreading model on social network to describe different opinions evolving in the campus.

2.1. Network Description

Students usually discuss the events happened and diffuse the information or opinions through social platforms. We here use a multiplex network to describe the multiplex social network with two layers and , where and , respectively, represent the sizes of networks and and and , respectively, stand for the number of edges in of networks and . Nodes stand for the student users, and edges mean the relationships among those users. If two students and communicate with each other in network , an edge exists between them.

We assume there are students in the social network and and connections among them. Therefore, the average degrees of the social network are and , respectively. The two networks’ degree distribution is denoted as and , respectively. Nodes with large degrees mean that they have more communication with friends and always act as hubs. A node in two networks represents that it communicates with others through two platforms. To build the multiplex network, we first generate the intralayer edges using the uncorrelated configuration model proposed in reference [43] according the degree distributions and . Then, we build the interlayer edges randomly, that is to say, we randomly match two nodes in the two subnetworks.

2.2. Competing Spreading Dynamics

In social networks, students always express different opinions on the same event, and as a result, additional information is competing with each other. In reality, adopting an idea is risky. Therefore, social reinforcement is included. Social reinforcement means that the adoption of the information requires multiple verifications from friends. We here adopt a generalized susceptible-informed-recovered model to describe the two information spreading dynamics. The susceptible node means that it does not know the information and can adopt the information. An informed node stands for that adopted the information and willing to transmit it to friends. A node in the recovered state means that it has lost interest in the information.

The two competing information spreading dynamics are denoted as and . We assume information and , respectively, spread on networks and simultaneously. The competing between two dynamics means that a node cannot adopt two information simultaneously. The two competing spreading dynamics evolve as follows. Initially, we randomly select a node infected by and in two subnetworks. At every time step, every node informed by information () will transmit the information to its neighboring nodes in network with probability . To eliminate the information’s credibility, we include the social reinforcement and assume a node adopting information must receive pieces. A susceptible node , respectively, obtains and pieces of and information at time step , and it adopts information with probability and adopts information with probability . The informed nodes are recovered with probability . The two competing spreading dynamics evolve until there are no nodes in the informed state. The pseudocode of the competing spreading dynamics as shown in Algorithm 1. At each time step, the time complexity of the algorithm is , where is the number of edges in the network. The space complexity is .

(1)Input: network , and dynamical parameters , , , ;
(2)Output: spreading sizes of information and and denoted as and ;
(3);
(4)Randomly a seed for information and , and put them into queue ;
(5)while is not empty do
(6);
(7);
(8) Initialize to be empty;
(9);
(10)for to do
(11)  if Node adopt the information then Node transmits the information to every susceptible neighbor in subnetwork with probability ;
(12)   if Node received the information successfully then
(13)    ;
(14)   end if
(15)  end if
(16)  if Node adopt the information then
(17)   Node transmits the information to every susceptible neighbor in subnetwork with probability ;
(18)   if Node received the information successfully then
(19)    ;
(20)   end if
(21)  end if
(22)  if and then
(23)   Node adopts information with probability ;
(24)   if Node adopts information then
(25)    Adding node into queue ;
(26)   end if
(27)   Node adopts information with probability ;
(28)   if Node adopts information then
(29)    Adding node into queue ;
(30)   end if
(31)  end if
(32)  if and then
(33)   Node adopts information with probability 1;
(34)   Adding node into queue ;
(35)  end if
(36)  if and then
(37)   Node adopts information with probability 1;
(38)   Adding node into queue ;
(39)  end if
(40)end for
(41)for to do
(42)  Recovering node with probability ;
(43)  if Node recovers then
(44)   Delete node from queue ;
(45)  end if
(46)end for
(47)for to do
(48)  Adding node to queue ;
(49)  Deleting node from queue ;
(50)end for
(51)
(52)end while

3. Theoretical Analysis

Mathematically, the competing information spreading dynamics on social networks can be described by using the heterogeneous mean-field approach [25, 26, 44], which is widely used in investigating the dynamics on networks. There are some basic assumptions behind the heterogeneous mean-field approach as follows:(i)The network size is large enough.(ii)There are no degree-degree correlations in the network.(iii)The infection probability of informed nodes is independent, i.e., there are no dynamical correlations among nodes. Based on the above assumptions, we derive the evolution equations for competing information spreading dynamics.

To describe the evolution equations, we use the following mathematical notations , , and to represent the probability that a node with degree in subnetwork is in the susceptible, informed, and recovered state for information at time . Similarly, we use , , and to represent the probability that a node with degree in subnetwork is in the susceptible, informed, and recovered state for information at time . Since the symmetry of the two spreading dynamics, we only describe the evolution equations for one of them. At time step , the fraction of nodes in susceptible, informed, and recovered for information is , , and , respectively. In the following, we derive the evolutions of , , and .

The fraction of decreases once node with degree in subnetwork becomes informed. On the one hand, node must receive pieces of information in subnetwork at time . In uncorrelated network, the probability of node connecting to an informed node with degree is . Therefore, the probability of node transmits the information to node successfully which is . The pieces of information that node received at time are

According to descriptions of the model, we know node adopts information with probabilitywhere if ; otherwise, . Note that is the other type of information. Based the above discussions, we know the evolution of as

With a similar discussion, we know the evolutions of and asrespectively.

4. Results Analysis

In this section, we numerically study the competing information spreading on multiplex networks detailedly. We study the competing information spreading on two types of multiplex networks. The first one is ER-ER multiplex network, i.e., networks and are Erdös–Rényi (ER), in which the degree distribution follows a Poisson distribution. That is to say, the degree distribution of subnetwork is , where is the average degree of subnetwork . We assume there is no interlayer degree-degree correlations, and thus the joint degree distribution is . For ER-ER multiplex networks, the degree distribution is homogeneous. In reality, many real-data revealed that the degree distribution is heterogeneous. We therefore use the SF-SF multiplex networks, in which the two subnetworks are both scale-free networks. The degree distribution of subnetwork is , where . Similar with ER-ER multiplex networks, we assume there is no interlayer degree-degree correlations, and the joint degree distribution is . In numerical simulations, we set the network size and . The results presented in the paper are averaged over at least 1000 times.

4.1. ER-ER Multiplex Networks

Because the social reinforcement is included in the information spreading dynamics, we first study the effects of initial seed size in Figure 1. To study the fraction of new informed nodes by each information, we investigate and versus for different values of and . For small values of , e.g., , the fraction of new informed nodes by the seeds depends on the fraction of seeds shown in Figure 1(a). When the seed is small (e.g., ), any values of information transmission probability cannot trigger the global outbreak for information . Since the social reinforcement is included, a small fraction of seeds cannot further adopt the information. With the increase in , the global outbreak of becomes possible when is above the threshold point . Note that can be numerically determined by studying the variability of , and at the threshold point, exhibits a peak. We do not show in this paper. Note that the global information of cannot outbreak when , as shown in Figure 1(b). For larger values of , such as and 0.8, we find that a small fraction of seeds cannot trigger the global outbreak for the two types of information (see Figures 1(c)1(f)). For larger values of and 0.1, increases with , meanwhile decreases with . We note that two types of information may coexist, and the coexistence conditions are discussed detailedly in Figure 2.

In Figure 2, we investigate the fractions of new informed nodes by the two types of information triggered by the initial seeds on a plane . As shown in Figure 2(a), the global information can not trigger small values . Previous studies indicated that the global information outbreak threshold of information is when there is only information spreading on network . However, when the two competing information spreads on multiplex networks simultaneously, the global information cannot break when slightly is more extensive than . We know that the competing information suppresses the spreading dynamics, and the global information outbreak threshold is much larger than . With the two information spreading dynamics’ symmetry mechanisms, we find the similar phenomena of in Figure 2(b).

One important question is the domination and coexistence of regions. Domination means that one information spreads to most nodes, while the other only transmits to a few nodes. Coexistence stands for the two types of information spread to a large fraction of nodes. To this end, we study in Figure 2(c). For the case of , we know information dominates. When , the information dominates. When , the two types of information compete for the nodes. Combing Figures 2(a)2(c), we speculate information dominates when and . Similarly, dominates when and . The two types of information are competing with the nodes and can coexist when , , and .

In Figure 3, we further study the effects of the average degrees on the competing information spreading dynamics when . When the two types of information spreading on the two networks with the same average degree, i.e., , we find that and exhibit the similar behavior, that is to say, and increase monotonously with the increase in . However, first increases with and then decreases when , e.g., and . That is to say, there exists a peak at which reaches its maximum. Meanwhile, always increases with regardless the values of and .

4.2. SF-SF Multiplex Networks

Finally, we study the two competing information spreading on SF-SF multiplex networks in Figure 4 and investigate the effects of degree heterogeneity of the two networks on the spreading dynamics. Compared with the ER-ER networks, the SF-SF networks promote the spreading dynamics since the existence of some nodes with large degrees. In addition, we find that the degree heterogeneity does not qualitatively affect the phenomena presented in Figure 2. That is to say, the domination and coexistence regions are the same. The domination region of information is and . Similarly, information domination region is and . The coexistence region is , , and .

5. Discussion

Coevolution spreading dynamics aims at describing interacting dynamics and revealing the phenomena induced by distinct interacting mechanisms. In this paper, we proposed a competing complex information model to describe two types of information diffusion on multiplex networks. A heterogeneous mean-field theory is developed to describe the model. Through extensive numerical simulations, we revealed the conditions of different regions. Specifically, no global information outbreak region is and ; information dominant is and ; information dominant is and ; and two information coexistence region is , , and . When the two competing complex information spread on a multiplex with heterogeneous degree distribution, the phenomena do not qualitatively be affected. The results presented in this paper may shed some light into further investigating the competing spreading, such as design of effective strategy to make one spreading always win regardless of the other spreading dynamics and network structures.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partially supported by the Ministry of Education Science and Technology Development Center University Research Innovation Fund: “Qingtai Digital Intelligence Integration” Collaborative Innovation, no. 2020QT20.