Search

We solve the weak percolation problem for multiplex networks with overlapping edges. In weak percolation, a vertex belongs to a connected component if at least one of its neighbors in each of the layers is in this component. This is a weaker condition than for a mutually connected component in interdependent networks, in which any two vertices must be connected by a path within each of the layers. The effect of the overlaps on weak percolation turns out to be opposite to that on the giant mutually connected component. While for the giant mutually connected component, overlaps do not change the critical phenomena, our theory shows that in two layers any (nonzero) concentration of overlaps drives the weak percolation transition to the ordinary percolation universality class. In three layers, the phase diagram of the problem contains two lines -- of a continuous phase transition and of a discontinuous one -- connected in various ways depending on how the layers overlap. In the case of only doubled overlapped edges, two of the end points of these lines coincide, resulting in a tricritical point like that seen in heterogeneous $k$-core percolation., Comment: 21 pages, 8 figures

We consider the effect of a nonvanishing fraction of initially infected nodes (seeds) on the SIR epidemic model on random networks. This is relevant when, for example, the number of arriving infected individuals is large, but also to the modeling of a large number of infected individuals, but also to more general situations such as the spread of ideas in the presence of publicity campaigns. This model is frequently studied by mapping to a bond percolation problem, in which edges in the network are occupied with the probability, $p$, of eventual infection along an edge connecting an infected individual to a susceptible neighbor. This approach allows one to calculate the total final size of the infection and epidemic threshold in the limit of a vanishingly small seed fraction. We show, however, that when the initial infection occupies a nonvanishing fraction $f$ of the network, this method yields ambiguous results, as the correspondence between edge occupation and contagion transmission no longer holds. We propose instead to measure the giant component of recovered individuals within the original contact network. This has an unambiguous interpretation and correctly captures the dependence of the epidemic size on $f$. We give exact equations for the size of the epidemic and the epidemic threshold in the infinite size limit. We observe a second order phase transition as in the original formulation, however with an epidemic threshold which decreases with increasing $f$. When the seed fraction $f$ tends to zero we recover the standard results., Comment: 10 pages, 6 figures

Weak multiplex percolation generalizes percolation to multi-layer networks, represented as networks with a common set of nodes linked by multiple types (colors) of edges. We report a novel discontinuous phase transition in this problem. This anomalous transition occurs in networks of three or more layers without unconnected nodes, $P(0)=0$. Above a critical value of a control parameter, the removal of a tiny fraction $\Delta$ of nodes or edges triggers a failure cascade which ends either with the total collapse of the network, or a return to stability with the system essentially intact. The discontinuity is not accompanied by any singularity of the giant component, in contrast to the discontinuous hybrid transition which usually appears in such problems. The control parameter is the fraction of nodes in each layer with a single connection, $\Pi=P(1)$. We obtain asymptotic expressions for the collapse time and relaxation time, above and below the critical point $\Pi_c$, respectively. In the limit $\Delta\to0$ the total collapse for $\Pi>\Pi_\text{c}$ takes a time $T \propto 1/(\Pi-\Pi_\text{c})$, while there is an exponential relaxation below $\Pi_\text{c}$ with a relaxation time $\tau \propto 1/[\Pi_\text{c}-\Pi]$., Comment: 12 pages, 6 figures

The outcome of SIR epidemics with heterogeneous infective lifetimes, or heterogeneous susceptibilities, can be mapped onto a directed percolation process on the underlying contact network. In this paper we study SIR models where heterogeneity is a result of the degree dependence of disease transmission rates. We develop numerical methods to determine the epidemic threshold, the epidemic probability and epidemic size close to the threshold for configuration model contact networks with arbitrary degree distribution and an arbitrary matrix of transmission rates (dependent on transmitting and receiving node degree). For the special case of separable transmission rates we obtain analytical expressions for these quantities. We propose a categorization of spreading processes based on the ratio of the probability of an epidemic and the expected size of an epidemic, and demonstrate that this ratio has a complex dependence on the degree distribution and the degree-dependent transmission rates. For scale-free contact networks and transmission rates that are power functions of transmitting and receiving node degrees, the epidemic threshold may be finite even when the degree distribution powerlaw exponent is below $\gamma < 3$. We give an expression, in terms of the degree distribution and transmission rate exponents, for the limit at which the epidemic threshold vanishes. We find that the expected epidemic size and the probability of an epidemic may grow nonlinearly above the epidemic threshold, with exponents that depend not only on the degree distribution powerlaw exponent, but on the parameters of the transmission rate degree dependence functions, in contrast to ordinary directed percolation and previously studied variations of the SIR model., Comment: 32 pages, 9 figures

In many systems consisting of interacting subsystems, the complex interactions between elements can be represented using multilayer networks. However percolation, key to understanding connectivity and robustness, is not trivially generalised to multiple layers. We describe a generalisation of percolation to multilayer networks: weak multiplex percolation. A node belongs to a connected component if at least one of its neighbours in each layer is in this component. We fully describe the critical phenomena of this process. In particular, in two layers, with finite second moments of the degree distributions, an unusual continuous transition with quadratic growth above the threshold. When the second moments diverge, the singularity is determined by the asymptotics of the degree distributions, creating a rich set of critical behaviours. In three or more layers we find a discontinuous hybrid transition which persists even in highly heterogeneous degree distributions, becoming continuous only when the powerlaw exponent reaches $1+ 1/(M-1)$ for $M$ layers., Comment: 42 pages, 8 figures

We explored the statistics of filtering of simple patterns on a number of deterministic and random graphs as a tractable simple example of information processing in complex systems. In this problem, multiple inputs map to the same output, and the statistics of filtering is represented by the distribution of this degeneracy. For a few simple filter patterns on a ring we obtained an exact solution of the problem and described numerically more difficult filter setups. For each of the filter patterns and networks we found a few numbers essentially describing the statistics of filtering and compared them for different networks. Our results for networks with diverse architectures appear to be essentially determined by two factors: whether the graphs structure is deterministic or random, and the vertex degree. We find that filtering in random graphs produces a much richer statistics than in deterministic graphs. This statistical richness is reduced by increasing the graph's degree., Comment: 21 pages, 8 figures, 3 tables

We describe the critical behavior of weak multiplex percolation, a generalization of percolation to multiplex or interdependent networks. A node can determine its active or inactive status simply by referencing neighboring nodes. This is not the case for the more commonly studied generalization of percolation to multiplex networks, the mutually connected clusters, which requires an interconnecting path within each layer between any two vertices in the giant mutually connected component. We study the emergence of a giant connected component of active nodes under the weak percolation rule, finding several non-typical phenomena. In two layers, the giant component emerges with a continuos phase transition, but with quadratic growth above the critical threshold. In three or more layers, a discontinuous hybrid transition occurs, similar to that found in the giant mutually connected component. In networks with asymptotically powerlaw degree distributions, defined by the decay exponent $\gamma$, the discontinuity vanishes but at $\gamma=1.5$ in three layers, more generally at $\gamma = 1+ 1/(M-1)$ in $M$ layers., Comment: 11 pages, 6 figures

In filtering, each output is produced by a certain number of different inputs. We explore the statistics of this degeneracy in an explicitly treatable filtering problem in which filtering performs the maximal compression of relevant information contained in inputs (arrays of zeroes and ones). This problem serves as a reference model for the statistics of filtering and related sampling problems. The filter patterns in this problem conveniently allow a microscopic, combinatorial consideration. This allows us to find the statistics of outputs, namely the exact distribution of output degeneracies, for arbitrary input sizes. We observe that the resulting degeneracy distribution of outputs decays as $e^{-c\log^\alpha \!d}$ with degeneracy $d$, where $c$ is a constant and exponent $\alpha>1$, i.e. faster than a power law. Importantly, its form essentially depends on the size of the input data set, appearing to be closer to a power-law dependence for small data set sizes than for large ones. We demonstrate that for sufficiently small input data set sizes typical for empirical studies, this distribution could be easily perceived as a power law. We extend our results to filter patterns of various sizes and demonstrate that the shortest filter pattern provides the maximum informative representations of the inputs., Comment: 17 pages, 8 figures, 1 table, Supplementary Material

The giant mutually connected component (GMCC) of an interdependent or multiplex network collapses with a discontinuous hybrid transition under random damage to the network. If the nodes to be damaged are selected in a targeted way, the collapse of the GMCC may occur significantly sooner. Finding the minimal damage set which destroys the largest mutually connected component of a given interdependent network is a computationally prohibitive simultaneous optimization problem. We introduce a simple heuristic strategy -- Effective Multiplex Degree -- for targeted attack on interdependent networks that leverages the indirect damage inherent in multiplex networks to achieve a damage set smaller than that found by any other non computationally intensive algorithm. We show that the intuition from single layer networks that decycling (damage of the $2$-core) is the most effective way to destroy the giant component, does not carry over to interdependent networks, and in fact such approaches are worse than simply removing the highest degree nodes., Comment: 9 pages, 9 figures

Many real complex systems cannot be represented by a single network, but due to multiple sub-systems and types of interactions, must be represented as a multiplex network. This is a set of nodes which exist in several layers, with each layer having its own kind of edges, represented by different colours. An important fundamental structural feature of networks is their resilience to damage, the percolation transition. Generalisation of these concepts to multiplex networks requires careful definition of what we mean by connected clusters. We consider two different definitions. One, a rigorous generalisation of the single-layer definition leads to a strong non-local rule, and results in a dramatic change in the response of the system to damage. The giant component collapses discontinuously in a hybrid transition characterised by avalanches of diverging mean size. We also consider another definition, which imposes weaker conditions on percolation and allows local calculation, and also leads to different sized giant components depending on whether we consider an activation or pruning process. This 'weak' process exhibits both continuous and discontinuous transitions., Comment: arXiv admin note: text overlap with arXiv:1312.3814

We present the theory of the k-core pruning process (progressive removal of nodes with degree less than k) in uncorrelated random networks. We derive exact equations describing this process and the evolution of the network structure, and solve them numerically and, in the critical regime of the process, analytically. We show that the pruning process exhibits three different behaviors depending on whether the mean degree of the initial network is above, equal to, or below the threshold _c corresponding to the emergence of the giant k-core. We find that above the threshold the network relaxes exponentially to the k-core. The system manifests the phenomenon known as "critical slowing down", as the relaxation time diverges when tends to _c. At the threshold, the dynamics become critical characterized by a power-law relaxation (1/t^2). Below the threshold, a long-lasting transient process (a "plateau" stage) occurs. This transient process ends with a collapse in which the entire network disappears completely. The duration of the process diverges when tends to _c. We show that the critical dynamics of the pruning are determined by branching processes of spreading damage. Clusters of nodes of degree exactly k are the evolving substrate for these branching processes. Our theory completely describes this branching cascade of damage in uncorrelated networks by providing the time dependent distribution function of branching. These theoretical results are supported by our simulations of the $k$-core pruning in Erdos-Renyi graphs., Comment: 12 pages, 10 figures

We investigate the dynamics of a broad class of stochastic copying processes on a network that includes examples from population genetics (spatially-structured Wright-Fisher models), ecology (Hubbell-type models), linguistics (the utterance selection model) and opinion dynamics (the voter model) as special cases. These models all have absorbing states of fixation where all the nodes are in the same state. Earlier studies of these models showed that the mean time when this occurs can be made to grow as different powers of the network size by varying the the degree distribution of the network. Here we demonstrate that this effect can also arise if one varies the asymmetry of the copying dynamics whilst holding the degree distribution constant. In particular, we show that the mean time to fixation can be accelerated even on homogeneous networks when certain nodes are very much more likely to be copied from than copied to. We further show that there is a complex interplay between degree distribution and asymmetry when they may co-vary; and that the results are robust to correlations in the network or the initial condition., Comment: 11 pages, 5 figures

We reveal the nature of the avalanche collapse of the giant viable component in multiplex networks under perturbations such as random damage. Specifically, we identify latent critical clusters associated with the avalanches of random damage. Divergence of their mean size signals the approach to the hybrid phase transition from one side, while there are no critical precursors on the other side. We find that this discontinuous transition occurs in scale-free multiplex networks whenever the mean degree of at least one of the interdependent networks does not diverge., Comment: 4 pages, 5 figures

We introduce time variation in the flip-rates of the Voter Model. This type of generalisation is relevant to models of ageing in language change, allowing the representation of changes in speakers' learning rates over their lifetime and may be applied to any other similar model in which interaction rates at the microscopic level change with time. The mean time taken to reach consensus varies in a nontrivial way with the rate of change of the flip-rates, varying between bounds given by the mean consensus times for static homogeneous (the original Voter Model) and static heterogeneous flip-rates. By considering the mean time between interactions for each agent, we derive excellent estimates of the mean consensus times and exit probabilities for any time scale of flip-rate variation. The scaling of consensus times with population size on complex networks is correctly predicted, and is as would be expected for the ordinary voter model. Heterogeneity in the initial distribution of opinions has a strong effect, considerably reducing the mean time to consensus, while increasing the probability of survival of the opinion which initially occupies the most slowly changing agents. The mean times to reach consensus for different states are very different. An opinion originally held by the fastest changing agents has a smaller chance to succeed, and takes much longer to do so than an evenly distributed opinion., Comment: 16 pages, 6 figures

We introduce the heterogeneous-$k$-core, which generalizes the $k$-core, and contrast it with bootstrap percolation. Vertices have a threshold $k_i$ which may be different at each vertex. If a vertex has less than $k_i$ neighbors it is pruned from the network. The heterogeneous-$k$-core is the sub-graph remaining after no further vertices can be pruned. If the thresholds $k_i$ are $1$ with probability $f$ or $k \geq 3$ with probability $(1-f)$, the process forms one branch of an activation-pruning process which demonstrates hysteresis. The other branch is formed by ordinary bootstrap percolation. We show that there are two types of transitions in this heterogeneous-$k$-core process: the giant heterogeneous-$k$-core may appear with a continuous transition and there may be a second, discontinuous, hybrid transition. We compare critical phenomena, critical clusters and avalanches at the heterogeneous-$k$-core and bootstrap percolation transitions. We also show that network structure has a crucial effect on these processes, with the giant heterogeneous-$k$-core appearing immediately at a finite value for any $f > 0$ when the degree distribution tends to a power law $P(q) \sim q^{-\gamma}$ with $\gamma < 3$., Comment: 10 pages, 4 figures

We consider bootstrap percolation on uncorrelated complex networks. We obtain the phase diagram for this process with respect to two parameters: $f$, the fraction of vertices initially activated, and $p$, the fraction of undamaged vertices in the graph. We observe two transitions: the giant active component appears continuously at a first threshold. There may also be a second, discontinuous, hybrid transition at a higher threshold. Avalanches of activations increase in size as this second critical point is approached, finally diverging at this threshold. We describe the existence of a special critical point at which this second transition first appears. In networks with degree distributions whose second moment diverges (but whose first moment does not), we find a qualitatively different behavior. In this case the giant active component appears for any $f>0$ and $p>0$, and the discontinuous transition is absent. This means that the giant active component is robust to damage, and also is very easily activated. We also formulate a generalized bootstrap process in which each vertex can have an arbitrary threshold., Comment: 9 pages, 3 figures

Dependencies between types in object-oriented software can be viewed as directed graphs, with types as nodes and dependencies as edges. The in-degree and out-degree distributions of such graphs have quite different forms, with the former resembling a power-law distribution and the latter an exponential distribution. This effect appears to be independent of application or type relationship. A simple generative model is proposed to explore the proposition that the difference arises because the programmer is aware of the out-degree of a type but not of its in-degree. The model reproduces the two distributions, and compares reasonably well to those observed in 14 different type relationships across 12 different Java applications., Comment: 9 pages, 8 figures

We investigate a set of stochastic models of biodiversity, population genetics, language evolution and opinion dynamics on a network within a common framework. Each node has a state, 0 < x_i < 1, with interactions specified by strengths m_{ij}. For any set of m_{ij} we derive an approximate expression for the mean time to reach fixation or consensus (all x_i=0 or 1). Remarkably in a case relevant to language change this time is independent of the network structure., Comment: 4+epsilon pages, two-column, RevTeX4, 3 eps figures; version accepted by Phys. Rev. Lett

We present a mathematical formulation of a theory of language change. The theory is evolutionary in nature and has close analogies with theories of population genetics. The mathematical structure we construct similarly has correspondences with the Fisher-Wright model of population genetics, but there are significant differences. The continuous time formulation of the model is expressed in terms of a Fokker-Planck equation. This equation is exactly soluble in the case of a single speaker and can be investigated analytically in the case of multiple speakers who communicate equally with all other speakers and give their utterances equal weight. Whilst the stationary properties of this system have much in common with the single-speaker case, time-dependent properties are richer. In the particular case where linguistic forms can become extinct, we find that the presence of many speakers causes a two-stage relaxation, the first being a common marginal distribution that persists for a long time as a consequence of ultimate extinction being due to rare fluctuations., Comment: 21 pages, 17 figures