Spectra of complex networks

Networks with Many Structural Scales: A Renormalization Group Perspective

Scale invariance profoundly influences the dynamics and structure of complex systems, spanning from critical phenomena to network architecture. Here, we propose a precise definition of scale-invariant networks by leveraging the concept of a constant entropy-loss rate across scales in a renormalization-group coarse-graining setting. This framework enables us to differentiate between scale-free and scale-invariant networks, revealing distinct characteristics within each class. Furthermore, we offer a comprehensive inventory of genuinely scale-invariant networks, both natural and artificially constructed, demonstrating, e.g., that the human connectome exhibits notable features of scale invariance. Our findings open new avenues for exploring the scale-invariant structural properties crucial in biological and sociotechnological systems.

Effects of clustering heterogeneity on the spectral density of sparse networks

We derive exact equations for the spectral density of sparse networks with an arbitrary distribution of the number of single edges and triangles per node. These equations enable a systematic investigation of the effects of clustering on the spectral properties of the network adjacency matrix. In the case of heterogeneous networks, we demonstrate that the spectral density becomes more symmetric as the fluctuations in the triangle-degree sequence increase. This phenomenon is explained by the small clustering coefficient of networks with a large variance of the triangle-degree distribution. In the homogeneous case of regular clustered networks, we find that both perturbative and nonperturbative approximations fail to predict the spectral density in the high-connectivity limit. This suggests that traditional large-degree approximations may be ineffective in studying the spectral properties of networks with more complex motifs. Our theoretical results are fully confirmed by numerical diagonalizations of finite adjacency matrices.

Analytical solution of susceptible-infected-recovered models on homogeneous networks

The ability to actually implement epidemic models is a crucial stake for public institutions, as they may be overtaken by the increasing complexity of current models and sometimes tend to revert to less elaborate models such as the susceptible-infected-recovered (SIR) model. In our work, we study a simple epidemic propagation model, called SIR-k, which is based on a homogeneous network of degree k, where each individual has the same number k of neighbors. This model represents a refined version of the basic SIR which assumes a completely homogeneous population. We show that nevertheless, analytical expressions, simpler and richer than the ones existing for the SIR model, can be derived for this SIR-k model. In particular, we obtain an exact implicit analytical solution for any k, from which quantities such as the epidemic threshold or the total number of agents infected during the epidemic can be obtained. We furthermore obtain simple exact explicit solutions for small ks, and in the large k limit we find a new formulation of the analytical solution of the basic SIR model, which comes with new insights.

Heterogeneous message passing for heterogeneous networks

Message passing (MP) is a computational technique used to find approximate solutions to a variety of problems defined on networks. MP approximations are generally accurate in locally treelike networks but require corrections to maintain their accuracy level in networks rich with short cycles. However, MP may already be computationally challenging on very large networks and additional costs incurred by correcting for cycles could be prohibitive. We show how the issue can be addressed. By allowing each node in the network to have its own level of approximation, one can focus on improving the accuracy of MP approaches in a targeted manner. We perform a systematic analysis of 109 real-world networks and show that our node-based MP approximation is able to increase both the accuracy and speed of traditional MP approaches. We find that, compared to conventional MP, a heterogeneous approach based on a simple heuristic is more accurate in 81% of tested networks, faster in 64% of cases, and both more accurate and faster in 49% of cases.

Random sparse generators of Markovian evolution and their spectral properties

The evolution of a complex multistate system is often interpreted as a continuous-time Markovian process. To model the relaxation dynamics of such systems, we introduce an ensemble of random sparse matrices which can be used as generators of Markovian evolution. The sparsity is controlled by a parameter φ, which is the number of nonzero elements per row and column in the generator matrix. Thus, a member of the ensemble is characterized by the Laplacian of a directed regular graph with D vertices (number of system states) and 2φD edges with randomly distributed weights. We study the effects of sparsity on the spectrum of the generator. Sparsity is shown to close the large spectral gap that is characteristic of nonsparse random generators. We show that the first moment of the eigenvalue distribution scales as ∼φ, while its variance is ∼sqrt[φ]. By using extreme value theory, we demonstrate how the shape of the spectral edges is determined by the tails of the corresponding weight distributions and clarify the behavior of the spectral gap as a function of D. Finally, we analyze complex spacing ratio statistics of ultrasparse generators, φ=const, and find that starting already at φ⩾2, spectra of the generators exhibit universal properties typical of Ginibre's orthogonal ensemble.

A 2D Lévy-flight model for the complex dynamics of real-life financial markets

We report on the emergence of scaling laws in the temporal evolution of the daily closing values of the S&P 500 index prices and its modeling based on the Lévy flights in two dimensions (2D). The efficacy of our proposed model is verified and validated by using the extreme value statistics in the random matrix theory. We find that the random evolution of each pair of stocks in a 2D price space is a scale-invariant complex trajectory whose tortuosity is governed by a 2/3 geometric law between the gyration radius R_g(t) and the total length ℓ(t) of the path, i.e., R_g(t)∼ℓ(t)^2/3. We construct a Wishart matrix containing all stocks up to a specific variable period and look at its spectral properties for over 30 years. In contrast to the standard random matrix theory, we find that the distribution of eigenvalues has a power-law tail with a decreasing exponent over time-a quantitative indicator of the temporal correlations. We find that the time evolution of the distance of 2D Lévy flights with index α=3/2 from origin generates the same empirical spectral properties. The statistics of the largest eigenvalues of the model and the observations are in perfect agreement.

Breakdown of random matrix universality in Markov models

Biological systems need to react to stimuli over a broad spectrum of timescales. If and how this ability can emerge without external fine-tuning is a puzzle. We consider here this problem in discrete Markovian systems, where we can leverage results from random matrix theory. Indeed, generic large transition matrices are governed by universal results, which predict the absence of long timescales unless fine-tuned. We consider an ensemble of transition matrices and motivate a temperature-like variable that controls the dynamic range of matrix elements, which we show plays a crucial role in the applicability of the large matrix limit: as the dynamic range increases, a phase transition occurs whereby the random matrix theory result is avoided, and long relaxation times ensue, in the entire "ordered" phase. We furthermore show that this phase transition is accompanied by a drop in the entropy rate and a peak in complexity, as measured by predictive information [Bialek, Nemenman, and Tishby Neural Comput. 13, 2409 (2001)NEUCEB0899-766710.1162/089976601753195969]. Extending the Markov model to a Hidden Markov model (HMM), we show that observable sequences inherit properties of the hidden sequences, allowing HMMs to be understood in terms of more accessible Markov models. We then apply our findings to fMRI data from 820 human subjects scanned at wakeful rest. We show that the data can be quantitatively understood in terms of the random model, and that brain activity lies close to the phase transition when engaged in unconstrained, task-free cognition-supporting the brain criticality hypothesis in this context.

Granger Causality among Graphs and Application to Functional Brain Connectivity in Autism Spectrum Disorder

Graphs/networks have become a powerful analytical approach for data modeling. Besides, with the advances in sensor technology, dynamic time-evolving data have become more common. In this context, one point of interest is a better understanding of the information flow within and between networks. Thus, we aim to infer Granger causality (G-causality) between networks' time series. In this case, the straightforward application of the well-established vector autoregressive model is not feasible. Consequently, we require a theoretical framework for modeling time-varying graphs. One possibility would be to consider a mathematical graph model with time-varying parameters (assumed to be random variables) that generates the network. Suppose we identify G-causality between the graph models' parameters. In that case, we could use it to define a G-causality between graphs. Here, we show that even if the model is unknown, the spectral radius is a reasonable estimate of some random graph model parameters. We illustrate our proposal's application to study the relationship between brain hemispheres of controls and children diagnosed with Autism Spectrum Disorder (ASD). We show that the G-causality intensity from the brain's right to the left hemisphere is different between ASD and controls.

Joint Embedding of Graphs

Feature extraction and dimension reduction for networks is critical in a wide variety of domains. Efficiently and accurately learning features for multiple graphs has important applications in statistical inference on graphs. We propose a method to jointly embed multiple undirected graphs. Given a set of graphs, the joint embedding method identifies a linear subspace spanned by rank one symmetric matrices and projects adjacency matrices of graphs into this subspace. The projection coefficients can be treated as features of the graphs, while the embedding components can represent vertex features. We also propose a random graph model for multiple graphs that generalizes other classical models for graphs. We show through theory and numerical experiments that under the model, the joint embedding method produces estimates of parameters with small errors. Via simulation experiments, we demonstrate that the joint embedding method produces features which lead to state of the art performance in classifying graphs. Applying the joint embedding method to human brain graphs, we find it extracts interpretable features with good prediction accuracy in different tasks.

Localization and Universality of Eigenvectors in Directed Random Graphs

Although the spectral properties of random graphs have been a long-standing focus of network theory, the properties of right eigenvectors of directed graphs have so far eluded an exact analytic treatment. We present a general theory for the statistics of the right eigenvector components in directed random graphs with a prescribed degree distribution and with randomly weighted links. We obtain exact analytic expressions for the inverse participation ratio and show that right eigenvectors of directed random graphs with a small average degree are localized. Remarkably, if the fourth moment of the degree distribution is finite, then the critical mean degree of the localization transition is independent of the degree fluctuations, which is different from localization in undirected graphs that is governed by degree fluctuations. We also show that in the high connectivity limit the distribution of the right eigenvector components is solely determined by the degree distribution. For delocalized eigenvectors, we recover in this limit the universal results from standard random matrix theory that are independent of the degree distribution, while for localized eigenvectors the eigenvector distribution depends on the degree distribution.

Characterization of Specific Ion Effects on PI(4,5)P2 Clustering: Molecular Dynamics Simulations and Graph-Theoretic Analysis

Numerous cellular functions mediated by phosphatidylinositol (4,5)-bisphosphate (PI(4,5)P2; PIP2) involve clustering of the lipid as well as colocalization with other lipids. Although the cation-mediated electrostatic interaction is regarded as the primary clustering mechanism, the ion-specific nature of the intermolecular network formation makes it challenging to characterize the clusters. Here we use all-atom molecular dynamics (MD) simulations of PIP2 monolayers and graph-theoretic analysis to gain insight into the phenomenon. MD simulations reveal that the intermolecular interactions preferentially occur between specific cations and phosphate groups (P1, P4, and P5) of the inositol headgroup with better-matched kosmotropic/chaotropic characters consistent with the law of matching water affinities (LMWA). Ca2+ is strongly attracted to P4/P5, while K+ preferentially binds to P1; Na+ interacts with both P4/P5 and P1. These specific interactions lead to the characteristic clustering patterns. Specificially, the size distributions and structures of PIP2 clusters generated by kosmotropic cations Ca2+ and Na+ are bimodal, with a combination of small and large clusters, while there is little clustering in the presence of only chaotropic K+; the largest clusters are obtained in systems with all three cations. The small-world network (a model with both local and long-range connections) best characterizes the clusters, followed by the random and the scale-free networks. More generally, the present results interpreted within the LMWA are consistent with the relative eukaryotic intracellular concentrations Ca2+ ≪ Na+ < Mg2+ < K+; that is, concentrations of Ca2+ and Na+ must be low to prevent damaging aggregation of lipids, DNA, RNA and phosphate-containing proteins.

The Maximum Eigenvalue of the Brain Functional Network Adjacency Matrix: Meaning and Application in Mental Fatigue Evaluation

The maximum eigenvalue of the adjacency matrix (AM) has been supposed to contain rich information about the corresponding network. An experimental study focused on revealing the meaning and application of the maximum eigenvalue is missing. To this end, AM was constructed using mutual information (MI) to determine the functional connectivity with electroencephalogram (EEG) data recorded with a mental fatigue model, and then was converted into both binary and weighted brain functional network (BFN) and corresponding random networks (RNs). Both maximum eigenvalue and corresponding network characters in BFNs and RNs were considered to explore the changes during the formation of mental fatigue. The results indicated that large maximum eigenvalue means more edges in the corresponding network, along with a high degree and a short characteristic path length both in weighted and binary BFNs. Interestingly, the maximum eigenvalue of AM was always a little larger than that of the corresponding random matrix (RM), and had an obvious linearity with the sum of the AM elements, indicating that the maximum eigenvalue can be able to distinguish the network structures which have the same mean degree. What is more, the maximum eigenvalue, which increased with the deepening of mental fatigue, can become a good indicator for mental fatigue estimation.

Eigenvalue repulsion and eigenvector localization in sparse non-Hermitian random matrices

Complex networks with directed, local interactions are ubiquitous in nature and often occur with probabilistic connections due to both intrinsic stochasticity and disordered environments. Sparse non-Hermitian random matrices arise naturally in this context and are key to describing statistical properties of the nonequilibrium dynamics that emerges from interactions within the network structure. Here we study one-dimensional (1D) spatial structures and focus on sparse non-Hermitian random matrices in the spirit of tight-binding models in solid state physics. We first investigate two-point eigenvalue correlations in the complex plane for sparse non-Hermitian random matrices using methods developed for the statistical mechanics of inhomogeneous two-dimensional interacting particles. We find that eigenvalue repulsion in the complex plane directly correlates with eigenvector delocalization. In addition, for 1D chains and rings with both disordered nearest-neighbor connections and self-interactions, the self-interaction disorder tends to decorrelate eigenvalues and localize eigenvectors more than simple hopping disorder. However, remarkable resistance to eigenvector localization by disorder is provided by large cycles, such as those embodied in 1D periodic boundary conditions under strong directional bias. The directional bias also spatially separates the left and right eigenvectors, leading to interesting dynamics in excitation and response. These phenomena have important implications for asymmetric random networks and highlight a need for mathematical tools to describe and understand them analytically.

Spectral analysis for weighted iterated q-triangulation networks

Deterministic weighted networks have been widely used to model real-world complex systems. In this paper, we study the weighted iterated q-triangulation networks, which are generated by iteration operation F(⋅). We add q(q∈N₊) new nodes on each old edge and connect them with two endpoints of the old edge. At the same time, the newly linked edges are given weight factor r(0<r≤1). From the construction of the network, we obtain all the eigenvalues and their multiplicities of its normalized Laplacian matrix from the two successive generations of the weighted iterated q-triangulation network. Further, as applications of spectra of the normalized Laplacian matrix, we study the Kemeny constant, the multiplicative degree-Kirchhoff index, and the number of weighted spanning trees and derive their exact closed-form expressions for weighted iterated q-triangulation networks.

Universal hypotrochoidic law for random matrices with cyclic correlations

The celebrated elliptic law describes the distribution of eigenvalues of random matrices with correlations between off-diagonal pairs of elements, having applications to a wide range of physical and biological systems. Here, we investigate the generalization of this law to random matrices exhibiting higher-order cyclic correlations between k tuples of matrix entries. We show that the eigenvalue spectrum in this ensemble is bounded by a hypotrochoid curve with k-fold rotational symmetry. This hypotrochoid law applies to full matrices as well as sparse ones, and thereby holds with remarkable universality. We further extend our analysis to matrices and graphs with competing cycle motifs, which are described more generally by polytrochoid spectral boundaries.

Spectra of networks containing short loops

The spectrum of the adjacency matrix plays several important roles in the mathematical theory of networks and network data analysis, for example in percolation theory, community detection, centrality measures, and the theory of dynamical systems on networks. A number of methods have been developed for the analytic computation of network spectra, but they typically assume that networks are locally treelike, meaning that the local neighborhood of any node takes the form of a tree, free of short loops. Empirically observed networks, by contrast, often have many short loops. Here we develop an approach for calculating the spectra of networks with short loops using a message passing method. We give example applications to some previously studied classes of networks and find that the presence of loops induces substantial qualitative changes in the shape of network spectra, generating asymmetries, multiple spectral bands, and other features.

Spectral properties of hyperbolic nanonetworks with tunable aggregation of simplexes

Cooperative self-assembly is a ubiquitous phenomenon found in natural systems which is used for designing nanostructured materials with new functional features. Its origin and mechanisms, leading to improved functionality of the assembly, have attracted much attention from researchers in many branches of science and engineering. These complex structures often come with hyperbolic geometry; however, the relation between the hyperbolicity and their spectral and dynamical properties remains unclear. Using the model of aggregation of simplexes introduced by Šuvakov et al. [Sci. Rep. 8, 1987 (2018)2045-232210.1038/s41598-018-20398-x], here we study topological and spectral properties of a large class of self-assembled structures or nanonetworks consisting of monodisperse building blocks (cliques of size n=3,4,5,6) which self-assemble via sharing the geometrical shapes of a lower order. The size of the shared substructure is tuned by varying the chemical affinity ν such that for significant positive ν sharing the largest face is the most probable, while for ν<0, attaching via a single node dominates. Our results reveal that, while the parameter of hyperbolicity remains δ_{max}=1 across the assemblies, their structure and spectral dimension d_{s} vary with the size of cliques n and the affinity when ν≥0. In this range, we find that d_{s}>4 can be reached for n≥5 and sufficiently large ν. For the aggregates of triangles and tetrahedra, the spectral dimension remains in the range d_{s}∈[2,4), as well as for the higher cliques at vanishing affinity. On the other end, for ν<0, we find d_{s}≂1.57 independently on n. Moreover, the spectral distribution of the normalized Laplacian eigenvalues has a characteristic shape with peaks and a pronounced minimum, representing the hierarchical architecture of the simplicial complexes. These findings show how the structures compatible with complex dynamical properties can be assembled by controlling the higher-order connectivity among the building blocks.

Perturbation analysis and comparison of network synchronization methods

In many networked systems, synchronization is important and useful, and how to enhance synchronizability is an interesting problem. Based on the matrix perturbation theory, we analyze five methods of network synchronization enhancement, including the link removal, node removal, dividing hub node, pull control, and pinning control methods, and obtain explicit expressions for eigenvalue changes. By these comparisons, we find that, among all these methods, the pull control method is remarkable, as it can extend the synchronization (coupling strength) region from both the left and right sides, for any controlled node. Extensive simulation results are given to support the accuracy of the perturbation-based analysis.

Spectra of random networks with arbitrary degrees

We derive a message-passing method for computing the spectra of locally treelike networks and an approximation to it that allows us to compute closed-form expressions or fast numerical approximates for the spectral density of random graphs with arbitrary node degrees-the so-called configuration model. We find that the latter approximation works well for all but the sparsest of networks. We also derive bounds on the position of the band edges of the spectrum, which are important for identifying structural phase transitions in networks.

Network spectra for drug-target identification in complex diseases: new guns against old foes

The fundamental understanding of altered complex molecular interactions in a diseased condition is the key to its cure. The overall functioning of these molecules is kind of jugglers play in the cell orchestra and to anticipate these relationships among the molecules is one of the greatest challenges in modern biology and medicine. Network science turned out to be providing a successful and simple platform to understand complex interactions among healthy and diseased tissues. Furthermore, much information about the structure and dynamics of a network is concealed in the eigenvalues of its adjacency matrix. In this review, we illustrate rapid advancements in the field of network science in combination with spectral graph theory that enables us to uncover the complexities of various diseases. Interpretations laid by network science approach have solicited insights into molecular relationships and have reported novel drug targets and biomarkers in various complex diseases.

Spectral properties of complex networks

This review presents an account of the major works done on spectra of adjacency matrices drawn on networks and the basic understanding attained so far. We have divided the review under three sections: (a) extremal eigenvalues, (b) bulk part of the spectrum, and (c) degenerate eigenvalues, based on the intrinsic properties of eigenvalues and the phenomena they capture. We have reviewed the works done for spectra of various popular model networks, such as the Erdős-Rényi random networks, scale-free networks, 1-d lattice, small-world networks, and various different real-world networks. Additionally, potential applications of spectral properties for natural processes have been reviewed.

Theory for the conditioned spectral density of noninvariant random matrices

We develop a theoretical approach to compute the conditioned spectral density of N×N noninvariant random matrices in the limit N→∞. This large deviation observable, defined as the eigenvalue distribution conditioned to have a fixed fraction k of eigenvalues smaller than x∈R, provides the spectrum of random matrix samples that deviate atypically from the average behavior. We apply our theory to sparse random matrices and unveil strikingly different and generic properties, namely, (i) their conditioned spectral density has compact support, (ii) it does not experience any abrupt transition for k around its typical value, and (iii) its eigenvalues do not accumulate at x. Moreover, our work points towards other types of transitions in the conditioned spectral density for values of k away from its typical value. These properties follow from the weak or absent eigenvalue repulsion in sparse ensembles and they are in sharp contrast to those displayed by classic or rotationally invariant random matrices. The exactness of our theoretical findings are confirmed through numerical diagonalization of finite random matrices.

Random matrix approach to plasmon resonances in the random impedance network model of disordered nanocomposites

Random impedance networks are widely used as a model to describe plasmon resonances in disordered metal-dielectric and other two-component nanocomposites. In the present work, the spectral properties of resonances in random networks are studied within the framework of the random matrix theory. We have shown that the appropriate ensemble of random matrices for the considered problem is the Jacobi ensemble (the MANOVA ensemble). The obtained analytical expressions for the density of states in such resonant networks show a good agreement with the results of numerical simulations in a wide range of metal filling fractions 0<p<1. A correspondence with the effective medium approximation is observed.

Emergent spectral properties of river network topology: an optimal channel network approach

Characterization of river drainage networks has been a subject of research for many years. However, most previous studies have been limited to quantities which are loosely connected to the topological properties of these networks. In this work, through a graph-theoretic formulation of drainage river networks, we investigate the eigenvalue spectra of their adjacency matrix. First, we introduce a graph theory model for river networks and explore the properties of the network through its adjacency matrix. Next, we show that the eigenvalue spectra of such complex networks follow distinct patterns and exhibit striking features including a spectral gap in which no eigenvalue exists as well as a finite number of zero eigenvalues. We show that such spectral features are closely related to the branching topology of the associated river networks. In this regard, we find an empirical relation for the spectral gap and nullity in terms of the energy dissipation exponent of the drainage networks. In addition, the eigenvalue distribution is found to follow a finite-width probability density function with certain skewness which is related to the drainage pattern. Our results are based on optimal channel network simulations and validated through examples obtained from physical experiments on landscape evolution. These results suggest the potential of the spectral graph techniques in characterizing and modeling river networks.

Degree heterogeneity and stability of ecological networks

A classic measure of ecological stability describes the tendency of a community to return to equilibrium after small perturbations. While many advances show how the network architecture of these communities severely constrains such tendencies, one of the most fundamental properties of network structure, i.e. degree heterogeneity-the variability of the number of links associated with each species, deserves further study. Here we show that the effects of degree heterogeneity on stability vary with different types of interspecific interactions. Degree heterogeneity consistently destabilizes ecological networks with both competitive and mutualistic interactions, while its effects on networks of predator-prey interactions such as food webs depend on prey contiguity, i.e. the extent to which the species consume an unbroken sequence of prey in community niche space. Increasing degree heterogeneity tends to stabilize food webs except those with the highest prey contiguity. These findings help explain why food webs are highly but not completely interval and, more broadly, deepen our understanding of the stability of complex ecological networks.

Localization in random bipartite graphs: Numerical and empirical study

We investigate adjacency matrices of bipartite graphs with a power-law degree distribution. Motivation for this study is twofold: first, vibrational states in granular matter and jammed sphere packings; second, graphs encoding social interaction, especially electronic commerce. We establish the position of the mobility edge and show that it strongly depends on the power in the degree distribution and on the ratio of the sizes of the two parts of the bipartite graph. At the jamming threshold, where the two parts have the same size, localization vanishes. We found that the multifractal spectrum is nontrivial in the delocalized phase, but still near the mobility edge. We also study an empirical bipartite graph, namely, the Amazon reviewer-item network. We found that in this specific graph the mobility edge disappears, and we draw a conclusion from this fact regarding earlier empirical studies of the Amazon network.

EIGENVECTOR-BASED CENTRALITY MEASURES FOR TEMPORAL NETWORKS. MULTISCALE MODELING & SIMULATION : A SIAM INTERDISCIPLINARY JOURNAL 2017;15:537-574. [PMID: 29046619 PMCID: PMC5643020 DOI: 10.1137/16m1066142]

Numerous centrality measures have been developed to quantify the importances of nodes in time-independent networks, and many of them can be expressed as the leading eigenvector of some matrix. With the increasing availability of network data that changes in time, it is important to extend such eigenvector-based centrality measures to time-dependent networks. In this paper, we introduce a principled generalization of network centrality measures that is valid for any eigenvector-based centrality. We consider a temporal network with N nodes as a sequence of T layers that describe the network during different time windows, and we couple centrality matrices for the layers into a supra-centrality matrix of size NT × NT whose dominant eigenvector gives the centrality of each node i at each time t. We refer to this eigenvector and its components as a joint centrality, as it reflects the importances of both the node i and the time layer t. We also introduce the concepts of marginal and conditional centralities, which facilitate the study of centrality trajectories over time. We find that the strength of coupling between layers is important for determining multiscale properties of centrality, such as localization phenomena and the time scale of centrality changes. In the strong-coupling regime, we derive expressions for time-averaged centralities, which are given by the zeroth-order terms of a singular perturbation expansion. We also study first-order terms to obtain first-order-mover scores, which concisely describe the magnitude of nodes' centrality changes over time. As examples, we apply our method to three empirical temporal networks: the United States Ph.D. exchange in mathematics, costarring relationships among top-billed actors during the Golden Age of Hollywood, and citations of decisions from the United States Supreme Court.

Eigenvalue tunneling and decay of quenched random network

We consider the canonical ensemble of N-vertex Erdős-Rényi (ER) random topological graphs with quenched vertex degree, and with fugacity μ for each closed triple of bonds. We claim complete defragmentation of large-N graphs into the collection of [p^{-1}] almost full subgraphs (cliques) above critical fugacity, μ_{c}, where p is the ER bond formation probability. Evolution of the spectral density, ρ(λ), of the adjacency matrix with increasing μ leads to the formation of a multizonal support for μ>μ_{c}. Eigenvalue tunneling from the central zone to the side one means formation of a new clique in the defragmentation process. The adjacency matrix of the network ground state has a block-diagonal form, where the number of vertices in blocks fluctuates around the mean value Np. The spectral density of the whole network in this regime has triangular shape. We interpret the phenomena from the viewpoint of the conventional random matrix model and speculate about possible physical applications.

Integrating cross-frequency and within band functional networks in resting-state MEG: A multi-layer network approach

Neuronal oscillations exist across a broad frequency spectrum, and are thought to provide a mechanism of interaction between spatially separated brain regions. Since ongoing mental activity necessitates the simultaneous formation of multiple networks, it seems likely that the brain employs interactions within multiple frequency bands, as well as cross-frequency coupling, to support such networks. Here, we propose a multi-layer network framework that elucidates this pan-spectral picture of network interactions. Our network consists of multiple layers (frequency-band specific networks) that influence each other via inter-layer (cross-frequency) coupling. Applying this model to MEG resting-state data and using envelope correlations as connectivity metric, we demonstrate strong dependency between within layer structure and inter-layer coupling, indicating that networks obtained in different frequency bands do not act as independent entities. More specifically, our results suggest that frequency band specific networks are characterised by a common structure seen across all layers, superimposed by layer specific connectivity, and inter-layer coupling is most strongly associated with this common mode. Finally, using a biophysical model, we demonstrate that there are two regimes of multi-layer network behaviour; one in which different layers are independent and a second in which they operate highly dependent. Results suggest that the healthy human brain operates at the transition point between these regimes, allowing for integration and segregation between layers. Overall, our observations show that a complete picture of global brain network connectivity requires integration of connectivity patterns across the full frequency spectrum.

Disentangling giant component and finite cluster contributions in sparse random matrix spectra

We describe a method for disentangling giant component and finite cluster contributions to sparse random matrix spectra, using sparse symmetric random matrices defined on Erdős-Rényi graphs as an example and test bed. Our methods apply to sparse matrices defined in terms of arbitrary graphs in the configuration model class, as long as they have finite mean degree.

Low-dimensional dynamics of structured random networks

Using a generalized random recurrent neural network model, and by extending our recently developed mean-field approach [J. Aljadeff, M. Stern, and T. Sharpee, Phys. Rev. Lett. 114, 088101 (2015)], we study the relationship between the network connectivity structure and its low-dimensional dynamics. Each connection in the network is a random number with mean 0 and variance that depends on pre- and postsynaptic neurons through a sufficiently smooth function g of their identities. We find that these networks undergo a phase transition from a silent to a chaotic state at a critical point we derive as a function of g. Above the critical point, although unit activation levels are chaotic, their autocorrelation functions are restricted to a low-dimensional subspace. This provides a direct link between the network's structure and some of its functional characteristics. We discuss example applications of the general results to neuroscience where we derive the support of the spectrum of connectivity matrices with heterogeneous and possibly correlated degree distributions, and to ecology where we study the stability of the cascade model for food web structure.

Correlation between weighted spectral distribution and average path length in evolving networks

The weighted spectral distribution (WSD) is a metric defined on the normalized Laplacian spectrum. In this study, synchronic random graphs are first used to rigorously analyze the metric's scaling feature, which indicates that the metric grows sublinearly as the network size increases, and the metric's scaling feature is demonstrated to be common in networks with Gaussian, exponential, and power-law degree distributions. Furthermore, a deterministic model of diachronic graphs is developed to illustrate the correlation between the slope coefficient of the metric's asymptotic line and the average path length, and the similarities and differences between synchronic and diachronic random graphs are investigated to better understand the correlation. Finally, numerical analysis is presented based on simulated and real-world data of evolving networks, which shows that the ratio of the WSD to the network size is a good indicator of the average path length.

Distinct types of eigenvector localization in networks

The spectral properties of the adjacency matrix provide a trove of information about the structure and function of complex networks. In particular, the largest eigenvalue and its associated principal eigenvector are crucial in the understanding of nodes’ centrality and the unfolding of dynamical processes. Here we show that two distinct types of localization of the principal eigenvector may occur in heterogeneous networks. For synthetic networks with degree distribution P(q) ~ q^−γ, localization occurs on the largest hub if γ > 5/2; for γ < 5/2 a new type of localization arises on a mesoscopic subgraph associated with the shell with the largest index in the K-core decomposition. Similar evidence for the existence of distinct localization modes is found in the analysis of real-world networks. Our results open a new perspective on dynamical processes on networks and on a recently proposed alternative measure of node centrality based on the non-backtracking matrix.

Neural Network Spectral Robustness under Perturbations of the Underlying Graph

Recent studies have been using graph-theoretical approaches to model complex networks (such as social, infrastructural, or biological networks) and how their hardwired circuitry relates to their dynamic evolution in time. Understanding how configuration reflects on the coupled behavior in a system of dynamic nodes can be of great importance, for example, in the context of how the brain connectome is affecting brain function. However, the effect of connectivity patterns on network dynamics is far from being fully understood. We study the connections between edge configuration and dynamics in a simple oriented network composed of two interconnected cliques (representative of brain feedback regulatory circuitry). In this article our main goal is to study the spectra of the graph adjacency and Laplacian matrices, with a focus on three aspects in particular: (1) the sensitivity and robustness of the spectrum in response to varying the intra- and intermodular edge density, (2) the effects on the spectrum of perturbing the edge configuration while keeping the densities fixed, and (3) the effects of increasing the network size. We study some tractable aspects analytically, then simulate more general results numerically, thus aiming to motivate and explain our further work on the effect of these patterns on the network temporal dynamics and phase transitions. We discuss the implications of such results to modeling brain connectomics. We suggest potential applications to understanding synaptic restructuring in learning networks and the effects of network configuration on function of regulatory neural circuits.

Spectra of weighted scale-free networks

Much information about the structure and dynamics of a network is encoded in the eigenvalues of its transition matrix. In this paper, we present a first study on the transition matrix of a family of weight driven networks, whose degree, strength, and edge weight obey power-law distributions, as observed in diverse real networks. We analytically obtain all the eigenvalues, as well as their multiplicities. We then apply the obtained eigenvalues to derive a closed-form expression for the random target access time for biased random walks occurring on the studied weighted networks. Moreover, using the connection between the eigenvalues of the transition matrix of a network and its weighted spanning trees, we validate the obtained eigenvalues and their multiplicities. We show that the power-law weight distribution has a strong effect on the behavior of random walks.

Complex network approach to fractional time series

In order to extract correlation information inherited in stochastic time series, the visibility graph algorithm has been recently proposed, by which a time series can be mapped onto a complex network. We demonstrate that the visibility algorithm is not an appropriate one to study the correlation aspects of a time series. We then employ the horizontal visibility algorithm, as a much simpler one, to map fractional processes onto complex networks. The degree distributions are shown to have parabolic exponential forms with Hurst dependent fitting parameter. Further, we take into account other topological properties such as maximum eigenvalue of the adjacency matrix and the degree assortativity, and show that such topological quantities can also be used to predict the Hurst exponent, with an exception for anti-persistent fractional Gaussian noises. To solve this problem, we take into account the Spearman correlation coefficient between nodes' degrees and their corresponding data values in the original time series.

Eigenvalues for the transition matrix of a small-world scale-free network: Explicit expressions and applications

The eigenvalues of the transition matrix for random walks on a network play a significant role in the structural and dynamical aspects of the network. Nevertheless, it is still not well understood how the eigenvalues behave in small-world and scale-free networks, which describe a large variety of real systems. In this paper, we study the eigenvalues for the transition matrix of a network that is simultaneously scale-free, small-world, and clustered. We derive explicit simple expressions for all eigenvalues and their multiplicities, with the spectral density exhibiting a power-law form. We then apply the obtained eigenvalues to determine the mixing time and random target access time for random walks, both of which exhibit unusual behaviors compared with those for other networks, signaling discernible effects of topologies on spectral features. Finally, we use the eigenvalues to count spanning trees in the network.

Ground-state phase-space structures of two-dimensional ±J spin glasses: A network approach

We illustrate a complex-network approach to study the phase spaces of spin glasses. By mapping the whole ground-state phase spaces of two-dimensional Edwards-Anderson bimodal (±J) spin glasses exactly into networks for analysis, we discovered various phase-space properties. The Gaussian connectivity distribution of the phase-space networks demonstrates that both the number of free spins and the visiting frequency of all microstates follow the Gaussian distribution. The spectra of phase-space networks are Gaussian, which is proven to be exact when the system is infinitely large. The phase-space networks exhibit community structures. By coarse graining to the community level, we constructed a network representing the entropy landscape of the ground state and discovered its scale-free property. The phase-space networks exhibit fractal structures, as a result of the rugged entropy landscape. Moreover, we show that the connectivity distribution, community structures, and fractal structures change drastically at the ferromagnetic-to-glass phase transition. These quantitative measurements of the ground states provide new insight into the study of spin glasses. The phase-space networks of spin glasses share a number of common features with those of lattice gases and geometrically frustrated spin systems and form a new class of complex networks with unique topology.

Linear processes in high dimensions: Phase space and critical properties

In this work we investigate the generic properties of a stochastic linear model in the regime of high dimensionality. We consider in particular the vector autoregressive (VAR) model and the multivariate Hawkes process. We analyze both deterministic and random versions of these models, showing the existence of a stable phase and an unstable phase. We find that along the transition region separating the two regimes the correlations of the process decay slowly, and we characterize the conditions under which these slow correlations are expected to become power laws. We check our findings with numerical simulations showing remarkable agreement with our predictions. We finally argue that real systems with a strong degree of self-interaction are naturally characterized by this type of slow relaxation of the correlations.

Origin and implications of zero degeneracy in networks spectra

The spectra of many real world networks exhibit properties which are different from those of random networks generated using various models. One such property is the existence of a very high degeneracy at the zero eigenvalue. In this work, we provide all the possible reasons behind the occurrence of the zero degeneracy in the network spectra, namely, the complete and partial duplications, as well as their implications. The power-law degree sequence and the preferential attachment are the properties which enhances the occurrence of such duplications and hence leading to the zero degeneracy. A comparison of the zero degeneracy in protein-protein interaction networks of six different species and in their corresponding model networks indicates importance of the degree sequences and the power-law exponent for the occurrence of zero degeneracy.

Universality in the spectral and eigenfunction properties of random networks

By the use of extensive numerical simulations, we show that the nearest-neighbor energy-level spacing distribution P(s) and the entropic eigenfunction localization length of the adjacency matrices of Erdős-Rényi (ER) fully random networks are universal for fixed average degree ξ≡αN (α and N being the average network connectivity and the network size, respectively). We also demonstrate that the Brody distribution characterizes well P(s) in the transition from α=0, when the vertices in the network are isolated, to α=1, when the network is fully connected. Moreover, we explore the validity of our findings when relaxing the randomness of our network model and show that, in contrast to standard ER networks, ER networks with diagonal disorder also show universality. Finally, we also discuss the spectral and eigenfunction properties of small-world networks.

Assortative and disassortative mixing investigated using the spectra of graphs

We investigate the impact of degree-degree correlations on the spectra of networks. Even though density distributions exhibit drastic changes depending on the (dis)assortative mixing and the network architecture, the short-range correlations in eigenvalues exhibit universal random matrix theory predictions. The long-range correlations turn out to be a measure of randomness in (dis)assortative networks. The analysis further provides insight into the origin of high degeneracy at the zero eigenvalue displayed by a majority of biological networks.

Localization transition, Lifschitz tails, and rare-region effects in network models

Effects of heterogeneity in the suspected-infected-susceptible model on networks are investigated using quenched mean-field theory. The emergence of localization is described by the distributions of the inverse participation ratio and compared with the rare-region effects appearing in simulations and in the Lifschitz tails. The latter, in the linear approximation, is related to the spectral density of the Laplacian matrix and to the time dependent order parameter. I show that these approximations indicate correctly Griffiths phases both on regular one-dimensional lattices and on small-world networks exhibiting purely topological disorder. I discuss the localization transition that occurs on scale-free networks at γ=3 degree exponent.

Spectra of random graphs with community structure and arbitrary degrees

Using methods from random matrix theory researchers have recently calculated the full spectra of random networks with arbitrary degrees and with community structure. Both reveal interesting spectral features, including deviations from the Wigner semicircle distribution and phase transitions in the spectra of community structured networks. In this paper we generalize both calculations, giving a prescription for calculating the spectrum of a network with both community structure and an arbitrary degree distribution. In general the spectrum has two parts, a continuous spectral band, which can depart strongly from the classic semicircle form, and a set of outlying eigenvalues that indicate the presence of communities.

Uncovering randomness and success in society

An understanding of how individuals shape and impact the evolution of society is vastly limited due to the unavailability of large-scale reliable datasets that can simultaneously capture information regarding individual movements and social interactions. We believe that the popular Indian film industry, “Bollywood”, can provide a social network apt for such a study. Bollywood provides massive amounts of real, unbiased data that spans more than 100 years, and hence this network has been used as a model for the present paper. The nodes which maintain a moderate degree or widely cooperate with the other nodes of the network tend to be more fit (measured as the success of the node in the industry) in comparison to the other nodes. The analysis carried forth in the current work, using a conjoined framework of complex network theory and random matrix theory, aims to quantify the elements that determine the fitness of an individual node and the factors that contribute to the robustness of a network. The authors of this paper believe that the method of study used in the current paper can be extended to study various other industries and organizations.