BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Peter Grassberger

Peter Grassberger appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.dis-nnphysics.soc-phphysics.data-anSocial and Information Networkscond-mat.softmath.PRPopulations and EvolutionGenomicsInformation Theorymath.ITmath.STNetworking and Internet Architecturephysics.chem-phphysics.comp-phphysics.geo-phq-fin.STQuantitative MethodsStatistics Theory

Catalog footprint

What is connected

41works
19topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

Kardar-Parisi-Zhang type dynamics with periodic tilt dependence of the propagation velocity in 1+1 dimensions

We consider the evolution of interfaces with a diffusive term and a generalized Kardar-Parisi-Zhang (KPZ) non-linearity, which results in a propagation velocity that depends periodically on the tilt of the interface. Using large scale simulations of a model class with these properties in 1+1 dimensions, we show that the fluctuations are in general still in the KPZ universality class, but a new universality class seems to appear in the limit of weak non-linearity. We argue that this is the typical behavior of any interface model with periodic tilt dependence.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Revisiting a Low-Dimensional Model with Short Range Interactions and Mean Field Critical Behavior

In all local low-dimensional models, scaling at critical points deviates from mean field behavior -- with one possible exception. This exceptional model with ``ordinary" behavior is an inherently non-equilibrium model studied some time ago by H.-M. Broker and myself. In simulations, its 2-dimensional version suggested that two critical exponents were mean-field, while a third one showed very small deviations. Moreover, the numerics agreed almost perfectly with an explicit mean field model. In the present paper we present simulations with much higher statistics, both for 2d and 3d. In both cases we find that the deviations of all critical exponents from their mean field values are non-leading corrections, and that the scaling is {\it precisely} of mean field type. As in the original paper, we propose that the mechanism for this is ``confusion", a strong randomization of the phases of feed-backs that can occur in non-equilibrium systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Immunization and targeted destruction of networks using explosive percolation

A new method (`explosive immunization' (EI)) is proposed for immunization and targeted destruction of networks. It combines the explosive percolation (EP) paradigm with the idea of maintaining a fragmented distribution of clusters. The ability of each node to block the spread of an infection (or to prevent the existence of a large cluster of connected nodes) is estimated by a score. The algorithm proceeds by first identifying low score nodes that should not be vaccinated/destroyed, analogously to the links selected in EP if they do not lead to large clusters. As in EP, this is done by selecting the worst node (weakest blocker) from a finite set of randomly chosen `candidates'. Tests on several real-world and model networks suggest that the method is more efficient and faster than any existing immunization strategy. Due to the latter property it can deal with very large networks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Percolation transitions in the survival of interdependent agents on multiplex networks, catastrophic cascades, and SOS

The "SOS" in the title does not refer to the international distress signal, but to "solid-on-solid" (SOS) surface growth. The catastrophic cascades are those observed by Buldyrev {\it et al.} in interdependent networks, which we re-interpret as multiplex networks with agents that can only survive if they mutually support each other, and whose survival struggle we map onto an SOS type growth model. This mapping not only reveals non-trivial structures in the phase space of the model, but also leads to a new and extremely efficient simulation algorithm. We use this algorithm to study interdependent agents on duplex Erdös-Rényi (ER) networks and on lattices with dimensions 2, 3, 4, and 5. We obtain new and surprising results in all these cases, and we correct statements in the literature for ER networks and for 2-d lattices. In particular, we find that $d=4$ is the upper critical dimension, that the percolation transition is continuous for $d\leq 4$ but -- at least for $d\neq 3$ -- not in the universality class of ordinary percolation. For ER networks we verify that the cluster statistics is exactly described by mean field theory, but find evidence that the cascade process is not. For $d=5$ we find a first order transition as for ER networks, but we find also that small clusters have a nontrivial mass distribution that scales at the transition point. Finally, for $d=2$ with intermediate range dependency links we propose a scenario different from that proposed in W. Li {\it et al.}, PRL {\bf 108}, 228702 (2012).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Oslo model, hyperuniformity, and the quenched Edwards-Wilkinson model

We present simulations of the 1-dimensional Oslo rice pile model in which the critical height at each site is randomly reset after each toppling. We use the fact that the stationary state of this sandpile model is hyperuniform to reach system of sizes $> 10^7$. Most previous simulations were seriously flawed by important finite size corrections. We find that all critical exponents have values consistent with simple rationals: $ν=4/3$ for the correlation length exponent, $D =9/4$ for the fractal dimension of avalanche clusters, and $z=10/7 $ for the dynamical exponent. In addition we relate the hyperuniformity exponent to the correlation length exponent $ν$. Finally we discuss the relationship with the quenched Edwards-Wilkinson (qEW) model, where we find in particular that the local roughness exponent is $α_{\rm loc} = 1$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymmetry of cross correlations between intra-day and overnight volatilities

We point out a stunning time asymmetry in the short time cross correlations between intra-day and overnight volatilities (absolute values of log-returns of stock prices). While overnight volatility is significantly (and positively) correlated with the intra-day volatility during the \textit{following} day (allowing thus non-trivial predictions), it is much less correlated with the intra-day volatility during the \textit{preceding} day. While the effect is not unexpected in view of previous observations, its robustness and extreme simplicity are remarkable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase Transitions in Cooperative Coinfections: Simulation Results for Networks and Lattices

We study the spreading of two mutually cooperative diseases on different network topologies, and with two microscopic realizations, both of which are stochastic versions of an SIR type model studied by us recently in mean field approximation. There it had been found that cooperativity can lead to first-order spreading/extinction transitions. However, due to the rapid mixing implied by the mean field assumption, first order transitions required non-zero initial densities of sick individuals. For the stochastic model studied here the results depend strongly on the underlying network. First order transitions are found when there are few short but many long loops: (i) No first order transitions exist on trees and on 2-d lattices with local contacts (ii) They do exist on Erdos-Renyi (ER) networks, on d-dimensional lattices with d >= 4, and on 2-d lattices with sufficiently long-ranged contacts; (iii) On 3-d lattices with local contacts the results depend on the microscopic details of the implementation; (iv) While single infected seeds can always lead to infinite epidemics on regular lattices, on ER networks one sometimes needs finite initial densities of infected nodes; (v) In all cases the first order transitions are actually "hybrid", i.e. they display also power law scaling usually associated with second order transitions. On regular lattices, our model can also be interpreted as the growth of an interface due to cooperative attachment of two species of particles. Critically pinned interfaces in this model seem to be in different universality classes than standard critically pinned interfaces in models with forbidden overhangs. Finally, the detailed results mentioned above hold only when both diseases propagate along the same network of links. If they use different links, results can be rather different in detail, but are similar overall.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Continuum Time Limit of Reaction-Diffusion Systems

The parity conserving branching-annihilating random walk (pc-BARW) model is a reaction-diffusion system on a lattice where particles can branch into $m$ offsprings with even $m$ and hop to neighboring sites. If two or more particles land on the same site, they immediately annihilate pairwise. In this way the number of particles is preserved modulo two. It is well known that the pc-BARW with $m=2$ in 1 spatial dimension has no phase transition (it is always subcritical), if the hopping is described by a continuous time random walk. In contrast, the $m=2$ 1-d pc-BARW has a phase transition when formulated in discrete time, but we show that the continuous time limit is non-trivial: When the time step $δt\to 0$, the branching and hopping probabilities at the critical point scale with different powers of $δt$. These powers are different for different microscopic realizations. Although this phenomenon is not observed in some other reaction-diffusion systems like, e.g. the contact process, we argue that it should be generic and not restricted to the 1-d pc-BARW model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Outbreaks of coinfections: the critical role of cooperativity

Modeling epidemic dynamics plays an important role in studying how diseases spread, predicting their future course, and designing strategies to control them. In this letter, we introduce a model of SIR (susceptible-infected-removed) type which explicitly incorporates the effect of {\it cooperative coinfection}. More precisely, each individual can get infected by two different diseases, and an individual already infected with one disease has an increased probability to get infected by the other. Depending on the amount of this increase, we observe different threshold scenarios. Apart from the standard continuous phase transition for single disease outbreaks, we observe continuous transitions where both diseases must coexist, but also discontinuous transitions are observed, where a finite fraction of the population is already affected by both diseases at the threshold. All our results are obtained in a mean field model using rate equations, but we argue that they should hold also in more general frameworks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Polymer collapse and crystallization in bond fluctuation models

While the $Θ$-collapse of single long polymers in bad solvents is usually a continuous (tri-critical) phase transition, there are exceptions where it is preempted by a discontinuous crystallization (liquid $\leftrightarrow$ solid) transition. For a version of the bond-fluctuation model (a model where monomers are represented as $2\times 2\times 2$ cubes, and bonds can have lengths between 2 and $\sqrt{10}$) it was recently shown by F. Rampf {\it et al.} that there exist distinct collapse and crystallization transitions for long but {\it finite} chains. But as the chain length goes to infinity, both transition temperatures converge to the same $T^*$, i.e. infinitely long polymers collapse immediately into a solid state. We explain this by the observation that polymers crystallize in the Rampf {\it et al.} model into a non-trivial cubic crystal structure (the `A15' or `Cr$_3$Si' Frank-Kasper structure) which has many degenerate ground states and, as a consequence, Bloch walls. If one controlls the polymer growth such that only one ground state is populated and Bloch walls are completely avoided, the liquid-solid transition is a smooth cross-over without any sharp transition at all.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SIR epidemics with long range infection in one dimension

We study epidemic processes with immunization on very large 1-dimensional lattices, where at least some of the infections are non-local, with rates decaying as power laws p(x) ~ x^{-sigma-1} for large distances x. When starting with a single infected site, the cluster of infected sites stays always bounded if $σ>1$ (and dies with probability 1, of its size is allowed to fluctuate down to zero), but the process can lead to an infinite epidemic for sigma <1. For sigma <0 the behavior is essentially of mean field type, but for 0 < sigma <= 1 the behavior is non-trivial, both for the critical and for supercritical cases. For critical epidemics we confirm a previous prediction that the critical exponents controlling the correlation time and the correlation length are simply related to each other, and we verify detailed field theoretic predictions for sigma --> 1/3. For sigma = 1 we find generic power laws with continuously varying exponents even in the supercritical case, and confirm in detail the predicted Kosterlitz-Thouless nature of the transition. Finally, the mass N(t) of supercritical clusters seems to grow for 0 < sigma < 1 like a stretched exponential. The latter implies that networks embedded in 1-d space with power-behaved link distributions have infinite intrinsic dimension (based on the graph distance), but are not small world.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Two-dimensional SIR epidemics with long range infection

We extend a recent study of susceptible-infected-removed epidemic processes with long range infection (referred to as I in the following) from 1-dimensional lattices to lattices in two dimensions. As in I we use hashing to simulate very large lattices for which finite size effects can be neglected, in spite of the assumed power law $p({\bf x})\sim |{\bf x}|^{-σ-2}$ for the probability that a site can infect another site a distance vector ${\bf x}$ apart. As in I we present detailed results for the critical case, for the supercritical case with $σ= 2$, and for the supercritical case with $0< σ< 2$. For the latter we verify the stretched exponential growth of the infected cluster with time predicted by M. Biskup. For $σ=2$ we find generic power laws with $σ-$dependent exponents in the supercritical phase, but no Kosterlitz-Thouless (KT) like critical point as in 1-d. Instead of diverging exponentially with the distance from the critical point, the correlation length increases with an inverse power, as in an ordinary critical point. Finally we study the dependence of the critical exponents on $σ$ in the regime $0<σ<2$, and compare with field theoretic predictions. In particular we discuss in detail whether the critical behavior for $σ$ slightly less than 2 is in the short range universality class, as conjectured recently by F. Linder {\it et al.}. As in I we also consider a modified version of the model where only some of the contacts are long range, the others being between nearest neighbors. If the number of the latter reaches the percolation threshold, the critical behavior is changed but the supercritical behavior stays qualitatively the same.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Agglomerative Percolation on Bipartite Networks: A Novel Type of Spontaneous Symmetry Breaking

Ordinary bond percolation (OP) can be viewed as a process where clusters grow by joining them pairwise, by adding links chosen randomly one by one from a set of predefined `virtual' links. In contrast, in agglomerative percolation (AP) clusters grow by choosing randomly a `target cluster' and joining it with all its neighbors, as defined by the same set of virtual links. Previous studies showed that AP is in different universality classes from OP for several types of (virtual) networks (linear chains, trees, Erdos-Renyi networks), but most surprising were the results for 2-d lattices: While AP on the triangular lattice was found to be in the OP universality class, it behaved completely differently on the square lattice. In the present paper we explain this striking violation of universality by invoking bipartivity. While the square lattice is a bipartite graph, the triangular lattice is not. In conformity with this we show that AP on the honeycomb and simple cubic (3-d) lattices -- both of which are bipartite -- are also not in the OP universality classes. More precisely, we claim that this violation of universality is basically due to a Z_2 symmetry that is spontaneously broken at the percolation threshold. We also discuss AP on bipartite random networks and suitable generalizations of AP on k-partite graphs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Discontinuous Percolation Transitions in Epidemic Processes, Surface Depinning in Random Media and Hamiltonian Random Graphs

Discontinuous percolation transitions and the associated tricritical points are manifest in a wide range of both equilibrium and non-equilibrium cooperative phenomena. To demonstrate this, we present and relate the continuous and first order behaviors in two different classes of models: The first are generalized epidemic processes (GEP) that describe in their spatially embedded version - either on or off a regular lattice - compact or fractal cluster growth in random media at zero temperature. A random graph version of GEP is mapped onto a model previously proposed for complex social contagion. We compute detailed phase diagrams and compare our numerical results at the tricritical point in d = 3 with field theory predictions of Janssen et al. [Phys. Rev. E 70, 026114 (2004)]. The second class consists of exponential ("Hamiltonian", or formally equilibrium) random graph models and includes the Strauss and the 2-star model, where 'chemical potentials' control the densities of links, triangles or 2-stars. When the chemical potentials in either graph model are O(logN), the percolation transition can coincide with a first order phase transition in the density of links, making the former also discontinuous. Hysteresis loops can then be of mixed order, with second order behavior for decreasing link fugacity, and a jump (first order) when it increases.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Information theoretic aspects of the two-dimensional Ising model

We present numerical results for various information theoretic properties of the square lattice Ising model. First, using a bond propagation algorithm, we find the difference $2H_L(w) - H_{2L}(w)$ between entropies on cylinders of finite lengths $L$ and 2L with open end cap boundaries, in the limit $L\to\infty$. This essentially quantifies how the finite length correction for the entropy scales with the cylinder circumference $w$. Secondly, using the transfer matrix, we obtain precise estimates for the information needed to specify the spin state on a ring encircling an infinite long cylinder. Combining both results we obtain the mutual information between the two halves of a cylinder (the "excess entropy" for the cylinder), where we confirm with higher precision but for smaller systems results recently obtained by Wilms et al. -- and we show that the mutual information between the two halves of the ring diverges at the critical point logarithmically with $w$. Finally we use the second result together with Monte Carlo simulations to show that also the excess entropy of a straight line of $n$ spins in an infinite lattice diverges at criticality logarithmically with $n$. We conjecture that such logarithmic divergence happens generically for any one-dimensional subset of sites at any 2-dimensional second order phase transition. Comparing straight lines on square and triangular lattices with square loops and with lines of thickness 2, we discuss questions of universality.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

PageRank and rank-reversal dependence on the damping factor

PageRank (PR) is an algorithm originally developed by Google to evaluate the importance of web pages. Considering how deeply rooted Google's PR algorithm is to gathering relevant information or to the success of modern businesses, the question of rank-stability and choice of the damping factor (a parameter in the algorithm) is clearly important. We investigate PR as a function of the damping factor d on a network obtained from a domain of the World Wide Web, finding that rank-reversal happens frequently over a broad range of PR (and of d). We use three different correlation measures, Pearson, Spearman, and Kendall, to study rank-reversal as d changes, and show that the correlation of PR vectors drops rapidly as d changes from its frequently cited value, $d_0=0.85$. Rank-reversal is also observed by measuring the Spearman and Kendall rank correlation, which evaluate relative ranks rather than absolute PR. Rank-reversal happens not only in directed networks containing rank-sinks but also in a single strongly connected component, which by definition does not contain any sinks. We relate rank-reversals to rank-pockets and bottlenecks in the directed network structure. For the network studied, the relative rank is more stable by our measures around $d=0.65$ than at $d=d_0$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Randomness, Information, and Complexity

We review possible measures of complexity which might in particular be applicable to situations where the complexity seems to arise spontaneously. We point out that not all of them correspond to the intuitive (or "naive") notion, and that one should not expect a unique observable of complexity. One of the main problems is to distinguish complex from disordered systems. This and the fact that complexity is closely related to information requires that we also give a review of information measures. We finally concentrate on quantities which measure in some way or other the difficulty of classifying and forecasting sequences of discrete symbols, and study them in simple examples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sampling properties of directed networks

For many real-world networks only a small "sampled" version of the original network may be investigated; those results are then used to draw conclusions about the actual system. Variants of breadth-first search (BFS) sampling, which are based on epidemic processes, are widely used. Although it is well established that BFS sampling fails, in most cases, to capture the IN-component(s) of directed networks, a description of the effects of BFS sampling on other topological properties are all but absent from the literature. To systematically study the effects of sampling biases on directed networks, we compare BFS sampling to random sampling on complete large-scale directed networks. We present new results and a thorough analysis of the topological properties of seven different complete directed networks (prior to sampling), including three versions of Wikipedia, three different sources of sampled World Wide Web data, and an Internet-based social network. We detail the differences that sampling method and coverage can make to the structural properties of sampled versions of these seven networks. Most notably, we find that sampling method and coverage affect both the bow-tie structure, as well as the number and structure of strongly connected components in sampled networks. In addition, at low sampling coverage (i.e. less than 40%), the values of average degree, variance of out-degree, degree auto-correlation, and link reciprocity are overestimated by 30% or more in BFS-sampled networks, and only attain values within 10% of the corresponding values in the complete networks when sampling coverage is in excess of 65%. These results may cause us to rethink what we know about the structure, function, and evolution of real-world directed networks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A review of Monte Carlo simulations of polymers with PERM

In this review, we describe applications of the pruned-enriched Rosenbluth method (PERM), a sequential Monte Carlo algorithm with resampling, to various problems in polymer physics. PERM produces samples according to any given prescribed weight distribution, by growing configurations step by step with controlled bias, and correcting "bad" configurations by "population control". The latter is implemented, in contrast to other population based algorithms like e.g. genetic algorithms, by depth-first recursion which avoids storing all members of the population at the same time in computer memory. The problems we discuss all concern single polymers (with one exception), but under various conditions: Homopolymers in good solvents and at the $Θ$ point, semi-stiff polymers, polymers in confining geometries, stretched polymers undergoing a forced globule-linear transition, star polymers, bottle brushes, lattice animals as a model for randomly branched polymers, DNA melting, and finally -- as the only system at low temperatures, lattice heteropolymers as simple models for protein folding. PERM is for some of these problems the method of choice, but it can also fail. We discuss how to recognize when a result is reliable, and we discuss also some types of bias that can be crucial in guiding the growth into the right directions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Agglomerative Percolation in Two Dimensions

We study a process termed "agglomerative percolation" (AP) in two dimensions. Instead of adding sites or bonds at random, in AP randomly chosen clusters are linked to all their neighbors. As a result the growth process involves a diverging length scale near a critical point. Picking target clusters with probability proportional to their mass leads to a runaway compact cluster. Choosing all clusters equally leads to a continuous transition in a new universality class for the square lattice, while the transition on the triangular lattice has the same critical exponents as ordinary percolation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Are Percolation Transitions always Sharpened by Making Networks Interdependent?

We study a model for coupled networks introduced recently by Buldyrev et al., Nature 464, 1025 (2010), where each node has to be connected to others via two types of links to be viable. Removing a critical fraction of nodes leads to a percolation transition that has been claimed to be more abrupt than that for uncoupled networks. Indeed, it was found to be discontinuous in all cases studied. Using an efficient new algorithm we verify that the transition is discontinuous for coupled Erdos-Renyi networks, but find it to be continuous for fully interdependent diluted lattices. In 2 and 3 dimension, the order parameter exponent $β$ is larger than in ordinary percolation, showing that the transition is less sharp, i.e. further from discontinuity, than for isolated networks. Possible consequences for spatially embedded networks are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Clustering Drives Assortativity and Community Structure in Ensembles of Networks

Clustering, assortativity, and communities are key features of complex networks. We probe dependencies between these attributes and find that ensembles with strong clustering display both high assortativity by degree and prominent community structure, while ensembles with high assortativity are much less biased towards clustering or community structure. Further, clustered networks can amplify small homophilic bias for trait assortativity. This marked asymmetry suggests that transitivity, rather than homophily, drives the standard nonsocial/social network dichotomy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exact solutions for mass-dependent irreversible aggregations

We consider the mass-dependent aggregation process (k+1)X -> X, given a fixed number of unit mass particles in the initial state. One cluster is chosen proportional to its mass and is merged into one either with k-neighbors in one dimension, or -- in the well-mixed case -- with k other clusters picked randomly. We find the same combinatorial exact solutions for the probability to find any given configuration of particles on a ring or line, and in the well-mixed case. The mass distribution of a single cluster exhibits scaling laws and the finite size scaling form is given. The relation to the classical sum kernel of irreversible aggregation is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Explosive Percolation is Continuous, but with Unusual Finite Size Behavior

We study four Achlioptas type processes with "explosive" percolation transitions. All transitions are clearly continuous, but their finite size scaling functions are not entire holomorphic. The distributions of the order parameter, the relative size $s_{\rm max}/N$ of the largest cluster, are double-humped. But -- in contrast to first order phase transitions -- the distance between the two peaks decreases with system size $N$ as $N^{-η}$ with $η> 0$. We find different positive values of $β$ (defined via $< s_{\rm max}/N > \sim (p-p_c)^β$ for infinite systems) for each model, showing that they are all in different universality classes. In contrast, the exponent $Θ$ (defined such that observables are homogeneous functions of $(p-p_c)N^Θ$) is close to -- or even equal to -- 1/2 for all models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Irreversible Aggregation and Network Renormalization

Irreversible aggregation is revisited in view of recent work on renormalization of complex networks. Its scaling laws and phase transitions are related to percolation transitions seen in the latter. We illustrate our points by giving the complete solution for the probability to find any given state in an aggregation process $(k+1)X\to X$, given a fixed number of unit mass particles in the initial state. Exactly the same probability distributions and scaling are found in one dimensional systems (a trivial network) and well-mixed solutions. This reveals that scaling laws found in renormalization of complex networks do not prove that they are self-similar.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Percolation Theory on Interdependent Networks Based on Epidemic Spreading

We consider percolation on interdependent locally treelike networks, recently introduced by Buldyrev et al., Nature 464, 1025 (2010), and demonstrate that the problem can be simplified conceptually by deleting all references to cascades of failures. Such cascades do exist, but their explicit treatment just complicates the theory -- which is a straightforward extension of the usual epidemic spreading theory on a single network. Our method has the added benefits that it is directly formulated in terms of an order parameter and its modular structure can be easily extended to other problems, e.g. to any number of interdependent networks, or to networks with dependency links.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Random Sequential Renormalization and Agglomerative Percolation in Networks: Application to Erd"os-R'enyi and Scale-free Graphs

We study the statistical behavior under random sequential renormalization(RSR) of several network models including Erd"os R'enyi (ER) graphs, scale-free networks and an annealed model (AM) related to ER graphs. In RSR the network is locally coarse grained by choosing at each renormalization step a node at random and joining it to all its neighbors. Compared to previous (quasi-)parallel renormalization methods [C.Song et.al], RSR allows a more fine-grained analysis of the renormalization group (RG) flow, and unravels new features, that were not discussed in the previous analyses. In particular we find that all networks exhibit a second order transition in their RG flow. This phase transition is associated with the emergence of a giant hub and can be viewed as a new variant of percolation, called agglomerative percolation. We claim that this transition exists also in previous graph renormalization schemes and explains some of the scaling laws seen there. For critical trees it happens as N/N0 -> 0 in the limit of large systems (where N0 is the initial size of the graph and N its size at a given RSR step). In contrast, it happens at finite N/N0 in sparse ER graphs and in the annealed model, while it happens for N/N0 -> 1 on scale-free networks. Critical exponents seem to depend on the type of the graph but not on the average degree and obey usual scaling relations for percolation phenomena. For the annealed model they agree with the exponents obtained from a mean-field theory. At late times, the networks exhibit a star-like structure in agreement with the results of Radicchi et. al. While degree distributions are of main interest when regarding the scheme as network renormalization, mass distributions (which are more relevant when considering 'supernodes' as clusters) are much easier to study using the fast Newman-Ziff algorithm for percolation, allowing us to obtain very high statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Random Sequential Renormalization of Networks I: Application to Critical Trees

We introduce the concept of Random Sequential Renormalization (RSR) for arbitrary networks. RSR is a graph renormalization procedure that locally aggregates nodes to produce a coarse grained network. It is analogous to the (quasi-)parallel renormalization schemes introduced by C. Song {\it et al.} (Nature {\bf 433}, 392 (2005)) and studied more recently by F. Radicchi {\it et al.} (Phys. Rev. Lett. {\bf 101}, 148701 (2008)), but much simpler and easier to implement. In this first paper we apply RSR to critical trees and derive analytical results consistent with numerical simulations. Critical trees exhibit three regimes in their evolution under RSR: (i) An initial regime $N_0^ν\lesssim N<N_0$, where $N$ is the number of nodes at some step in the renormalization and $N_0$ is the initial size. RSR in this regime is described by a mean field theory and fluctuations from one realization to another are small. The exponent $ν=1/2$ is derived using random walk arguments. The degree distribution becomes broader under successive renormalization -- reaching a power law, $p_k\sim 1/k^γ$ with $γ=2$ and a variance that diverges as $N_0^{1/2}$ at the end of this regime. Both of these results are derived based on a scaling theory. (ii) An intermediate regime for $N_0^{1/4}\lesssim N \lesssim N_0^{1/2}$, in which hubs develop, and fluctuations between different realizations of the RSR are large. Crossover functions exhibiting finite size scaling, in the critical region $N\sim N_0^{1/2} \to \infty$, connect the behaviors in the first two regimes. (iii) The last regime, for $1 \ll N\lesssim N_0^{1/4}$, is characterized by the appearance of star configurations with a central hub surrounded by many leaves. The distribution of sizes where stars first form is found numerically to be a power law up to a cutoff that scales as $N_0^{ν_{star}}$ with $ν_{star}\approx 1/4$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Edge direction and the structure of networks

Directed networks are ubiquitous and are necessary to represent complex systems with asymmetric interactions---from food webs to the World Wide Web. Despite the importance of edge direction for detecting local and community structure, it has been disregarded in studying a basic type of global diversity in networks: the tendency of nodes with similar numbers of edges to connect. This tendency, called assortativity, affects crucial structural and dynamic properties of real-world networks, such as error tolerance or epidemic spreading. Here we demonstrate that edge direction has profound effects on assortativity. We define a set of four directed assortativity measures and assign statistical significance by comparison to randomized networks. We apply these measures to three network classes---online/social networks, food webs, and word-adjacency networks. Our measures (i) reveal patterns common to each class, (ii) separate networks that have been previously classified together, and (iii) expose limitations of several existing theoretical models. We reject the standard classification of directed networks as purely assortative or disassortative. Many display a class-specific mixture, likely reflecting functional or historical constraints, contingencies, and forces guiding the system's evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lower Bounds on Mutual Information

We correct claims about lower bounds on mutual information (MI) between real-valued random variables made in A. Kraskov {\it et al.}, Phys. Rev. E {\bf 69}, 066138 (2004). We show that non-trivial lower bounds on MI in terms of linear correlations depend on the marginal (single variable) distributions. This is so in spite of the invariance of MI under reparametrizations, because linear correlations are not invariant under them. The simplest bounds are obtained for Gaussians, but the most interesting ones for practical purposes are obtained for uniform marginal distributions. The latter can be enforced in general by using the ranks of the individual variables instead of their actual values, in which case one obtains bounds on MI in terms of Spearman correlation coefficients. We show with gene expression data that these bounds are in general non-trivial, and the degree of their (non-)saturation yields valuable insight.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sequence alignment, mutual information, and dissimilarity measures for constructing phylogenies

Existing sequence alignment algorithms use heuristic scoring schemes which cannot be used as objective distance metrics. Therefore one relies on measures like the p- or log-det distances, or makes explicit, and often simplistic, assumptions about sequence evolution. Information theory provides an alternative, in the form of mutual information (MI) which is, in principle, an objective and model independent similarity measure. MI can be estimated by concatenating and zipping sequences, yielding thereby the "normalized compression distance". So far this has produced promising results, but with uncontrolled errors. We describe a simple approach to get robust estimates of MI from global pairwise alignments. Using standard alignment algorithms, this gives for animal mitochondrial DNA estimates that are strikingly close to estimates obtained from the alignment free methods mentioned above. Our main result uses algorithmic (Kolmogorov) information theory, but we show that similar results can also be obtained from Shannon theory. Due to the fact that it is not additive, normalized compression distance is not an optimal metric for phylogenetics, but we propose a simple modification that overcomes the issue of additivity. We test several versions of our MI based distance measures on a large number of randomly chosen quartets and demonstrate that they all perform better than traditional measures like the Kimura or log-det (resp. paralinear) distances. Even a simplified version based on single letter Shannon entropies, which can be easily incorporated in existing software packages, gave superior results throughout the entire animal kingdom. But we see the main virtue of our approach in a more general way. For example, it can also help to judge the relative merits of different alignment algorithms, by estimating the significance of specific alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Interacting Branching Process as a Simple Model of Innovation

We describe innovation in terms of a generalized branching process. Each new invention pairs with any existing one to produce a number of offspring, which is Poisson distributed with mean p. Existing inventions die with probability p/τat each generation. In contrast to mean field results, no phase transition occurs; the chance for survival is finite for all p > 0. For τ= \infty, surviving processes exhibit a bottleneck before exploding super-exponentially - a growth consistent with a law of accelerating returns. This behavior persists for finite τ. We analyze, in detail, the asymptotic behavior as p \to 0.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Clustering Phase Transitions and Hysteresis: Pitfalls in Constructing Network Ensembles

Ensembles of networks are used as null models in many applications. However, simple null models often show much less clustering than their real-world counterparts. In this paper, we study a model where clustering is enhanced by means of a fugacity term as in the Strauss (or "triangle") model, but where the degree sequence is strictly preserved -- thus maintaining the quenched heterogeneity of nodes found in the original degree sequence. Similar models had been proposed previously in [R. Milo et al., Science 298, 824 (2002)]. We find that our model exhibits phase transitions as the fugacity is changed. For regular graphs (identical degrees for all nodes) with degree k > 2 we find a single first order transition. For all non-regular networks that we studied (including Erdos - Renyi and scale-free networks) we find multiple jumps resembling first order transitions, together with strong hysteresis. The latter transitions are driven by the sudden emergence of "cluster cores": groups of highly interconnected nodes with higher than average degrees. To study these cluster cores visually, we introduce q-clique adjacency plots. We find that these cluster cores constitute distinct communities which emerge spontaneously from the triangle generating process. Finally, we point out that cluster cores produce pitfalls when using the present (and similar) models as null models for strongly clustered networks, due to the very strong hysteresis which effectively leads to broken ergodicity on realistic time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Local persistence in directed percolation

We reconsider the problem of local persistence in directed site percolation. We present improved estimates of the persistence exponent in all dimensions from 1+1 to 7+1, obtained by new algorithms and by improved implementations of existing ones. We verify the strong corrections to scaling for 2+1 and 3+1 dimensions found in previous analyses, but we show that scaling is much better satisfied for very large and very small dimensions. For d > 4 (d is the spatial dimension), the persistence exponent depends non-trivially on d, in qualitative agreement with the non-universal values calculated recently by Fuchs {\it et al.} (J. Stat. Mech.: Theor. Exp. P04015 (2008)). These results are mainly based on efficient simulations of clusters evolving under the time reversed dynamics with a permanently active site and a particular survival condition discussed in Fuchs {\it et al.}. These simulations suggest also a new critical exponent $ζ$ which describes the growth of these clusters conditioned on survival, and which turns out to be the same as the exponent, η+δin standard notation, of surviving clusters under the standard DP evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Logarithmic corrections in (4+1)-dimensional directed percolation

We simulate directed site percolation on two lattices with 4 spatial and 1 time-like dimensions (simple and body-centered hypercubic in space) with the standard single cluster spreading scheme. For efficiency, the code uses the same ingredients (hashing, histogram re-weighing, and improved estimators) as described in Phys. Rev. {\bf E 67}, 036101 (2003). Apart from providing the most precise estimates for $p_c$ on these lattices, we provide a detailed comparison with the logarithmic corrections calculated by Janssen and Stenull [Phys. Rev. {\bf E 69}, 016125 (2004)]. Fits with the leading logarithmic terms alone would give estimates of the powers of these logarithms which are too big by typically 50%. When the next-to-leading terms are included, each of the measured quantities (the average number of sites wetted at time $t$, their average distance from the seed, and the probability of cluster survival) can be fitted nearly perfectly. But these fits would not be mutually consistent. With a consistent set of fit parameters, one obtains still much improvement over the leading log - approximation. In particular we show that there is one combination of these three observables which seems completely free of logarithmic terms.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Scaling of loop-erased walks in 2 to 4 dimensions

We simulate loop-erased random walks on simple (hyper-)cubic lattices of dimensions 2,3, and 4. These simulations were mainly motivated to test recent two loop renormalization group predictions for logarithmic corrections in $d=4$, simulations in lower dimensions were done for completeness and in order to test the algorithm. In $d=2$, we verify with high precision the prediction $D=5/4$, where the number of steps $n$ after erasure scales with the number $N$ of steps before erasure as $n\sim N^{D/2}$. In $d=3$ we again find a power law, but with an exponent different from the one found in the most precise previous simulations: $D = 1.6236\pm 0.0004$. Finally, we see clear deviations from the naive scaling $n\sim N$ in $d=4$. While they agree only qualitatively with the leading logarithmic corrections predicted by several authors, their agreement with the two-loop prediction is nearly perfect.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Comment on "Central limit behavior in deterministic dynamical systems"

We check claims for a generalized central limit theorem holding at the Feigenbaum (infinite bifurcation) point of the logistic map, made recently by U. Tirnakli, C. Beck, and C. Tsallis (Phys. Rev. {\bf 75}, 040106(R) (2007)). We show that there is no obvious way that these claims can be made consistent with high statistics simulations. We also refute more recent claims by the same authors that extend the claims made in the above reference.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Earthquake recurrence as a record breaking process

Extending the central concept of recurrence times for a point process to recurrent events in space-time allows us to characterize seismicity as a record breaking process using only spatiotemporal relations among events. Linking record breaking events with edges between nodes in a graph generates a complex dynamical network isolated from any length, time or magnitude scales set by the observer. For Southern California, the network of recurrences reveals new statistical features of seismicity with robust scaling laws. The rupture length and its scaling with magnitude emerges as a generic measure for distance between recurrent events. Further, the relative separations for subsequent records in space (or time) form a hierarchy with unexpected scaling properties.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Monte Carlo Algorithm for Least Dependent Non-Negative Mixture Decomposition

We propose a simulated annealing algorithm (called SNICA for "stochastic non-negative independent component analysis") for blind decomposition of linear mixtures of non-negative sources with non-negative coefficients. The de-mixing is based on a Metropolis type Monte Carlo search for least dependent components, with the mutual information between recovered components as a cost function and their non-negativity as a hard constraint. Elementary moves are shears in two-dimensional subspaces and rotations in three-dimensional subspaces. The algorithm is geared at decomposing signals whose probability densities peak at zero, the case typical in analytical spectroscopy and multivariate curve resolution. The decomposition performance on large samples of synthetic mixtures and experimental data is much better than that of traditional blind source separation methods based on principal component analysis (MILCA, FastICA, RADICAL) and chemometrics techniques (SIMPLISMA, ALS, BTEM) The source codes of SNICA, MILCA and the MI estimator are freely available online at http://www.fz-juelich.de/nic/cs/software

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Polymers grafted to porous membranes

We study a single flexible chain molecule grafted to a membrane which has pores of size slightly larger than the monomer size. On both sides of the membrane there is the same solvent. When this solvent is good, i.e. when the polymer is described by a self avoiding walk, it can fairly easily penetrate the membrane, so that the average number of membrane crossings tends, for chain length $N\to\infty$, to a positive constant. The average numbers of monomers on either side of the membrane diverges in this limit, although their ratio becomes infinite. For a poor solvent, in contrast, the entire polymer is located, for large $N$, on one side of the membrane. For good and for theta solvents (ideal polymers) we find scaling laws, whose exponents can in the latter case be easily understood from the behaviour of random walks.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Slow Logarithmic Decay of Magnetization in the Zero Temperature Dynamics of an Ising Spin Chain: Analogy to Granular Compaction

We study the zero temperature coarsening dynamics in an Ising chain in presence of a dynamically induced field that favors locally the `-' phase compared to the `+' phase. At late times, while the `+' domains still coarsen as $t^{1/2}$, the `-' domains coarsen slightly faster as $t^{1/2}\log (t)$. As a result, at late times, the magnetization decays slowly as, $m(t)=-1 +{\rm const.}/{\log (t)}$. We establish this behavior both analytically within an independent interval approximation (IIA) and numerically. In the zero volume fraction limit of the `+' phase, we argue that the IIA becomes asymptotically exact. Our model can be alternately viewed as a simple Ising model for granular compaction. At late times in our model, the system decays into a fully compact state (where all spins are `-') in a slow logarithmic manner $\sim 1/{\log (t)}$, a fact that has been observed in recent experiments on granular systems.

0 citations0 reviewsNo code linkedOpen access