Abstract

Information spreading dynamics on temporal networks have attracted significant attention in the field of network science. Extensive real-data analyses revealed that network memory widely exists in the temporal network. This paper proposes a mathematical model to describe the information spreading dynamics with the network memory effect. We develop a Markovian approach to describe the model. Using the Monte Carlo simulation method, we find that network memory may suppress and promote the information spreading dynamics, which depends on the degree heterogeneity and fraction of bigots. The network memory effect suppresses the information spreading for small information transmission probability. The opposite situation happens for large value of information transmission probability. Moreover, network memory effect may benefit the information spreading, which depends on the degree heterogeneity of the activity-driven network. Our results presented in this paper help us understand the spreading dynamics on temporal networks.

1. Introduction

Extensive real-data analyses revealed that social network exhibits strong temporal properties [13], i.e., the edges and nodes do not always exist at any time, and may vary with time. For instance, in scientist collaborative networks, two researchers may collaborate to publish a paper but rarely collaborate at every time step [49]. Besides, two researchers may build their first collaboration. Another example is that in the transportation network, two cities may build expressways and high-speed trains. Thus, a new edge is added. The emergence and disappearance of edges and nodes widely exist for the online social network due to the login or log out of the online platforms. Therefore, the temporal network is a widely used method to describe the social network, in which nodes represent individuals and edges stand for their relationships.

For the information spreading on temporal networks, researchers from different disciplines made great contributions [1019]. In what follows, we first review the progress of information spreading on social network. Different from the static networks, i.e., the network topology does not change with time, researchers found some important results [15, 2023]. When the information is spreading on the static networks, scholars found that the existence of some hubs may eliminate the threshold point [24, 25]. Specifically, any values of information transmission probability can trigger the information spreading on social networks. Based on these results, we can understand why information can always spread on social platforms. Researchers further revealed that the network community, clustering, and degree-degree correlations could alter the spreading dynamics of information [26, 27]. In reality, sharing a piece of information is risky, and thus affirming its reality and reliability is fatal. Researchers used the threshold-based model to include this factor in the spreading dynamics, such as the Watts threshold model and other generalized models. For that threshold-based information spreading model, the phase transition of the dynamical system is always discontinuous, i.e., first-order phase transition [28]. Wang et al. [29, 30] proposed nonredundant information spreading dynamics and revealed a transition between the continuous and discontinuous transition in the system.

When researchers studied the information spreading dynamics on temporal networks, scholars found that network temporality can suppress or promote the information spreading [1, 20, 21, 3137]. Xue et al. proposed a mathematical model to describe a heterogeneous population, where a fraction of nodes adopts a complex contagion. They used a Markovian approach to describe the spreading dynamics. They found that the promotion or suppression of the network temporality is determined by the heterogeneities of population and degree distribution. Scholtes et al. [38] found that the non-Markovian temporal network may speed up or slow down the information spreading. Wang et al. [39] proposed a heuristic immunization strategy for information spreading and demonstrated the effectiveness of this strategy in real-world data.

For the temporal network, an essential factor is the network memory [38, 4042], which means that edges that existed in the current time step may already occur in previous time steps. Sun et al. [43] revealed that network memory inhibits the spreading process for SIR models, in which the epidemic threshold is enlarged while the spreading size decreases. How memory affects the information spreading dynamics is an important question. To address this, we propose a mathematical model on temporal networks with memory. Then, we develop a Markovian theory for the dynamical model. Through extensive Monte Carlo simulations, we systematically investigate the dynamics. Finally, we conclude.

2. Information Spreading Model

In this section, we introduce the information spreading dynamics on temporal networks with network memory.

2.1. Activity-Driven Network with Memory

Mathematically, the temporal network can be described as , where represents the temporal network at time step . For network , we use the temporal adjacency matrix to represent the topology of . If , there is an edge between nodes and at time . Otherwise, .

To build the activity-driven network with memory effect, we generalize the activity-driven model proposed in [4446]. We build memory activity-driven network as follows.(i)Assign value for network size and potential activity according to a given distribution . In this paper, we assume follows a power-law distribution. That is to say, , where is the potential activity distribution exponent, , , and . Mathematically, the larger the value of , the more homogeneous the degree distribution of the temporal network. Therefore, we can change to investigate the degree heterogeneity of the temporal network.(ii)Generating temporal network : for each node , there are two possible ways to build edges. If node becomes active with probability , where is a parameter, node forwardly connects to randomly selected nodes. If node is inactive, it can only receive the connections from other active nodes.(iii)At the end of the time step , we delete every edge with probability . Therefore, the memory effect is induced. The higher the value of , the more substantial the memory effect of the temporal network. For the case of , the temporal network is memoryless. When , the network is static.

According to the above steps, we know that the average degree of network is .

2.2. Information Spreading Model

We here adopt the information spreading model proposed in [47]. This model uses a generalized susceptible-infected-susceptible (SIS) model to describe the dynamics of information spreading. The susceptible nodes mean that they do not receive any information but may receive the information. The infected nodes represent that they have obtained the information and willing to share it with neighbors. In this model, we assume that there are two types of nodes, i.e., activists and bigots. The activists are willing to share the information with friends. Thus, we assume that they have a smaller adoption threshold . The bigots are less likely to accept the information. Therefore, we set those bigots with a higher adoption threshold . In this model, we randomly select a fraction of nodes as the bigots and the remaining nodes as activists.

The information spreading dynamics evolve as follows:(i)Randomly selecting fraction of nodes to receive the information and setting the rest as susceptible nodes.(ii)At every time step, every infected node tries to transmit the information to its very susceptible neighbor with probability . If node succeeds in receiving a piece of information, it would be represented by the mark +1 and the like. In this situation, we consider the state of node . If node is an activist, it becomes infected. If node is a bigot, it becomes infected only when its received information is larger than .(iii)Every infected node becomes a susceptible state with probability .

When the infected nodes no longer exist or the spreading dynamics has run 10000 times, the dynamics ends.

3. Theoretical Analyses

To obtain mathematical analysis results on the spreading dynamics, we use a generalized discrete Markovian approach, which is inspired by Refs. [39, 4850]. In theory, we assume that there are no dynamical correlations among the state of neighbors. That is to say, the infection probabilities of a susceptible node from informed neighbors are dependent. At any time step , node can only be susceptible or infected state. Define when node i is in the infected state and when in susceptible state. The probability of node in the susceptible and infected states is and , respectively. For the sake of simplicity, we denote and as and , respectively.

For the evolution of , we should consider two situations: the evolutions of bigots and activists, which are denoted as and , respectively. Therefore, we have

In the following, we study the evolutions of and . For an activist, such as node , the evolution of includes two situations. On the one hand, node is in the infected state at time step and does not recover to susceptible state with probability . On the other hand, node is susceptible at time step and gets infected by neighbors at this time step with probability , whererepresents the probability of node remaining susceptible at time step . The evolution of is

For the bigots, the evolution of is more complex. The bigots will only be infected if the information received is greater than the threshold . Therefore, we should first compute the number of information received by node i and denote as n pieces of information received by node at time step . Using the results presented in Ref. [47], we have

Until now, we have presented the expression of the evolution equations of the information spreading dynamics. In the steady state, i.e., , we have . We denote the fraction of nodes in the infected state in the steady state as , which can be computed as

We use as the order parameter of the system in the following numerical studies.

4. Results

In this section, we use the Monte Carlo method to study the information spreading dynamics on temporal network with memory effect. When generating the temporal network, we set the average degree of every temporal network with average degree . In the dynamical system, we set the order parameter aswhere and . We set . To numerically locate the steady state of the dynamical system, we compute the average values of . When the average value of variance is smaller than , the steady state of the system is reached.

In Figure 1, we first investigate the information spreading dynamics without bigots (i.e., ). For different strengths of the memory effect, we find two distinct regions. When the information transmission probability is small, a strong memory effect suppresses the information spreading dynamics (see and ). That is to say, the temporal network is not beneficial in information transmission. However, for larger values of , the temporal network is beneficial for information spreading. For static networks, the network topology does change with time. The information can be accessed by more susceptible nodes for small values of since the giant connected cluster is larger. However, more different nodes will be connected to the temporal giant connected cluster for temporal networks and promote the information spreading for larger . Once the network memory effect is not strong enough, compared with the information spreading on static and temporal networks, the information is suppressed regardless of the information transmission probability. Thus, we know that the middle memory effect can suppress the information spreading.

We further investigate the information spreading on a given strength of memory effect for different values of bigots in Figure 2. We note that both the fraction of bigots and memory affect the information spreading. For a given network topology, we find that increases with since bigots are not beneficial in transmitting the information. Besides, we note that the increase pattern of versus is different for distinct values of . When , increases continuously with . However, discontinuously increases with when is large (e.g., in Figures 2(b)2(d) and 2(f)).

The network memory effect may benefit the information spreading, which depends on the degree heterogeneity of the activity-driven network. When the potential exponent is and 3.5, the network memory effect suppresses the information spreading (Figures 2(a), 2(c), 2(d), and 2(f)). However, when the potential exponent is , the network memory effect promotes the information spreading, as shown in Figures 2(b) and 2(e).

Finally, we investigate the phase transition of the system with different initial seed sizes in Figure 3. Generally speaking, the phenomena are the same as those stated in Figure 2. We find that the system always has a hysteresis loop no matter what value and is assigned. In more detail, the final information spreading size depends on . The larger , the higher in the hysteresis loop region. We also note that the larger , the smaller . That is to say, the bigots hinder the information spreading dynamics.

5. Conclusions

In this paper, we proposed a mathematical model to investigate the effects of network memory on the information spreading dynamics on temporal networks. We first proposed a dynamical model in which the temporal network has a memory effect. Specifically, the current temporal network can remember the previous interconnections. Furthermore, an information diffusion model is developed on this type of temporal network. Then, we used a Markovian approach to describe the information spreading dynamics. Finally, we used the Monte Carlo simulation method to study the information spreading model numerically and found that network memory may promote and suppress the dynamics. The effects depend on the heterogeneous degree distribution and the fraction of bigots in the populations.

The results presented in this paper may shed some light into investigating the dynamics of information on temporal networks. On the one hand, the Markovian theory may be used to study other dynamics on temporal networks. On the other hand, the memory effect of temporal network should be included when studying other dynamics. Finally, some further studies about memory of temporal network should be investigated, for instance, developing more accurate theory and designing more realistic models to describe the spreading dynamics.

Data Availability

The datasets used in the present study are available from the first author upon reasonable request ([email protected]).

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 (no. U1733203), the Safety Foundation of CAAC (no. AQ20200019), and the Foundation of CAFUC (no. J2021-072).