Abstract

Information spreading on multiplex networks has been investigated widely. For multiplex networks, the relations of each layer possess different extents of intimacy, which can be described as weighted multiplex networks. Nevertheless, the effect of weighted multiplex network structures on information spreading has not been analyzed comprehensively. We herein propose an information spreading model on a weighted multiplex network. Then, we develop an edge-weight-based compartmental theory to describe the spreading dynamics. We discover that under any adoption threshold of two subnetworks, reducing weight distribution heterogeneity does not alter the growth pattern of the final adoption size versus information transmission probability while accelerating information spreading. For fixed weight distribution, the growth pattern changes with the heterogeneous of degree distribution. There is a critical initial seed size, below which no global information outbreak can occur. Extensive numerical simulations affirm that the theoretical predictions agree well with the numerical results.

1. Introduction

The Internet [1, 2], web 2.0 [3], and intelligent mobile terminals [4] have completely transformed people’s daily lives; these infrastructures have assisted people in constructing virtual social networks [58], which map people’s connections in real life into virtual cyberspaces and extend people’s social scope significantly, not only with schoolmates, colleagues, and family members, but also with strangers in different countries. Previously, various types of social networks such as Facebook [9], Twitter [10], Weibo [11], and Wechat [12] have been used. Social networks are important in information transmission between users [13], such as commodity recommendation [14], news spreading [15, 16], advertisement diffusion [17, 18], health behavior adoption [19], discovery of new friends [20], and social innovation propagation [21]. In addition, understanding the mechanism of information spreading in social networks can facilitate in the design of antivirus strategies [2224], control rumors [25], avoid economic risks [26, 27], restrain social unrest [28], etc. Because of important values in economic activities and engineering applications, social information spreading has attracted significant attention [2832].

Based on complex network theory [33, 34], regarding a user as a node and a connective relation as an edge [35], researchers have mined various underlying factors to investigate the information spreading mechanism and unveil the fundamental laws. Wide-ranging factors have been investigated, such as the initial seed size [36], clustering coefficient [37], community structure [38], temporal network [39], role of synergy [40, 41], heterogeneous adoption thresholds [42], and limited imitation [4345]. Additionally, information spreading on a multiplex network [4648] has aroused the interest of researchers who perform numerous investigations about communication channel alternations on multiplex networks [49], heterogeneous behavioral adoption on multiplex networks [16], opinion competition on multiplex networks [50], effect of multiple social networks on user awareness [51], social support for suppressing epidemics on multiplex networks [52], the interplay of social influence in multiplex networks [53], replicator dynamics on multiplex networks [54], etc.

Moreover, in social networks, because of multiple confirmations of the credibility and legitimacy of information, information spreading exhibits the effect of social reinforcement [55, 56]. Repeated transmission of the same information from one neighbor cannot enhance the credibility nor facilitate adoption. Therefore, social reinforcement derived from the memory of nonredundant information has been investigated [57], which should be considered in information spreading. Apart from the abovementioned factors, edges represent the connective relations between users and are always modeled in a binary state [58]. In real social networks, different types of individuals, such as colleagues, schoolmates, family members, and strangers may have intimate or distant relations. Even for the same type of individuals, they may exhibit different extents of intimacy. Therefore, similar to epidemic spreading on a weighted network [59], the relations in a social network should be reasonably modeled as weighted edges [60, 61]. Noticeably, compared with a one-layered complex network, research on information spreading in multiplex networks is more complicated and challenging. Apart from the intraheterogeneity of weights between edges in the same layer, interheterogeneity of weight distributions between two layers exists as well. Unfortunately, studies performed hitherto lack profound analysis regarding information spreading on weighted multiplex networks with social reinforcement derived from the memory of nonredundant information.

Herein, in terms of social reinforcement effect generated from the memory of nonredundant information and heterogeneous weight distributions in two subnetworks, we propose an information spreading model on a weighted multiplex network. To further understand the intrinsic mechanism of information spreading and quantitatively predict the effects of spreading in time evolution and final steady state, we develop an edge-weight-based compartmental theory. Extensive simulations based on our proposed numerical method suggest that theoretical predictions agree well with numerical simulations. Combining the theoretical predictions and numerical simulations, we discovered that reducing weight distribution heterogeneity does not alter the growth pattern of the final adoption size but can accelerate information spreading and promote the outbreak of information spreading. Furthermore, we discovered that a critical seed size determines the global information outbreak, above which reducing the weight distribution heterogeneity materially affects the facilitation of information spreading. Additionally, reducing the degree distribution heterogeneity can alter the growth pattern of the final adoption size .

The remainder of the paper is organized as follows: In Section 2, we describe the information spreading model on a weighted multiplex network. In Section 3, we propose an edge-weight-based compartmental theory to analyze information spreading on a weighted multiplex network. In Section 4, the numerical simulation method and parameter settings are introduced. In Section 5, we compare the theoretical predictions with the numerical simulation results and unveil the effect mechanism of the weighted multiplex network structure on information spreading. In Section 6, conclusions are presented.

2. Model Descriptions

For investigating information spreading on a weighted multiplex social network, we constructed an information spreading model based on an N-node two-layer weighted multiplex network, in which layers A and B represent two different independent social networks (illustrated in Figure 1). Correlated nodes in two layers indicate that they are the same individual, and edges between nodes in each subnetwork represent interindividual social relations. To eliminate the unnecessary intradegree-degree correlation, the uncorrelated configuration model [62] was manipulated to construct the multiplex network following two independent degree distributions and , in which and indicate the degrees of node i in layer A and B, respectively. The degree vector of node i is set as . In this paper, we assume that there is no correlation between the degree and weight distributions and the joint degree distribution is . Self-loop and multiple edges are forbidden. Regarding social reinforcement, we assume that each individual has different adoption willingness in different social subnetworks. Consequently, node i, representing an arbitrary individual, has adoption thresholds and in layers A and B, respectively, and a smaller adoption threshold suggests more considerable information adoption willingness. Moreover, the intimate relations between individuals typically hold different strengths in the same social subnetwork and in two distinctive subnetworks. Therefore, we further construct the weighted multiplex network with two dependent edge-weight distributions, and . Subsequently, if an edge with weight exists between nodes i and j in the layer , individual i in the subnetwork X transmits the information to j with probabilitywhere is the unit transmission probability in the subnetwork . Given , enlarging can monotonically increase , which coincides with the positive correlation between relationship weight and transmission preference.

In real social networks, individuals always trustfully adopt information after receiving a certain number of it from distinctive neighbors, owing to social reinforcement originating from the memory of nonredundant information transmission [57]. Accordingly, we applied the generalized susceptible-adopted-recovered (SAR) model to describe the state of an individual at each time step in the information spreading process. At the same time step, an individual can only remain in one state: S-state (susceptible), A-state (adopted), or R-state (recovered). The S-state implies that an individual has not adopted the information. The A-state indicates that an individual has adopted the information and can transmit the information to neighbors. The R-state means that an individual has lost interest in the information and stops diffusion.

Furthermore, the detailed information spreading process is as follows: Initially, randomly selecting a fraction of nodes as initial seeds (in the A-state). At each time step, each A-state node i attempts to transmit the information to its S-state neighbor j simultaneously through an edge of weight with probability in layer A and through an edge of weight with probability in layer B according to equation (1). Because of nonredundant information transmission, once an edge between nodes i and j in layer A or B successfully delivers the information from i to j, such an edge is not allowed for repeated information delivery. An adopted node, however, can perform many attempts until it turns into the R-state. If the S-state node j in the layer successfully receives one piece of information from the A-state neighbor node j, its accumulative number of information pieces in the layer will increase by one, i.e., . Then, node j will compare the accumulative number with its adoption threshold . In reality, for confirming the legitimacy and credibility of the information, an individual can adopt the information only if he/she corroborates with the information in both subnetworks. Therefore, the susceptible node j can enter the A-state when and co-occur in both layers A and B; otherwise, it remains in the S-state. Regarding social reinforcement, because the accumulative number of information pieces determines the information adoption, the information spreading process is characterized by the non-Markovian effect. At the same time step, after information transmission, the A-state node j will lose interest in the information and enter the R-state with probability γ, i.e., it will stop diffusing information in both layers A and B. Once the node changes into the R-state, it will retreat from further information spreading. Finally, the disappearance of all A-state nodes signals the termination of the spreading dynamics.

3. Edge-Weight-Based Compartmental Theory

To include the social reinforcement effect, the mean-filed theory [63] and an advanced approach [64] can be used. Inspired by References [57, 65, 66], an edge-weight-based compartmental theory is conceived to analyze the spreading process with theoretically strong correlation among the three states of nodes in a multiplex the network featured by strong weight heterogeneity. The theory assumes that the information spreading is on large, sparse, uncorrelated, and local tree-like weighted multiplex networks.

In the edge-weight-based compartmental theory, a node in the cavity state [67] implies that it can only receive information from but cannot transmit information to its neighbors. We denote as the probability that the information has not been passed to the S-state neighbors alongside a randomly selected edge in layer X by time t. Furthermore, in a weight multiplex network, we denote as the probability that a node has not been informed through a randomly selected edge with weight in layer X by time t. Following weight distribution in layer A and in layer B, we can obtain

An arbitrary S-state node i with degree vector receives and pieces of information by time t with the probabilitiesrespectively.

As mentioned in Section 2, an S-state node with degree vector can adopt the information only if the cumulative numbers and occur simultaneously; inversely, the S-state node having received and pieces of information by time t will remain in the S-state with probability

We define the probability that a susceptible node with any degree in the layer remaining in the S-state if it receives pieces of information by time t as

Furthermore, the fraction of nodes in the S-state at time t can be expressed as

If we wish to obtain , we should solve , which depends on in a weighted network. Therefore, we focus on at first. According to [57], in layer X, the neighbor connected by an edge with weight to a node in the cavity state can be in the S-state, A-state, or R-state; therefore, comprises the following three parts:where , , and indicate the probabilities that a neighbor connected by an edge with weight to a node in the cavity state remains in the S-, A-, and R-states, respectively, and has not transmitted information by time t.

If neighbor j with degree vector is linked by an edge of weight to the node i in the cavity state in layer , then neighbor j can only receive information from other edges in layer . Therefore, neighbor j receives pieces of information in layer A with probability and, with no constraint in layer B, receives pieces of information in layer B with probability . Subsequently, neighbor j remains in the S-state in layer A with probability

Analogously, neighbor j is connected with node i by an edge with weight in layer B and remains in the S-state with probability

Accordingly, given the joint degree distribution , in layer X, an edge with weight connects to an S-state neighbor j with probabilitywhere the expression denotes the probability an edge connecting to a neighbor with degree in layer X.

Next, we proceed with the evolution of and in the layer . If, through an edge with weight , the information is successfully transmitted to the neighbor in layer X with probability , the change in is calculated as

If neighbor j in the A-state has not transmitted information to node i through an edge of weight with probability in layer X but enters into the R-state with probability γ at time t, the value of will increase, and thus the dynamical differential value of can be obtained by

The initial conditions are and ; therefore, combining equation (12) with (13), we can obtain

After inserting equations (11) and (14) into equation (8), we can obtain

In equation (15), substituting with according to equation (8), we further deduce

When , , and we have the steady-state value as

Inserting into equations (2) and (3), we thus obtain

For a convenient expression, we rewrite asand as

A discontinuous growth pattern will occur [68] when equation (19) is tangent to equation (20) with and . Moreover, critical conditions can be acquired from the following equation:

From the theoretical analysis above, the general solutions for adoption thresholds and can be derived, but the theoretical analysis itself may not be easily comprehended. Therefore, an intuitive illustration exemplifies the general theory with a particular case, i.e., . In this case, equations (9) and (10) are rewritten as

Owing to independence, we obtain and substitute (23) and (24) into equation (17) as follows:where represents the generation function of degree distribution , and represents the generation function of excess degree distribution . Additionally, the generation functions of degree distribution are expressed as follows:where . Subsequently, substituting equations (23) and (24) into equation (21), we can deduce the critical conditions.

4. Numerical Method

Based on Erdös–Rényi (ER) [69] and the scale-free (SF) [62] random network model, extensive numerical simulations on artificial two-layer networks were performed to verify the theoretical analysis above. Without loss of generality, we set the network size and mean-degree . For convenience, the information transmission unit probability in the two layers is set as  =   = β. Additionally, the weight distribution of the network layer follows with and mean weight . At the steady state in which no node in the A-state exists, we compute the spreading range by measuring the fraction of nodes in the R-state . The simulation results are obtained by averaging over at least independent dynamical realizations on 30 artificial networks.

To obtain the critical condition from the simulations, the relative variance χ [70] is applied in our numerical verification, using the following method:where is the ensemble average, and and are the first and secondary moments of , respectively. The peaks of the curve of relative variance χ versus the information transmission probability correspond to the critical points.

5. Results and Discussion

In reality, from the perspective of a degree distribution model, layers in some multiplex social networks exhibit the homogeneous creation mechanism, such as MSN versus ICQ, but layers in other multiplex social networks may exhibit heterogeneous creation mechanisms, such as ICQ versus Facebook. Therefore, from the aspects of homogeneous and heterogeneous degree distributions of two layers, the effects of weight distribution heterogeneity on information spreading in multiplex social networks is discussed in this section in terms of the parameters in Section 4.

5.1. ER-ER Weighted Multiplex Networks

First, we investigate information spreading on a two-layered weighted ER-ER network with Poisson degree distribution in layer A and in layer B.

We first explore the time evolutions of node densities in three states (S-state, A-state, and R-state) under different weight distribution exponents and , as shown in Figure 2. Under the same information transmission unit probability β but different adoption thresholds and seed sizes, Figure 2 shows the plots of curves with in the top panel and in the bottom panel. Regardless of the top or the bottom panel, with time t, the size of decreases gradually from to 0, changes from to 0 eventually, and indicating the final adoption size increases gradually from 0 to 1 at the end. With the increase in the weight distribution exponent from in Figures 2(a) and 2(e), in Figures 2(b) and 2(f), and in Figures 2(c) and 2(g), the evolution time consumed decreases from 19 to 13 in the top panel and 11 to 9 in the bottom panel. Additionally, the final adoption size augments at the same time t (see the comparison in Figures 2(d) and 2(h)). The phenomena occurring during the time evolution suggest that increasing the weight distribution exponent (reducing weight distribution heterogeneity) can facilitate information spreading in a weighted multiplex network.

Next, given the adoption thresholds , , and seed size , we investigate the growth pattern of the final adoption size versus information transmission unit probability β. In Figure 3(a), with , , and , increasing the weight distribution exponent such as from to does not change the growth pattern of the final adoption size versus β but reduces the critical information transmission unit probability , which can be numerically captured by the peaks of relative variance χ in Figure 3(e) with the same settings of Figure 3(a). Figure 3(b), with different , exhibits the same pattern as Figure 3(a), with the peaks of relative variance χ in Figure 3(f). Even enlarging the seed size () in Figure 3(c) with the same adoption threshold of Figure 3(b), the final adoption size versus β increases in the same pattern regardless of the increase in weight distribution. Besides, keeping the seed size as the same of Figure 3(c) (), Figure 3(d), with further increase of the adoption threshold , also shows similar phenomena as those in Figures 3(a)3(c), verified by the peaks of χ in Figure 3(h). The above phenomena suggest that increasing the weight distribution exponent (decreasing weight distribution heterogeneity) will not alter the growth pattern of versus β in any pair of adoption thresholds but can promote the outbreak of information spreading. Noticeably, the theoretical solutions (lines) coincide well with the numerical simulations (symbols).

As shown in Figure 3, the initial seed size affects the information spreading dynamics. We investigate the effect of seed size on information spreading on a weighted multiplex social network, as shown in Figure 4. In Figure 4(a) (), increasing the seed size from to not only changes the growth pattern from discontinuous to continuous, but also decreases the critical probability to promote the outbreak of information spreading. Subsequently, if the weight distribution exponent is enlarged, such as in Figure 4(d) and in Figure 4(g), increasing the seed size still advances information propagation. Comparing Figure 4(a), 4(d), and 4(g) vertically, altering the initial fraction can change the growth pattern of versus β but cannot change the fact that increasing the weight distribution exponent (reducing weight distribution heterogeneity) retains the growth pattern while facilitating the outbreak of information spreading. Moreover, in a horizontal comparison, although enlarging the adoption thresholds, such as in Figures 4(b), 4(e), and 4(h), or in Figure 4(c), 4(f), and 4(i) renders information adoption more difficult, it is clear that increasing still transforms the growth pattern of from discontinuous to continuous phase transition and promotes the outbreak of information spreading, as in the case of . Additionally, changes in adoption thresholds and seed size do not qualitatively alter phenomena.

For a more comprehensive study, we obtained the results of the final adoption size on the parameter plane, as illustrated in Figure 5. Here, we set the subgraphs in column 1 with ; in column 2 with ; in column 3 with ; in row 1 with ; in row 2 with ; and in row 3 with . In the following, we discuss the findings based on Figure 5 in detail. In Figure 5(a) (), a critical exists, below which at any information transmission unit probability β no outbreak of information spreading will occur and above which outbreak and global adoption will occur. In a vertical comparison, we retain the adoption thresholds in Figures 5(a), 5(d), and 5(h) and further increase the weight distribution exponents in Figure 5(d) () and Figure 5(h) (). The results indicate that the three subgraphs have the same critical seed size , below which increasing the weight distribution exponent (reducing the weight distribution heterogeneity) still cannot induce the outbreak of information spreading, but above which increasing the weight distribution exponent can accelerate the outbreak at a smaller information transmission unit probability β. Moreover, in a horizontal comparison, we retain the weight distribution exponent in Figures 5(a), 5(b), and 5(c) with and increase the adoption thresholds in Figure 5(b) () and Figure 5(c) (). Under these circumstances, we discovered that the critical seed size is enlarged from 0.003 to 0.06, which implies increasing the difficulty in information adoption. Analogously, other subgraphs exhibit similar findings. The phenomena in Figure 5 suggest that a critical seed size exists, above which outbreak of information spreading and global adoption will occur, and the weight distribution exponent (reducing the weight distribution heterogeneity) is essential for facilitating information spreading.

5.2. ER-SF Weighted Multiplex Networks

For cases of heterogeneous degree distributions, we investigate information spreading on a two-layered weighted ER-SF network with Poisson degree distribution in layer A and with power-law degree distribution in layer B, where and parameter denotes the degree exponent of layer B. Additionally, we set the minimum degree and maximum degree .

On the weighted ER-SF networks of distinct degree distribution exponents, the growth pattern is shown in Figure 6, with the common setting in all subgraphs. In Figure 6(a) with degree distribution exponent , we discovered that enlarging the weight distribution exponent, such as from to , can increase the final adoption size at the same β and decrease the critical probability to promote the outbreak of information spreading. Moreover, we proceed to increase the degree distribution exponent (reduce the degree distribution heterogeneity) in Figure 6(b) () and Figure 6(c) (). Consequently, from Figure 6(a) to 6(c), given and , the growth pattern changes from continuous to the discontinuous growth pattern. However, at any given , changing the weight distribution exponent will not alter the growth pattern of . The phenomena in Figure 6 suggest that changing the degree distribution heterogeneity will not qualitatively alter the fact that reducing weight distribution heterogeneity can promote information spreading.

6. Conclusions

With the development of the Internet and mobile devices, social networks have prevailed in people’s lives. Typically, an individual does not only use one social network, but multiple social networks simultaneously. Moreover, the information does not only diffuse in one social network but also spreads simultaneously on several social networks, e.g., information spreads on Twitter and Facebook simultaneously. Because of different extents of intimacy among friends, colleagues, and family members, connections in social networks should be modeled with different weights. As such, a multiplex social network should be a multilayer weighted network. Therefore, information spreading on a multiplex social network should be investigated on a multilayer weighted network. Unfortunately, the mechanism of information spreading on a weighted multiplex network has not been investigated comprehensively.

Besides, social information spreading possesses the characteristics of social reinforcement derived from the memory of nonredundant information. Herein, in terms of social reinforcement and weight distribution, we first proposed a general information spreading model on a weighted multiplex network to depict the dynamics of information spreading on a two-layered weighted network. Subsequently, we conceived a generalized edge-weight-based compartmental theory to analyze the mechanism of information spreading on a weighted multiplex network with the effect of the memory of nonredundant information. Based on the numerical method, extensive simulations confirmed that the theoretical predictions coincided well with the simulation results.

Combining numerical simulations with theoretical analyses, we discovered that under any adoption threshold of two layers, reducing weight distribution heterogeneity did not alter the growth pattern of the final adoption size but could accelerate information spreading and promote the outbreak of information spreading. Furthermore, we discovered that a critical seed size existed, below which no outbreak of information spreading occurred and above which reducing weight distribution heterogeneity was meaningful to the facilitation of information spreading. Additionally, we discovered that changing degree distribution heterogeneity could transform the growth pattern but did not qualitatively affect the promotion of reducing weight distribution heterogeneity to information spreading. Above all, increasing the weight distribution exponent can decrease critical transmission probability, without the influence of the degree distribution exponent.

In this study, we unveiled the effect of weighted multiplex network structures on information spreading. This study can inspire further studies related to the design of optimized strategies to control information spreading on multiplex social networks. Additionally, more accurate theories to depict the social contagion on real two-layer coupled networks should be further explored.

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 supported by the National Natural Science Foundation of China (nos. 61602048, and 61673086) and the Fundamental Research Funds of the Central Universities.