An Evolution Model of Rumor Spreading Based on WeChat Social Circle

Lubang Wang* and Yue Guo**

Abstract

Abstract: With the rapid development of the Internet and the Mobile Internet, social communication based on the network has become a life style for many people. WeChat is an online social platform, for about one billion users, therefore, it is meaningful to study the spreading and evolution mechanism of the rumor on the WeChat social circle. The Rumor was injected into the WeChat social circle by certain individuals, and the communication and the evolution occur among the nodes within the circle; after the refuting-rumor-information injected into the circle, subsequently, the density of four types of nodes, including the Susceptible, the Latent, the Infective, and the Recovery changes, which results in evolving the WeChat social circle system. In the study, the evolution characteristics of the four node types are analyzed, through construction of the evolution equation. The evolution process of the rumor injection and the refuting-rumor-information injection is simulated through the structure of the virtual social network, and the evolution laws of the four states are depicted by figures. The significant results from this study suggest that the spreading and evolving of the rumors are closely related to the nodes degree on the WeChat social circle.

Keywords: Evolution Equation , Rumor Spreading , Social Networking , WeChat Social Circle

1. Introduction

Network social, as a kind of virtual social, is considered as the expansion of the real social networking in the internet space. Multimedia-on-demand system in the network has increased the social people’s links by the virtual network [1] and the behavior of people on the network is a kind of mapping of the real behavior. As a part of their social activities, people play different roles and accept, inject and spread information in different social networking circles, by joining different social networking platforms. The current literature suggest that the social network structure has the small world effect [2], or the scale-free characteristic on the degree distribution [3]. The topology association between individuals is formed, due to information dissemination between the individuals, and the topology of different social networks has been studied by many scholars. The previous research [4-7] studied the structure of social network from the perspective of complex network. The topological structure of social network, including the property of degree distribution, the clustering coefficient, the vertex degree, and the correlation coefficient, have been subjected to study by many researchers.

The spread of rumors on social networks has become one of the research hotspots, during the past ten years. The studies of rumors in the fields of communication, psychology, sociology and other fields mainly analyze the connotation of the rumor [8], the influence factors of the rumor transmission [9], and the psychological motivation [10], by means of investigation or experiment. Daleya and Kenda [11] and Maki and Thompson [12] studied the spread of rumors from the perspective of epidemic-like spreading models. The DK model proposed by Daley and Kendal [11] divides the population into three different categories including people who do not know the rumor (Susceptible), people who know and disseminate the rumor (Infective), and people who know but do not spread the rumor (Recovery). The number of people who eventually heard rumors is almost unrelated to the number of first unknowns in the spread of rumor, that is, there is no critical value of the spread of rumors [11]. Maki and Thompson [12] suggested that rumors are spread through the direct contact between the communicator and the other individuals. The MT model proposed by these researchers suggests that only the initial communicator will become recovery when one communicator encounters another communicator [12]. However, both research on the DK model and the MK model lacks the deep consideration of the structure of the rumor spreading network. The later scholars have performed a great deal of expansion to address this issue. Doerr et al. [13] studied the impact of the network topology on the spread of rumors. Zanette [14, 15] proposed the rumor spreading model, in the small world network and suggested the threshold of the rumor spreading. Moreno et al. [16] studied the dynamic mechanism of the rumor spreading in the scale-free networks. Based on the classical rumor spreading model, Pan et al. [17] studied the rumor spreading in the scale-free network with power law distribution and variable clustering coefficient. The findings from their research show the higher the clustering coefficient, the better the suppression of rumor spread [17]. Thus, the research of the rumor spreading is transformed from study of the propagation characteristics based on the model itself, to the research on the propagation oriented network topology structure.

The study of the rumor communication based on various indicators (i.e., the degree, the clustering coefficient, the distance, etc.) of network structure has a realistic basis, but the platform for the spreading rumors is the social medium, such as Twitter, Sina microblog, WeChat and other similar platforms. The communication law of information in the virtual communities and blogs is analyzed by name of authors [18]. The environment of the blog communication preserve some of the characteristics of a social network, therefore, we can see the blog network as a social network based on every blogger node. The literature [19] extracted characteristic factors from four aspects, including the publishing users, the receiving user, the user intimacy, and the timeliness of information, which affect the user’s forwarding behavior. On the basis of the dynamic model of infectious disease, the user behavior analysis and contact nodes were introduced, subsequently, a SCIR model was proposed, based on the user behavior analysis. The proposed model was used to analyze the association and evolution among nodes, in micro-blog social network. A two stage A-SC1 C2IR model was constructed by Song et al. [20] and introduced the dynamic equation of the model. The simulation results suggested that the prematurely intervention of authoritative information can effectively control the spread of micro-blog rumor. Lu et al. [21] studied the crowd of WeChat rumor spread and constructed the S1S2IR double S rumor spreading model, consequently, they provided the steady state value of the mean field equation of the model. Luo [22] compared the causes and governance schemes of the rumor spread between the micro-blog and the WeChat. Therefore, after the propagation research oriented to network topology structure, some literatures began to pay attention to the rumor spreading research, based on the real social background or the social platform.

The above literatures have carried out a more in-depth study on the related fields of the rumor, such as the social network, the communication model, and the rumor communication platform, as the main factors in this field. Due to the huge number of users and the huge information flow of WeChat, the research value of the rumor spread based on the WeChat platform is particularly prominent. In the study of WeChat rumors in [21], the two states of S1 S2 match the channels of rumors (i.e., the public number and the common node). However, in reality, most of the information published by the public number is the content of the common node transmission and the WeChat’s Circle, therefore, the display of this WeChat’s Circle is an important channel for the spread of rumors. Because common node transmission is a one to one transmission, and the spread of WeChat’s Circle is a display for all friends. This display is generally persistent and affects the friends’ nodes, repeatedly. Therefore, taking into account the particularity of WeChat’s Circle, this research creates a SLIR (Susceptible, Latent, Infective, Recovery) model for the rumor spreading. Subsequently, the rumor spreading law is deduced and the real data simulation is carried out based on the model.

2. Rumor Spreading Model SLIR

Tencent’s WeChat mobile dating software is the most popular and influential universal social application software. In this study, the rumor evolution model in the social network is based on the user functions and privileges of the WeChat’s dating software. Of course, the functions and user permissions of the WeChat and other social software will change throughout the time, but the basic functions of information sharing and information acquisition remains relatively unchanged. At present, the WeChat basically share and obtain the information through sending WeChat’s Circles, and information are released from friends group, and information transmission between the friends.

Based on the reality of WeChat social network, it constructs an undirected and unweighted attribute graph , in which V represents the set of nodes of the graph [TeX:] $$\mathrm{V}=\left\{V_{i} | i \in N\right\}$$. And E is a set of the edge of [TeX:] $$\mathrm{E}=\left\{\overline{{V}_{i} V_{j}}| i, j \in N \wedge i \neq j\right\}$$, where [TeX:] $$\overline{V_{i} V_{j}}$$ indicates that [TeX:] $$V_{i}$$ and [TeX:] $$V_{j}$$ have an edge connection, in other words, the information between the two nodes is connected ; and A is a collection of attributes of all nodes of the graph, and a set of attributes is defined as the four states in [TeX:] $$\mathrm{A}=\{S, L, I, R\}$$, which is corresponded to: S is a healthy and susceptible state[[TeX:] $$Susceptible$$], L is a latent state[[TeX:] $$Latent$$], [TeX:] $$I$$ is the state of spreading[[TeX:] $$Infective$$], and R is an immune state[[TeX:] $$Recovery$$], respectively.

According to the spreading law of rumor information in social network, the basic situation of the four kinds of nodes of [TeX:] $$S, L, I,$$ and [TeX:] $$R$$ in the network are as follows: susceptible node S indicate that no information about rumors is received in a social network; The infective node I indicates that the node accepts the rumor information from its neighbor nodes and implements activities to propagate the information, such as sending information in WeChat’s Circles. The latent node L indicates that the rumor information has been received or seen, but the nodes themselves do not share or disseminate information, however, they may share information or disseminate information at any time (now or later); recovery node R indicates that the node has accepted the message of refuting the rumor, which comes from the neighbor nodes (Refute the rumor by asking or looking at WeChat’s Circles).

The transfer between nodes in susceptible state, latent state, infective state, and recovery state, which depends not only on the state of the node itself, but also on the state of its neighbor nodes. On this basis, the following spreading rules are defined:

(1) If a susceptible Node S is contacted with an infective node I, the node S will convert to node L with the probability [TeX:] $$P_{1}$$. The latent node L will be converted to the infective node I with the probability [TeX:] $$P_{3}$$, and it is necessary to explain that the probability [TeX:] $$P_{3} : P_{3}$$ is no related to topology structure between nodes, the conversion [TeX:] $$\mathrm{V}_{\mathrm{j}}(L) \rightarrow \mathrm{V}_{\mathrm{j}}(I)$$ of the node [TeX:] $$\mathrm{V}_{\mathrm{j}}$$ state L to the I is completely subject to the subjective factors of the node itself (psychological factors);

(2) When the infective node [TeX:] $$I$$ meets the recovery node R, it converts to the node R by probability [TeX:] $$P_{2}$$; when the latent node L meets the recovery node R it will be converted to the node R by probability [TeX:] $$P_{2}$$;

(3) The rumor comes from outside and is injected into a point [TeX:] $$\mathrm{V}_{\mathrm{j}}$$ that makes [TeX:] $$\mathrm{V}_{\mathrm{j}}(S) \rightarrow \mathrm{V}_{\mathrm{j}}(I)$$. In other words, if there is no external injection, the closed social network has no so-called rumors. In the same way, the information of refuting the rumor also comes from the external injection.

Fig. 1 represents the rules of communication. It needs to be explained that the probability [TeX:] $$P_{3}$$ reflects the subjective state of the node itself. Some people see the rumor and believe it, hence they spread the rumor, while some people believe the rumor, but they reject to spread it. So, the probability [TeX:] $$P_{3}$$ is generally less than 1 in the real-world scenario. As for the “external injection” in Rule 3 is that the rumor does not come from the “social network circle”, but it comes from the acquisition of the individual outside the “social network circle”, perhaps, Injected into “the social network circle” by somebody after being coined.

Fig. 1.
The relation figure of node state conversion.

In the transformation rule of the spreading model (1), “if a susceptible node S is in contact with a infective node [TeX:] $$I$$, node S converts to the latent node [TeX:] $$L$$ with the probability [TeX:] $$P_{3}$$.” It can be understood as, [TeX:] $$\mathrm{V}_{\mathrm{j}}(\mathrm{S}) \stackrel{\mathrm{p}_{1}}{\longrightarrow} \mathrm{V}_{\mathrm{j}}(\mathrm{L})$$, and the converted result [TeX:] $$\mathrm{V}_{\mathrm{j}}(L)$$ will be converted to infective node [TeX:] $$I$$ with the probability [TeX:] $$P_{3}$$, and still is a latent node [TeX:] $$L$$ with the probability [TeX:] $$1 - P_{3}$$, as shown in Fig. 2.

Fig. 2.
Explanatory drawing of node S, L, I state conversion illustration.

The reality is that the conversion of [TeX:] $$\mathrm{V}_{\mathrm{j}}(\mathrm{L}) \stackrel{\mathrm{p}_{3}}{\longrightarrow} \mathrm{V}_{\mathrm{j}}(\mathrm{I})$$ is usually happened while [TeX:] $$\mathrm{V}_{\mathrm{j}}(\mathrm{S}) \stackrel{\mathrm{p}_{1}}{\longrightarrow} \mathrm{V}_{\mathrm{j}}(\mathrm{L})$$ changes. Since, a node [TeX:] $$S$$ is contacted with rumor mode [TeX:] $$I$$, if node [TeX:] $$S$$ trusts the rumor, and the rumor is forwarded in no time (node [TeX:] $$S$$ becomes [TeX:] $$I$$); if node [TeX:] $$S$$ does not trust the rumor, and the rumor is not forwarded (node [TeX:] $$S$$ will be latent node [TeX:] $$L$$), at that time and afterwards. So we can think of the approximation of the spreading model, as shown in Fig. 3 (the revised rule provides convenience for the subsequent equation description).

Fig. 3.
Revised drawing of node state transformation relationship.

In this way, the spreading model transformation rules (1) can also be expressed as: If a susceptible node S is contacted with a infective node [TeX:] $$I$$, the susceptible node will convert to the infective node [TeX:] $$I$$ with the probability [TeX:] $$p_{1} * \alpha$$, and will convert to the latent node L with the probability [TeX:] $$p_{1} * \beta$$ (among them, [TeX:] $$\beta=1-\alpha$$ ).

3. Dynamic Analysis of Spreading Model

3.1 Average Transfer Probability of Node State

Assume that a node [TeX:] $$V^{j}$$ is a susceptible node, which is expressed as [TeX:] $$V^{j}(S) \cdot P^{j}_{S S}$$ indicates the probability that node [TeX:] $$V^{j}(S)$$ is still in a susceptible state in the time period of [TeX:] $$[t, t+\Delta t]$$. Accordingly, [TeX:] $$P_{S L}^{J}$$ represents the probability of the transition of node [TeX:] $$V^{j}(S)$$ into the latent node [TeX:] $$V^{j}(L)$$ in the time period of [TeX:] $$[t, t+\Delta t] \cdot P^{j}_{S I}$$ represents the probability that the node [TeX:] $$V^{j}(S)$$ is transformed into the infective node [TeX:] $$V^{j}(I)$$ in the time period of [TeX:] $$[t, t+\Delta t]$$. In similar tags, the probability expression of the node state invariance is also [TeX:] $$P^{j}_{L L} P^{j}_{I I}$$, and the probability expression of the node state change is [TeX:] $$P^{j}_{L I}, P^{j}_{L R}, P^{j}_{I R},$$ with the similar meaning to the above explanation. Then,

(1)
[TeX:] $$P^{j}_{S S}=(1-\Delta t * p 1)^{G}$$

Among them, [TeX:] $$G=G(t)$$ indicates the number of infective nodes in the neighbor of the [TeX:] $$V^{j}(S)$$ node at the t moment. It is assumed that the node [TeX:] $$V^{j}(S)$$ contains K edges, and [TeX:] $$G(t)$$ is a random variable with two distribution, as follows:

(2)
[TeX:] $$\prod(G, t)=\left(\begin{array}{l} {k} \\ {G} \end{array}\right) \omega(k, t)^{G}(1-\omega(k, t))^{(k-G)}$$

where, [TeX:] $$\omega(k, t)$$ represents the probability of connecting a susceptible node with a K edge to an infective node at [TeX:] $$t$$ time. [TeX:] $$\omega(k, t)$$ can be written as:

(3)
[TeX:] $$\omega(k, t)=\sum_{k^{\prime}} p\left(k^{\prime} | k\right) p\left(I_{k^{\prime}} | S_{k}\right) \approx \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)$$

where, [TeX:] $$\sum_{k} p\left(k^{\prime} | k\right)$$ is a degree correlation function, which represents the conditional probability of a node with a degree of [TeX:] $$k$$ and a node with a degree [TeX:] $$k^{\prime} ; p\left(I_{k^{\prime}} | S_{k}\right)$$ represents the probability of being in a infective state under the condition of node with a [TeX:] $$k^{\prime}$$ edges that are connected to a susceptible node with a degree of [TeX:] $$k \cdot \rho^{I}\left(k^{\prime}, t\right)$$ represents the density of the infective nodes of the scale of [TeX:] $$k^{\prime}$$ at the time of [TeX:] $$t$$.

Therefore, the average transfer probability [TeX:] $$\bar{p}_{s s}(k, t)$$ of a susceptible node with a degree of K, during the [TeX:] $$[t, t+\Delta t]$$ preiod can be calculated as follows:

(4)
[TeX:] $$ \begin{aligned} \overline{p}_{s s}(k, t) &=\sum_{G=0}^{k}\left(\begin{array}{c} {k} \\ {G} \end{array}\right) \omega(k, t)^{G}(1-\omega(k, t))^{(k-G)}\left(1-\Delta t * p_{1}\right)^{G} \\ &=\sum_{G=0}^{k}\left(\begin{array}{c} {k} \\ {G} \end{array}\right)\left(\left(1-\Delta t * p_{1}\right) \omega(k, t)\right)^{G}(1-\omega(k, t))^{(k-G)} \\ &=\left(\left(1-\Delta t * p_{1}\right) \omega(k, t)+(1-\omega(k, t))\right)^{k}=\left(1-\Delta t * p_{1} \omega(k, t)\right)^{k} \\ &=\left(1-\Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)\right)^{k} \end{aligned} $$

In the same way, it can be concluded that the conversion probability of the latent node [TeX:] $$V^{j}(L)$$, neighbored with the recovery [TeX:] $$V(R)$$, is [TeX:] $$p_{L L}^{j}=\left(1-\Delta t * p_{2}\right)^{G}$$. Meanwhile, the infective node [TeX:] $$V^{j}(I)$$ has the same state change probability [TeX:] $$p_{II}^{j}=\left(1-\Delta t * p_{2}\right)^{G}$$ with the recovery node [TeX:] $V(R)$$ around it.

As a result, it can be concluded that the average transfer probability of a node with a degree of k in the latent state or infective state within a period of time [TeX:] $$[t, t+\Delta t]$$ is:

(5)
[TeX:] $$ \begin{aligned} \overline{p}_{L L}(k, t) &=\sum_{G=0}^{k}\left(\begin{array}{c} {k} \\ {G} \end{array}\right) \omega(k, t)^{G}(1-\omega(k, t))^{(k-G)}\left(1-\Delta t * p_{2}\right) \\ &=\sum_{G=0}^{k}\left(\begin{array}{l} {k} \\ {G} \end{array}\right)\left(\left(1-\Delta t * p_{2}\right) \omega(k, t)\right)^{G}(1-\omega(k, t))^{(k-G)} \\ &=\left(\left(1-\Delta t * p_{2}\right) \omega(k, t)+(1-\omega(k, t))\right)^{k}=\left(1-\Delta t * p_{2} \omega(k, t)\right)^{k} \\ &=\left(1-\Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right)\right)^{k} \end{aligned} $$

Therefore, the average transfer probability of a node with a degree of k in the infective state within a period of time [TeX:] $$[t, t+\Delta t]$$ is:

(6)
[TeX:] $$ \begin{aligned} \overline{p}_{II}(k, t) &=\left(1-\Delta t * p_{2} \omega(k, t)\right)^{k} \\ &=\left(1-\Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right)\right)^{k} \end{aligned} $$

3.2 The Density Change Rate of Each State

[TeX:] $$N(k, t)$$ is used to express the total number of the nodes in the network with the degree of k, and [TeX:] $S(k, t), L(k, t), I(k, t), R(k, t)$$ represent the number of the susceptible nodes, the latent nodes, the infective nodes, and the recovery nodes, respectively. The total number of the nodes in the network is as follows:

[TeX:] $$S(k, t)+L(k, t)+I(k, t)+R(k, t)=N(k, t)$$

The changes of the number of the susceptible node with a degree of K in the network in the period of time [TeX:] $$[t, t+\Delta t]$$ are:

(7)
[TeX:] $$ \begin{aligned} S(k, t+\Delta t) &=S(k, t)-S(k, t)\left(1-\overline{p}_{s s}(k, t)\right) \\ &=S(k, t)-S(k, t)\left(1-\left(1-\Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)\right)^{k}\right) \end{aligned} $$

There are:

[TeX:] $$ S(k, t+\Delta t)-S(k, t)=-S(k, t)\left(1-\left(1-\Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)\right)^{k}\right) $$

There are:

(8)
[TeX:] $$ \frac{S(k, t+\Delta t)-S(k, t)}{N(k, t)}=-\frac{S(k, t)}{N(k, t)}\left(1-\left(1-\Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)\right)^{k}\right) $$

By the binomial theorem:

[TeX:] $$ \left(1-\Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)\right)^{k}=\sum_{i=0}^{k}\left(\begin{array}{c} {k} \\ {i} \end{array}\right)\left[-\Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)\right]^{i} $$
[TeX:] $$ \approx 1-k \Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) $$

The higher order power is skipped here. Then:

(9)
[TeX:] $$ \frac{S(k, t+\Delta t)-S(k, t)}{N(k, t)}=-\frac{S(k, t)}{N(k, t)} k \Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) $$

Then:

(10)
[TeX:] $$ \frac{\left(\frac{S(k, t+\Delta t)}{N(k, t)}-\frac{S(k, t)}{N(k, t)}\right)}{\Delta t}=-\frac{S(k, t)}{N(k, t)} k * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) $$

In Eq. (10), when [TeX:] $$\Delta t \rightarrow 0$$, the density change immediate rate of the susceptible node is obtained as follows:

(11)
[TeX:] $$ \frac{\partial \rho^{S}(k, t)}{\partial t}=-\rho^{S}(k, t) k * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) $$

The changes in the number of latent nodes with a degree of k in the network during the period of time [TeX:] $$[t, t+\Delta t]$$ are:

(12)
[TeX:] $$ L(k, t+\Delta t)=L(k, t)+S(k, t)\left(1-\overline{p}_{s s}(k, t)\right)-N(L \rightarrow I)-L(k, t)\left(1-\overline{p}_{L L}(k, t)\right) $$

The changes in the number of infective nodes with a degree of k in the network during the period of time [TeX:] $$[t, t+\Delta t]$$$$ are:

(13)
[TeX:] $$ I(k, t+\Delta t)=I(k, t)+S(k, t)\left(1-\overline{p}_{s s}(k, t)\right)+N(L \rightarrow I)-I(k, t)\left(1-\overline{p}_{I I}(k, t)\right) $$

According to the description and interpretation of the model rules (1) in Figs. 2 and 3, the equations (12) and (13) are further analyzed as follows:

It needs to be explained: [TeX:] $$N(L \rightarrow I)=S(k, t)\left(1-\overline{p}_{S S}(k, t)\right) * \alpha$$

Among them: [TeX:] $$\alpha \in(0,1)$$

[Note: [TeX:] $$N(L \rightarrow I)$$ is the number of latent nodes changing to the number of infective nodes] then: [TeX:] $$L(k, t+\Delta t)-L(k, t)=S(k, t)\left(1-\overline{p}_{s s}(k, t)\right) \beta-L(k, t)\left(1-\overline{p}_{L L}(k, t)\right),(\beta=1-\alpha)$$

Among them :

[TeX:] $$ \begin{aligned} &\left(1-\overline{p}_{s s}(k, t)\right) \approx k \Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)\\ &\left(1-\overline{p}_{L L}(k, t)\right) \approx k \Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) \end{aligned} $$

There are:

(14)
[TeX:] $$ L(k, t+\Delta t)-L(k, t)=S(k, t) k \Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) o^{I}\left(k^{\prime}, t\right) \beta-L(k, t) k \Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) $$

Then:

(15)
[TeX:] $$ \frac{\partial \rho^{L}(k, t)}{\partial t}=\rho^{s}(k, t) k * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) \beta-\rho^{L}(k, t) k * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) $$

Likewise, with the conversion of Eq. (13), we can achieve the following equation:

[TeX:] $$ I(k, t+\Delta t)=I(k, t)+S(k, t)\left(1-\overline{p}_{s s}(k, t)\right)+N(L \rightarrow I)-I(k, t)\left(1-\overline{p}_{II}(k, t)\right) $$

There are:

[TeX:] $$ I(k, t+\Delta t)=I(k, t)+S(k, t)\left(1-\overline{p}_{S S}(k, t)\right) \alpha-I(k, t)\left(1-\overline{p}_{II}(k, t)\right) $$

There are:

(16)
[TeX:] $$ I(k, t+\Delta t)-I(k, t)=S(k, t) k \Delta t * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) \alpha-I(k, t) k \Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) $$

Then:

(17)
[TeX:] $$ \frac{\partial \rho^{\mathrm{I}}(k, t)}{\partial t}=\rho^{s}(k, t) k * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) \alpha-\rho^{I}(k, t) k * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) $$

Finally, the density change immediate rate of the R nodes:

[TeX:] $$ \begin{aligned} R(k, t+\Delta t) &=R(k, t)+L(k, t)\left(1-\overline{p}_{L L}(k, t)\right)+I(k, t)\left(1-\overline{p}_{II}(k, t)\right) \\ &\left(1-\overline{p}_{L L}(k, t)\right) \approx k \Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) \end{aligned} $$

and,

[TeX:] $$ \left(1-\overline{p}_{I I}(k, t)\right) \approx k \Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) $$

Then:

(18)
[TeX:] $$ R(k, t+\Delta t)-R(k, t)=[L(k, t)+I(k, t)] k \Delta t * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) $$

The final results are as follows:

(19)
[TeX:] $$ \frac{\partial \rho^{R}(k, t)}{\partial t}=\left[\rho^{L}(k, t)+\rho^{I}(k, t)\right] k * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right) $$

3.3 Analysis of Dynamic Equation

From the above formulas (11), (15), (17), (19), the dynamic equations set are obtained:

[TeX:] $$ \left\{\begin{array}{l} {\frac{\partial \rho^{S}(k, t)}{\partial t}=-\rho^{s}(k, t) k * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)}{\qquad(11)}\\ {\frac{\partial \rho^{L}(k, t)}{\partial t}=\rho^{s}(k, t) k * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) \beta-\rho^{L}(k, t) k * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right)}\qquad(15) \\ {\frac{\partial \rho^{I}(k, t)}{\partial t}=\rho^{s}(k, t) k * p_{1} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right) \alpha-\rho^{I}(k, t) k * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right)}\qquad(17) \\ {\frac{\partial \rho^{R}(k, t)}{\partial t}=\left[\rho^{L}(k, t)+\rho^{I}(k, t)\right) {k} * p_{2} \sum_{k^{\prime}} p\left(k^{\prime} | k )\rho^{R}\left(k^{\prime}, t\right)\right.}\qquad(19) \end{array}\right. $$

From Eq. (11), we know that the density change rate of susceptible nodes [TeX:] $$\frac{\partial \rho^{s}(k, t)}{\partial t}$$ is positively related to the amount of the susceptible nodes (density) [TeX:] $$\rho^{S}(k, t)$$, the relative amount of infective nodes [TeX:] $$\sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)$$ and k and [TeX:] $$p_{1}$$ in closed system.

On the one hand, the composite parameter [TeX:] $$\sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)$$ reflects the total quantity (density) of the infective nodes, and on the other hand it is related to the topology of the susceptible node. In general, the more the [TeX:] $$\mathrm{V}^{\mathrm{j}}(\mathrm{S})$$ connection nodes, the greater the probability of touching the rumor. The most extreme case is when [TeX:] $$\mathrm{V}^{\mathrm{j}}(\mathrm{S})$$ is an isolated node, it can only be infected by the rumors outside the social system, but not the social network. Hence, the parameter [TeX:] $$p_{1}$$ reaction is the sensitivity of the system node to the rumor, which has a lot of sociality, for instance, it is difficult for the knowledgeable people to easily believe the rumor.

It is known from Eq. (19) that the density change rate of the recovery nodes [TeX:] $$\frac{\partial \rho^{R}(k, t)}{\partial t}$$ is positively related to the amount of the sum of density of the latent nodes [TeX:] $$\rho^{L}(k, t)$$ and the density of the infective nodes [TeX:] $$\rho^{I}(k, t)$$. Also, [TeX:] $$\frac{\partial \rho^{R}(k, t)}{\partial t}$$ is postively related to the total population density of the recovery nodes [TeX:] $$\sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right)$$ and the k and [TeX:] $$p_{2}$$. Also the measure of [TeX:] $$\frac{\partial \rho^{R}(k, t)}{\partial t}$$ is related to the topology of [TeX:] $$\mathrm{V}^{\mathrm{j}}(L)$$ and [TeX:] $$\mathrm{V}^{\mathrm{j}}(I)$$, and is also related to the sensitivity of the rumor and refuting rumors of the social groups' nodes.

The analyses of Eqs. (15) and (17) can also confirm that two density conversion rates are all related to the topological structure of social network interior points and the sensitivity of social network nodes to the rumors and the refuting information. It is necessary to add that the parameters [TeX:] $$\alpha$$ and [TeX:] $$\beta$$ in the system have basic conditions: [TeX:] $$\alpha+\beta=1$$, [TeX:] $$\alpha$$ and [TeX:] $$\beta$$ are rational numbers greater than 0, and less than 1. The value of [TeX:] $$\alpha$$ and [TeX:] $$\beta$$ is subject to a complex psychological factor, subjected to futher study, in future.

4. Data Simulation Analysis

4.1 Network Construction

In this study, a virtual network is built based on WeChat social circles, using simulation. This network is constructed by building node connections through the WeChat social circle in real life, subsequently, four states (i.e., the SLIR) of some nodes are set, the evolution of states are analyzed, accordingly. The details are as follows:

The WeChat is a social network tool which has become popular during recent years. According to Tencent’s mid-2016 performance report, the number of active users in the WeChat is 806 million. Therefore, it is of practical significance to use the social network built by the WeChat users to analyze the evolution simulation of the rumors and the refuting rumors. However, the related information between the WeChat users involves the business confidentiality of the Tencent, the acquisition of one hand data and the experiment violate of the relevant laws and regulations.

Considering the feasibility of the experiment, a small social circle of about 1,000 people are considered to do the experiment and create a virtual social network, through sample data and fictitious data. On this basis, the experiment forms a social network by a moderate fictitious association of nodes, based on a modest fictitious node association based on the ZZ street WeChat social network in the YY town of XX city. The constructed social network for simulation experiments is about 1000 nodes. Finding out some institutions or companies from the ZZ street, through questionnaires and visits, we can basically identify which units have a group of employees on the scale of the WeChat, and which unit employees have frequent WeChat contacts. On the basis, based on the data from more than 700 people, a number of freelancers are added to the social network, to identify the connections between this group and the main group (i.e., 700 people)—relatives, neighbors, trade, etc. The person is numbered by a unit if he belongs to a unit, or, he is numbered by a freelance organization if he does not belong to any unit. On this basis, a WeChat community with 1,002 members is formed.

The survey revealed that this 1,002-member community has a typical small world network characteristics, for example, the employees of the same unit are closely related (large clustering coefficient), hence, the randomly connected freelancer is enough to make the flow of information fast (short path). Following the literature [23-25], the 1,002 network nodes allowed to have a small world characteristics. On the other hand, based on the survey, there are a few people which were especially “get along fine”, and were able to chat with many people (individual nodes have a large degree). In contrast, some extremely introverted personality defects, were often unable to even contact their office colleague (individual nodes have extremely small degrees), which are identified as the characteristics of scale free network.

Therefore, following the theory proposed in the literature [3], in combination with the investigation of several unusual nodes, the degree of association of individual nodes (which is the few “most active individuals” in the survey) is appropriately increased. In contrast, the degree of association of individual nodes (which is the few “extremely introverted individuals” in the survey) is reduced. On this basis, some characteristics of the scale-free degree are formed within the WeChat community (the WeChat social circle). In addition to the above requirements, ensuring that each node is not an isolated node (if it is isolated, it is deleted), the virtual network with 1,002 nodes is finally formed. The minimum nodal degree is 11, the maximum is 390, the average is 80.56, and the median of the nodal degree is 71. The general situation of node association is presented in Fig. 4.

Fig. 4.
The virtual social network of 1,002 nodes of the construction.
4.2 Simulation Data Analysis

The image of the degree distribution function is drawn by the statistics of the correlation degree of the 1,002 nodes, as shown in Fig. 5. As it can be observed in Fig. 5, the nodes and associations based on reality have a certain degree of scale-free characteristics on the basis of a modest fictitious basis.

Subsequently, the rumor is injected into the system, through three different ways: (1) a class of nodes with smaller degrees (the degree is near the minimum); (2) a class of nodes with larger degrees (the degree is near the maximum); and (3) a class of nodes in the vicinity of the average (the degree is near the median). In the next stage, 5 points taken from each of these three types of nodes which were infused the rumors. Consequently, the change of the density of the infective nodes and the latent nodes are obtained, as shown in Fig. 6 (The value of [TeX:] $$P_{3}$$ is 0.6, The value of [TeX:] $$P_{1}$$ is 0.2, The value of [TeX:] $$P_{2}$$ is 0.1).

Fig. 5.
degree distribution of virtual social network nodes.
Fig. 6.
Density variation of the three classes of infective nodes.

The rumors are injected into the system in three categories, when the number of infective nodes reaches the peak. Subsequently, the information of refuting the rumor was injected into three types of nodes. Fig. 7 represents the changes in the density of recovery nodes after the three types of nodes are injected.

Fig. 7.
Density variation of the three classes of recovery nodes in a virtual social network

The above two types of experiments, regardless of being the injection of rumors or the injection of refuting-rumor-information, suggest that the evolution of the system after injection is significantly associated with the degree distribution of the injected node. The rumor information is injected from some nodes with node degree greater than 350, and the density growth of infective nodes is much larger. While, the injection from other nodes with node degree less than 20, hence, the density growth of infective nodes is much slower. In addition, the refuting-rumor-information is injected from some nodes with larger value degree, and the density growth of recovery nodes is much faster. While, the injection from other nodes with smaller value degree, hence, the density growth of recovery nodes is much slower.

The parameter [TeX:] $$\sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{I}\left(k^{\prime}, t\right)$$ of Eq. (11) explains the relationship between the evolution of the system and the topology of the node, after the rumor injection. The same applies to the parameter [TeX:] $$\sum_{k^{\prime}} p\left(k^{\prime} | k\right) \rho^{R}\left(k^{\prime}, t\right)$$ in the Eqs. (15), (17), and (19). A clear conclusion is that the selection of information injection nodes of the rumor and the refuting-rumor-information affect the subsequent evolution of the various nodes in the system.

The situation can be more complicated, when the rumor and the refuting-rumor-information are injected into the system, at the same time. The other possibility is that to select a number of nodes near the median degree to infuse the rumor and the refuting-rumor-information, and select the average of the 20 trials by many experiments. The density variations of the four types of nodes in SLIR are presented in Fig. 8 (the value of [TeX:] $$P_{3}$$ is 0.5, the value of [TeX:] $$P_{1}$$ is 0.2, the value of [TeX:] $$P_{2}$$ is 0.1).

The result of the experiment identify that the influence of the early spread of the rumor is obviously higher than the influence of the refuting-rumor-information, in case of injecting the rumor and the refuting-rumor-information at the same time (the total density of the infective nodes and the latent nodes is obviously larger than the density of the recovery nodes, at first). With the expansion of the transmission range and the increase of the total density of the affected nodes, the influence of the recovery nodes is bigger and bigger. Finally, the recovery nodes occupy the absolute superiority of the whole system. Surprisingly, always a part of the susceptible nodes in the system remains uninfected. This is probably because the rumors and the refuting-rumor-information are almost simultaneously injected, some nodes with small degrees are earlier to be immunized (being recovery nodes), and some other nodes that are closely related to it, may have no chance of being infected by rumors (never being infective nodes).

Fig. 8.
The average density changes after the injection of rumor and refuting-rumor-information.

5. Summary and Prospect

In this research, an evolution model of rumor spreading is subjected to study, based on WeChat social circle. The evolution of rumor spreading is depicted through the propagation dynamics equation, where, the evolution of the four nodes in the SLIR model is subject to the topology of the nodes, the density of the four kinds of nodes, the information sensitivity of the nodes ([TeX:] $$P_{1}, P_{2}$$), and some psychological factors ([TeX:] $$P_{3}$$). The relative factors, about the evolution of the four states, are obtained through different equations. Subsequently, the four states in the form of mutual cause and effect, form a dynamic evolution process. This study constructs a social network through a realistic WeChat community, and simulates the spread of the rumors in the WeChat’s social circle, based on this social network. The simulation results show that the spreading efficiency of the nodes with various node degrees is significantly different, which also confirms the early deduction from the dynamic equation. In addition, the node degree affects the efficiency of evolution. However, the final density of I nodes or the R nodes are not significantly affected. This result can be related to the total number of nodes in the simulation network and the topology of the whole network.

In the era of Internet and mobile Internet, the dissemination of information in social networks, including the rumors, has become a norm. Current literature review suggest significant studies are performed in regards to the rumor spreading, However, this research stands out from the the existing literature, due to providing an evolution model, based on the reality oriented simulation. The construction of the rumor evolution model SLIR and the construction and deduction of propagation dynamics equation in this study provide a significant reference for the future studies on the rumor spreading model. The simulation realization facing the real background, especially for the injection of the rumors or the refuting-rumor-information, provides a valuable method for tracking and controlling of the rumors in the real-world scenarios. The background of this article is mainly the social network constructed by the WeChat social circle, which is obviously different from the social circle based on the portal platform, as well as the micro-blog dating for mutual attention and information display. Effective equation deduction solves the problem of logical trust and provides the certainty in research. In addition, simulation of realistic social networks is more believable for authenticity. This information interaction on the WeChat platform, in a larger scale, is considered as a limitation for this kind of research, due to the privacy of citizens’ personal information and the highly confidential information of the Tencent Inc. Subsequent research can be carried out on other platforms, such as Facebook or a new social platform in China, called “Douyin”. Furthermore, some parameters of this model, such as [TeX:] $$p_{1}, p_{2}, \alpha$$ and [TeX:] $$\beta$$,can be studied in depth. In particular, studying the psychological factors, [TeX:] $$\alpha$$ and [TeX:] $$\beta$$, helps to identify the influence of various factors on them, through comprehensive methods of statistics and social psychology, such as the individual age, the education, the work experience, and the work nature influence.

Acknowledgement

This research was supported by Zhejiang Province philosophy and social science planning project, China (No. 16NDJC274YB) and Key Research Institute of Philosophy and Social Sciences of Zhejiang Province—Modern Port Service Industry and Creative Culture Research Center.

Biography

Lubang Wang
https://orcid.org/0000-0001-9104-7258

He was born in Anhui Province, China, February 16, 1974. Currently, he is an associate professor in Department of Electronic Commerce, College of Logistics and E-commerce, Zhejiang Wanli University, China. He received M.S. degree in School of Information from Sun Yat-sen University in 2004, and in July 2011, he received Ph.D. in School of Management from University of Shanghai for Science and Technology. His current research interests include knowledge management and public opinion spreading.

Biography

Yue Guo
https://orcid.org/0000-0002-9472-9171

He received Ph.D. degree in School of Management from Wuhan University of Technology in 2011. He is a professor in the School of Economics and Management, Ningbo University of Technology, China. His current research interests include information management and enterprise management.

References

  • 1 W. C. Hou, Y. Pan, D. Che, "Utilizing fragmented bandwidth in a staggered striping multimedia system," Journal of Information Processing Systems, vol. 4, no. 1, pp. 1-8, 2008.doi:[[[10.3745/JIPS.2008.4.1.001]]]
  • 2 D. J. Watts, S. H. Strogatz, "Collective dynamics of small-world networks," Nature, vol. 393, no. 6684, pp. 440-442, 1998.custom:[[[-]]]
  • 3 A. L. Barabasi, R. Albert, "Emergence of scaling in random networks," Science, vol. 286, no. 5439, pp. 509-512, 1999.doi:[[[10.1515/9781400841356.349]]]
  • 4 R. Kumar, J. Novak, A. Tomkins, "Structure and evolution of online social networks," in Proceeding of the 12th ACM SIGKDD International Conference on Knowledge Discovery And Data Mining, Philadelphia, PA, 2006;pp. 611-617. custom:[[[-]]]
  • 5 H. Hu, X. Wang, "Evolution of a large online social network," Physics Letters A, vol. 373, no. 12-13, pp. 1p2p3p1p2p   1105-1110, 2009.doi:[[[10.1016/j.physleta.2009.02.004]]]
  • 6 A. Mislove, M. Marcon, K. P. Gummadi, P. Druschel, B. Bhattacharjee, "Measurement and analysis of online social networks," in Proceedings of the 7th ACM SIGCOMM Conference on Internet Measurement, San Diego, CA, 2007;pp. 29-42. custom:[[[-]]]
  • 7 H. Chun, H. Kwak, Y. Eom, Y. Ahn, S. Moon, H. Jeong, "Comparison of online social relations in terms of volume vs. interaction: a case study of cyworld," in Proceedings of the 8th ACM SIGCOMM Conference on Internet Measurement, V ouliagmeni, Greece, 2008;pp. 57-70. custom:[[[-]]]
  • 8 J. N. Kapeerer, Rumor: The W orld's Oldest Media, China: Shanghai People’s Publishing House, Shanghai, 2008.custom:[[[-]]]
  • 9 R. K. Garrett, "Troubling consequences of online political rumoring," Human Communication Research, vol. 37, no. 2, pp. 255-274, 2011.doi:[[[10.1111/j.1468-2958.2010.01401.x]]]
  • 10 G. L. Cohen, J. Aronson, C. M. Steele, "When beliefs yield to evidence: reducing biased evaluation by affirming the self," Personality and Social Psychology Bulletin, vol. 26, no. 9, pp. 1151-1164, 2000.doi:[[[10.1177/01461672002611011]]]
  • 11 D. J. Daley, D. G. Kendall, "Epidemics and rumours," Nature, vol. 204, no. 4963, pp. 1118-1118, 1964.custom:[[[-]]]
  • 12 D. Maki, M. Thomson, Mathematical Models and Applications, NJ: Prentice-Hall, Englewood Cliffs, 1973.custom:[[[-]]]
  • 13 B. Doerr, M. Fouz, T. Friedrich, "Social networks spread rumors in sublogarithmic time," Electronic Notes in Discrete Mathematics, vol. 38, pp. 303-308, 2011.doi:[[[10.1016/j.endm.2011.09.050]]]
  • 14 D. H. Zanette, "Critical behavior of propagation on small-world networks," Physical Review E, vol. 64, no. 5, 2001.custom:[[[-]]]
  • 15 D. H. Zanette, "Dynamics of rumor propagation on small-world networks," Physical Review E, vol. 65, no. 4, 2002.doi:[[[10.1103/physreve.65.041908]]]
  • 16 Y. Moreno, M. Nekovee, A. Pacheco, "Dynamics of rumor spreading in complex networks," Physical Review E, vol. 69, no. 6, 2004.doi:[[[10.1103/PhysRevE.69.066130]]]
  • 17 Z. F. Pan, X. F. Wang, X. Li, "Simulation investigation on rumor spreading on scale-free network with tunable clustering," Journal of System Simulation, vol. 18, no. 8, pp. 2346-2348, 2006.custom:[[[-]]]
  • 18 A. Java, P. Kolari, T. Finin, and T. Oates, 2006;, https://ebiquity.umbc.edu/paper/abstract/id/300/Modeling-the-Spread-of-Influence-on-the-Blogosphere
  • 19 Z. Y. Zheng, F. Guo, Z. F. Wang, D. Li, "Study on microblog propagation model based on analysis of user behavior," Computer Science, vol. 43, no. 12, pp. 41-46, 2016.custom:[[[-]]]
  • 20 Z. J. Song, R. Shi, J. Wang, "Influence of authoritative information release on rumors microblog spreading in emergency events," Journal of Intelligence, vol. 35, no. 12, pp. 41-47, 2016.custom:[[[-]]]
  • 21 D. Lu, J. L. Guo, H. Y. He, "Based on the double S model research in WeChat rumors spreading," Mathematics in Practice and Theory, vol. 2017, no. 16, pp. 157-163, 2017.custom:[[[-]]]
  • 22 Y. X. Luo, "A comparative study of spread of Weibo and WeChat rumors from the perspective of public opinion," Journalism Knowledge, vol. 2014, no. 4, pp. 6-8, 2014.custom:[[[-]]]
  • 23 W. Dou, Y. Jia, H. M. Wang, W. Q. Song, P. Zou, "A P2P approach for global computing," in Proceedings International Parallel and Distributed Processing Symposium, Nice, France, 2003;custom:[[[-]]]
  • 24 Z. Liu, W. Don, W. M. Zhang, P. Zou, "Paradropper: a general-purpose global computing environment built on peer-to-peer overlay network," in Proceedings of the 23rd International Conference on Distributed Computing Systems W orkshops, Providence, RI, 2003;pp. 954-957. custom:[[[-]]]
  • 25 W. Dou, H. M. Wang, Y. Jia, P. Zou, "A rumor-spreading analog on unstructured P2P broadcast mechanism," Journal of Computer Research and Development, vol. 2004, no. 9, pp. 1460-1465, 2004.custom:[[[-]]]
The relation figure of node state conversion.
Explanatory drawing of node S, L, I state conversion illustration.
Revised drawing of node state transformation relationship.
The virtual social network of 1,002 nodes of the construction.
degree distribution of virtual social network nodes.
Density variation of the three classes of infective nodes.
Density variation of the three classes of recovery nodes in a virtual social network
The average density changes after the injection of rumor and refuting-rumor-information.