Slow dynamics and rare-region effects in the contact process on weighted tree networks

Droplet finite-size scaling of the majority-vote model on scale-free networks

We discuss the majority vote model coupled with scale-free networks and investigate its critical behavior. Previous studies point to a nonuniversal behavior of the majority vote model, where the critical exponents depend on the connectivity. At the same time, the effective dimension D_{eff} is unity for a degree distribution exponent 5/2<γ<7/2. We introduce a finite-size theory of the majority vote model for uncorrelated networks and present generalized scaling relations with good agreement with Monte Carlo simulation results. Our finite-size approach has two sources of size dependence: an external field representing the influence of the mass media on consensus formation and the scale-free network cutoff. The critical exponents are nonuniversal, dependent on the degree distribution exponent, precisely when 5/2<γ<7/2. For γ≥7/2, the model is in the same universality class as the majority vote model on Erdős-Rényi random graphs. However, for γ=7/2, the critical behavior includes additional logarithmic corrections.

Impact of the distribution of recovery rates on disease spreading in complex networks

We study a general epidemic model with arbitrary recovery rate distributions. This simple deviation from the standard setup is sufficient to prove that heterogeneity in the dynamical parameters can be as important as the more studied structural heterogeneity. Our analytical solution is able to predict the shift in the critical properties induced by heterogeneous recovery rates. We find that the critical value of infectivity tends to be smaller than the one predicted by quenched mean-field approaches in the homogeneous case and that it can be linked to the variance of the recovery rates. Our findings also illustrate the role of dynamical-structural correlations, where we allow a power-law network to dynamically behave as a homogeneous structure by an appropriate tuning of its recovery rates. Overall, our results demonstrate that heterogeneity in the recovery rates, eventually in all dynamical parameters, is as important as the structural heterogeneity.

Robustness of Griffiths effects in homeostatic connectome models

I provide numerical evidence for the robustness of the Griffiths phase (GP) reported previously in dynamical threshold model simulations on a large human brain network with N=836733 connected nodes. The model, with equalized network sensitivity, is extended in two ways: introduction of refractory states or by randomized time-dependent thresholds. The nonuniversal power-law dynamics in an extended control parameter region survives these modifications for a short refractory state and weak disorder. In case of temporal disorder the GP shrinks and for stronger heterogeneity disappears, leaving behind a mean-field type of critical transition. Activity avalanche size distributions below the critical point decay faster than in the original model, but the addition of inhibitory interactions sets it back to the range of experimental values.

Griffiths phases in infinite-dimensional, non-hierarchical modular networks

Griffiths phases (GPs), generated by the heterogeneities on modular networks, have recently been suggested to provide a mechanism, rid of fine parameter tuning, to explain the critical behavior of complex systems. One conjectured requirement for systems with modular structures was that the network of modules must be hierarchically organized and possess finite dimension. We investigate the dynamical behavior of an activity spreading model, evolving on heterogeneous random networks with highly modular structure and organized non-hierarchically. We observe that loosely coupled modules act as effective rare-regions, slowing down the extinction of activation. As a consequence, we find extended control parameter regions with continuously changing dynamical exponents for single network realizations, preserved after finite size analyses, as in a real GP. The avalanche size distributions of spreading events exhibit robust power-law tails. Our findings relax the requirement of hierarchical organization of the modular structure, which can help to rationalize the criticality of modular systems in the framework of GPs.

Griffiths phase on hierarchical modular networks with small-world edges

The Griffiths phase has been proposed to induce a stretched critical regime that facilitates self-organizing of brain networks for optimal function. This phase stems from the intrinsic structural heterogeneity of brain networks, i.e., the hierarchical modular structure. In this work, the Griffiths phase is studied in modified hierarchical networks with small-world connections based on the 3-regular Hanoi network. Through extensive simulations, the hierarchical level-dependent inter-module wiring probabilities are identified to determine the emergence of the Griffiths phase. Numerical results and the complementary spectral analysis of the relevant networks can be helpful for a deeper understanding of the essential structural characteristics of finite-dimensional networks to support the Griffiths phase.

Critical dynamics on a large human Open Connectome network

Extended numerical simulations of threshold models have been performed on a human brain network with N=836733 connected nodes available from the Open Connectome Project. While in the case of simple threshold models a sharp discontinuous phase transition without any critical dynamics arises, variable threshold models exhibit extended power-law scaling regions. This is attributed to fact that Griffiths effects, stemming from the topological or interaction heterogeneity of the network, can become relevant if the input sensitivity of nodes is equalized. I have studied the effects of link directness, as well as the consequence of inhibitory connections. Nonuniversal power-law avalanche size and time distributions have been found with exponents agreeing with the values obtained in electrode experiments of the human brain. The dynamical critical region occurs in an extended control parameter space without the assumption of self-organized criticality.

Sampling methods for the quasistationary regime of epidemic processes on regular and complex networks

A major hurdle in the simulation of the steady state of epidemic processes is that the system will unavoidably visit an absorbing, disease-free state at sufficiently long times due to the finite size of the networks where epidemics evolves. In the present work, we compare different quasistationary (QS) simulation methods where the absorbing states are suitably handled and the thermodynamical limit of the original dynamics can be achieved. We analyze the standard QS (SQS) method, where the sampling is constrained to active configurations, the reflecting boundary condition (RBC), where the dynamics returns to the pre-absorbing configuration, and hub reactivation (HR), where the most connected vertex of the network is reactivated after a visit to an absorbing state. We apply the methods to the contact process (CP) and susceptible-infected-susceptible (SIS) models on regular and scale free networks. The investigated methods yield the same epidemic threshold for both models. For CP, that undergoes a standard collective phase transition, the methods are equivalent. For SIS, whose phase transition is ruled by the hub mutual reactivation, the SQS and HR methods are able to capture localized epidemic phases while RBC is not. We also apply the autocorrelation time as a tool to characterize the phase transition and observe that this analysis provides the same finite-size scaling exponents for the critical relaxation time for the investigated methods. Finally, we verify the equivalence between RBC method and a weak external field for epidemics on networks.

The topology of large Open Connectome networks for the human brain

The structural human connectome (i.e. the network of fiber connections in the brain) can be analyzed at ever finer spatial resolution thanks to advances in neuroimaging. Here we analyze several large data sets for the human brain network made available by the Open Connectome Project. We apply statistical model selection to characterize the degree distributions of graphs containing up to nodes and edges. A three-parameter generalized Weibull (also known as a stretched exponential) distribution is a good fit to most of the observed degree distributions. For almost all networks, simple power laws cannot fit the data, but in some cases there is statistical support for power laws with an exponential cutoff. We also calculate the topological (graph) dimension D and the small-world coefficient σ of these networks. While σ suggests a small-world topology, we found that D < 4 showing that long-distance connections provide only a small correction to the topology of the embedding three-dimensional space.

Griffiths effects of the susceptible-infected-susceptible epidemic model on random power-law networks

We provide numerical evidence for slow dynamics of the susceptible-infected-susceptible model evolving on finite-size random networks with power-law degree distributions. Extensive simulations were done by averaging the activity density over many realizations of networks. We investigated the effects of outliers in both highly fluctuating (natural cutoff) and nonfluctuating (hard cutoff) most connected vertices. Logarithmic and power-law decays in time were found for natural and hard cutoffs, respectively. This happens in extended regions of the control parameter space λ(1)<λ<λ(2), suggesting Griffiths effects, induced by the topological inhomogeneities. Optimal fluctuation theory considering sample-to-sample fluctuations of the pseudothresholds is presented to explain the observed slow dynamics. A quasistationary analysis shows that response functions remain bounded at λ(2). We argue these to be signals of a smeared transition. However, in the thermodynamic limit the Griffiths effects loose their relevancy and have a conventional critical point at λ(c)=0. Since many real networks are composed by heterogeneous and weakly connected modules, the slow dynamics found in our analysis of independent and finite networks can play an important role for the deeper understanding of such systems.

Griffiths phases and localization in hierarchical modular networks

We study variants of hierarchical modular network models suggested by Kaiser and Hilgetag [ Front. in Neuroinform., 4 (2010) 8] to model functional brain connectivity, using extensive simulations and quenched mean-field theory (QMF), focusing on structures with a connection probability that decays exponentially with the level index. Such networks can be embedded in two-dimensional Euclidean space. We explore the dynamic behavior of the contact process (CP) and threshold models on networks of this kind, including hierarchical trees. While in the small-world networks originally proposed to model brain connectivity, the topological heterogeneities are not strong enough to induce deviations from mean-field behavior, we show that a Griffiths phase can emerge under reduced connection probabilities, approaching the percolation threshold. In this case the topological dimension of the networks is finite, and extended regions of bursty, power-law dynamics are observed. Localization in the steady state is also shown via QMF. We investigate the effects of link asymmetry and coupling disorder, and show that localization can occur even in small-world networks with high connectivity in case of link disorder.

Localization transition, Lifschitz tails, and rare-region effects in network models

Effects of heterogeneity in the suspected-infected-susceptible model on networks are investigated using quenched mean-field theory. The emergence of localization is described by the distributions of the inverse participation ratio and compared with the rare-region effects appearing in simulations and in the Lifschitz tails. The latter, in the linear approximation, is related to the spectral density of the Laplacian matrix and to the time dependent order parameter. I show that these approximations indicate correctly Griffiths phases both on regular one-dimensional lattices and on small-world networks exhibiting purely topological disorder. I discuss the localization transition that occurs on scale-free networks at γ=3 degree exponent.

Slow, bursty dynamics as a consequence of quenched network topologies

Bursty dynamics of agents is shown to appear at criticality or in extended Griffiths phases, even in case of Poisson processes. I provide numerical evidence for a power-law type of intercommunication time distributions by simulating the contact process and the susceptible-infected-susceptible model. This observation suggests that in the case of nonstationary bursty systems, the observed non-Poissonian behavior can emerge as a consequence of an underlying hidden Poissonian network process, which is either critical or exhibits strong rare-region effects. On the contrary, in time-varying networks, rare-region effects do not cause deviation from the mean-field behavior, and heterogeneity-induced burstyness is absent.

Spectral analysis and slow spreading dynamics on complex networks

The susceptible-infected-susceptible (SIS) model is one of the simplest memoryless systems for describing information or epidemic spreading phenomena with competing creation and spontaneous annihilation reactions. The effect of quenched disorder on the dynamical behavior has recently been compared to quenched mean-field (QMF) approximations in scale-free networks. QMF can take into account topological heterogeneity and clustering effects of the activity in the steady state by spectral decomposition analysis of the adjacency matrix. Therefore, it can provide predictions on possible rare-region effects, thus on the occurrence of slow dynamics. I compare QMF results of SIS with simulations on various large dimensional graphs. In particular, I show that for Erdős-Rényi graphs this method predicts correctly the occurrence of rare-region effects. It also provides a good estimate for the epidemic threshold in case of percolating graphs. Griffiths Phases emerge if the graph is fragmented or if we apply a strong, exponentially suppressing weighting scheme on the edges. The latter model describes the connection time distributions in the face-to-face experiments. In case of a generalized Barabási-Albert type of network with aging connections, strong rare-region effects and numerical evidence for Griffiths Phase dynamics are shown. The dynamical simulation results agree well with the predictions of the spectral analysis applied for the weighted adjacency matrices.

Slow epidemic extinction in populations with heterogeneous infection rates

We explore how heterogeneity in the intensity of interactions between people affects epidemic spreading. For that, we study the susceptible-infected-susceptible model on a complex network, where a link connecting individuals i and j is endowed with an infection rate β(ij)=λw(ij) proportional to the intensity of their contact w(ij), with a distribution P(w(ij)) taken from face-to-face experiments analyzed in Cattuto et al. [PLoS ONE 5, e11596 (2010)]. We find an extremely slow decay of the fraction of infected individuals, for a wide range of the control parameter λ. Using a distribution of width a we identify two large regions in the a-λ space with anomalous behaviors, which are reminiscent of rare region effects (Griffiths phases) found in models with quenched disorder. We show that the slow approach to extinction is caused by isolated small groups of highly interacting individuals, which keep epidemics alive for very long times. A mean-field approximation and a percolation approach capture with very good accuracy the absorbing-active transition line for weak (small a) and strong (large a) disorder, respectively.

Rare regions of the susceptible-infected-susceptible model on Barabási-Albert networks

I extend a previous work to susceptible-infected-susceptible (SIS) models on weighted Barabási-Albert scale-free networks. Numerical evidence is provided that phases with slow, power-law dynamics emerge as the consequence of quenched disorder and tree topologies studied previously with the contact process. I compare simulation results with spectral analysis of the networks and show that the quenched mean-field (QMF) approximation provides a reliable, relatively fast method to explore activity clustering. This suggests that QMF can be used for describing rare-region effects due to network inhomogeneities. Finite-size study of the QMF shows the expected disappearance of the epidemic threshold λ(c) in the thermodynamic limit and an inverse participation ratio ~0.25, meaning localization in case of disassortative weight scheme. Contrarily, for the multiplicative weights and the unweighted trees, this value vanishes in the thermodynamic limit, suggesting only weak rare-region effects in agreement with the dynamical simulations. Strong corrections to the mean-field behavior in case of disassortative weights explains the concave shape of the order parameter ρ(λ) at the transition point. Application of this method to other models may reveal interesting rare-region effects, Griffiths phases as the consequence of quenched topological heterogeneities.