Delocalization transition for the Google matrix

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.

Multifractality in random networks with power-law decaying bond strengths

In this paper we demonstrate numerically that random networks whose adjacency matrices A are represented by a diluted version of the power-law banded random matrix (PBRM) model have multifractal eigenfunctions. The PBRM model describes one-dimensional samples with random long-range bonds. The bond strengths of the model, which decay as a power-law, are tuned by the parameter μ as A_{mn}∝|m-n|^{-μ}; while the sparsity is driven by the average network connectivity α: for α=0 the vertices in the network are isolated and for α=1 the network is fully connected and the PBRM model is recovered. Though it is known that the PBRM model has multifractal eigenfunctions at the critical value μ=μ_{c}=1, we clearly show [from the scaling of the relative fluctuation of the participation number I_{2} as well as the scaling of the probability distribution functions P(lnI_{2})] the existence of the critical value μ_{c}≡μ_{c}(α) for α<1. Moreover, we characterize the multifractality of the eigenfunctions of our random network model by the use of the corresponding multifractal dimensions D_{q}, that we compute from the finite network-size scaling of the typical eigenfunction participation numbers exp〈lnI_{q}〉.

Information Entropy of Tight-Binding Random Networks with Losses and Gain: Scaling and Universality

We study the localization properties of the eigenvectors, characterized by their information entropy, of tight-binding random networks with balanced losses and gain. The random network model, which is based on Erdős–Rényi (ER) graphs, is defined by three parameters: the network size N, the network connectivity α, and the losses-and-gain strength γ. Here, N and α are the standard parameters of ER graphs, while we introduce losses and gain by including complex self-loops on all vertices with the imaginary amplitude iγ with random balanced signs, thus breaking the Hermiticity of the corresponding adjacency matrices and inducing complex spectra. By the use of extensive numerical simulations, we define a scaling parameter ξ≡ξ(N,α,γ) that fixes the localization properties of the eigenvectors of our random network model; such that, when ξ<0.1 (10<ξ), the eigenvectors are localized (extended), while the localization-to-delocalization transition occurs for 0.1<ξ<10. Moreover, to extend the applicability of our findings, we demonstrate that for fixed ξ, the spectral properties (characterized by the position of the eigenvalues on the complex plane) of our network model are also universal; i.e., they do not depend on the specific values of the network parameters.

How the site degree influences quantum probability on inhomogeneous substrates

We investigate the effect of the node degree and energy E on the electronic wave function for regular and irregular structures, namely, regular lattices, disordered percolation clusters, and complex networks. We evaluate the dependency of the quantum probability for each site on its degree. For a class of biregular structures formed by two disjoint subsets of sites sharing the same degree, the probability P_{k}(E) of finding the electron on any site with k neighbors is independent of E≠0, a consequence of an exact analytical result that we prove for any bipartite lattice. For more general nonbipartite structures, P_{k}(E) may depend on E as illustrated by an exact evaluation of a one-dimensional semiregular lattice: P_{k}(E) is large for small values of E when k is also small, and its maximum values shift towards large values of |E| with increasing k. Numerical evaluations of P_{k}(E) for two different types of percolation clusters and the Apollonian network suggest that this observed feature might be generally valid.

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.

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.

Google matrix of the citation network of Physical Review

We study the statistical properties of spectrum and eigenstates of the Google matrix of the citation network of Physical Review for the period 1893-2009. The main fraction of complex eigenvalues with largest modulus is determined numerically by different methods based on high-precision computations with up to p = 16384 binary digits that allow us to resolve hard numerical problems for small eigenvalues. The nearly nilpotent matrix structure allows us to obtain a semianalytical computation of eigenvalues. We find that the spectrum is characterized by the fractal Weyl law with a fractal dimension d(f) ≈ 1. It is found that the majority of eigenvectors are located in a localized phase. The statistical distribution of articles in the PageRank-CheiRank plane is established providing a better understanding of information flows on the network. The concept of ImpactRank is proposed to determine an influence domain of a given article. We also discuss the properties of random matrix models of Perron-Frobenius operators.

Scattering and transport properties of tight-binding random networks

We study numerically scattering and transport statistical properties of tight-binding random networks characterized by the number of nodes N and the average connectivity α. We use a scattering approach to electronic transport and concentrate on the case of a small number of single-channel attached leads. We observe a smooth crossover from insulating to metallic behavior in the average scattering matrix elements <|S(mn)|(2)>, the conductance probability distribution w(T), the average conductance <T>, the shot noise power P, and the elastic enhancement factor F by varying α from small (α→0) to large (α→1) values. We also show that all these quantities are invariant for fixed ξ=αN. Moreover, we proposes a heuristic and universal relation between <|S(mn)|(2)>, <T>, and P and the disorder parameter ξ.

Google matrix analysis of DNA sequences

For DNA sequences of various species we construct the Google matrix [Formula: see text] of Markov transitions between nearby words composed of several letters. The statistical distribution of matrix elements of this matrix is shown to be described by a power law with the exponent being close to those of outgoing links in such scale-free networks as the World Wide Web (WWW). At the same time the sum of ingoing matrix elements is characterized by the exponent being significantly larger than those typical for WWW networks. This results in a slow algebraic decay of the PageRank probability determined by the distribution of ingoing elements. The spectrum of [Formula: see text] is characterized by a large gap leading to a rapid relaxation process on the DNA sequence networks. We introduce the PageRank proximity correlator between different species which determines their statistical similarity from the view point of Markov chains. The properties of other eigenstates of the Google matrix are also discussed. Our results establish scale-free features of DNA sequence networks showing their similarities and distinctions with the WWW and linguistic networks.

Adiabatic quantum algorithm for search engine ranking

We propose an adiabatic quantum algorithm for generating a quantum pure state encoding of the PageRank vector, the most widely used tool in ranking the relative importance of internet pages. We present extensive numerical simulations which provide evidence that this algorithm can prepare the quantum PageRank state in a time which, on average, scales polylogarithmically in the number of web pages. We argue that the main topological feature of the underlying web graph allowing for such a scaling is the out-degree distribution. The top-ranked log(n) entries of the quantum PageRank state can then be estimated with a polynomial quantum speed-up. Moreover, the quantum PageRank state can be used in "q-sampling" protocols for testing properties of distributions, which require exponentially fewer measurements than all classical schemes designed for the same task. This can be used to decide whether to run a classical update of the PageRank.

Spectral properties of directed random networks with modular structure

We study spectra of directed networks with inhibitory and excitatory couplings. We investigate in particular eigenvector localization properties of various model networks for different values of correlation among their entries. Spectra of random networks with completely uncorrelated entries show a circular distribution with delocalized eigenvectors, whereas networks with correlated entries have localized eigenvectors. In order to understand the origin of localization we track the spectra as a function of connection probability and directionality. As connections are made directed, eigenstates start occurring in complex-conjugate pairs and the eigenvalue distribution combined with the localization measure shows a rich pattern. Moreover, for a very well distinguished community structure, the whole spectrum is localized except few eigenstates at the boundary of the circular distribution. As the network deviates from the community structure there is a sudden change in the localization property for a very small value of deformation from the perfect community structure. We search for this effect for the whole range of correlation strengths and for different community configurations. Furthermore, we investigate spectral properties of a metabolic network of zebrafish and compare them with those of the model networks.

Spectral properties of the Google matrix of the World Wide Web and other directed networks

We study numerically the spectrum and eigenstate properties of the Google matrix of various examples of directed networks such as vocabulary networks of dictionaries and university World Wide Web networks. The spectra have gapless structure in the vicinity of the maximal eigenvalue for Google damping parameter α equal to unity. The vocabulary networks have relatively homogeneous spectral density, while university networks have pronounced spectral structures which change from one university to another, reflecting specific properties of the networks. We also determine specific properties of eigenstates of the Google matrix, including the PageRank. The fidelity of the PageRank is proposed as a characterization of its stability.

Google matrix and Ulam networks of intermittency maps

We study the properties of the Google matrix of an Ulam network generated by intermittency maps. This network is created by the Ulam method which gives a matrix approximant for the Perron-Frobenius operator of dynamical map. The spectral properties of eigenvalues and eigenvectors of this matrix are analyzed. We show that the PageRank of the system is characterized by a power law decay with the exponent beta dependent on map parameters and the Google damping factor alpha . Under certain conditions the PageRank is completely delocalized so that the Google search in such a situation becomes inefficient.

Google matrix, dynamical attractors, and Ulam networks

We study the properties of the Google matrix generated by a coarse-grained Perron-Frobenius operator of the Chirikov typical map with dissipation. The finite-size matrix approximant of this operator is constructed by the Ulam method. This method applied to the simple dynamical model generates directed Ulam networks with approximate scale-free scaling and characteristics being in certain features similar to those of the world wide web with approximate scale-free degree distributions as well as two characteristics similar to the web: a power-law decay in PageRank that mirrors the decay of PageRank on the world wide web and a sensitivity to the value alpha in PageRank. The simple dynamical attractors play here the role of popular websites with a strong concentration of PageRank. A variation in the Google parameter alpha or other parameters of the dynamical map can drive the PageRank of the Google matrix to a delocalized phase with a strange attractor where the Google search becomes inefficient.