key: cord-0043386-8nwhxa2r
authors: Tizzani, Michele; Lenti, Simone; Ubaldi, Enrico; Vezzani, Alessandro; Castellano, Claudio; Burioni, Raffaella
title: Epidemic spreading and aging in temporal networks with memory
date: 2018-12-18
journal: nan
DOI: 10.1103/physreve.98.062315
sha: 4dd0c46c0f436b67eae701caafc66dc121dd58fc
doc_id: 43386
cord_uid: 8nwhxa2r

Time-varying network topologies can deeply influence dynamical processes mediated by them. Memory effects in the pattern of interactions among individuals are also known to affect how diffusive and spreading phenomena take place. In this paper we analyze the combined effect of these two ingredients on epidemic dynamics on networks. We study the susceptible-infected-susceptible (SIS) and the susceptible-infected-recovered (SIR) models on the recently introduced activity-driven networks with memory. By means of an activity-based mean-field approach, we derive, in the long-time limit, analytical predictions for the epidemic threshold as a function of the parameters describing the distribution of activities and the strength of the memory effects. Our results show that memory reduces the threshold, which is the same for SIS and SIR dynamics, therefore favoring epidemic spreading. The theoretical approach perfectly agrees with numerical simulations in the long-time asymptotic regime. Strong aging effects are present in the preasymptotic regime and the epidemic threshold is deeply affected by the starting time of the epidemics. We discuss in detail the origin of the model-dependent preasymptotic corrections, whose understanding could potentially allow for epidemic control on correlated temporal networks.

Time-varying network topologies can deeply influence dynamical processes mediated by them. Memory effects in the pattern of interactions among individuals are also known to affect how diffusive and spreading phenomena take place. In this paper we analyze the combined effect of these two ingredients on epidemic dynamics on networks. We study the susceptible-infected-susceptible (SIS) and the susceptible-infected-recovered (SIR) models on the recently introduced activity-driven networks with memory. By means of an activity-based mean-field approach, we derive, in the long-time limit, analytical predictions for the epidemic threshold as a function of the parameters describing the distribution of activities and the strength of the memory effects. Our results show that memory reduces the threshold, which is the same for SIS and SIR dynamics, therefore favoring epidemic spreading. The theoretical approach perfectly agrees with numerical simulations in the long-time asymptotic regime. Strong aging effects are present in the preasymptotic regime and the epidemic threshold is deeply affected by the starting time of the epidemics. We discuss in detail the origin of the model-dependent preasymptotic corrections, whose understanding could potentially allow for epidemic control on correlated temporal networks. DOI: 10.1103/PhysRevE.98.062315

In many social and information systems, the timescales for the evolution of the interaction network are often comparable to the timescales of the dynamical processes taking place on top of them [1] . The diffusion of online information or the spreading of transmitted diseases in a population are typical examples of such processes, for which a focus on a static representation of the network is not able to capture the very influence of the rapidly varying topology [2] [3] [4] [5] [6] [7] . In addition, recent advances in technology have allowed one to measure and monitor the evolution of interactions with an unprecedented time resolution [8] , calling for new theories to understand the effect of time-varying topologies on dynamical processes.

Interactions and the creation of links are generated by the agents' activity, a quantity that can be easily measured from available large-scale and time-resolved data sets [9] . An interesting line of modeling has developed, aiming at including explicitly the effect of activity distributions on network dynamics: activity-driven networks [10] . In activity-driven models, each agent is endowed with a degree of freedom that encodes the propensity of the individual to engage in a social event, establishing a link with another agent in the system. Notably, measured activities are typically highly heterogeneous and this has strong effects on network evolution.

When links are randomly established among agents, activity-driven models have been studied in detail [10] [11] [12] [13] , uncovering the effects of heterogeneous activity distributions on network topology and on dynamical processes, such as random walks and epidemic processes. However, in general, agents do not connect randomly to their peers [14] [15] [16] . During their activity, individuals remember their friends and their social circles and they are more inclined to interact with already known pals, establishing strong and weak ties with their peers [17, 18] . Recent works have tackled this problem by applying a data-driven approach. A tie allocation mechanism in real systems has been measured, introducing a memory process on top of activity-driven models [19, 20] . As reasonably expected, social interactions are not randomly established but they are rather concentrated towards already contacted nodes, with a reinforcement process encoded in a single measurable memory parameter. The memory process tunes the network evolution, which can be predicted at large times [21] [22] [23] , and it is also expected to influence dynamical processes. Non-Markovian dynamics can indeed change the spreading rate in a diffusion process, slowing it down in some cases and speeding it up in others [24] [25] [26] [27] [28] [29] [30] . Disparate effects have been shown to occur also in epidemic spreading on activity-driven networks, where memory can lower or increase the epidemic threshold in the susceptible-infected-susceptible (SIS) or susceptible-infected-recovered (SIR) model, respectively [31] . This happens when the epidemic process and the network evolution start at the same time. However, the network evolution in the presence of non-Markovian effects could introduce aging in the process [32] , as usually observed in processes with memory also in other fields [33] , and this could further influence the spreading dynamics.

In this paper we analyze in detail SIS and SIR epidemic processes on activity-driven time-varying networks in the presence of memory. Introducing an activity-based mean-field (ABMF) approximation, we derive analytically, as a function of the activity distribution and of the parameter tuning the memory, a prediction for the threshold, holding for both the SIS and the SIR models. The result shows that memory overall reinforces the effects of activity fluctuations, leading to a lower value for the epidemic threshold. We prove that the ABMF approximation is equivalent to an epidemic model defined on an effective static network, which we also investigate. Numerical simulations evidence that the ABMF approach provides exact results when the epidemics start after the network has evolved for a long time. In this regime indeed the mean field holds since agents have been connected to a large number of pairs, while the creation of new links becomes negligible.

We show, however, that strong aging effects are present in the preasymptotic regime and the epidemic threshold is deeply affected by the starting time of the epidemics. In particular, at weak memory the creation of new links cannot be discarded and this increases the epidemic threshold according to the memoryless predictions. On the other hand, for strong memory, at short times the dynamics displays correlations among infection probabilities of nodes which have already been in contact. The correlations give rise to backtracking effects that cannot be ignored. In this case, typically the threshold of the SIS and SIR models are respectively smaller and larger with respect to the mean-field prediction. We explain in detail the origin of such deviations, opening new perspectives for epidemic control of disease and information spreading on temporal networks with high correlations.

The paper is organized as follows. In Sec. II we first summarize the activity-driven model for network topology in the presence of memory and then define epidemic dynamics on top of it. In Sec. III, after a brief recapitulation of the analytical approach to epidemic dynamics on memoryless activity-driven networks, we describe in detail how the approach is modified to deal with networks with memory, deriving predictions for the epidemic threshold. In Sec. IV we compare analytical predictions with numerical results, obtained by considering both an effective static network and the full time evolution of the topology. Section V presents some concluding remarks and perspectives for future work.

In activity-driven models [10] , each node v i (i = 1, . . . , N) of the graph G t has an activity a i assigned randomly according to a given distribution F (a). The dynamics occurs over discrete temporal steps of length t. At each step, with probability a i t, the vertex v i becomes active and gets linked to m other vertices. Connections last for a temporal interval t. At the next time step t + t all existing edges are deleted and the procedure is iterated. Note that the activity a has the dimensions of a probability per unit time. Real data observations indicate that human interactions are very often characterized by skewed and long-tailed activity distributions [10] , so F (a) is typically assumed to be a power law F (a) = Ba −(ν+1) with ε a i A. Since in our simulations we will keep the time interval t = 1, the upper cutoff is naturally set to A = 1.

In order to take into account the tendency of individuals to persist in their social connections, a reinforcement mechanism was recently introduced in activity-driven models [19, 21] . The nodes are endowed with a memory of their previous contacts and they have a propensity to establish contacts preferably with individuals belonging to their social circle. For an active node v i , which has already contacted k i (t ) different nodes at time t, this process is described by assuming that the node connects with a new node with probability

while it establishes a connection with a previously contacted node with complementary probability 1 − p[k i (t )]. The probability depends on the degree of the integrated network at time t, k i (t ), i.e., the number of nodes that v i has contacted up to time t. We will call A ij (t ) the adjacency matrix of this integrated network. The parameter β i > 0 tunes the memory process. For β i ≈ 0 the probability p[k i (t )] ≈ 1 weakly depends on the growing degree k i (t ), while at large β i , p[k i (t )] rapidly decays with k i (t ). The constant c i sets an intrinsic value for the number of connections that node v i is able to engage in before memory effects become relevant. Empirical measures on several data sets [21] are compatible with constant values of β i = β and c i = c, so the function p(·) turns out to be independent of i. In particular, the exponent β ranges from β ≈ 0.15 in the citation networks to β ≈ 0.5 in Twitter mentions and to β ≈ 1.2 in the mobile phone calls. In this paper we study the dependence of the epidemic threshold on the value of the exponent β, while we set c = 1 and m = 1, with no loss of generality. As shown in Ref. [21] , the asymptotic form of the degree distribution for the integrated network can be derived analytically. In particular, in the regime 1 k N the degree of nodes of activity a is narrowly distributed around the average valuek

i.e., the degree of each node increases sublinearly in time, with a prefactor depending on its activity. The prefactor C(a) is determined by the condition

In the memoryless case β = 0, where an active node connects always with a randomly chosen vertex, Eq. (3) gives C(a) = a + a , recovering the result of [34] . Hereafter we will define, in general, by g = da F (a)g(a) the average of a function of the activity g(a) over the network.

We now turn to the spreading of infectious diseases on activity-driven temporal networks with memory. We start by considering the standard SIS model, the simplest description of a disease not conferring immunity. In the SIS model nodes can be in two states, either susceptible (S) or infectious (I). An infected node can turn spontaneously susceptible with rate μ, while an infected node transmits the infection over an edge to a susceptible neighbor with rate λ. The two elementary events are therefore

In the SIR model the disease confers immunity, so the nodes can be in three states, susceptible (S), infectious (I), and recovered (R), which are immune to a new infection. The dynamics is described by the following reaction scheme:

The epidemic process on activity-driven networks is implemented by iterating discrete time steps of duration t.

(i) At the beginning of each time step there are N disconnected vertices.

(ii) With probability a i t a vertex v i becomes active and connects to a previously linked node with probability 1 − p(k i ) or with a new node v j with probability p(k i ); in this second case k i (t ), k j (t ), and A ij (t ) are increased by one unit.

(iii) If one of the nodes connected by the link is infected and the other one is susceptible, the susceptible becomes infected with probability λ.

(iv) A vertex v j , if infected, becomes susceptible (SIS) or recovers (SIR) with probability μ t.

In activity-driven models, λ is a pure number, i.e., the probability that in a single contact the infection is actually transmitted, while μ is still the rate of recovery for a single individual. Ignoring the inhomogeneity in the activities, one can estimate the total rate for the infection process per node as λ k , where k = 2 a is the average degree per unit time; this is the quantity to be compared with the recovery rate per node μ.

The epidemic spreading for the memoryless case β = 0 has been studied in [10] by adopting an ABMF approach. The epidemic state of a node, when averaged over all possible dynamical evolutions, only depends on the value of its activity a i . In particular, one can define the probability ρ(a i , t ) that a node with activity a i is infected at time t. The corresponding evolution equation is

The first term on the right-hand side is due to recovery events; the second term takes into account the event that a susceptible node of class a i becomes active and contracts the disease by connecting to an infected individual, while the third term is the analogous term for the case of a susceptible node that, independently of one's own activity, is contacted by an infected active individual. The description in terms of quantities that only depend on the activity is conceptually analogous to the heterogeneous mean-field approach for dynamical processes on static networks [35] . In that case, one assumes that the only property determining the epidemic state of a node is the degree k and then derives equations for the probabilities ρ k . An important difference must however be stressed. Assuming the epidemic state to depend only on the degree is an approximation for static networks, because it neglects the quenched nature of the network structure that makes properties of nodes, with the same degree but embedded in different local environments, different. In practice, this assumption is equivalent to replacing the actual adjacency matrix of the network (A ij equal to 0 or 1 depending on the presence of the connection between v i and v j ) with an annealed adjacency matrix P ij = k i k j /( k N ) [36] , expressing the probability that vertices v i and v j with degree k i and k j are connected. The annealed approach is an approximation for static networks, while it is exact for networks where connections are continuously reshuffled at each time step of the dynamics, since the reshuffling process destroys local correlations. Since in memoryless activitydriven networks connections are extracted anew at each time step, the ABMF approach provides exact results in this case.

Equation (6) can be analyzed by means of a linear stability analysis, yielding, for large N , the threshold [10]

The same result can be derived for the SIR case.

In the presence of memory, interactions occur preferably with a subset of the other nodes (the social circle) creating correlations. Therefore, we implement a different, individualbased, mean-field approach, explicitly keeping track of the evolution of social contacts (i.e., of the memory). Let us first consider the SIS model. The observable of interest is the probability ρ i (t ) that node v i is infected at time t. Its evolution can be written as

Here j i indicates the sum over the nodes j not yet connected to i and N − k j (t ) − 1 is their number. The quantity A ij (t ) is the adjacency matrix of the time-integrated network at time t, i.e., it is equal to 1 if v i and v j have been in contact at least once in the past and 0 otherwise. In Eq. (8), the only approximation made is that the dynamical state of every node is considered to be independent of the state of the partner in the interaction; in other words, we neglect the existence of dynamical correlations among nodes which are created by the partially quenched nature of the interaction pattern due to memory. It is exactly the same approximation that is involved by the individual-based mean-field approach for static networks [37] . The first term on the right-hand side of Eq. (8) is the recovery rate of ρ i (t ). The second term, describing the infection process, is the product of λ times the probability for v i to be susceptible and, in curly brackets, the fraction of infected nodes contacted by v i per unit time. In the curly brackets, the first and the second term describe the case where v i is active and connects to the infected node v j taking into account that the link can be an old or a new one, respectively. In the same way, the third and the fourth term represent the probabilities that v i is contacted by an infected and active node v j , and the link is old or new, respectively.

Since both A ij (t ) and k i (t ) depend on the evolution time t, the behavior of the epidemics can strongly depend on the starting time of the outbreak, giving rise to aging effects that will be investigated in numerical simulations. When the epidemic starts at very large times, an analytic approach can be considered. In this regime, with 1 k i (t ) N, we expect that the creation of new contacts can be ignored and that the dynamical correlations are asymptotically negligible, since the connectivity of the integrated network becomes large. If the epidemic starts at very large times, therefore, we can apply an heterogeneous mean-field approximation for A ij (t ), allowing for an analytical solution of the problem which we expect to be asymptotically exact.

Let us consider the regime of large times, where 1 k i (t ) N for all nodes. This means that each node has already had a large number of contacts, but that number is not so large that the integrated network cannot be considered as a complete graph, i.e., it is still sparse. In the limit of large N there is clearly a large temporal interval such that this condition is fulfilled. The condition 1 k i (t ) N allows us to replace in Eq.

Considering only leading terms, Eq. (8) becomes

To proceed further we perform the equivalent of the heterogeneous mean-field approximation for static networks, i.e., we replace the time-integrated adjacency matrix A ij (t ) with its annealed form P ij (t ), i.e., the probability that v i and v j have been in contact in the past. The evolution of P ij (t ) is described by the master equation In the temporal interval of interest we can use again the relations holding for large times p(k i ) ≈ k −β i and N − k j − 1 ≈ N . Moreover, for large times, the degree of a node of activity a i can be estimated by its average valuek(a i , t ), given by Eq. (2). So we obtain

where we have defined

Equation (11) can be readily solved, yielding

In the regime t 1/(1+β ) N , P ij (t ) becomes

Notice that P ij (t ) is a topological feature of the activity-driven network, independent of the epidemic process.

We now introduce into Eq. (9) the annealed expression for the integrated adjacency matrix A ij (t ) ≈ P ij (t ) = P (a i , a j , t ) and for the connectivity k i (t ) =k(a i , t ). In this way the equations depend on the nodes v i and v j only through their activities a i and a j . The equation for the probability ρ(a, t ) that a generic node of activity a is infected at time t is therefore

where we have replaced the sums over nodes with integrals over the activities 1/N j → da F (a ) and used Eq. (3), which can be rewritten as

Equation (15) is effectively an ABMF approach, since all the information on the behavior of the node v i depends on its activity a i . Note that, although Eqs. (8) and (9) underlying Eq. (10) has transformed the approach into an ABMF one, conceptually analogous to the heterogeneous mean-field approximation on static networks, where all the information on node v i is encoded in its degree k i .

It is important to remark that in Eq. (9) the time dependences of P (a i , a j , t ) ∝ t 1/(1+β ) and of the average degreek(a i , t ) ∝ t 1/(1+β ) cancel out, so the right-hand side of Eq. (15) does not depend explicitly on time. This suggests that in this temporal regime the epidemic can be seen as an activity-driven process taking place on an effective static graph, where the probability for nodes v i and v j to be linked is given by Eq. (14) and the quantity t 1/(1+β ) /N is a fixed quantity τ whose value only determines the average degree of the network. Performing simulations over an ensemble of these effective static networks and averaging the results, one should then reproduce the predictions of the ABMF approach, Eq. (15).

From Eq. (15), by performing a linear stability analysis around the absorbing state ρ(a, t ) = 0 (see the Appendix), it is possible to compute analytically the epidemic threshold (λ/μ) c for any value of the reinforcement parameter β and of the exponent of the analytical distribution ν. Since for large times the node degrees diverge and correlations can be neglected, we expect the linear stability analysis to provide the correct estimate of the epidemic threshold when the epidemics start at very long times, i.e., when the degrees k i (t ) have already become very large.

The results of the linear stability analysis are presented in Fig. 1 , showing that the thresholds are smaller than in the memoryless case. This lower value is a consequence of the fact that memory reinforces the activity fluctuations, and in these models fluctuations clearly reduce the epidemic threshold, as shown by Eq. (7) . The effect can be simply understood since nodes with large activity have also a large degree, therefore they are easily involved in epidemic contacts not only because they are frequently activated but also because they are frequently contacted by other nodes. In this way memory reinforces the effect of activity fluctuations. In this framework, Fig. 1 also shows that at large ν, i.e., for increasingly smaller fluctuations, the difference with the memoryless model vanishes. In particular, for F (a) = δ(a − a 0 ), i.e., when the activity does not fluctuate, one obtains from Eq.

, which is the same equation as in the memoryless case. This also explains the quite surprising observation that the threshold is a growing function of β, converging to the memoryless case as β → ∞. Indeed, the tail of the degree distribution decays at large k as k −[(1+β )ν+1] [21] , therefore at large β we get a faster decay and smaller degree fluctuations. For the same reason, in the limit β → 0 the difference with the memoryless case is maximal, since degree inhomogeneities are stronger in this case.

We remark that in Eq. (15) , as in the memoryless case, dynamical correlations are ignored. However, we expect that at finite times, due to the finite connectivity of the integrated graph, the effect of correlations becomes important. The memory process leads to the formation of small clusters of mutually connected high activity vertices, which become reservoirs of the disease in the SIS model. The high frequency of mutual contacts allows for reinfection, favoring the overall survival of the epidemic spreading in the system. In this way, social circles with high activity play a role analogous to that played by the maximum K-core or the hub and its immediate neighbors for SIS epidemics in static networks [38, 39] . To clarify the effect of dynamical correlations at finite time, in the next section we compare the analytical predictions with the results of numerical simulations. As a final remark, we note that, in the asymptotic ABMF approach, the linear stability analysis also holds for the SIR model, implying that the epidemic threshold is the same of the SIS model. However, in the SIR model reinfection is not allowed, so the initial presence of small clusters of mutually connected high activity vertices effectively inhibits the spread of the disease. For this reason, we expect that finite connectivity (i.e., finite time) increases the epidemic threshold with respect to the mean-field result, as we will check in numerical simulations.

As discussed above, Eq. (15) can be interpreted as a heterogeneous mean-field equation for a SIS epidemic on an effective static network where the probability that v i and v j are connected is

Here τ 1 is a constant, g(a) = a/[C(a)] β , and C(a) is a function that can be evaluated numerically for β > 0, while for β = 0 it takes the simple form C(a) = a + a . The constant τ can be tuned in order to set the average degree of the network, because 

We now study the SIS epidemic evolution on the effective static network. Given the activity of each node, extracted according to the distribution F (a), for each of the possible pairs of nodes we place an edge with probability given by Eq. (17) . On top of this quenched topology we run a memoryless activity-driven SIS dynamics, starting with 10% of the nodes in the infected state, until the stationary state state is reached and we record the fraction of infected nodes. We repeat the procedure many times for each value of λ, while μ is fixed to 0.015. We determine the threshold as the position of the maximum of the susceptibility [40] χ = N (ρ 2 − ρ 2 )/ρ, where the overbar denotes the average over dynamical realizations at fixed topology. We repeat all this for several networks obtained using different sequences of activities and different samplings of P ij and we average the epidemic thresholds found for each of them.

We first check that, as long as 1 k N , the results are independent of the exact value of k , as predicted by the theory. Figure 2 shows, for β = 1, that the effective threshold initially grows with k but then reaches a plateau, in accordance with the expectations.

In Fig. 3 we report the dependence of the effective epidemic threshold as a function of β for three values of the average degree k , compared with the predictions of the ABMF theory with and without memory. We observe that, as k grows, numerical results tend to coincide with theoretical predictions. Nice agreement with the ABMF results is obtained for k ≈ 100, which is a large but realistic number of connections in a social system. We remark that k ≈ 100 is large enough for observing mean-field behavior but also it is much smaller than the total number of sites N = 5 × 10 4 , so the system is not fully connected and degree fluctuations are important. On the other hand, for small values of k the value of the threshold is smaller than the one predicted theoretically. Indeed, on effective static networks with small connectivity we expect the presence of clusters of mutually interconnected nodes to be relevant, as they are able to reinfect each other several times. It is well known that for the SIS model these backtracking effects tend to lower the epidemic threshold since social circles with high activity favor the overall survival of the epidemic.

Let us now consider simulations of the epidemic spreading on the full time-evolving network. We consider a graph of size N = 5 × 10 4 with activity distributed according to F (a) = Ba −(ν+1) (ν = 2.4) and a cutoff ε = 10 −2 . To extract the activities of individual nodes we perform an importance sampling so that, even in the finite-size system, the moments a and a 2 coincide with their expected values.

We first start the temporal evolution of the network and at a later time t 0 we let the epidemic begin. We evaluate at t 0 the average connectivity of the nodes k 0 which measures the evolution of the network at the starting time. In both the SIS and SIR models, we use two different initial conditions. The first is to randomly select the node (RA) to infect at time t 0 , while the second is to infect the most active node (MA) at time t 0 . We keep the recovery rate fixed at μ = 1.5 × 10 −2 and vary the probability of infection λ to study the dependence of its critical value on the memory parameter β.

In the SIS model, we determine the epidemic threshold using the lifespan method [41] . We plot (see Fig. 4 ), as a function of the parameter λ, the average lifespan of simulations ending before the coverage (i.e., the the fraction of distinct sites ever infected) reaches a preset value that we take equal to 1/2. The threshold is estimated as the value of λ for which the lifespan has a peak.

The epidemic thresholds of numerical simulations are compared with theoretical predictions in Figs. 5 (RA case) and 6 (MA case). Numerical results converge toward the analytical prediction as k 0 becomes larger, while there are strong deviations for small k 0 .

For small k 0 two competing effects are at work. First, infections are mediated by an effective static network with small connectivity, therefore we expect backtracking effects to enhance epidemic spreading and to lower the threshold. Second, for small k 0 , moves connecting new partners are also possible. In these moves nodes are chosen randomly in the whole system according to a memoryless dynamics, which displays a higher epidemic threshold. So there exists a competition between backtracking correlations and memoryless moves which reduce and increase the threshold, respectively. Clearly, for large k 0 both effects become negligible and the ABMF result is recovered. Figures 5 and 6 show that at large β backtracking effects (leading to small thresholds) are strong when the evolution starts from the most active site, while they are negligible with random initial conditions. The most active node indeed has the largest degree and one forms a cluster of highly activated nodes where the high frequency of mutual contacts allows for reinfections and positive correlations. Conversely, the average site has a small connectivity and can activate new links with high probability giving rise essentially to a memoryless epidemic dynamics.

The case β = 0 coincides with the memoryless case (ML) and the dynamics never occurs on the effective static network. Figures 5 and 6 show that at very small β, even for the largest considered value k 0 = 120, the creation of new links dominates the dynamics increasing the epidemic threshold towards the memoryless case. However, we expect that for large enough k 0 , at any β > 0, the dynamics is dominated by memory and infections spread on the effective static network recovering the ABMF result. This originates a singular behavior where the limits β → 0 and k 0 → ∞ do not commute.

The results of simulations of the SIR process are displayed in Figs. 7 and 8 for the RA and MA cases, respectively. The threshold is estimated from the peak of the variability = √ N 2 R − N R 2 / N R , i.e., the standard deviation of the number of recovered nodes N R at the end of the simulation [42] . As for SIS, at large k 0 dynamical correlations can be neglected and simulations recover the ABMF result. Simulations clearly show that now correlations at small k 0 inhibit the epidemic spreading and the critical threshold becomes larger. As in the SIS model, at small β the memory is negligible and the dynamics is driven by the creation of new links, so the threshold values are close to the memoryless case, even if for any β > 0 we expect for large enough k 0 a convergence to the ABMF prediction. On the other hand, for larger β, the epidemics evolve on the integrated network, dynamical correlations become important, and the thresholds grow even larger than in the memoryless case.

The analytical and numerical results presented in this paper provide a complete understanding of the interplay between the temporal evolution of the activity-driven network with memory and the epidemic process occurring on top of it. The time when the epidemic process begins has crucial consequences. In the long-time limit the reinforcement mechanism of the topological evolution completely inhibits the formation of new connections. If the activity-driven epidemic dynamics starts at this stage, it takes place on a topology which is in practice static. All nodes have a very large number of connections and this makes mean-field theory asymptotically exact. The epidemic threshold, which is the same for SIS and SIR dynamics, is reduced with respect to the memoryless case, because memory enhances the effect of activity fluctuations.

If instead the epidemic process starts before memory has completely taken over, interesting model-dependent preasymptotic effects are observed. The fundamental observation is that at this stage nodes with large activity tend to interact with their social circles, while less active nodes still tend to explore the system creating new connections. The first type of interaction tends to enhance SIS spreading, while the second tends to suppress it. This leads to positive or negative corrections to the asymptotic value of the threshold, depending on the initial conditions and on the reinforcement parameter β. In the SIR case instead, since reinfection is not possible, the interaction within social circles is strongly detrimental for the epidemic propagation, so the asymptotic value of the threshold is always reached from above. Our results allow us to fully understand the contrasting effects of strong ties on SIS and SIR dynamics observed in Ref. [31] and it offers different possibilities in control of epidemic spreading on temporal networks with high correlations.

Several possible extensions of the model considered here are possible to make it more realistic, in terms of both the topological evolution and the spreading process. Among them, probably the most interesting would be the inclusion of burstiness in the agents' activity. The combined effect on activity-driven models of tie reinforcement and nonexponentially distributed interevent times has been studied in Refs. [22, 43] . The inclusion of a spreading dynamics in this framework is a promising and challenging avenue for future research.

The dynamical process is described by the ABMF equation (15) , which we rewrite as λ Ag − μ λ A λ g λ λ Ag 2 λ Ag − μ λ g 2 λ g λ A 2 g λ A 2 λ Ag − μ λ A λ A 2 g 2 λ A 2 g λ Ag 2 λ Ag − μ ⎞ ⎟ ⎟ ⎠ .

(A4)

The state ρ(a) = 0 is stable provided all eigenvalues of J are negative, hence the epidemic threshold is given by the value (λ/μ) c such that largest eigenvalue of the Jacobian matrix is zero. Numerical evaluation of the matrix J and of its eigenvalues can be obtained by first solving numerically Eq. (3) to get C(a) and g(a) and then calculating the averages defining J .

Proc. Natl. Acad. Sci. USA

Dynamical Processes on Complex Networks

Proc. Natl. Acad. Sci. USA

Ageing and Dynamical Scaling Far from Equilibrium