Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
15topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 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.

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.

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.

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.

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.

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.

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$.

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.

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.

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.

preprint2011arXiv

Explosive Percolation is Continuous, but with Unusual Finite Size Behavior

We study four Achlioptas type processes with &#34;explosive&#34; 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.

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.

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.

preprint2011arXiv

Random Sequential Renormalization and Agglomerative Percolation in Networks: Application to Erd&#34;os-R&#39;enyi and Scale-free Graphs

We study the statistical behavior under random sequential renormalization(RSR) of several network models including Erd&#34;os R&#39;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 &#39;supernodes&#39; as clusters) are much easier to study using the fast Newman-Ziff algorithm for percolation, allowing us to obtain very high statistics.

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$.

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&#39;s evolution.

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.

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 &#34;triangle&#34;) 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 &#34;cluster cores&#34;: 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.

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.

preprint2008arXiv

Comment on &#34;Central limit behavior in deterministic dynamical systems&#34;

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.

preprint2006arXiv

Monte Carlo Algorithm for Least Dependent Non-Negative Mixture Decomposition

We propose a simulated annealing algorithm (called SNICA for &#34;stochastic non-negative independent component analysis&#34;) 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