Search

In the analysis of large-scale network data, a fundamental operation is the comparison of observed phenomena to the predictions provided by null models: when we find an interesting structure in a family of real networks, it is important to ask whether this structure is also likely to arise in random networks with similar characteristics to the real ones. A long-standing challenge in network analysis has been the relative scarcity of reasonable null models for networks; arguably the most common such model has been the configuration model, which starts with a graph $G$ and produces a random graph with the same node degrees as $G$. This leads to a very weak form of null model, since fixing the node degrees does not preserve many of the crucial properties of the network, including the structure of its subgraphs. Guided by this challenge, we propose a new family of network null models that operate on the $k$-core decomposition. For a graph $G$, the $k$-core is its maximal subgraph of minimum degree $k$; and the core number of a node $v$ in $G$ is the largest $k$ such that $v$ belongs to the $k$-core of $G$. We provide the first efficient sampling algorithm to solve the following basic combinatorial problem: given a graph $G$, produce a random graph sampled nearly uniformly from among all graphs with the same sequence of core numbers as $G$. This opens the opportunity to compare observed networks $G$ with random graphs that exhibit the same core numbers, a comparison that preserves aspects of the structure of $G$ that are not captured by more local measures like the degree sequence. We illustrate the power of this core-based null model on some fundamental tasks in network analysis, including the enumeration of networks motifs.

The performance of wireless networks is fundamentally limited by the aggregate interference, which depends on the spatial distributions of the interferers, channel conditions, and user traffic patterns (or queueing dynamics). These factors usually exhibit spatial and temporal correlations and thus make the performance of large-scale networks environment-dependent (i.e., dependent on network topology, locations of the blockages, etc.). The correlation can be exploited in protocol designs (e.g., spectrum-, load-, location-, energy-aware resource allocations) to provide efficient wireless services. For this, accurate system-level performance characterization and evaluation with spatial-temporal correlation are required. In this context, stochastic geometry models and random graph techniques have been used to develop analytical frameworks to capture the spatial-temporal interference correlation in large-scale wireless networks. The objective of this article is to provide a tutorial on the stochastic geometry analysis of large-scale wireless networks that captures the spatial-temporal interference correlation (and hence the signal-to-interference ratio (SIR) correlation). We first discuss the importance of spatial-temporal performance analysis, different parameters affecting the spatial-temporal correlation in the SIR, and the different performance metrics for spatial-temporal analysis. Then we describe the methodologies to characterize spatial-temporal SIR correlations for different network configurations (independent, attractive, repulsive configurations), shadowing scenarios, user locations, queueing behavior, relaying, retransmission, and mobility. We conclude by outlining future research directions in the context of spatial-temporal analysis of emerging wireless communications scenarios.

Consider a graph $G$, representing a social network. Assume that initially each node is colored either black or white, which corresponds to a positive or negative opinion regarding a consumer product or a technological innovation. In the majority model, in each round all nodes simultaneously update their color to the most frequent color among their connections. Experiments on the graph data from the real world social networks (SNs) suggest that if all nodes in an extremely small set of high-degree nodes, often referred to as the elites, agree on a color, that color becomes the dominant color at the end of the process. We propose two countermeasures that can be adopted by individual nodes relatively easily and guarantee that the elites will not have this disproportionate power to engineer the dominant output color. The first countermeasure essentially requires each node to make some new connections at random while the second one demands the nodes to be more reluctant towards changing their color (opinion). We verify their effectiveness and correctness both theoretically and experimentally. We also investigate the majority model and a variant of it when the initial coloring is random on the real world SNs and several random graph models. In particular, our results on the Erd\H{o}s-R\'{e}nyi and regular random graphs confirm or support several theoretical findings or conjectures by the prior work regarding the threshold behavior of the process. Finally, we provide theoretical and experimental evidence for the existence of a poly-logarithmic bound on the expected stabilization time of the majority model., Comment: Accepted to the International Joint Conference on Artificial Intelligence (IJCAI)

From social interactions to the human brain, higher-order networks are key to describe the underlying network geometry and topology of many complex systems. While it is well known that network structure strongly affects its function, the role that network topology and geometry has on the emerging dynamical properties of higher-order networks is yet to be clarified. In this perspective, the spectral dimension plays a key role since it determines the effective dimension for diffusion processes on a network. Despite its relevance, a theoretical understanding of which mechanisms lead to a finite spectral dimension, and how this can be controlled, represents nowadays still a challenge and is the object of intense research. Here we introduce two non-equilibrium models of hyperbolic higher-order networks and we characterize their network topology and geometry by investigating the interwined appearance of small-world behavior, $\delta$-hyperbolicity and community structure. We show that different topological moves determining the non-equilibrium growth of the higher-order hyperbolic network models induce tunable values of the spectral dimension, showing a rich phenomenology which is not displayed in random graph ensembles. In particular, we observe that, if the topological moves used to construct the higher-order network increase the area$/$volume ratio, the spectral dimension continuously decreases, while the opposite effect is observed if the topological moves decrease the area$/$volume ratio. Our work reveals a new link between the geometry of a network and its diffusion properties, contributing to a better understanding of the complex interplay between network structure and dynamics., Comment: 14 pages, 10 figures

By studying the commuting graphs of conjugacy classes of the sequence of Heisenberg groups $H_{2n+1}(p)$ and their limit $H_\infty(p)$ we find pseudo-random behavior (and the random graph in the limiting case). This makes a nice case study for transfer of information between finite and infinite objects. Some of this behavior transfers to the problem of understanding what makes understanding the character theory of the uni-upper-triangular group (mod p) “wild.” Our investigations in this paper may be seen as a meditation on the question: is randomness simple or is it complicated?

Here we present the entropic dynamics formalism for networks. That is, a framework for the dynamics of graphs meant to represent a network derived from the principle of maximum entropy and the rate of transition is obtained taking into account the natural information geometry of probability distributions. We apply this framework to the Gibbs distribution of random graphs obtained with constraints on the node connectivity. The information geometry for this graph ensemble is calculated and the dynamical process is obtained as a diffusion equation. We compare the steady state of this dynamics to degree distributions found on real-world networks.

Simple binary-state coordination models are widely used to study collective socio-economic phenomena such as the spread of innovations or the adoption of products on social networks. The common trait of these systems is the occurrence of large-scale coordination events taking place abruptly, in the form of a cascade process, as a consequence of small perturbations of an apparently stable state. The conditions for the occurrence of cascade instabilities have been largely analysed in the literature, however for the same coordination models no sufficient attention was given to the relation between structural properties of (Nash) equilibria and possible outcomes of dynamical equilibrium selection. Using methods from the statistical physics of disordered systems, the present work investigates both analytically and numerically, the statistical properties of such Nash equilibria on networks, focusing mostly on random graphs. We provide an accurate description of these properties, which is then exploited to shed light on the mechanisms behind the onset of coordination/miscoordination on large networks. This is done studying the most common processes of dynamical equilibrium selection, such as best response, bounded-rational dynamics and learning processes. In particular, we show that well beyond the instability region, full coordination is still globally stochastically stable, however equilibrium selection processes with low stochasticity (e.g. best response) or strong memory effects (e.g. reinforcement learning) can be prevented from achieving full coordination by being trapped into a large (exponentially in number of agents) set of locally stable Nash equilibria at low/medium coordination (inefficient equilibria). These results should be useful to allow a better understanding of general coordination problems on complex networks., Comment: Revtex style, 56 pages, 21 figures

The primate heteromodal cortex presents an evident functional modularity at a mesoscopic level, with physiological and anatomical evidence pointing to it as likely substrate of long-term memory. In order to investigate some of its properties, a model of multimodular autoassociator is studied. Each of the many modules represents a neocortical functional ensemble of recurrently connected neurons and operates as a Hebbian autoassociator, storing a number of local features which it can recall upon cue. The global memory patterns are made of combinations of features sparsely distributed across the modules. Intermodular connections are modelled as a finite-connectivity random graph. Any pair of features in any respective pair of modules is allowed to be involved in several memory patterns; the coarse-grained modular network dynamics is defined in such a way as to overcome the consequent ambiguity of associations. Effects of long-range homeostatic synaptic scaling on network performance are also assessed. The dynamical process of cued retrieval almost saturates a natural upper bound while producing negligible spurious activation. The extent of finite-size effects on storage capacity is quantitatively evaluated. In the limit of infinite size, the functional relationship between storage capacity and number of features per module reduces to that which other authors found by methods from equilibrium statistical mechanics, which suggests that the origin of the functional form is of a combinatorial nature. In contrast with its apparent inevitability at intramodular level, long-range synaptic scaling results to be of minor relevance to both retrieval and storage capacity, casting doubt on its existence in the neocortex. A conjecture is also posited about how statistical fluctuation of connectivity across the network may underpin spontaneous emergence of semantic hierarchies through learning., Comment: v2: 17 pages, 7+2 figures, pdflatex. Minor Revision: Further observations and references were included, the wording of some paragraphs was improved, and several typos were corrected

Epidemic processes on random graphs or networks are marked by localization of activity that can trap the dynamics into a metastable state, confined to a subextensive part of the network, before visiting an absorbing configuration. Quasistationary (QS) method is a technique to deal with absorbing states for finite sizes and has played a central role in the investigation of epidemic processes on heterogeneous networks where localization is a hallmark. The standard QS method possesses high computer and algorithmic complexity for large systems besides parameters whose choice are not systematic. However, simpler approaches, such as a reflecting boundary condition (RBC), are not able to capture the localization effects as the standard QS method does. In the present work, we propose a QS method that consists of reactivating nodes proportionally to the time they were active along the preceding simulation. The method is compared with the standard QS and RBC methods for the susceptible-infected-susceptible model on complex networks, which is a prototype of a dynamic process with strong and localization effects. We verified that the method performs as well the as standard QS in all investigated simulations, providing the same scaling exponents, epidemic thresholds, and localized phases, thus overcoming the limitations of other simpler approaches. We also report that the present method has significant reduction of the computer and algorithmic complexity than the standard QS methods. So, this method arises as a simpler and efficient tool to analyze localization on heterogeneous structures through QS simulations., Comment: 11 pages, 9 figures

While classical spin systems in random networks have been intensively studied, much less is known about quantum magnets in random graphs. Here, we investigate interacting quantum spins on small-world networks, building on mean-field theory and extensive quantum Monte Carlo simulations. Starting from one-dimensional (1D) rings, we consider two situations: all-to-all interacting and long-range interactions randomly added. The effective infinite dimension of the lattice leads to a magnetic ordering at finite temperature $T_\mathrm{c}$ with mean-field criticality. Nevertheless, in contrast to the classical case, we find two distinct power-law behaviors for $T_\mathrm{c}$ versus the average strength of the extra couplings. This is controlled by a competition between a characteristic length scale of the random graph and the thermal correlation length of the underlying 1D system, thus challenging mean-field theories. We also investigate the fate of a gapped 1D spin chain against the small-world effect., Comment: 12 pages, 8 figures

Exponential random graphs are important to model the structure of real-world complex networks. Here we solve the two-star model with degree-degree correlations in the sparse regime. The model constraints the average correlation between the degrees of adjacent nodes (nearest neighbors) and between the degrees at the end-points of two-stars (next nearest neighbors). We compute exactly the network free energy and show that this model undergoes a first-order transition to a condensed phase. For non-negative degree correlations between next nearest neighbors, the degree distribution inside the condensed phase has a single peak at the largest degree, while for negative degree correlations between next nearest neighbors the condensed phase is characterized by a bimodal degree distribution. We calculate the degree assortativities and show they are non-monotonic functions of the model parameters, with a discontinuous behavior at the first-order transition. The first-order critical line terminates at a second-order critical point, whose location in the phase diagram can be accurately determined. Our results can help to develop more detailed models of complex networks with correlated degrees., Comment: 13 pages, 11 figures

Subgraphs such as cliques, loops and stars form crucial connections in the topologies of real-world networks. Random graph models provide estimates for how often certain subgraphs appear, which in turn can be tested against real-world networks. These subgraph counts, however, crucially depend on the assumed degree distribution. Fitting a degree distribution to network data is challenging, in particular for scale-free networks with power-law degrees. In this paper we develop robust subgraph counts that do not depend on the entire degree distribution, but only on the mean and mean absolute deviation (MAD), summary statistics that are easy to obtain for most real-world networks. By solving an optimization problem, we provide tight (the sharpest possible) bounds for the subgraph counts, for all possible subgraphs, and for all networks with degree distributions that share the same mean and MAD. We identify the extremal random graph that attains the tight bounds as the graph with a specific three-point degree distribution. We leverage the bounds to obtain robust scaling laws for how the numbers of subgraphs grow as function of the network size. The scaling laws indicate that sparse power-law networks are not the most extreme networks in terms of subgraph counts, but dense power-law networks are. The robust bounds are also shown to hold for several real-world data sets., Comment: 14 pages

The separation dimension of a hypergraph $G$ is the smallest natural number $d$ for which there is an embedding of $G$ into $\mathbb{R}^d$, such that any pair of disjoint edges is separated by some hyperplane normal to one of the axes. The perfect separation dimension further requires that any pair of disjoint edges is separated by the same amount of such (pairwise nonparallel) hyperplanes. While it is known that for any fixed $r \ge 2$, the separation dimension of any $n$-vertex $r$-graph is $O(\log n)$, the perfect separation dimension is much larger. In fact, no polynomial upper-bound for the perfect separation dimension of $r$-uniform hypergraphs is known. In our first result we essentially resolve the case $r=2$, i.e. graphs. We prove that the perfect separation dimension of $K_n$ is linear in $n$, up to a small polylogarithmic factor. In fact, we prove it is at least $n/2-1$ and at most $n(\log n)^{1+o(1)}$. Our second result proves that the perfect separation dimension of almost all graphs is also linear in $n$, up to a logarithmic factor. This follows as a special case of a more general result showing that the perfect separation dimension of the random graph $G(n,p)$ is w.h.p. $\Omega(n p /\log n)$ for a wide range of values of $p$, including all constant $p$. Finally, we prove that significantly relaxing perfection to just requiring that any pair of disjoint edges of $K_n$ is separated the same number of times up to a difference of $c \log n$ for some absolute constant $c$, still requires the dimension to be $\Omega(n)$. This is perhaps surprising as it is known that if we allow a difference of $7\log_2 n$, then the dimension reduces to $O(\log n)$., Comment: 22 pages, to appear in J. Combin. Des

Combinatorial optimization algorithms for graph problems are usually designed afresh for each new problem with careful attention by an expert to the problem structure. In this work, we develop a new framework to solve any combinatorial optimization problem over graphs that can be formulated as a single player game defined by states, actions, and rewards, including minimum spanning tree, shortest paths, traveling salesman problem, and vehicle routing problem, without expert knowledge. Our method trains a graph neural network using reinforcement learning on an unlabeled training set of graphs. The trained network then outputs approximate solutions to new graph instances in linear running time. In contrast, previous approximation algorithms or heuristics tailored to NP-hard problems on graphs generally have at least quadratic running time. We demonstrate the applicability of our approach on both polynomial and NP-hard problems with optimality gaps close to 1, and show that our method is able to generalize well: (i) from training on small graphs to testing on large graphs; (ii) from training on random graphs of one type to testing on random graphs of another type; and (iii) from training on random graphs to running on real world graphs.

The quantum approximate optimization algorithm (QAOA) has become a cornerstone of contemporary quantum applications development. Here we show that the \emph{density} of problem constraints versus problem variables acts as a performance indicator. Density is found to correlate strongly with approximation inefficiency for fixed depth QAOA applied to random graph minimization problem instances. Further, the required depth for accurate QAOA solution to graph problem instances scales critically with density. Motivated by Google's recent experimental realization of QAOA, we preform a reanalysis of the reported data reproduced in an ideal noiseless setting. We found that the reported capabilities of instances addressed experimentally by Google, approach a rapid fall-off region in approximation quality experienced beyond intermediate-density. Our findings offer new insight into performance analysis of contemporary quantum optimization algorithms and contradict recent speculation regarding low-depth QAOA performance benefits., Comment: (feedback welcome) 11 pages; 5 composite figures

We develop random graph models where graphs are generated by connecting not only pairs of vertices by edges but also larger subsets of vertices by copies of small atomic subgraphs of arbitrary topology. This allows the for the generation of graphs with extensive numbers of triangles and other network motifs commonly observed in many real world networks. More specifically we focus on maximum entropy ensembles under constraints placed on the counts and distributions of atomic subgraphs and derive general expressions for the entropy of such models. We also present a procedure for combining distributions of multiple atomic subgraphs that enables the construction of models with fewer parameters. Expanding the model to include atoms with edge and vertex labels we obtain a general class of models that can be parametrized in terms of basic building blocks and their distributions that includes many widely used models as special cases. These models include random graphs with arbitrary distributions of subgraphs, random hypergraphs, bipartite models, stochastic block models, models of multilayer networks and their degree corrected and directed versions. We show that the entropy for all these models can be derived from a single expression that is characterized by the symmetry groups of atomic subgraphs., Comment: 15 pages, 2 figures

Many learning tasks require us to deal with graph data which contains rich relational information among elements, leading increasing graph neural network (GNN) models to be deployed in industrial products for improving the quality of service. However, they also raise challenges to model authentication. It is necessary to protect the ownership of the GNN models, which motivates us to present a watermarking method to GNN models in this paper. In the proposed method, an Erdos-Renyi (ER) random graph with random node feature vectors and labels is randomly generated as a trigger to train the GNN to be protected together with the normal samples. During model training, the secret watermark is embedded into the label predictions of the ER graph nodes. During model verification, by activating a marked GNN with the trigger ER graph, the watermark can be reconstructed from the output to verify the ownership. Since the ER graph was randomly generated, by feeding it to a non-marked GNN, the label predictions of the graph nodes are random, resulting in a low false alarm rate (of the proposed work). Experimental results have also shown that, the performance of a marked GNN on its original task will not be impaired. Moreover, it is robust against model compression and fine-tuning, which has shown the superiority and applicability., Comment: https://hzwu.github.io/

The organisation of a network in a maximal set of nodes having at least $k$ neighbours within the set, known as $k$-core decomposition, has been used for studying various phenomena. It has been shown that nodes in the innermost $k$-shells play a crucial role in contagion processes, emergence of consensus, and resilience of the system. It is known that the $k$-core decomposition of many empirical networks cannot be explained by the degree of each node alone, or equivalently, random graph models that preserve the degree of each node (i.e., configuration model). Here we study the $k$-core decomposition of some empirical networks as well as that of some randomised counterparts, and examine the extent to which the $k$-shell structure of the networks can be accounted for by the community structure. We find that preserving the community structure in the randomisation process is crucial for generating networks whose $k$-core decomposition is close to the empirical one. We also highlight the existence, in some networks, of a concentration of the nodes in the innermost $k$-shells into a small number of communities., Comment: 34 pages, 4 tables, and 15 figures. Main + supplementary. Final version accepted for publication

We consider the statistical inference problem of recovering an unknown perfect matching, hidden in a weighted random graph, by exploiting the information arising from the use of two different distributions for the weights on the edges inside and outside the planted matching. A recent work has demonstrated the existence of a phase transition, in the large size limit, between a full and a partial recovery phase for a specific form of the weights distribution on fully connected graphs. We generalize and extend this result in two directions: we obtain a criterion for the location of the phase transition for generic weights distributions and possibly sparse graphs, exploiting a technical connection with branching random walk processes, as well as a quantitatively more precise description of the critical regime around the phase transition., Comment: 19 pages, 8 figures

We develop a formalism to compute the statistics of the second largest eigenpair of weighted sparse graphs with $N\gg 1$ nodes, finite mean connectivity and bounded maximal degree, in cases where the top eigenpair statistics is known. The problem can be cast in terms of optimisation of a quadratic form on the sphere with a fictitious temperature, after a suitable deflation of the original matrix model. We use the cavity and replica methods to find the solution in terms of self-consistent equations for auxiliary probability density functions, which can be solved by an improved population dynamics algorithm enforcing eigenvector orthogonality on-the-fly. The analytical results are in perfect agreement with numerical diagonalisation of large (weighted) adjacency matrices, focussing on the cases of random regular and Erd\H{o}s-R\'enyi graphs. We further analyse the case of sparse Markov transition matrices for unbiased random walks, whose second largest eigenpair describes the non-equilibrium mode with the largest relaxation time. We also show that the population dynamics algorithm with population size $N_P$ does not actually capture the thermodynamic limit $N\to\infty$ as commonly assumed: the accuracy of the population dynamics algorithm has a strongly non-monotonic behaviour as a function of $N_P$, thus implying that an optimal size $N_P^\star=N_P^\star(N)$ must be chosen to best reproduce the results from numerical diagonalisation of graphs of finite size $N$., Comment: 51 pages, 7 figures. Section 3.1 expanded and clarified w.r.t. the published version

We consider evolutionary games on a population whose underlying topology of interactions is determined by a binomial random graph $G(n,p)$. Our focus is on 2-player symmetric games with 2 strategies played between the incident members of such a population. Players update their strategies synchronously. At each round, each player selects the strategy that is the best response to the current set of strategies its neighbours play. We show that such a system reduces to generalised majority and minority dynamics. We show rapid convergence to unanimity for $p$ in a range that depends on a certain characteristic of the payoff matrix. In the presence of a bias among the pure Nash equilibria of the game, we determine a sharp threshold on $p$ above which the largest connected component reaches unanimity with high probability. For $p$ below this critical value, where this does not happen, we identify those substructures inside the largest component that remain discordant throughout the evolution of the system., Comment: 42 pages, 1 figure

In this article, we discuss the random graph, Barab\'asi-Albert (BA) model, and lattice networks from a unified view point, with the parameter $\omega$ with values $1,0,-1$ characterizing these networks, respectively. The parameter is related to the preferential attachment of nodes in the networks and has different weights for the incoming and outgoing links. In addition, we discuss the correspondence between quantum statistics and the networks. Positive and negative $\omega$ correspond to Bose and Fermi-like statistics, respectively, and we obtain the distribution that connects the two. When $\omega$ is positive, it is related to the threshold of Bose-Einstein condensation (BEC). As $\omega$ decreases, the area of the BEC phase is narrowed, and disappears in the limit $\omega=0$. When $\omega$ is negative, nodes have limits in the number of attachments for newly added nodes (outgoing links), which corresponds to Fermi statistics. We also observe the Fermi degeneracy of the network. When $\omega=-1$, a standard Fermion-like network is observed. Fermion networks are realized in the cryptocurrency network "Tangle.", Comment: 19 pages, 6 figures

In this paper we develop a framework for multivariate functional approximation by a suitable Gaussian process via an exchangeable pairs coupling that satisfies a suitable approximate linear regression property, thereby building on work by Barbour (1990) and Kasprzak (2020). We demonstrate the applicability of our results by applying it to joint subgraph counts in an Erd\H{o}s-Renyi random graph model on the one hand and to vectors of weighted, degenerate $U$-processes on the other hand. As a concrete instance of the latter class of examples, we provide a bound for the functional approximation of a vector of success runs of different lengths by a suitable Gaussian process which, even in the situation of just a single run, would be outside the scope of the existing theory.

We investigate spatial random graphs defined on the points of a Poisson process in d-dimensional space, which combine scale-free degree distributions and long-range effects. Every Poisson point is assigned an independent weight. Given the weight and position of the points, we form an edge between any pair of points independently with a probability depending on the two weights of the points and their distance. Preference is given to short edges and connections to vertices with large weights. We characterize the parameter regime where there is a non-trivial percolation phase transition and show that it depends not only on the power-law exponent of the degree distribution but also on a geometric model parameter. We apply this result to characterize robustness of age-based spatial preferential attachment networks.

We present a modified diffusive epidemic process (DEP) that has a finite threshold on scale-free graphs, motivated by the COVID-19 pandemic. The DEP describes the epidemic spreading of a disease in a non-sedentary population, which can describe the spreading of a real disease. Our main modification is to use the Gillespie algorithm with a reaction time t max, exponentially distributed with mean inversely proportional to the node population in order to model the individuals’ interactions. Our simulation results of the modified model on Barabasi–Albert networks are compatible with a continuous absorbing-active phase transition when increasing the average concentration. The transition obeys the mean-field critical exponents β = 1, γ′ = 0 and ν ⊥ = 1/2. In addition, the system presents logarithmic corrections with pseudo-exponents β ̂ = γ ̂ ′ = − 3 / 2 on the order parameter and its fluctuations, respectively. The most evident implication of our simulation results is if the individuals avoid social interactions in order to not spread a disease, this leads the system to have a finite threshold in scale-free graphs.

We establish rapid mixing of the random-cluster Glauber dynamics on random $\Delta$-regular graphs for all $q\ge 1$ and $p2$ the Glauber dynamics on random $\Delta$-regular graphs undergoes an exponential slowdown at $p_u(q,\Delta)$. More precisely, we show that for every $q\ge 1$, $\Delta\ge 3$, and $p, Comment: 35 pages, 6 figures

Resolving a conjecture of K\"uhn and Osthus from 2012, we show that $p= 1/\sqrt{n}$ is the threshold for the random graph $G_{n,p}$ to contain the square of a Hamilton cycle., Comment: 6 pages

In this paper, we investigate the distributed convex optimization problem over a multi-agent system with Markovian switching communication networks. The objective function is the sum of each agent's local objective function, which cannot be known by other agents. The communication network is assumed to switch over a set of weight-balanced directed graphs with a Markovian property.We propose a consensus sub-gradient algorithm with two time-scale step-sizes to handle the uncertainty due to the Markovian switching topologies and the absence of global gradient information. With a proper selection of step-sizes, we prove the almost sure convergence of all agents' local estimates to the same optimal solution when the union graph of the Markovian network' states is strongly connected and the Markovian network is irreducible. Simulations are given for illustration of the results., Comment: 10 pages, 2 figures

We present a unified exact tensor network approach to compute the ground state energy, identify the optimal configuration, and count the number of solutions for spin glasses. The method is based on tensor networks with the Tropical Algebra defined on the semiring. Contracting the tropical tensor network gives the ground state energy; differentiating through the tensor network contraction gives the ground state configuration; mixing the tropical algebra and the ordinary algebra counts the ground state degeneracy. The approach brings together the concepts from graphical models, tensor networks, differentiable programming, and quantum circuit simulation, and easily utilizes the computational power of graphical processing units (GPUs). For applications, we compute the exact ground state energy of Ising spin glasses on square lattice up to 1024 spins, on cubic lattice up to 216 spins, and on 3 regular random graphs up to 220 spins, on a single GPU; We obtain exact ground state energy of (+/-)J Ising spin glass on the chimera graph of D-Wave quantum annealer of 512 qubits in less than 100 seconds and investigate the exact value of the residual entropy of (+/-)J spin glasses on the chimera graph; Finally, we investigate ground-state energy and entropy of 3-state Potts glasses on square lattices up to size 18 x 18. Our approach provides baselines and benchmarks for exact algorithms for spin glasses and combinatorial optimization problems, and for evaluating heuristic algorithms and mean-field theories., Comment: code: https://github.com/TensorBFS/TropicalTensors.jl

Maximum entropy null models of networks come in different flavors that depend on the type of constraints under which entropy is maximized. If the constraints are on degree sequences or distributions, we are dealing with configuration models. If the degree sequence is constrained exactly, the corresponding microcanonical ensemble of random graphs with a given degree sequence is the configuration model per se. If the degree sequence is constrained only on average, the corresponding grand-canonical ensemble of random graphs with a given expected degree sequence is the soft configuration model. If the degree sequence is not fixed at all but randomly drawn from a fixed distribution, the corresponding hypercanonical ensemble of random graphs with a given degree distribution is the hypersoft configuration model, a more adequate description of dynamic real-world networks in which degree sequences are never fixed but degree distributions often stay stable. Here, we introduce the hypersoft configuration model of weighted networks. The main contribution is a particular version of the model with power-law degree and strength distributions, and superlinear scaling of strengths with degrees, mimicking the properties of some real-world networks. As a byproduct, we generalize the notions of sparse graphons and their entropy to weighted networks., Comment: 26 pages, 10 figures

In this paper we study the influence of additive noise in randomized consensus algorithms. Assuming that the update matrices are symmetric, we derive a closed form expression for the mean square error induced by the noise, together with upper and lower bounds that are simpler to evaluate. Motivated by the study of Open Multi-Agent Systems, we concentrate on Randomly Induced Discretized Laplacians, a family of update matrices that are generated by sampling subgraphs of a large undirected graph. For these matrices, we express the bounds by using the eigenvalues of the Laplacian matrix of the underlying graph or the graph's average effective resistance, thereby proving their tightness. Finally, we derive expressions for the bounds on some examples of graphs and numerically evaluate them., Comment: 6 pages, submitted to IEEE Control Systems Letters

An important challenge in the field of exponential random graphs (ERGs) is the fitting of non-trivial ERGs on large graphs. By utilizing fast matrix block-approximation techniques, we propose an approximative framework to such non-trivial ERGs that result in dyadic independence (i.e., edge independent) distributions, while being able to meaningfully model both local information of the graph (e.g., degrees) as well as global information (e.g., clustering coefficient, assortativity, etc.) if desired. This allows one to efficiently generate random networks with similar properties as an observed network, and the models can be used for several downstream tasks such as link prediction. Our methods are scalable to sparse graphs consisting of millions of nodes. Empirical evaluation demonstrates competitiveness in terms of both speed and accuracy with state-of-the-art methods -- which are typically based on embedding the graph into some low-dimensional space -- for link prediction, showcasing the potential of a more direct and interpretable probabalistic model for this task., Comment: Accepted for DSAA 2020 conference

Consider an infinite planar graph with uniform polynomial growth of degree d > 2. Many examples of such graphs exhibit similar geometric and spectral properties, and it has been conjectured that this is necessary. We present a family of counterexamples. In particular, we show that for every rational d > 2, there is a planar graph with uniform polynomial growth of degree d on which the random walk is transient, disproving a conjecture of Benjamini (2011). By a well-known theorem of Benjamini and Schramm, such a graph cannot be a unimodular random graph. We also give examples of unimodular random planar graphs of uniform polynomial growth with unexpected properties. For instance, graphs of (almost sure) uniform polynomial growth of every rational degree d > 2 for which the speed exponent of the walk is larger than 1/d, and in which the complements of all balls are connected. This resolves negatively two questions of Benjamini and Papasoglou (2011)., 26 pages, 6 figures

The eternal graph colouring problem, recently introduced by Klostermeyer and Mendoza, is a version of the graph colouring game, where two players take turns properly colouring a graph. In this note, we study the eternal game chromatic number of random graphs. We show that with high probability $\chi_{g}^{\infty}(G_{n,p}) = (\frac{p}{2} + o(1))n$ for odd $n$, and also for even $n$ when $p=\frac{1}{k}$ for some $k \in \mathbb{N}$. The upper bound applies for even $n$ and any other value of $p$ as well, but we conjecture in this case this upper bound is not sharp. Finally, we answer a question posed by Klostermeyer and Mendoza., Comment: 16 pages

We study the detection and the reconstruction of a large very dense subgraph in a social graph with $n$ nodes and $m$ edges given as a stream of edges, when the graph follows a power law degree distribution, in the regime when $m=O(n. \log n)$. A subgraph $S$ is very dense if it has $\Omega(|S|^2)$ edges. We uniformly sample the edges with a Reservoir of size $k=O(\sqrt{n}.\log n)$. Our detection algorithm checks whether the Reservoir has a giant component. We show that if the graph contains a very dense subgraph of size $\Omega(\sqrt{n})$, then the detection algorithm is almost surely correct. On the other hand, a random graph that follows a power law degree distribution almost surely has no large very dense subgraph, and the detection algorithm is almost surely correct. We define a new model of random graphs which follow a power law degree distribution and have large very dense subgraphs. We then show that on this class of random graphs we can reconstruct a good approximation of the very dense subgraph with high probability. We generalize these results to dynamic graphs defined by sliding windows in a stream of edges., Comment: 32 pages

We propose stochastic control policies to cope with uncertain and variable gas extractions in natural gas networks. Given historical gas extraction data, these policies are optimized to produce the real-time control inputs for nodal gas injections and for pressure regulation rates by compressors and valves. We describe the random network state as a function of control inputs, which enables a chance-constrained optimization of these policies for arbitrary network topologies. This optimization ensures the real-time gas flow feasibility and a minimal variation in the network state up to specified feasibility and variance criteria. Furthermore, the chance-constrained optimization provides the foundation of a stochastic pricing scheme for natural gas networks, which improves on a deterministic market settlement by offering the compensations to network assets for their contribution to uncertainty and variance control. We analyze the economic properties, including efficiency, revenue adequacy and cost recovery, of the proposed pricing scheme and make them conditioned on the network design., Comment: for associated GitHub repository, see https://github.com/anubhavratha/ng_stochastic_control_and_pricing

This paper provides the first empirical network analysis of the Argentine interbank money market. Its main topological features are examined applying graph theory, focusing on the unsecured overnight loans settled from 2003 to 2017. The network, where banks are the nodes and the operations between them represent the links, exhibits low density, a higher reciprocity than comparable random graphs, short average distances and its clustering coefficient remains above that of a random network of equal size. Furthermore, the network is prominently disassortative. Its structural metrics experienced significant volatility, in correlation with the economic activity fluctuations and regulatory shifts. Signs of nodes' random-like behavior are detected during contractions. The degree distributions fit better to a Lognormal distribution than to a Poisson or a Power Law. Additionally, different node centrality measures are computed. It is found that a higher centrality enables a node to settle more convenient bilateral interest rates compared to the average market rate, identifying a statistical and economically significant effect by means of a regression analysis. These results constitute a relevant input for systemic risk assessment and provide solid empirical foundations for future theoretical modelling and shock simulations, especially in the context of underdeveloped financial systems., Comment: 24 pages, 8 figures, 4 tables

Let L be subset of { 3 , 4 , … } and let X n , M ( L ) be the number of cycles belonging to unicyclic components whose length is in L in the random graph G ( n , M ) . We find the limiting distribution of X n , M ( L ) in the subcritical regime M = c n with c 1 / 2 and the critical regime M = n 2 ( 1 + μ n − 1 / 3 ) with μ = O ( 1 ) . Depending on the regime and a condition involving the series ∑ l ∈ L z l / ( 2 l ) , we obtain in the limit either a Poisson or a normal distribution as n → ∞ .

Tree-child networks are a recently-described class of directed acyclic graphs that have risen to prominence in phylogenetics (the study of evolutionary trees and networks). Although these networks have a number of attractive mathematical properties, many combinatorial questions concerning them remain intractable. In this paper, we show that endowing these networks with a biologically relevant ranking structure yields mathematically tractable objects, which we term ranked tree-child networks (RTCNs). We explain how to derive exact and explicit combinatorial results concerning the enumeration and generation of these networks. We also explore probabilistic questions concerning the properties of RTCNs when they are sampled uniformly at random. These questions include the lengths of random walks between the root and leaves (both from the root to the leaves and from a leaf to the root); the distribution of the number of cherries in the network; and sampling RTCNs conditional on displaying a given tree. We also formulate a conjecture regarding the scaling limit of the process that counts the number of lineages in the ancestry of a leaf. The main idea in this paper, namely using ranking as a way to achieve combinatorial tractability, may also extend to other classes of networks.

We study the problem of robust global synchronization of pulse-coupled oscillators (PCOs) over directed graphs. It is known that when the digraphs are strongly connected, global synchronization can be achieved by using a class of deterministic binary set-valued resetting controllers (Poveda and Teel, 2019). However, for large-scale networks, these algorithms are not scalable because some of their tuning parameters have upper bounds of the order O 1 N , where N is the number of agents. This paper resolves this scalability issue by presenting several new results about global synchronization of PCOs with more general network topologies using the frameworks of deterministic and stochastic hybrid dynamical systems. First, we establish that similar deterministic binary resetting algorithms can achieve robust global and fixed-time synchronization in any rooted acyclic digraph. Moreover, in this case, we show that the synchronization dynamics are now scalable as the tuning parameters of the algorithm are network independent, i.e., of order O ( 1 ) . However, the algorithms cannot be further extended to all rooted digraphs. We establish this new impossibility result by introducing a counter-example with a particular rooted digraph for which global synchronization cannot be achieved, irrespective of the tuning parameters. Nevertheless, we show that if the binary resetting algorithms are modified by accommodating an Erdos–Renýi type random graph model, then the resulting stochastic resetting dynamics will guarantee global synchronization almost surely for all rooted digraphs and, moreover, the tuning parameters of the dynamics are network independent. Stability and robustness properties of the resetting algorithms are studied using the tools from set-valued hybrid dynamical systems. Numerical simulations are provided at the end of the paper for demonstration of the main results.

We analyze the fluctuations of incomplete $U$-statistics over a triangular array of independent random variables. We give criteria for a Central Limit Theorem (CLT, for short) to hold in the sense that we prove that an appropriately scaled and centered version of the U-statistic converges to a normal random variable. Our method of proof relies on a martingale CLT. A possible application -- a CLT for the hitting time for random walk on random graphs -- will be presented in \cite{LoTe20b}

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\H{o}s-R\'{e}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: (1) we prove that the random intersection graph converges in total variation to $\mathcal{G}(n, p)$ when $d = \tilde{\omega}(n^3)$, and does not if $d = o(n^3)$, resolving an open problem in Fill et al. (2000), Rybarczyk (2011) and Kim et al. (2018); (2) we provide conditions under which the matrix of intersection sizes of random family of sets converges in total variation to a symmetric matrix with independent Poisson entries, yielding the first total variation convergence result for $\tau$-random intersection graphs to $\mathcal{G}(n, p)$; and (3) we show that the random geometric graph on $\mathbb{S}^{d - 1}$ with edge density $p$ converges in total variation to $\mathcal{G}(n, p)$ when $d = \tilde{\omega}\left(\min\{ pn^3, p^2 n^{7/2} \} \right)$, yielding the first progress towards a conjecture of Bubeck et al. (2016). The first of these three results was obtained simultaneously and independently by Bubeck, Racz and Richey., Comment: 62 pages

Motivated by the evidence that real-world networks evolve in time and may exhibit non-stationary features, we propose an extension of the Exponential Random Graph Models (ERGMs) accommodating the time variation of network parameters. Within the ERGM framework, a network realization is sampled from a static probability distribution defined parametrically in terms of network statistics. Inspired by the fast growing literature on Dynamic Conditional Score-driven models, in our approach, each parameter evolves according to an updating rule driven by the score of the conditional distribution. We demonstrate the flexibility of the score-driven ERGMs, both as data generating processes and as filters, and we prove the advantages of the dynamic version with respect to the static one. Our method captures dynamical network dependencies, that emerge from the data, and allows for a test discriminating between static or time-varying parameters. Finally, we corroborate our findings with the application to networks from real financial and political systems exhibiting non stationary dynamics.

We numerically study the Kuramoto model's synchronization consisting of the two groups of conformist-contrarian and excitatory-inhibitory phase oscillators with equal intrinsic frequency. We consider random and small-world (SW) topologies for the connectivity network of the oscillators. In random networks, regardless of the contrarian to conformist connection strength ratio, we found a crossover from the $\pi$-state to the blurred $\pi$-state and then a continuous transition to the incoherent state by increasing the fraction of contrarians. However, for the excitatory-inhibitory model in a random network, we found that for all the values of the fraction of inhibitors, the two groups remain in-phase, and the transition point of fully synchronized to incoherent state reduced by strengthening the ratio of inhibitory to excitatory links. In the SW networks, we found that the order parameters for both models have not monotonic behavior in terms of the fraction of contrarians and inhibitors. Up to the optimal fraction of contrarians inhibitors, the synchronization rises by introducing the contrarians and inhibitors and then falls. We discuss that this non-monotonic behavior in synchronization is due to the weakening of the defects formed already in conformists and excitatory agents in SW networks., Comment: 8 pages, f figures

A novel distributed algorithm for estimating the maximum of the node initial state values in a network, in the presence of additive communication noise is proposed. Conventionally, the maximum is estimated locally at each node by updating the node state value with the largest received measurements in every iteration. However, due to the additive channel noise, the estimate of the maximum at each node drifts at each iteration and this results in nodes diverging from the true max value. Max-plus algebra is used as a tool to study this ergodic process. The subadditive ergodic theorem is invoked to establish a constant growth rate for the state values due to noise, which is studied by analyzing the max-plus Lyapunov exponent of the product of noise matrices in a max-plus semiring. The growth rate of the state values is upper bounded by a constant which depends on the spectral radius of the network and the noise variance. Upper and lower bounds are derived for both fixed and random graphs. Finally, a two-run algorithm robust to additive noise in the network is proposed and its variance is analyzed using concentration inequalities. Simulation results supporting the theory are also presented., Comment: Accepted by IEEE Transactions on Signal and Information Processing over Networks

In a complete graph $K_n$ with edge weights drawn independently from a uniform distribution $U(0,1)$ (or alternatively an exponential distribution $\operatorname{Exp}(1)$), let $T_1$ be the MST (the spanning tree of minimum weight) and let $T_k$ be the MST after deletion of the edges of all previous trees $T_i$, $i

We analyze a dynamic random undirected graph in which newly added vertices are connected to those already present in the graph either using, with probability p , an anti-preferential attachment mechanism or, with probability 1 − p , a preferential attachment mechanism. We derive the asymptotic degree distribution in the general case and study the asymptotic behavior of the expected degree process in the general and that of the degree process in the pure anti-preferential attachment case. Degree homogenization mainly affects convergence rates for the former case and also the limiting degree distribution in the latter. Lastly, we perform a simulative study of a variation of the introduced model allowing for anti-preferential attachment probabilities given in terms of the current maximum degree of the graph.

Random graphs with a given degree sequence are often constructed using the configuration model, which yields a random multigraph. We may adjust this multigraph by a sequence of switchings, eventually yielding a simple graph. We show that, assuming essentially a bounded second moment of the degree distribution, this construction with the simplest types of switchings yields a simple random graph with an almost uniform distribution, in the sense that the total variation distance is $o(1)$. This construction can be used to transfer results on distributional convergence from the configuration model multigraph to the uniform random simple graph with the given vertex degrees. As examples, we give a few applications to asymptotic normality. We show also a weaker result yielding contiguity when the maximum degree is too large for the main theorem to hold., Comment: 29 pages. New references in v2

This paper examines bounds on upper tails for cycle counts in $G_{n,p}$. For a fixed graph $H$ define $\xi_H= \xi_H^{n,p}$ to be the number of copies of $H$ in $G_{n,p}$. It is a much studied and surprisingly difficult problem to understand the upper tail of the distribution of $\xi_H$, for example, to estimate \begin{equation*} \mathbb{P}(\xi_H > 2 \mathbb{E}\xi_H). \end{equation*} The best known result for general $H$ and $p$ is due to Janson, Oleszkiewicz, and Ruci\'nski, who, in 2004, proved \begin{align}\label{a:JOR} \exp[-O_{H, \eta}(M_H(n,p) \ln(1/p))]& (1+\eta)\mathbb{E} \xi_H)\\& \frac{\ln^{1/(l-2)}n}{n}$.

In this paper, we propose a new averaging model for modeling the competitive influence of $K$ candidates among $n$ voters in an election process. For such an influence propagation model, we address the question of how many seeded voters a candidate needs to place among undecided voters in order to win an election. We show that for a random network generated from the stochastic block model, there exists a percolation threshold for a candidate to win the election if the number of seeded voters placed by the candidate exceeds the threshold. By conducting extensive experiments, we show that our theoretical percolation thresholds are very close to those obtained from simulations for random networks and the errors are within $10\%$ for a real-world network., Comment: 11 pages, 9 figures, this article is the complete version (with proofs) of the IEEE Global Communications Conference 2019 review paper