Google matrix

Construction of Google matrix G

Then the final Google matrix G can be expressed via S as:

${\displaystyle G_{ij}=\alpha S_{ij}+(1-\alpha )\frac{1}{N}(1)}$

By the construction the sum of all non-negative elements inside each matrix column is equal to unity. The numerical coefficient ${\displaystyle \alpha }$ is known as a damping factor.

Usually S is a sparse matrix and for modern directed networks it has only about ten nonzero elements in a line or column, thus only about 10N multiplications are needed to multiply a vector by matrix G.

Examples of Google matrix

An example of the matrix ${\displaystyle S}$ construction via Eq.(1) within a simple network is given in the article CheiRank.

For the actual matrix, Google uses a damping factor ${\displaystyle \alpha }$ around 0.85. The term ${\displaystyle (1-\alpha )}$ gives a surfer probability to jump randomly on any page. The matrix ${\displaystyle G}$ belongs to the class of Perron-Frobenius operators of Markov chains. The examples of Google matrix structure are shown in Fig.1 for Wikipedia articles hyperlink network in 2009 at small scale and in Fig.2 for University of Cambridge network in 2006 at large scale.

Spectrum and eigenstates of G matrix

For ${\displaystyle 0<\alpha <1}$ there is only one maximal eigenvalue ${\displaystyle \lambda =1}$ with the corresponding right eigenvector which has non-negative elements ${\displaystyle P_i}$ which can be viewed as stationary probability distribution. These probabilities ordered by their decreasing values give the PageRank vector ${\displaystyle P_i}$ with the PageRank ${\displaystyle K_i}$ used by Google search to rank webpages. Usually one has for the World Wide Web that ${\displaystyle P\propto 1/K^{\beta }}$ with ${\displaystyle \beta \approx 0.9}$. The number of nodes with a given PageRank value scales as ${\displaystyle N_P\propto 1/P^{\nu }}$ with the exponent ${\displaystyle \nu =1+1/\beta \approx 2.1}$. The left eigenvector at ${\displaystyle \lambda =1}$ has constant matrix elements. With ${\displaystyle 0<\alpha }$ all eigenvalues move as ${\displaystyle {\lambda }_i\to \alpha {\lambda }_i}$ except the maximal eigenvalue ${\displaystyle \lambda =1}$, which remains unchanged. The PageRank vector varies with ${\displaystyle \alpha }$ but other eigenvectors with ${\displaystyle {\lambda }_i<1}$ remain unchanged due to their orthogonality to the constant left vector at ${\displaystyle \lambda =1}$. The gap between ${\displaystyle \lambda =1}$ and other eigenvalue being ${\displaystyle 1-\alpha \approx 0.15}$ gives a rapid convergence of a random initial vector to the PageRank approximately after 50 multiplications on ${\displaystyle G}$ matrix.

At ${\displaystyle \alpha =1}$ the matrix ${\displaystyle G}$ has generally many degenerate eigenvalues ${\displaystyle \lambda =1}$ (see e.g. [6]). Examples of the eigenvalue spectrum of the Google matrix of various directed networks is shown in Fig.3 from and Fig.4 from.

The Google matrix can be also constructed for the Ulam networks generated by the Ulam method [8] for dynamical maps. The spectral properties of such matrices are discussed in [9,10,11,12,13,15]. In a number of cases the spectrum is described by the fractal Weyl law [10,12].

The Google matrix can be constructed also for other directed networks, e.g. for the procedure call network of the Linux Kernel software introduced in [15]. In this case the spectrum of ${\displaystyle \lambda }$ is described by the fractal Weyl law with the fractal dimension ${\displaystyle d\approx 1.3}$ (see Fig.5 from ). Numerical analysis shows that the eigenstates of matrix ${\displaystyle G}$ are localized (see Fig.6 from ). Arnoldi iteration method allows to compute many eigenvalues and eigenvectors for matrices of rather large size [13].

Other examples of ${\displaystyle G}$ matrix include the Google matrix of brain [17] and business process management [18], see also. Applications of Google matrix analysis to DNA sequences is described in [20]. Such a Google matrix approach allows also to analyze entanglement of cultures via ranking of multilingual Wikipedia articles abouts persons [21]

Historical notes

The Google matrix with damping factor was described by Sergey Brin and Larry Page in 1998 [22], see also articles on PageRank history [23], [24].