Ferromagnetic ordering in graphs with arbitrary degree distribution

Understanding the relationship between the heterogeneous structure of complex networks and cooperative phenomena occurring on them remains a key problem in network science. Mean-field theories of spin models on networks constitute a fundamental tool to tackle this problem and a cornerstone of statistical physics, with an impressive number of applications in condensed matter, biology, and computer science. In this work we derive the mean-field equations for the equilibrium behavior of vector spin models on high-connectivity random networks with an arbitrary degree distribution and with randomly weighted links. We demonstrate that the high-connectivity limit of spin models on networks is not universal in that it depends on the full degree distribution. Such nonuniversal behavior is akin to a remarkable mechanism that leads to the breakdown of the central limit theorem when applied to the distribution of effective local fields. Traditional mean-field theories on fully-connected models, such as the Curie–Weiss, the Kuramoto, and the Sherrington–Kirkpatrick model, are only valid if the network degree distribution is highly concentrated around its mean degree. We obtain a series of results that highlight the importance of degree fluctuations to the phase diagram of mean-field spin models by focusing on the Kuramoto model of synchronization and on the Sherrington–Kirkpatrick model of spin-glasses. Numerical simulations corroborate our theoretical findings and provide compelling evidence that the present mean-field theory describes an intermediate regime of connectivity, in which the average degree c scales as a power c ∝ N b (b < 1) of the total number N ≫ 1 of spins. Our findings put forward a novel class of spin models that incorporate the effects of degree fluctuations and, at the same time, are amenable to exact analytic solutions.

We analyze maximum entropy random graph ensembles with constrained degrees, drawn from arbitrary degree distributions, and a tuneable number of three-cycles (triangles). We find that such ensembles generally exhibit two transitions, a clustering and a shattering transition, separating three distinct regimes. At the clustering transition, the graphs change from typically having only isolated cycles to forming cycle clusters. At the shattering transition the graphs break up into many small cliques to achieve the desired three-cycle density. The locations of both transitions depend nontrivially on the system size. We derive a general formula for the three-cycle density in the regime of isolated cycles, for graphs with degree distributions that have finite first and second moments. For bounded degree distributions we present further analytical results on cycle densities and phase transition locations, which, while non-rigorous, are all validated via MCMC sampling simulations. We show that the shattering transition is of an entropic nature, occurring for all three-cycle density values, provided the system is large enough.

We generalize the poissonian evolving random graph model of Bauer and Bernard to deal with arbitrary degree distributions. The motivation comes from biological networks, which are well-known to exhibit non poissonian degree distribution. A node is added at each time step and is connected to the rest of the graph by oriented edges emerging from older nodes. This leads to a statistical asymmetry between incoming and outgoing edges. The law for the number of new edges at each time step is fixed but arbitrary. Thermodynamical behavior is expected when this law has a large time limit. Although (by construction) the incoming degree distributions depend on this law, this is not the case for most qualitative features concerning the size distribution of connected components, as long as the law has a finite variance. As the variance grows above 1/4, the average being <1/2, a giant component emerges, which connects a finite fraction of the vertices. Below this threshold, the distribution of component sizes decreases algebraically with a continuously varying exponent. The transition is of infinite order, in sharp contrast with the case of static graphs. The local-in-time profiles for the components of finite size allow to give a refined description of the system.

Large real-world networks typically follow a power-law degree distribution. To study such networks, numerous random graph models have been proposed. However, real-world networks are not drawn at random. Therefore, Brach et al. (27th symposium on discrete algorithms (SODA), pp 1306–1325, 2016) introduced two natural deterministic conditions: (1) a power-law upper bound on the degree distribution (PLB-U) and (2) power-law neighborhoods, that is, the degree distribution of neighbors of each vertex is also upper bounded by a power law (PLB-N). They showed that many real-world networks satisfy both properties and exploit them to design faster algorithms for a number of classical graph problems. We complement their work by showing that some well-studied random graph models exhibit both of the mentioned PLB properties. PLB-U and PLB-N hold with high probability for Chung–Lu Random Graphs and Geometric Inhomogeneous Random Graphs and almost surely for Hyperbolic Random Graphs. As a consequence, all results of Brach et al. also hold with high probability or almost surely for those random graph classes. In the second part we study three classical textsf {NP}-hard optimization problems on PLB networks. It is known that on general graphs with maximum degree Delta, a greedy algorithm, which chooses nodes in the order of their degree, only achieves a Omega (ln Delta )-approximation for Minimum Vertex Cover and Minimum Dominating Set, and a Omega (Delta )-approximation for Maximum Independent Set. We prove that the PLB-U property with beta >2 suffices for the greedy approach to achieve a constant-factor approximation for all three problems. We also show that these problems are APX-hard even if PLB-U, PLB-N, and an additional power-law lower bound on the degree distribution hold. Hence, a PTAS cannot be expected unless P = NP. Furthermore, we prove that all three problems are in MAX SNP if the PLB-U property holds.

The capacity (or maximum flow) of an unicast network is known to be equal to the minimum s-t cut capacity due to the max-flow min-cut theorem. If the topology of a network (or link capacities) is dynamically changing or unknown, it is not so trivial to predict statistical properties on the maximum flow of the network. In this paper, we present a probabilistic analysis for evaluating the accumulate distribution of the minimum s-t cut capacity on random graphs. The graph ensemble treated in this paper consists of weighted graphs with arbitrary specified degree distribution. The main contribution of our work is a lower bound for the accumulate distribution of the minimum s-t cut capacity. From some computer experiments, it is observed that the lower bound derived here reflects the actual statistical behavior of the minimum s-t cut capacity of random graphs with specified degrees.

Social networks generally display a positively skewed degree distribution and higher values for clustering coefficient and degree assortativity than would be expected from the degree sequence. For some types of simulation studies, these properties need to be varied in the artificial networks over which simulations are to be conducted. Various algorithms to generate networks have been described in the literature but their ability to control all three of these network properties is limited. We introduce a spatially constructed algorithm that generates networks with constrained but arbitrary degree distribution, clustering coefficient and assortativity. Both a general approach and specific implementation are presented. The specific implementation is validated and used to generate networks with a constrained but broad range of property values.

The “classical” random graph models, in particular G(n,p), are “homogeneous,” in the sense that the degrees (for example) tend to be concentrated around a typical value. Many graphs arising in the real world do not have this property, having, for example, power-law degree distributions. Thus there has been a lot of recent interest in defining and studying “inhomogeneous” random graph models. One of the most studied properties of these new models is their “robustness”, or, equivalently, the “phase transition” as an edge density parameter is varied. For G(n,p), p = c/n, the phase transition at c = 1 has been a central topic in the study of random graphs for well over 40 years. Many of the new inhomogeneous models are rather complicated; although there are exceptions, in most cases precise questions such as determining exactly the critical point of the phase transition are approachable only when there is independence between the edges. Fortunately, some models studied have this property already, and others can be approximated by models with independence. Here we introduce a very general model of an inhomogeneous random graph with (conditional) independence between the edges, which scales so that the number of edges is linear in the number of vertices. This scaling corresponds to the p = c/n scaling for G(n,p) used to study the phase transition; also, it seems to be a property of many large real-world graphs. Our model includes as special cases many models previously studied. We show that, under one very weak assumption (that the expected number of edges is “what it should be”), many properties of the model can be determined, in particular the critical point of the phase transition, and the size of the giant component above the transition. We do this by relating our random graphs to branching processes, which are much easier to analyze. We also consider other properties of the model, showing, for example, that when there is a giant component, it is “stable”: for a typical random graph, no matter how we add or delete o(n) edges, the size of the giant component does not change by more than o(n). © 2007 Wiley Periodicals, Inc. Random Struct. Alg., 31, 3–122, 2007

Advances in exponential random graph ( p*) models applied to a large social network

In this paper we introduce new models of random graphs, arising from Latin squares which include random Cayley graphs as a special case. We investigate some properties of these graphs including their clique, independence and chromatic numbers, their expansion properties as well as their connectivity and Hamiltonicity. The results obtained are compared with other models of random graphs and several similarities and differences are pointed out. For many properties our results for the general case are as strong as the known results for random Cayley graphs and sometimes improve the previously best results for the Cayley case. © 2011 Wiley Periodicals, Inc. Random Struct. Alg., 2011

Entropy of the ensemble of acyclic graphs

In this talk, we discuss some abstract models of sensor networks and also some basic algorithmic problems. We first go through the basic architecture and communication capabilities of a single, ultra small sized, sensor. We distinguish broadcast capability (at a certain radius r around the sensor) and directional broadcasting at an angle (which can be reduced to almost zero in case of optical communication). Based on such features, we first go through the (well known) model of Random Geometric Graphs and its threshold properties for connectivity, average degree, chromatic number etc. Then, we define the model of Random Sector Graphs, in order to capture unidirectional sending at an angle. We provide thresholds for connectivity, clique and chromatic number there also and point out its differences from the Random Geometric Graphs model. We also discuss other models related to ways of configuring sensor nets, such as “for each point (sensor) connect to its k nearest neighbours” or “for each sensor, connect to k points within its communication range r”. We also refer to the Combinatorial model of Random Intersection Graphs (which abstracts the geometry out). It is interesting to see that some properties (e.g., sensors density necessary for connectivity) are “invariant” in all these models. We then go through some fundamental algorithmic issues and their existing solutions. We first discuss the problem of finding efficient local protocols aiming in propagating a local event E to a sink. We discuss several efficiency measures like hop count, energy spent, short paths. Randomised local protocols are then sketched, useful for robust and energy efficient local event reporting in case of a sensor net with failures at some nodes (or with energy depleted nodes). In this framework we indicate some energy-time tradeoffs for propagation of local information and notice their similarity with tradeoffs between area and time in VLSI circuits. The problem of propagation of local information in areas with obstacles (lakes) is also discussed. We then examine efficient routing policies and especially refer to greedy techniques in the case of existence of virtual coordinates. The topology control problem is the next topic. We present several local protocols for configuring the net and establishing global properties like connectivity or short paths. Finally, we discuss problems related to scenarios where a sensor net is “embedded” into a larger fixed-connectivity network. There, we provide the abstraction of a “sensors cloud” and discuss the possible enhancement of properties like maximum flow or connectivity or diameter of the resulting hybrid network. Our talk aims in showing the emergence of a new algorithmic subfield, useful to the pragmatic considerations in the actual design and control of such networks.

Random graph models are used to describe the complex structure of real-world networks in diverse fields of knowledge. Studying their behaviour and fitting properties are still critical challenges that, in general, require model-specific techniques. An important line of research is to develop generic methods able to fit and select the best model among a collection. Approaches based on spectral density (i.e. distribution of the graph adjacency matrix eigenvalues) appeal to that purpose: they apply to different random graph models. Also, they can benefit from the theoretical background of random matrix theory. This work investigates the convergence properties of model fitting procedures based on the graph spectral density and the corresponding cumulative distribution function. We also review the convergence of the spectral density for the most widely used random graph models. Moreover, we explore through simulations the limits of these graph spectral density convergence results, particularly in the case of the block model, where only partial results have been established. random graphs, spectral density, model fitting, model selection, convergence.

A key challenge in the social network community is the problem of network generation -- that is, how can we create synthetic networks that match characteristics traditionally found in most real world networks? Important characteristics that are present in social networks include a power law degree distribution, small diameter, and large amounts of clustering. However, most current network generators, such as the Chung Lu and Kronecker models, largely ignore the clustering present in a graph and focus on preserving other network statistics, such as the power law distribution. Models such as the exponential random graph model have a transitivity parameter that can capture clustering, but they are computationally difficult to learn, making scaling to large real world networks intractable. In this work, we propose an extension to the Chung Lu random graph model, the Transitive Chung Lu (TCL) model, which incorporates the notion transitive edges. Specifically, it combines the standard Chung Lu model with edges that are formed through transitive closure (e.g., by connecting a 'friend of a friend'). We prove TCL's expected degree distribution is equal to the degree distribution of the original input graph, while still providing the ability to capture the clustering in the network. The single parameter required by our model can be learned in seconds on graphs with millions of edges, networks can be generated in time that is linear in the number of edges. We demonstrate the performance of TCL on four real-world social networks, including an email dataset with hundreds of thousands of nodes and millions of edges, showing TCL generates graphs that match the degree distribution, clustering coefficients and hop plots of the original networks.

Random graphs with latent geometric structure are popular models of social and biological networks, with applications ranging from network user profiling to circuit design. These graphs are also of purely theoretical interest within computer science, probability and statistics. A fundamental initial question regarding these models is: when are these random graphs affected by their latent geometry and when are they indistinguishable from simpler models without latent structure, such as the Erdős–Rényi graph $${\mathcal {G}}(n, p)$$ ? We address this question for two of the most well-studied models of random graphs with latent geometry—the random intersection and random geometric graph. Our results are as follows: (a) The random intersection graph is defined by sampling n random sets $$S_1, S_2, \ldots , S_n$$ by including each element of a set of size d in each $$S_i$$ independently with probability $$\delta $$ , and including the edge $$\{i, j\}$$ if $$S_i \cap S_j \ne \emptyset $$ . We prove that the random intersection graph converges in total variation to an Erdős–Rényi graph if $$d = {\tilde{\omega }}(n^3)$$ , and does not if $$d = o(n^3)$$ , for both dense and sparse edge densities p. This resolves an open problem in Fill et al. (Random Struct Algorithms 16(2):156–176, 2000), Rybarczyk (Random Struct Algorithms 38(1–2):205–234, 2011) and Kim et al. (Random Struct Algorithms 52(4):662–679, 2018). The same result was obtained simultaneously and independently by Bubeck et al. (When random intersection graphs lose geometry. Manuscript, 2019). (b) We strengthen the preceding argument to show that the matrix of random intersection sizes $$|S_i \cap S_j|$$ converges in total variation to a symmetric matrix with independent Poisson entries. This yields the first total variation convergence result for $$\tau $$ -random intersection graphs, where the edge $$\{i, j\}$$ is included if $$|S_i \cap S_j| \ge \tau $$ . More precisely, our results imply that, if p is bounded away from 1, then the $$\tau $$ -random intersection graph with edge density p converges to $${\mathcal {G}}(n, p)$$ if $$d = \omega (\tau ^3 n^3)$$ . (c) The random geometric graph on $${\mathbb {S}}^{d - 1}$$ is defined by sampling $$X_1, X_2, \dots , X_n$$ uniformly at random from $${\mathbb {S}}^{d - 1}$$ and including the edge $$\{i, j\}$$ if $$\Vert X_i - X_j \Vert _2 \le \tau $$ . A result of Bubeck et al. (Random Struct Algorithms 49:503–532, 2016) showed that this model converges to $${\mathcal {G}}(n, p)$$ in total variation, where p is chosen so that the models have matching edge densities, as long as $$d = \omega (n^3)$$ . It was conjectured in Bubeck et al. (2016) that this threshold decreases drastically for p small. We make the first progress towards this conjecture by showing convergence if $$d = {\tilde{\omega }}\left( \min \{ pn^3, p^2 n^{7/2} \} \right) $$ . Our proofs are a hybrid of combinatorial arguments, direct couplings and applications of information inequalities. Previous upper bounds on the total variation distance between random graphs with latent geometry and $${\mathcal {G}}(n, p)$$ have typically not been both combinatorial and information-theoretic, while this interplay is essential to the sharpness of our bounds.

Ad-regular graph has largest or rst (adjacency matrix) eigenvalue 1 = d. Consider for an even d 4, a random d-regular graph model formed from d=2 uniform, independent permutations on f1;:::;ng. We shall show that for any >0 we have all eigenvalues aside from 1 = d are bounded by 2 p d 1+ with probability 1 O(n ), where =d p d 1+ 1 =2e 1. We also show that this probability is at most 1 c=n 0 , for a constant c and a 0 that is either or +1 (\more often than + 1). We prove related theorems for other models of random graphs, including models with d odd. These theorems resolve the conjecture of Alon, that says that for any > 0a ndd, the second largest eigenvalue of \most random dregular graphs are at most 2 p d 1+ (Alon did not specify precisely what \most should mean or what model of random graph one should take).