Quantum Google in a Complex Network
We investigate the behavior of the recently proposed quantum Google algorithm, or quantum PageRank, in large complex networks. Applying the quantum algorithm to a part of the real World Wide Web, we find that the algorithm is able to univocally reveal the underlying scale-free topology of the network and to clearly identify and order the most relevant nodes (hubs) of the graph according to their importance in the network structure. Moreover, our results show that the quantum PageRank algorithm generically leads to changes in the hierarchy of nodes. In addition, as compared to its classical counterpart, the quantum algorithm is capable to clearly highlight the structure of secondary hubs of the network, and to partially resolve the degeneracy in importance of the low lying part of the list of rankings, which represents a typical shortcoming of the classical PageRank algorithm. Complementary to this study, our analysis shows that the algorithm is able to clearly distinguish scale-free networks from other widespread and important classes of complex networks, such as Erdős-Rényi networks and hierarchical graphs. We show that the ranking capabilities of the quantum PageRank algorithm are related to an increased stability with respect to a variation of the damping parameter that appears in the Google algorithm, and to a more clearly pronounced power-law behavior in the distribution of importance among the nodes, as compared to the classical algorithm. Finally, we study to which extent the increased sensitivity of the quantum algorithm persists under coordinated attacks of the most important nodes in scale-free and Erdős-Rényi random graphs.
pacs:03.67.Ac, 03.67.Hk, 89.20.Hh, 05.40.Fb
It is of great interest to explore and classify the large amount of information that is stored in huge complex networks like the World Wide Web (WWW). A central problem of bringing order to classical information stored in networks such as the WWW amounts to rank nodes containing such information according to their relevance. A highly successful and nowadays widespread tool for this purpose has been the PageRank algorithm Brin and Page (1998); Page et al. (1999), which lies at the core of Google’s ranking engine. In the foreseeable future where large-scale quantum networks have become a reality, classifying the quantum information stored therein will become a priority. It is in this sense that the recently introduced quantum PageRank algorithm Paparo and Martin-Delgado (2012) is an important achievement as it constitutes a quantization of the classical PageRank protocol. This new quantum algorithm has shown, applied to small networks, a striking behavior with respect to its classical counterpart, such as producing a different hierarchy of nodes together, paired with a better performance. In this paper we investigate the properties of the quantum algorithm for networks which model large real-world complex systems. We also test the algorithm on real-world data stemming from a part of the WWW.
Complex networks are more and more pervasive and essential in our everyday’s life, and thus the importance of network science. Consequently, considerable research is devoted to analyze and understand networks like the World Wide Web, the Internet, networks associated to transportation and communication systems and even biological and social networks. Starting with the seminal papers by Watts and Strogatz on small-world networks Watts and Strogatz (1998) and by Barabási and Albert on scale-free networks Barabási and Albert (1999), researchers realized that most relevant networks belong to a class known as small-world scale-free networks, as they exhibit both strong local clustering (nodes have many mutual neighbors) and a small average path length while sharing another important characteristic: the number of links of nodes usually obeys a power-law distribution (the network is scale-free). Moreover, it has been found that many real networks, including the WWW, are also self-similar, see Song et al. (2005). Such properties can often be related to a modular and hierarchical structure and organization which is essential for their communication and dynamical processes Ravasz and Barabási (2003); Barabási and Oltvai (2004); Barabási et al. (2003). On the other hand, this hierarchical structure could explain the existence of nodes with a relatively large number of links (or hubs), which play a critical role in the information flow of the system. Hubs are also associated with a low average distance in the network. Several general reviews and books on complex networks are now available, to which we refer the reader interested in more information on this topic Newman (2003, 2010); Cohen and Havlin (2010).
Networks considered in this paper are modeled by three classes of graphs: The first type are Erdős-R\a’enyi random graphs Erdős and Rényi (1959) . These graphs are constructed by connecting a given set of nodes with directed edges, each one added according to a certain fixed probability. The second class of graphs are scale-free graphs which were introduced by Barabási and Albert to model the WWW Barabási and Albert (1999); Barabási (2013). A graph is dynamically formed by a continuous addition of new vertices which are connected preferentially to vertices which already have a large degree. We consider here a version for directed graphs published in Bollobás et al. (2003). The third family are hierarchical graphs which are also scale-free but their clustering and degree distributions are negatively correlated as hubs have a smaller clustering coefficient than nodes with a lower degree. It has been found that hierarchical graphs constitute also a good model for the WWW Ravasz and Barabási (2003); Barabási and Oltvai (2004).
In the next sections we apply the quantum Google algorithm proposed in ref. Paparo and Martin-Delgado (2012) to representatives of the three classes of graphs. Our focus lies on directed scale-free networks, and also on hierarchical graphs as they are good models for the WWW, but we in addition also consider Erdős -Rényi random networks as a reference in order to contrast the results which we find for the former models.
In particular, two main fundamental questions will be addressed:
1/ Does the quantum PageRank algorithm preserve the structure of a scale-free network? In other words, we ask whether the ranking distribution obtained by the quantum algorithm also follows the same pattern of node-importance as the underlying scale-free network, and does not get mixed up with the distribution corresponding to a random Erdős-R\a’enyi network. Similarly, for a Erdős-R\a’enyi network we study whether the behavior of the quantum Page Rank algorithm is intrinsic to these networks and to which extent it differs from the one for scale-free networks.
2/ Can one improve on the information gained on some of the properties of a scale-free network as far as the quantum Google algorithm is concerned? Specifically, scale-free networks are robust against random disturbances which are unavoidable in a noisy environment. However, scale-free networks are vulnerable to coordinated attacks.
These questions are outlined in more detail throughout the paper and clear answers are given. We hereby summarize briefly some of our main results:
i/ We find that the quantum PageRank algorithm detects each of the representatives classes of complex networks, in the sense that the classification results show that the ranking of the nodes follow a similar law than the degrees of the network. In particular, for a scale-free network this distribution remains scale-free under the application of the quantum PageRank algorithm. This implies that the quantum PageRank algorithm is expected to be robust with respect to random external noise.
ii/ Remarkably enough, the detection of the hubs for a network in the class of scale-free networks with the QG algorithm is clearly enhanced with respect to the classical PageRank algorithm. In particular, the quantum PageRank algorithm is able to more clearly reveal the existence of secondary hubs in scale-free networks.
iii/ We find that the quantum PageRank algorithm and the resulting rankings are more stable than the classical PageRank protocol with respect to the variation of the damping parameter that appears in the Google algorithm. The low dependence on this parameter, which inevitably had to be arbitrarily tuned in the classical algorithm, provides the quantum algorithm with a higher objectivity of the importance rankings.
iv/ Our study shows that the quantum PageRank algorithm displays for scale-free graphs a power law scaling behavior of the importances of the nodes. Furthermore, this power law behavior is more favorable than the one the classical algorithm exhibits: Indeed, a smoother behavior is connected to a more harmonious distribution of importance among the nodes, which enables the algorithm to better uncover the structure of hubs in the underlying scale-free graphs.
v/ The enhanced sensitivity of the quantum PageRank algorithm to structural details of the networks is related to an increased sensitivity of the quantum algorithm under coordinated attacks of the most important nodes in scale-free networks.
i.1 Conceptual Setting of the Quantum PageRank Algorithm
Motivated by the fact that in a near-future scenario a certain class of quantum network will be operative Elliott (2006); Poppe et al. (2008); Sasaki et al. (2011); Swi (2009); Lancho et al. (2010); Länger and Lenhart (2009), but not yet a scalable quantum computer, in Paparo and Martin-Delgado (2012) a class of quantum algorithms to rank the nodes in a quantum network was put forward. The algorithms in this class must be compatible with the classical one. Indeed, existing projects for large-scale quantum networks contemplate using the backbone of existing communication networks, upgrading them to include the quantum hardware in order to store and manipulate quantum information. In particular, the directed structure of the graph must be preserved, a feature which is crucial to measure a node’s authority. An instance of this class of algorithms was explicitly constructed and analyzed for graphs of small size and it was found to serve as a valid quantum counterpart of Google’s PageRank algorithm (for details see ref. Paparo and Martin-Delgado (2012)).
To perform such a quantum task, it is important that the ranking algorithm must in principle incorporate some of the quantum properties of the network, like quantum fluctuations, and be objective. Indeed, the latter property earned much of Google’s PageRank’s success and was achieved embedding in the algorithm the random walk of a surfer exploring the WWW based on simple sensible rules. Within the same line of reasoning, we set up a quantum walk based algorithm that mimics the exploration of nodes in a quantum network, in the setting when these are represented by states of a Hilbert space. The simple rules are encoded in the quantum dynamics, and in doing so the quantum nature of the networks and the information stored is properly taken into account. In the setting where a fully fledged large-scale quantum computer is not yet available, a key property of the quantum algorithm is that it be efficiently simulatable on a classical computer, that is, it must belong to the computational complexity class P. Our algorithm is based on a single particle quantum walk and thus efficiently simulatable.
Furthermore, the fact that the quantum algorithm contains a quantum walk at its heart paves the way to the analysis of its dynamics from a purely physical perspective: In our work, we analyze extensively the localization properties of the quantum walk contained in the algorithm applied to several classes of networks (see sect. IV) showing its effects on ranking. To do so, we introduce the new concept of a quantum Inverse Participation Ratio (IPR). This is a generalization of the classical IPR, which is the main quantity that has been used extensively to probe localization properties in the study of classical random walks.
Besides its application to future quantum networks, the Quantum PageRank algorithm can indeed also be regarded as a valuable quantum tool, that can be efficiently run on a a classical computer, to perform a ”classical” task, namely to rank nodes in a classical networks. Indeed, embedding nontrivially the network connectivity structure in the quantum dynamics, our protocol turns out to show several features that improve those present in the classical algorithm.
Here we do not concentrate on possible quantum speedups and a detailed resource analysis for the quantum algorithm, which lies in the computational complexity P of efficiently simulatable algorithms. Instead, in this work we focus on the advantages that ranking nodes in classical networks using quantum algorithms displays, such as an increased resolution in the structure analysis of scale-free graphs or an increased stability with respect to the variation of the damping parameter.
i.2 Operational Summary of the Quantum PageRank Algorithm
In this section we briefly review the quantum PageRank algorithm from an operational point of view, whereas the reader can find the details of the construction in Paparo and Martin-Delgado (2012). This quantum PageRank algorithm satisfies all the properties of this class and represents a valid quantization of Google’s PageRank algorithm. A step-by-step illustration of the quantum algorithm is presented in Fig. 1.
In Google PageRank’s algorithm the ranking is performed setting up a random walk on the network. The walk uses as a transition matrix, known as the Google matrix . The Google matrix is the weighted sum of two transition matrices. The first walk is driven by a modified connectivity matrix where outgoing links to all other nodes have been added to every node that has no outgoing link. The second walk is a simple random hopping that connects every node to any other node. Accordingly, the Google matrix associated to a given graph is defined as follows:
where is a matrix with entries all set to and is the number of nodes. The parameter is known as the damping parameter.
The ranking of the nodes is performed measuring the probability to find the walker on each node when the stationary distribution has been reached i.e. when .
In the quantum PageRank algorithm the idea is to set up a quantum walk on the nodes of the network and to perform the ranking of the nodes measuring the probability of finding the quantum walker on such nodes.
The quantum walk is a quantization of the Markov chain underlying the classical PageRank algorithm and is set up using Szegedy’s method Szegedy (2004) which allows one to take into account explicitly the connectivity structure, and the directedness of the network in particular. One can quantize a Markov chain on a -vertex graph that has as transition matrix, the Google matrix . This is performed by introducing a discrete-time quantum walk embedding the stochastic matrix on the same graph.
The Hilbert space is the span of all vectors representing the (directed) edges of the graphs i.e. , where the order of the spaces in the tensor product is meaningful here because we are dealing with a directed graph.
For each vertex let us define the quantum state vector,
that is a superposition of the quantum states representing the edges outgoing from the vertex. The weights are given by the (square root of the) Google matrix .
The idea of the quantum PageRank satisfying all the properties that define the class of possible ranking algorithms, is to set up a quantum walk starting from the initial state and whose dynamics is governed by the quantum evolution operator where is the swap operator i.e. and
Although the Hilbert space is clearly dimensional it can be shown that the dynamics takes place in an invariant subspace which is at most dimensional. This allows one to numerically treat networks with a larger number of nodes.
The rankings of the nodes in the quantum network is calculated using the instantaneous quantum PageRank:
Another quantity, called the average quantum PageRank can be defined:
Whereas this quantity can be shown to converge for sufficiently large, the Instantaneous Quantum PageRank does not converge in time Aharonov et al. (2001).
In Paparo and Martin-Delgado (2012) nontrivial features of the quantum PageRank were uncovered, such as the instantaneous outperformance of the algorithm or the violation of the hierarchy as predicted by the classical algorithm.
These properties motivate us to investigate the persistence of these novel effects on larger complex networks. This paper is organized as follows: In section II we present an analysis on graphs of the scale-free type for networks with hundreds of nodes, and in particular for a real-world network. In III) we have extended the analysis of the quantum PageRank algorithm on Erdös-Rényi and hierarchical graphs. From this we conclude that the behavior of the quantum PageRank is characteristic of each type of complex networks, and moreover, it is different from the classical PR algorithm: the classical ranking changes by quantum fluctuations. The new quantum dynamics incorporated in the task of ranking also raises other important questions on the properties of the quantum walk embedded in the algorithm such as the localization phenomena of the walker on the network is analyzed in section IV, while the stability of the ranking with respect to noise (or damping parameter) is studied in section V. In section VI we address the question of whether the power law behavior displayed by the classical PageRanks is preserved by the quantized algorithm and compute the scaling exponent. In section VII se investigate a very practical situation and study how sensitive the quantum PageRank algorithm is under coordinated attacks in scale-free graphs. Section VIII is devoted to conclusions and outlook, and an appendix provides some technical details.
Ii Quantum PageRank on Scale-Free Networks
We will analyze the predictions of the quantum PageRank algorithm on complex networks. We will focus on random scale-free networks because of their widespread appearance and relevance in real-world applications. In the next section we will also deal with the important cases of random (Erdös - Rényi) and hierarchical networks for completeness in order to check whether the quantum PageRank algorithm preserves the characteristics of different classes of complex networks. Moreover, the study of random and hierarchical networks results to be useful to confront its features with scale-free networks and to draw interesting conclusions.
Random scale-free graphs Albert and Barabási (2002); Boccaletti et al. (2006) are ubiquitous in nature. They appear as good models of the World Wide Web Barabási et al. (2000), airline networks Barrat et al. (2004) or metabolic networks Jeong et al. (2001, 2000), just to name a few. These are networks that display a small fraction of hubs, i.e. nodes with a high connectivity, a property that follows from the degree distribution that shows a scale-free behavior, . Scale-free networks exhibit intriguing properties which have been studied extensively, such as robustness against uncoordinated attacks Albert et al. (2000); Callaway et al. (2000); Vázquez and Moreno (2003), good navigability Boguna et al. (2008); Carmi et al. (2009); Lee and Holme (2012) and controllability Liu et al. (2011); Nepusz and Vicsek (2012); Nicosia et al. (2012).
One of the first models proposed to describe scale-free networks is the preferential attachment model Barabási and Albert (1999); Barabási (2013). In this model links are preferentially formed to already highly connected nodes. A random directed scale-free model for the WWW was also introduced in Reka et al. (1999) and a generalization appeared in Bollobás et al. (2003). To produce the characteristic power-law degree distribution of degrees, the models consider two main mechanisms: growth and preferential attachment. A graph is dynamically formed by a continuous addition of new vertices and each new vertex is joined to several existing vertices selected proportionally to their in and out degrees. The generalized model allows also the introduction of directed edges between two already existing nodes. In this work we use graphs created with this model as implemented in NetworkX Hagberg et al. (2008).
We next discuss results originating from the application of the quantum PageRank to complex networks of the scale-free type with hundreds of nodes (). To this end, we have performed a numerical analysis using the quantum PageRank algorithm. We find that the algorithm clearly identifies that the networks are of scale-free type. The algorithm is able to point out the most important hubs. This is a task already well performed by the classical PageRank. However, the quantum PageRank algorithm has improved ranking capabilities, not concentrating all the importance on these few nodes. Indeed, it is capable of unveiling the structure of the graph highlighting also the secondary hubs (see figures 3 , 4, 5 and captions therein).
Furthermore, we find that the hierarchy as predicted by the classical PageRank is not preserved. This is a property already found in Paparo and Martin-Delgado (2012) for smaller networks. From our study we are able to clearly conclude that far from being an artefact of choosing small networks, this results to be a generic feature of the quantum PageRank algorithm. We also found that the quantum PageRank is able to lift the degeneracy of the nodes that have a lower importance. This feature can be seen clearly in figure 2 where we analyze a subgraph of the WWW obtained by exploring pages linking to www.epa.gov and available from Pajek Batagelj and Mrvar (2006).
Iii Quantum PageRank on Erdős-Rényi Networks and Hierarchical Networks
iii.1 Erdős-Rényi Networks
In this section we will briefly introduce the Erdős-R\a’enyi class of random networks and analyze the performance of the quantum PageRank algorithm applied to them.
This class of graphs was introduced by Paul Erdős and Alfred Rényi more than fifty years ago Erdős and Rényi (1959, 1960, 1961), and is of particular importance in the context of Graph Theory. There are different equivalent methods to describe this family. To allow an easy computer implementation we use the following: a graph of order can be constructed by connecting vertices randomly by adding edges with a given probability, which is independent from other edges. We use directed versions of the graphs created with NetworkX Hagberg et al. (2008).
The graphs falling into this class follow a Poissonian degree distribution, i.e., where is the average degree. Thus, most nodes have a degree not far from the average and therefore the graphs do not display relevant hubs when applying the PageRank algorithm.
We have performed a numerical study on graphs of the Erdős-R\a’eny type and found that the quantum PageRank algorithm displays a sharp change in hierarchy (see figure 6 b). We can draw two fundamental conclusions: i/ The behavior of the quantum PageRank algorithm is intrinsic and characteristic of the Erdős-R\a’enyi class of networks, in particular, the no-hub-behavior is reflected at the quantum level. ii/ However, the quantum PageRank algorithm changes the ranking of the classical algorithm, in this sense outperforming the classical PageRank hierarchy.
iii.2 Hierarchical Networks
Some relevant real-life networks which describe technological and biological systems, such as the WWW, some electronic circuits and protein or metabolic networks are usually scale-free but have also a modular structure Ravasz et al. (2002); Song et al. (2005). That is, they are composed of modules that group different sets of nodes. These modules can be distinguished by the fact that nodes belonging to the same module are usually strongly connected. On the other side, modules are relatively weakly connected among them. Thus, even when the networks are scale-free, their hubs use to have a low clustering as they joint different modules. Several authors claim that a signature for a hierarchical network is that, other than the small-world scale-free characteristics, the scaling of the clustering of the vertices of the graph with their degree follows Ravasz and Barabási (2003); Barabási and Oltvai (2004).
Hierarchical network models usually are constructed from recursive rules. For example, we can start from a complete graph and connect to a selected root node replicas of . Next, replicas of the new whole structure are added to this root. At this step the graph will have vertices. The process continues until we reach the desired graph order. There are many variations for these hierarchical networks, depending on the initial graph, the introduction of extra edges among the different copies of the complete subgraphs, etc. However, once the starting graph is given, these networks do not have adjustable parameters and their main characteristics are fixed.
In Barabási et al. (2001), Barabási et al. introduced a simple hierarchical family of networks and showed it had a small-world scale-free nature. The model was generalized in Ravasz and Barabási (2003) and further studied in Noh (2003). For our analysis we have designed a directed version based on these graphs, see figure 7b. In this case the starting point is a directed 3-cycle.
Another interesting family of hierarchical directed graphs has been obtained by giving directions to the edges of the construction published in Comellas and Miralles (2009a, b), see figure 7a. The graphs are in this case small-world, self-similar, unclustered and outerplanar (a planar graph is called outerplanar if it has an embedding where all vertices lie on the boundary of the exterior face). However, they are not scale-free, but follow an exponential distribution. It has been shown that many algorithms which are NP-complete for general graphs perform polynomial in outerplanar graphs Brandstädt et al. (1987).
We have performed a numerical study on hierarchical networks using the quantum PageRank algorithm. We analyzed two families of graphs (see figure 7 for the construction) and we find that the hierarchy (similarly to what was found in ref. Paparo and Martin-Delgado (2012) for the binary tree) is preserved by the average PageRanks. Interestingly, though, the quantum PageRank is able to highlight the connectivity structure of the nodes that belong to the same level in the hierarchical construction (see figures 8 and 9). We observe that the difference in importance between nodes belonging to the same hierarchical level but with different local connectivity is amplified when calculated using the quantum PageRank algorithm.
Iv Localization-Delocalization Transition
We have seen in the previous sections that the quantum PageRank algorithm is able to distinguish between networks of Erdös-Rényi and scale-free type. In particular, the quantum algorithm, as opposed to the classical one, is also able to highlight the secondary hubs in the case of scale-free networks. Furthermore, regarding the classical and quantum algorithms as walks on a network, a necessary condition that few nodes, the hubs, account for the majority of the importance is that the walker be localized, i.e. the number of nodes with a significant probability to find the walker is negligible with respect to the number of nodes in the network.
In this section we will, therefore, study the localization properties of the walker on different networks, of Erdös-Rényi and scale-free type. We will briefly review the case of the walker in the classical PageRank algorithm and explore which phases does the quantum walker choose in the case of the quantum PageRank as a function of the structural properties of the underlying network.
In order to measure the localization phenomenon we will employ the Inverse Participation Ratio (IPR). This concept was introduced in the context of condensed matter and, more specifically, to study the Anderson localization phenomenon in disordered systems (see for example Evers and Mirlin (2008)). It is particularly useful to study localization-delocalization transitions.
The IPR , in the case of a classical walker, is defined as :
where is an integer parameter which can be freely chosen and is fixed.
Let us consider a classical walker on a network. We can have two extreme behaviors. The first one is that the walker is completely delocalized i.e. the probability of finding it on a site is uniform. Therefore, if we introduce a random variable whose realizations are the sites of the lattice, we can write . The other limiting case is that the walker is localized only on one site, i.e. the probability of finding the walker is a Kronecker delta, that is , if the walker is localized on say, site . The IPR (eq. (5)) yields for the two limiting behaviors:
Thus the IPR, displaying respectively, a power law or a constant behavior as a function of the number of nodes, is a legitimate witness of the localization of the walker.
We can rewrite as: . In order to study the localization-delocalization transition, it is useful to introduce the normalized anomalous dimension :
which interpolates between the two phases if the system undergoes a transition from a localized regime (where ) to a delocalized one (where ).
In Georgeot et al. (2010) the localization-delocalization transition for a classical walker performing a random walk (with transition matrix given in (1) ) was characterized studying its dependence on the damping parameter where . This study is important to understand at a deeper level the classical PageRank algorithm. Indeed, as we have anticipated, in the case of a scale-free graph observing localization over a broad range of values for is indeed a necessary condition for the algorithm, quantum or classical, to perform well the task of ranking the nodes. The PageRank vector , is given by:
and represents the stationary probability distribution of the walker on the network.
In ref. Georgeot et al. (2010) it was found that delocalization is absent for a very large range of values of ranging from to .
This is quite natural because for the range of values of the damping factor stated above, the second term, corresponding to random hopping in the Google matrix (see eq. (1)) is suppressed. This yields the localization effect. On the other side, for close to , it was found that the walker is delocalized over the network. This is understood easily, indeed for close to only the second term in the Google matrix survives and the walker is only subject to random uniform hopping between any pair of nodes. This leads to a trivial phenomenon of delocalization of the walker over the whole network.
Let us now focus on the localization phenomenon in the case of the quantum walk on a network according to the quantum PageRank protocol. In order to carry out the analysis we need to generalize the definition of the IPR given above reinterpreting the notion of probability of finding a walker on a node when dealing with the quantum PageRank. In doing so, we will choose as a guiding principle the interpretation of the average quantum PageRank of a node (see eq. (4)) as the probability of finding the quantum walker on a particular node. Therefore, we will employ the definition given in eq. (4) that we will rewrite as:
where is the (strong) measurement operator on the second space (indexing the nodes where the edges point to, see subsect. I.2), i.e. .
We are now in a position to define the IPR in the case of a quantum walk:
Also in the case of a quantum walk we have and it is evident that one can extract the localization phase of a walker from the scaling exponent of the IPR as a function of the number of nodes . Indeed, from equation:
it is clear that the witness of the localization lies in the slope of the graph of the aforementioned log-log plot.
We consider two kinds of networks, of the scale-free and of the Erdös-Rényi type. In order to study the localization phenomenon we generated networks with different numbers of nodes belonging to the two aforementioned classes. We then calculate the IPR (in the following we will fix the parameter ) in order to understand whether for the quantum walker was localized or delocalized on the network.
We find that the IPR in the case of the class of scale-free networks does not vary appreciably (see fig. 10) signaling localization of the walker on the graph. Notice however that also in ref. Georgeot et al. (2010) for these values of a similar behavior was found.
We have analyzed also the graphs in the Erdös-Rényi class performing the same steps as above. Our study shows that, albeit this being random graphs and that for this value of the damping parameter the walk is strongly influenced by the topology of the network, the quantum walker is delocalized in this case. Indeed it can be seen from figure 10 that the behavior of the logarithm of the IPR is linear in the logarithm of the number of nodes. This is a clear witness of the delocalization phenomenon.
This behavior is remarkable for two reasons. The first reason is that in the classical random walk case localization of the walker was found Georgeot et al. (2010). The second reason is that albeit randomness being expected to give rise to localization, in the case of graphs of Erdös-Rényi type, our study shows that the opposite is true.
We conclude that the scale-free graph seems to favor a localization phase in both walks, random and quantum, respectively for the classical and quantumPageRank algorithm. This result in the latter case was obtained with a value of . This is consistent with a good ranking of nodes in a network. Indeed, in order to unveil the main hubs the random or quantum walk must be able to localize the walkers on few important nodes. Interestingly, instead, the Erdös Rényi graphs seem to prefer a delocalized phase albeit the networks being grown randomly. Both for the classical and the quantum PageRank one finds delocalization. This can be correlated with the absence of a small number of main hubs in this class of networks.
As for the classical PageRank the localization-delocalization transition was characterized as a function of the damping parameter, it is important to understand how the quantum PageRank depends on the value of . We will study this dependence and characterize it in the next section.
V Stability of the Quantum Google Algorithm with respect to the Noise Parameter
In the previous section we analyzed the localization properties of the quantum PageRank algorithm. We studied this phenomenon having fixed the damping parameter using the value . The natural question is to study how the quantum PageRank varies with respect to the variation of .
The stability of the quantum PageRank is an important issue to consider also because the damping parameter is arbitrarily tuned to a specific value. Indeed there is no a priori argument to fix the value of and the value was originally chosen in the classical PageRank protocol to mimic the behavior of a surfer (or a random walker) that randomly hops to any other page once every seven times. Only a posteriori it turned out that this is indeed a sensible choice given that the network is small-world, being in fact a crucial ingredient for the PageRank algorithm to yield reasonable ranking results. In view of the ad-hoc choice of the precise of the parameter , it is a very desirable property that the output of the algorithm be stable, i.e. the ranking vary slowly with respect to the variation of the damping parameter. A question that was addressed in the computer science community has been to quantify how susceptible to changes in this parameter the classical PageRank algorithm is. It was found that the effect of this parameter on ranking is large and that two rankings obtained by running the algorithm using different values of this parameter can be very different Georgeot et al. (2010).
To tackle this problem in the quantum case we will make use of two quantities. The first one is related to the quantum fidelity (see e.g. Galindo and Martin-Delgado (2002); Nielsen and Chuang (2010)) that provides a way to measure the distance between two quantum states. The second quantity that we will use is the classical fidelity. It is employed for the same task when dealing with probability distributions.
In ref.Georgeot et al. (2010) the stability of the classical PageRank was studied as a function of the damping parameter. Since the PageRank vectors are classical probability distributions one can measure the distance between two PageRank vectors, calculated using different values of the damping parameter, with the classical fidelity. The latter can be written as:
It was found that there is a plateau around and that the fidelity is rather flat around thus implying that the classical PageRank is rather robust against perturbations (see ref. Georgeot et al. (2010)).
The quantum fidelity is a quantity that measures the distance of two quantum states:
which in the case of commuting density matrices reduces to the classical fidelity.
Another valid measure of the distance of two quantum states is the trace distance:
where denotes the square root of the (positive) operator .
The fidelity and the trace distance turn out to be equivalent measures of distance. Indeed if the fidelity of quantum states is near to one then their trace distance is close to zero and viceversa
Let us rewrite explicitly the definition of the averaged quantum PageRank adding the dependence on (that enters in the initial state and in the evolution operator of the walk) as:
where and for bookkeeping purposes it has been made explicit to which spaces the density matrix refers to.
Let us discuss how to apply the concept of trace distance in our case. We might measure the instantaneous distance , thus measuring the distance of the quantum states. However, it is more significant to take the time-average of the distance between the partial traces of the states, because we are interested in the quantum PageRank as an observable, rather than the state itself. Not being interested in the distance of states, we can make use of a less refined measure of distance (that is, one that can appreciate less the difference of two states). This can be written in the following form (see the appendix for the derivation):
To summarize, we will study the stability making use of the classical fidelity and the quantity in (16), that is a simpler measure descending from the trace distance (which turns out to be equivalent to the quantum fidelity).
We perform an analysis of how the quantum PageRank varies with respect to the rankings when the value of goes from to . We analyze a scale free with 128 nodes that was generated with networkx, a python module. The results show clearly that the quantum PageRanks vary very little when the quantum walk underlying the quantum PageRank has a different damping parameter . Indeed it can be seen from figure 11 that the minimum fidelity between two values of and does not go below the value of . One should compare with the analysis of the classical PageRank Georgeot et al. (2010), where the fidelity between different values of and can be approximately . It can thus be inferred that the ranking is rather robust when it is performed with the quantum PageRank.
We have also investigated the behavior around the value of the damping parameter . From the classical fidelity between the quantum PageRanks at and at ranging from to one can see that there is a plateau around the value of (see figure 11) extending especially for smaller values of . There is a dip for which is due to the fact that the ranking is very sensitive to changes in the damping parameter when this approaches . For this value only the second term in the Google matrix , giving random hopping, survives.
The analysis was made more precise and the conclusions more cogent by considering also the measure of distance of rankings originating from the trace distance (see (16)). Also in this case the overall robustness of the ranking performed with the quantum PageRank is evident. One can see from figure 11 that the maximum value of this measure of distance is comparing any two values of and ranging in the aforementioned interval. It can be observed also that the region where this ranking is more robust is where . Indeed from figure 11 the blue region is wider (and correspondingly the ranking more robust with respect to perturbation of the value of ).
We find also in this case that figure 12 the trace distance is rather smooth for . One can see that also in this case there is a curious peak for for similarly to the previous case.
We conclude that the quantum PageRank as measured by the classical fidelity or by the trace distance seems to vary mildly with respect to the variation of the damping parameter . Indeed the minimum fidelity between any two distributions of importance arising from the quantum PageRank is . The maximum of the distance obtained by using the trace distance between any two states (with different ) is . This means that the quantum PageRank is very robust with respect to variation of the parameter that controls the fraction of random hopping. It is much more robust than in the classical case. From our analysis of the stability of the classical PageRank we find that for extremal values of , the value of the fidelity is less than (see fig. 13). In ref. Georgeot et al. (2010) the minimum fidelity was found to be very close to between PageRanks’ rankings corresponding to extremal values of the damping parameter .
Vi Power Law Behavior for Quantum PageRank
In section II we have found that the quantum PageRank is able to highlight the structure of secondary hubs.
In addition to this finding we have shown in section IV that this strength of the quantum algorithm is associated with the fact that the quantum walker on networks of scale-free type is in a localized phase, i.e. the nodes with a significant average quantum PageRank are a negligible fraction of the nodes in the network. On the other side, for the classical algorithm, it has been shown Donato et al. (2004); Pandurangan et al. (2002) that for real networks the nodes’ classical PageRanks , sorted in descending order, behave following a power law. That is, the classical PageRanks decrease approximately like , where . This is a clear witness of the fact that the algorithm is able to identify the hubs. Furthermore, the scaling exponent measures the relative importance given to the hubs with respect to the other nodes.
A similar study is therefore important in the case of the quantum PageRank and will be carried out in this section. Indeed, it is desirable that the quantum PageRank display a power law scaling behavior. Indeed, such a behavior is distinctive of the fact that the algorithm be able to uncover the scale free nature of the network.
In order to point out the scaling behavior of the quantum PageRanks of the nodes it is clear from the conjectured form
that we can extract such behavior from the slope of the log-log plot of the quantum PageRanks versus the (sorted) index of the as can be seen from
Here we perform the analysis on scale free networks of 256 nodes. We calculate the classical and quantum PageRanks of the nodes and after having sorted the nodes in descending order we analyze the log-log plot of the classical and quantum PageRanks versus the nodes’ index.
Considering one instance of a graph in this class one can clearly see (cf. figure 14) that both rankings obtained using the classical and the quantum PageRank display a power law behavior. The plot displays three areas, corresponding respectively to the hubs, to the intermediate part of the list and the low part of the ranking (which in the classical case is degenerate in importance).
The fact that this behavior persists also in the case of the quantum PageRank is due to the effect, similar to what has been observed in Georgeot et al. (2010) for the classical case (and in section IV for the quantum Pagerank that the walker is in a localized phase in the case of scale free networks). Consequently the most highly ranked nodes tend to concentrate nearly the totality of the importance. This can be clearly seen in figure 14, area , where the hubs’ importances lie above the line.
Furthermore, it is clear that the scaling coefficients are different in the quantum and classical case. We find that the quantum PageRank has a smoother behavior, giving less relative importance to the nodes in the high part of the list. On the other hand it is also able to better rank in the low part of the list (where the classical PageRank gives highly degenerate values) lifting the degeneracy. (cf. the area in figure 14).
The quantum PageRank is therefore able to better distinguish the nodes’ importances in the lowest part of the list. This is because the difference between the importance of the nodes in the higher and lower part of the list is lower. Furthermore, the power law behavior interpolates a greater portion of the data with respect to the classical case as can be noted by inspection (the area in figure 14 extends much more in the case of the quantum algorithm).
To complete the analysis we consider an ensemble of scale free networks in order to display the ensemble’s properties rather than the particular instance’s. The ensemble consists of 29 scale free networks. It is clear that the properties found in the instance of the graph in figure 14 persist also when considering a mean property of the ensemble (see figure 14). We conclude that these properties are generic and not an artifact of considering only one instance of the class of scale-free networks.
Finally, we consider the real-world network, a subgraph of the WWW obtained by exploring pages linking to www.epa.gov Batagelj and Mrvar (2006), which we refer as EPA in the following (see figure 15). Also in this case we find a power law behavior with . The plot displays three areas as mentioned above for the other cases and our conclusions are valid also for this real-world network.
Vii Sensitivity of the quantum PageRank algorithm under coordinated attacks in scale-free graphs
In this section, we aim to study of the sensitivity of the quantum PageRank algorithm under coordinated attacks. More precisely, we ask how much the list of quantum PageRanks of an -node graph changes as a whole if the most important nodes (hubs) of a network fail (e.g. due to a hacker attack) and the quantum quantum PageRank algorithm is run on the remaining - node graph. Motivated by the fact that the real-world internet belongs to the class of scale-free networks, we focus in our study on the scalefree networks of mesoscopic size (graphs of 16 and 32 nodes).
Operationally, we proceed in our numerical study as follows:
(i) First, we determine for the initial -node graph the quantum PageRank values and the ordered list of nodes according to the quantum PageRank algorithm.
(ii) Next, we take out the most important node (main hub) from the graph by eliminating the node itself as well as all its in- and outgoing links from and to the other nodes. On this reduced graph (with corresponding modified connectivity matrix ) we carry out the quantum PageRank algorithm to determine of the modified list of quantum PageRanks, with a possibly modified order list of nodes. We note that on this modified graph the algorithm differs quantitatively from the one run on the complete -node network, as the Hilbert space, the initial state and the underlying coherent dynamics are different due to the modified connectivity matrix and resulting modified Google matrix (see Sec. I.2).
(iii) We compare the ordered list of () nodes according to the quantum PageRank values with the original list for the -node graph, where the most important node is taken out (see Fig. 16a). To quantify the overall difference between these two lists of elements, we employ Kendall’s coefficient Kendall and Smith (1939). This function returns for lists in which the order of all elements is the same (irrespective of the individual values associated to each element of the list), zero for lists whose order of elements is exactly the opposite, and values in between for lists where the order of elements partially differs.
This procedure of steps (i) to (iii) can be iterated to take out subsequently the most important nodes according to the initial quantum PageRank list of the -node graph. The resulting quantum PageRank list of the reduced - node graph is then compared to the initial list (not including the most important nodes).
To compare the sensitivity of the quantum PageRank algorithm with the classical PageRank algorithm, we perform the same type of coordinated attacks in the classical scenario, i.e., we analyze how the ordered list of classical PageRank values changes when the most important (according to the classical PageRank protocol) nodes of the graph are taken out and the classical PageRank algorithm is run to determine the importance of nodes in the reduced graphs.
The results are shown in Fig. 16b and c. We find that when the most important nodes (hubs) are attacked and fail, and the PageRank is recalculated for the reduced graphs, the order of the importance of the remaining nodes changes with respect to the list of the initial complete graph both in the classical and in the quantum case. In the quantum PageRank algorithm, attacks on hubs have a stronger effect than for the classical algorithm. This behavior can be understood by the fact that whereas for the classical case there is a large degeneracy of importance values of nodes of low PageRank, quantum fluctuations partially lift this degeneracy – see discussion in Sec. II and the insets in Fig. 2. Thus, when hubs of the network are attacked and fail, the order of less important nodes – whose importance values slightly differ, can truly change in the quantum case, whereas the degeneracy of a larger number of nodes persists in the classical case. In other works, the increased capability of the quantum algorithm to resolve more structural details of the directed graphs, comes at the cost of an increased sensitivity to structural changes of the network.
In this paper we have studied the behavior of the quantum PageRank algorithm, developed in Ref. Paparo and Martin-Delgado (2012), when applied to complex networks. We have found that the quantum algorithm is able to clearly distinguish the structure of the underlying network. More specifically, the quantum PageRank’s behavior is distinctive for the three classes of complex networks studied in this work: scale-free networks, graphs of Erdös-Rényi type and hierarchical networks. In particular, we have observed that the quantum algorithm when applied on scale-free networks is able to highlight the structure of the secondary hubs and to resolve the degeneracy in importance of the low lying part of the list of rankings, which represents a typical shortcoming of the classical PageRank algorithm. Although best suited for scale-free graphs the quantum PageRank is also able to univocally uncover whether graphs lie in the Erdös-Rényi class. Applied to hierarchical graphs the algorithm has the capability to better reveal the hierarchy of levels, of which the graph is composed, and to highlight the connectivity structure within every hierarchy layer better than its classical PageRank counterpart.
Regarding the quantum PageRank algorithm as a directed quantum walk, we have studied the localization properties of the quantum walker in the quantum protocol. By an analysis of the Inverse Participation Ratio (IPR), we have observed localization of the quantum walker in the case of the quantum PageRank applied to scale-free networks under standard conditions (damping parameter ). This finding is consistent the ability of the quantum algorithm to highlight hubs of the network. In contrast, for Erdös-Rényi graphs delocalization was found, which is in accordance with the absence of hubs for this class of networks.
Furthermore, we have analyzed the robustness of the quantum algorithm with respect to variations of the damping parameter. We find a very high degree of robustness as compared to the classical PageRank protocol, which indicates that the value of this parameter, whose choice is to some extent arbitrary, turns out to be not crucial for the quantum algorithm to work reliably.
Furthermore, we have found that the distribution of importance values of quantum PageRanks over the nodes of scale-free networks follow - as has been previously found for the PageRank - also a power law behavior. However, the corresponding scaling exponent is for the quantum protocol smaller than in the classical case, indicating a smoother ranking of nodes. In contrast to the classical algorithm, in the quantum protocol the hubs of the graphs do not concentrate the whole importance and the algorithm lifts the degeneracy of the large set of nodes with low importance values. This increased ranking capability comes at the cost of being more sensitive to structural changes to the network such as coordinate attacks on hubs.
Remarkably, the described characteristics of the quantum PageRank even persist if the algorithm is applied to real-world networks. We have studied and successfully tested the performance of the algorithm by applying it to a real-world network, originating from the hyperlink structure of www.epa.org Batagelj and Mrvar (2006), thereby showing that its intriguing properties are not restricted to an application of the algorithm to artificially, numerically grown networks.
The classical PageRank algorithm has been the subject of exact studies yielding analytical results, and other interesting studies de Lima Marquezino et al. (2013); Venegas-Andraca (2012); Whitfield (2012). It would be nice if the quantum algorithm can yield also exact analytical results and for this the class of hierarchical networks is a good candidate.
A related subject of current study is the entanglement properties of quantum complex networks Cuquet and Calsamiglia (2012); Acín et al. (2007); Cuquet and Calsamiglia (2009). A different line of studies has pursued the application of the quantum adiabatic algorithm to the classical PageRank algorithm Garnerone et al. (2012).
In future work it will be interesting to analyze in more detail the impact of random failures of nodes in large networks of differing topology. From an algorithmic point of view, it is an interesting task to develop a dissipative version of this algorithm and to understand its performance and robustness properties in such scenario. Dissipation has already been considered as an element with respect to some aspect of the algorithm Garnerone (2012); Sánchez-Burillo et al. (2012), but the development of a truly dissipative version in the spirit of dissipative quantum algorithms and computation Verstraete et al. (2009); Diehl et al. (2008) remains an open question. Furthermore the growing field of complex quantum networks would benefit from a version of the algorithm that is able to rank nodes in the more general case where qubits are located at the nodes of the network. An important question in this scenario is whether an algorithm based on a multi-particle quantum walk Childs et al. (2013); Childs (2009) is needed in this context, or if there exists for this task an efficiently simulatable algorithm that belonging to the computational complexity class .
Acknowledgements. The authors acknowledge the Centro de Supercomputacion y Visualizacion de Madrid (CeSViMa) for CPU time on the Magerit2 cluster. This work has been supported by the Spanish MINECO grants and the European Regional Development Fund under projects FIS2012-33152, MTM2011-28800-C02-01, CAM research consortium QUITEMAD S2009-ESP-1594, European Commission PICC: FP7 2007-2013, Grant No. 249958 and UCM-BS grant GICC-910758.
From eq. (15):
where for bookkeeping purposes in the derivations that follow it has been made explicit to which spaces the density matrix refers to. Using the fact that trace preserving quantum operations are contractive. That is:
and making use of another property, namely that the trace distance of two states is attainable as the maximum over all strong (i.e. PVM) measurement outcomes of the difference of the two states:
that, in our case, specializes to:
we obtain Finally it is useful to measure the average distance as we are interested in the stability of the average quantum PageRank:
eq. (22) can be rewritten as:
which in our case is clearly stated as:
Where we can take the absolute value because the distance between two states is a nonnegative number.
- This can be seen from a general formula that relates the two measures of distance (see e.g. Nielsen and Chuang (2010) chap. 9):
- S. Brin and L. Page, Computer networks and ISDN systems 30, 107 (1998).
- L. Page, S. Brin, R. Motwani, and T. Winograd, (1999).
- G. Paparo and M. Martin-Delgado, Scientific Reports 2 (2012).
- D. Watts and S. Strogatz, Nature 393, 440 (1998).
- A.-L. Barabási and R. Albert, Science 286, 509 (1999).
- C. Song, S. Havlin, and H. A. Makse, Nature 433, 392 (2005).
- E. Ravasz and A.-L. Barabási, Phys. Rev. E 67, 026112 (2003).
- A.-L. Barabási and Z. N. Oltvai, Nat. Rev. Genet. 5, 101 (2004).
- A.-L. Barabási, Z. Dezso, E. Ravasz, S.-H. Yook, and Z. Oltvai, in AIP Conference Proceedings, Vol. 661 (2003) p. 1.
- M. E. Newman, SIAM Review 45, 167 (2003).
- M. Newman, Networks: an introduction (Oxford University Press, Inc., 2010).
- R. Cohen and S. Havlin, Complex networks: structure, robustness and function (Cambridge University Press, 2010).
- P. Erdős and A. Rényi, Publ. Math. Debrecen 6, 290 (1959).
- A.-L. Barabási, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 371 (March 28, 2013), 10.1098/rsta.2012.0375.
- B. Bollobás, C. Borgs, J. Chayes, and O. Riordan, in Proceedings of the fourteenth annual ACM-SIAM symposium on Discrete algorithms (Society for Industrial and Applied Mathematics, 2003) pp. 132–139.
- C. Elliott, Quantum Communications and cryptography , 83 (2006).
- A. Poppe, M. Peev, and O. Maurhart, International Journal of Quantum Information 6, 209 (2008).
- M. Sasaki, M. Fujiwara, H. Ishizuka, W. Klaus, K. Wakui, M. Takeoka, S. Miki, T. Yamashita, Z. Wang, A. Tanaka, et al., Optics Express 19, 10387 (2011).
- http://swissquantum.idquantique.com/ (2009).
- D. Lancho, J. Martinez, D. Elkouss, M. Soto, and V. Martin, Quantum Communication and Quantum Networking , 142 (2010).
- T. Länger and G. Lenhart, New Journal of Physics 11, 055051 (2009).
- M. Szegedy, in Foundations of Computer Science, 2004. Proceedings. 45th Annual IEEE Symposium on (IEEE, 2004) pp. 32–41.
- D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani, in Proceedings of the thirty-third annual ACM symposium on Theory of computing (ACM, 2001) pp. 50–59.
- R. Albert and A.-L. Barabási, Rev. Mod. Phys. 74, 47 (2002).
- S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D.-U. Hwang, Phys. Rep. 424, 175 (2006).
- A.-L. Barabási, R. Albert, and H. Jeong, Physica A: Statistical Mechanics and its Applications 281, 69 (2000).
- A. Barrat, M. Barthelemy, R. Pastor-Satorras, and A. Vespignani, Proceedings of the National Academy of Sciences of the United States of America 101, 3747 (2004).
- H. Jeong, P. Mason, A.-L. Barabasi, and Z. Oltvai, Nature 411, 41 (2001).
- H. Jeong, B. Tombor, R. Albert, Z. N. Oltvai, and A.-L. Barabási, Nature 407, 651 (2000).
- R. Albert, H. Jeong, and A.-L. Barabási, Nature 406, 378 (2000).
- D. S. Callaway, M. E. Newman, S. H. Strogatz, and D. J. Watts, Phys. Rev. Lett. 85, 5468 (2000).
- A. Vázquez and Y. Moreno, Phys. Rev. E 67, 015101 (2003).
- M. Boguna, D. Krioukov, and K. Claffy, Nat. Phys. 5, 74 (2008).
- S. Carmi, S. Carter, J. Sun, and D. Ben-Avraham, Phys. Rev. Lett. 102, 238702 (2009).
- S. H. Lee and P. Holme, Phys. Rev. Lett. 108, 128701 (2012).
- Y.-Y. Liu, J.-J. Slotine, and A.-L. Barabási, Nature 473, 167 (2011).
- T. Nepusz and T. Vicsek, Nat. Phys. 8, 568 (2012).
- V. Nicosia, R. Criado, M. Romance, G. Russo, and V. Latora, Scientific Reports 2 (2012).
- A. Reka, H. Jeong, and A.-L. Barabasi, Nature 401, 130 (1999).
- A. Hagberg, P. Swart, and D. S Chult, Exploring network structure, dynamics, and function using NetworkX, Tech. Rep. (Los Alamos National Laboratory (LANL), 2008).
- V. Batagelj and A. Mrvar, Web page http://vlado. fmf. uni-lj. si/pub/networks/data (2006).
- P. Erdős and A. Rényi, Magyar Tud. Akad. Mat. Kutató Int. Közl 5, 17 (1960).
- P. Erdős and A. Rényi, Acta Mathematica Hungarica 12, 261 (1961).
- E. Ravasz, A. L. Somera, D. A. Mongru, Z. N. Oltvai, and A.-L. Barabási, Science 297, 1551 (2002).
- A.-L. Barabási, E. Ravasz, and T. Vicsek, Physica A: Statistical Mechanics and its Applications 299, 559 (2001).
- J. D. Noh, Phys. Rev. E 67, 045103 (2003).
- F. Comellas and A. Miralles, Physica A: Statistical Mechanics and its Applications 388, 2227 (2009a).
- F. Comellas and A. Miralles, Journal of Physics A: Mathematical and Theoretical 42, 425001 (2009b).
- A. Brandstädt, J. P. Spinrad, et al., Graph classes: a survey, Vol. 3 (Society for Industrial Mathematics, 1987).
- F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 (2008).
- B. Georgeot, O. Giraud, and D. L. Shepelyansky, Phys. Rev. E 81, 056109 (2010).
- A. Galindo and M. A. Martin-Delgado, Rev. Mod. Phys. 74, 347 (2002).
- M. A. Nielsen and I. L. Chuang, Quantum computation and quantum information (Cambridge university press, 2010).
- This can be seen from a general formula that relates the two measures of distance (see e.g. Nielsen and Chuang (2010) chap. 9): .
- D. Donato, L. Laura, S. Leonardi, and S. Millozzi, The European Physical Journal B-Condensed Matter and Complex Systems 38, 239 (2004).
- G. Pandurangan, P. Raghavan, and E. Upfal, Computing and Combinatorics , 1 (2002).
- M. G. Kendall and B. B. Smith, The annals of mathematical statistics , 275 (1939).
- F. de Lima Marquezino, R. Portugal, and S. Boettcher, Phys. Rev. A 87, 012329 (2013).
- S. E. Venegas-Andraca, Quantum Information Processing , 1 (2012).
- J. D. Whitfield, (2012).
- M. Cuquet and J. Calsamiglia, Phys. Rev. A 86, 042304 (2012).
- A. Acín, J. I. Cirac, and M. Lewenstein, Nat. Phys. 3, 256 (2007).
- M. Cuquet and J. Calsamiglia, Phys. Rev. Lett. 103, 240503 (2009).
- S. Garnerone, P. Zanardi, and D. A. Lidar, Phys. Rev. Lett. 108, 230506 (2012).
- S. Garnerone, Phys. Rev. A 86, 032342 (2012).
- E. Sánchez-Burillo, J. Duch, J. Gómez-Gardeñes, and D. Zueco, Scientific Reports 2 (2012).
- F. Verstraete, M. M. Wolf, and J. I. Cirac, Nat. Phys. 5, 633 (2009).
- S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. Büchler, and P. Zoller, Nat. Phys. 4, 878 (2008).
- A. M. Childs, D. Gosset, and Z. Webb, Science 339, 791 (2013).
- A. M. Childs, Phys. Rev. Lett. 102, 180501 (2009).