Complex Networks: from Classical to Quantum

Recent progress in applying complex network theory to problems faced in quantum information and computation has resulted in a beneficial crossover between two fields. Complex network methods have successfully been used to characterize quantum walk and transport models, entangled communication networks, graph theoretic models of emergent space-time and in detecting community structure in quantum systems. Information physics is setting the stage for a theory of complex and networked systems with quantum information-inspired methods appearing in complex network science, including information-theoretic distance and correlation measures for network characterization. Novel quantum induced effects have been predicted in random graphs---where edges represent entangled links---and quantum computer algorithms have recently been proposed to offer super-polynomial enhancement for several network and graph theoretic problems. Here we review the results at the cutting edge, pinpointing the similarities and reconciling the differences found in the series of results at the intersection of these two fields.

Quantum mechanics has long been predicted to help solve computational problems in physics [1], chemistry [2], and machine learning [3] and to offer quantum security enhancement in communications [4], including a quantum secure Internet [5].Network and graph theory fundamentally arises in nearly all aspects of quantum information and computation.Rapid experimental progress has pushed quantum computing and communication devices into truly data-intensive domains, where even the classical network describing a quantum system can exhibit complex features, giving rise to what appears as a paradigm shift needed to face a fundamental type of complexity [6][7][8][9][10][11][12][13].Methods originating in complex networks-traditionally based on statistical mechanics-are now being generalized to the quantum domain in order to address these new quantum complexity challenges.
Two types of quantum networks have been of primary focus in the series of pioneering results we review.The first consists of quantum systems whose connections are represented by entangled states [6,32,33].These quantum networks are used in secure quantum communication systems.The second area of focus consists of networks of quantum systems, such as atoms or superconducting quantum electronics, whose connections are physical [34][35][36][37][38][39].Such systems are used to develop quantum-enhanced algorithms or quantum information transport systems, both modeled by quantum walks on complex networks.At a fundamental level, the two types of quantum networks are described by quantum information theory, allowing one to extend the spectrum of network descriptors-such as ranking indicators, similarity and correlation measures-inside the quantum domain.
Interestingly, the same tools can then be appropriately modified to apply to traditional complex networks, suggesting the existence of a framework-network information theory-suitable for application to both classical and quantum networked systems [13,40,41].This bidirectional cross-over is carving out a coherent path forward built fundamentally on the intersection of these two fields (see Fig. 1).Fundamentally, several quantum effects are still outside of the predictive range of applicability of complex network theory.Future work should build on recent breakthroughs and head towards a new theory of complex networks which augments the current statistical mechanics approach to complex networks, with a theory built fundamentally on quantum mechanics.Such a unified path forward appears to be through the language of information theory.
Here, we make an effort to review some of the crucial steps towards the creation of a network theory based fundamentally on quantum effects.Therefore, we do arXiv:1702.08459v1[quant-ph] 27 Feb 2017 not cover several topics that, nevertheless, deserve to be mentioned as part of the field.These include, in no particular order, quantum gravity theories based on complex networks [42][43][44][45], synchronization in and on quantum networks [46], quantum random circuits [47,48], classical spin models and quantum statistics successfully used in complex network theory [49][50][51] (see [52] for a thorough review).

I. QUANTUM NETWORKS BASED ON ENTANGLED STATES
The reader not familiar with some basic concepts of network science might find useful a few definitions.A network is an abstract representation of relationships (encoded by edges) between units (encoded by nodes) of a complex system.Edges can be directed, i.e. they can represent information incoming to or outgoing from a node and, in general, they can be weighted by real numbers.The number of incoming, outgoing and total edges is known as incoming, outgoing and total degree of a node, respectively.The sum of the corresponding weights defines the incoming, outgoing and total strength of that node, respectively.Networks are often characterized by how node degree and strength are distributed and correlated.Systems modeled by uncorrelated networks with homogeneous degree distribution are known as Erdos-Renyi networks, whereas systems with power-law degree distribution are known as scale-free networks.We refer to [16,53] for reviews of network concepts and models.
To define quantum networks based on entangled states, let us start from the state of each i th qubit, written without loss of generality, as with |0 and |1 the preferred or 'computational' basis.The qubit is in a pure, coherent superposition of the two basis states and any measurement in this same basis will cause the state to collapse onto |0 or |1 , with probability cos 2 (α i ) and sin 2 (β i ), respectively.Let us consider a quantum system with two qubits, i.e. i = 1 and 2. The basis of this system is given by the so called, tensor product, of the two basis states: |00 , |01 , |10 and |11 .
If the two qubits are not entangled, i.e. their states are independent from each other, then the state of the overall system can be written as e.g.|ψ 12 = |ψ 1 ⊗ |ψ 2 , whereas this is not possible if the two qubits are entangled.A generalization of this description to the case of mixed states is obtained in terms of the non-negative density matrix ρ; a unit trace Hermitian operator representing the state of the system as an ensemble of (unknown) pure states.
Instead of distributing entanglement on regular graphs, such as uniform lattices typically studied in condensed matter physics, it has been shown that it is possible to tune the amount of entanglement between two nodes in such a way that it equals the probability to have a link in (classical) Erdos-Renyi graphs [6].Such random graphs can be defined by the family of networks G(N, p), where N is the number of nodes and p the probability to find a link between any two nodes.The probability scales with the size of the network following a power law p ∝ N −z , with z ≥ 0. In classical network theory, there exists a critical value for the probability p c (N ) for which, if p > p c (N ) a given subgraph of n nodes and l links has higher probability to be observed.The classical result is that this critical probability scales with N as p c (N ) ∝ N −n/l .Acin, Cirac, and Lewenstein [6] formulated an elegant extension of this picture to the quantum realm by replacing each link with an entangled pair of particles, where the probability p i,j = p that the link exists between nodes i and j is substituted by a quantum state ρ i,j := ρ of two qubits, one at each node.One can build a quantum network where each node consists of N − 1 qubits which are entangled, in pairs, with qubits of other nodes.However, in this case, although the connections are identical and pure they encode non-maximally entangled pairs.For pure states of qubits, ρ = |φ φ|, with Here, 0 ≤ p ≤ 1 quantifies the entanglement of links and the state of the overall quantum random graph can be denoted by |G(N, p) .If each link, i.e. each entangled pair, attempts to convert its state to the maximally entangled one (p = 1/2) through local operations and classical communication (LOCC), the optimal probability of successful conversion is exactly p.It follows that the fraction of existing entangled states converted to maximally entangled ones by LOCC corresponds to the probability of having a link between nodes in the corresponding classical random network [6].By varying the value of the parameter z, i.e. how the critical probability scales with system size, it is possible to control the number and type of subgraphs present in a quantum network of N nodes.This is useful to create special multipartite states, such as the Greenberger-Horne-Zeilinger state which exhibits nonclassical correlations [54].The striking result is that it is possible to obtain, with probability approaching unity, a quantum state with the topology of any finite subgraph for N approaching infinity and z = −2.This bridge between complex network theory and quantum theory provides a powerful tool to investigate the critical properties of a quantum system.For instance, in the case of regular lattices, it has been shown that the probability p opt to establish a perfect quantum channel between the nodes can be mapped to the probability of distributing links among each pair of nodes in the lattice [6], a scenario that can be studied using the well-established bond-percolation theory from statistical physics.This result allows one to calculate the critical probability above which the system will exhibit an infinite connected cluster and, in the case of qubits, it has been shown that the probability of having an entangled path with infi- FIG. 1. Cross-pollination between the fields of complex networks and quantum information science.Each of the shaded regions represent published findings that map out the field from theoretical, experimental and computational perspectives.The top area classifies the analytical tools inspired by quantum information for classical network analysis and vice versa-both covered in this review-as well as the algorithms, genuinely quantum, developed to address specific problems in network science.
The classification in the bottom area includes quantum networks based on entangled states and on physical connections-covered in this review-as well as quantum network models of space-time, random quantum circuits, random tensor networks and geometry.
nite length-i.e., an infinite sequence of entangled states connecting an infinite number of qubits-is unity.However, for product states this probability is zero, denoting the existence of a sharp transition between these two scenarios.However, local measurements based on this approach, called classical entanglement percolation (CEP), are not optimal, in general, to generate maximally entangled states: CEP is not even asymptotically optimal for two-dimensional lattices and new quantum protocols based on quantum entanglement percolation have to be used instead [6].A novel critical phenomenon, defining an entanglement phase transition, emerges from this new strategy, where the critical parameter is the degree of entanglement required to be distributed in order to establish a quantum channel with probability that does not decay exponentially with the size of the system, at variance with CEP.This type of enhancement with respect to the classical case has been reported for different network topologies, such as Erdos-Renyi, scale-free and small-world networks [32].
The static entangled states providing the network connectivity described here, will be replaced in the next section by dynamical processes on networked quantum systems.

FIG. 2. Entanglement connecting distant network nodes.
A pair of distant cavities in a quantum communication network are driven by a shared squeezed light source S. In the systems steady state the two atoms A and B entangle, forming a network edge [33].

II. QUANTUM NETWORKS BASED ON PHYSICAL CONNECTIVITY
Another wide area where network concepts have found applicability consists of quantum systems physically interconnected, such as atoms or superconducting quantum electronics [34][35][36][37][38][39].These types of systems provide fertile ground where quantum algorithms are tested [55][56][57][58][59] and quantum information transport systems are studied [60][61][62].Typical modeling approaches are based on so called 'quantum walks' on complex networks, with recent studies showing that quantum information tasks, typically designed for simple topologies, retain performance in very disordered structures [63].Stochastic (non-quantum) walks are also a central model in complex network theorysee the review [64].
Any quantum process can be viewed as a single particle walk on a graph.Single-particle quantum walks represent a universal model of quantum computation--meaning that any algorithm for a quantum computer can be translated into a quantum walk on a graph-and, additionally, quantum walks have been widely studied in the realm of quantum search on graphs, in both continuous and discrete time via coined walks (see e.g.[55][56][57][58]-in particular the graph optimality results [63]).The computational advantages of quantum versus stochastic random walk based algorithms has attracted wide interest with typical focus being on general graphs which consequently do not exhibit complex features.However, many works have compared properties of stochastic [65,66] and quantum random walks [60][61][62] on complex networks [67].
Network topology has further been shown as a means to direct transport by adding complex numbers-while maintaining Hermiticity-to the networks adjacency matrix in 'chiral quantum walks' [62,68,69] (note that chiral walks were realized experimentally in [12]).Open system walks which mix stochastic and quantum effects in 'open' evolutions [61] have aided in the study of quantum effects in biological exciton transport (again, modeled as a quantum walk) and developments in a quantum version of Google's PageRank [8][9][10] has been seen, providing a practical solution to overcome the degeneracy issues affecting the classical version and enhancing node ranking in large networks.Recently, Faccin et al. have analytically solved a model which shed light on some key differences between stochastic and quantum walks on complex networks [7].These differences push forward a general understanding which can lead to a theory explaining novel complex features in quantum systems.
Quantum walks on complex networks represent both a practical model of transport [67] as well as an interesting stage of comparison between the quantum and stochastic cases.As a closed quantum system exhibits fluctuations in the probabilities in time, typically a long time average is considered.Physically, this is the best approximation one can hope for, provided that there is no knowledge of when the walk started.In this case, the probability to find a quantum walker in the i-th node is given by where |0 is the initial state and U t = e −iQt is the unitary evolution operator defined by the quantum generator Q.
Interference between subspaces of different energy vanish in the long time average so we obtain an expression for the probability (P Q ) i in terms of the energy eigen-space projectors Π j of the Hamiltonian H Q , Here Π j = k |φ k j φ k j | projects onto the subspace spanned by the eigenvalues |φ k j of H Q corresponding to the same eigenvalue λ j .

Quantum-enhanced page-ranking.
The nonsymmetric adjacency matrix representing the directed connectivity of the World Wide Web, a.k.a. the Google matrix G, satisfies the Perron-Frobenius theorem [70] and hence there is a maximal eigenvalue corresponding to an eigenvector of positive entries Gp = p.The eigenvector p corresponds to the limiting distribution of occupation probabilities of a random web surfer-it represents a unique attractor for the dynamics independently of the initial state.The vector p is known as the Page-Rank.
Several recent studies embed G into a quantum system and consider quantum versions of Google's Page-Rank [8][9][10].Garnerone et al. [71] relied on an adiabatic quantum algorithm to compute the Page-Rank of a given directed network, whereas Burillo et al. [11] rely on a mixture of unitary and dissipative evolution to define a ranking that converges faster than classical PageRank.The page-ranking vector p is an eigenvector of I − G corresponding to the zero eigenvalue (the lowest).This fact leads to a definition of a Hermitian operator which can play the role of a Hamiltonian, defined as: though highly non-local, its ground state represents the target Page-Rank which could be found by adiabatic quantum annealing into the ground state.Using a quantum computer to accelerate the calculation of various network properties has been considered widely [55][56][57][58].
As Page-Rank relies on finding the vector corresponding to the lowest eigenvalue of the Google matrix, the adiabatic algorithm opens the door up to accelerate network calculations using quantum computers.
Directing transport by symmetry breaking in chiral walks.Chiral quantum walks, introduced by Zimborás et al. in [62] and realized experimentally in [12] append complex numbers to the adjacency matrix (playing the role of the system Hamiltonian) while still maintaining the Hermitian property [12,62,68,69].These complex phases in many cases do not effect transfer probabilities: the theory explaining this finding was developed in [12,62], without relying on approximations or averaging.The case of open systems has been investigated as well in [62].In the scenarios where the addition of complex phases affects transfer probabilities, the underlying system breaks time-reversal symmetry and, consequently, the probability flow into the quantum system is biased.This fact enables directed state transfer without requiring a biased (or non-local) distribution in the initial states, or coupling to an environment.
When the underlying graph is bipartite (e.g. a graph whose vertices can be divided into two disjoint sets such as a square lattice), time-reversal symmetry in the transport probabilities can not be broken.Transport suppression is indeed possible however [62].Bipartite graphs include trees, linear chains and generally, graphs with only even cycles.These results point to a subtle interplay between the topology of the underlying graph, giving rise to a new challenge for dynamical control of probability transfer when considering walks on complex networks [12,62,68,69].

Open quantum walks.
The area of open quantum systems [72] studies noise and its effects in quantum systems.The adiabatic version of Page-Rank [9] uses a quantum stochastic quantum walk as proposed by [61] (see also [73,74] for studies on open walks).
Quantum stochastic walks are defined by a quantum walk undergoing dissipative dynamics.The latter follows the quantum master equation in the Lindbladian form: where L k represents a jump operator while [•, •] and {•, •} are commutator and anti-commutator respectively.The network topology is embedded by choosing H equal to the adjacency matrix of the symmetrized network and In this picture, node ranking is defined by an activity vector α computed at the steady state ρ ss .
Paparo et al. [8,10] introduced a Szegedy type of Markov chain quantization [75] of the random walk.In order to quantize the Markov chain defined by the Google matrix G of N nodes, one introduces a Hilbert space H = span{|i 1 |j 2 , i, j ∈ [0, N ]} and the superposition of outgoing edges from node i: Each step of the quantum walk U is defined by a coin flip 2Π − 1 and a swap operation S which ensures unitarity [8] U = S(2Π − 1) the swap operator is S = ij |ij ji|.
In the case of quantum Page-Rank, this is set to the instantaneous probability P (i, t) of finding the walker at node i at the time-step t.To obtain a fixed value for the quantum Page-Rank a time average is calculated as long as with its variance as a measure for quantum fluctuations.
Another approach [11] involves defining a Markov quantum evolution similar to (7), with a tuning parameter α: where the Hamiltonian H is the symmeterized adjacency matrix and L k = L ij = G ij |i i| represent the jump operators which consider the directness of the network.With this definition, for values of α ∈ (0, 1], a stationary state is guaranteed.In this case, for α = 0 we revert to the unitary evolution while for α = 1 we revert to the stochastic case.
The authors [8,10,11] show how this definition of quantum Page-Rank resolves problems of degeneracy in the classical Page-Rank definition, enhances the importance of secondary hubs and, for certain values of α, the algorithm exhibits faster convergence.

Q u a n t u m
Sto c h a s t i c Spec t r a l S i m il ar FIG. 3. Known mappings between quantum and stochastic generators.Here G = (V, E) is a graph with adjacency matrix A, D is a diagonal matrix of node degrees.This yields the graph Laplacian L = D − A, and hence, the stochastic walk generator LS = LD −1 , from this a similarity transform results in LQ = D −1/2 LD −1/2 , which generates a valid quantum walk and exhibits several interesting connections to the classical case.The mapping L −→ L L preserves the lowest 0 energy ground state, opening the door for adiabatic quantum annealing which solves computational problems by evolving a system into its ground state.

III. TOWARDS UNIFIED ANALYSIS OF NETWORK COMPLEXITY
The interaction between network science and quantum information science has led to the development of theoretical and computational tools that benefitted from both fields.On the one hand, quantum-inspired tools, such as information entropies and quantum distance measures, have been successfully applied to practical problems concerning classical complex networks [13].On the other hand, classical network descriptors have been ported to the quantum realm to gain better insights about the structure and the dynamics of networked quantum systems [76].
The cross-pollination between the two briefly reviewed in this section, is still ongoing with vibrant future research opportunities.
Information entropy of classical networks.Historically, the concept of entropy has been successfully used to quantify the complexity of many systems [77,78].Recently, the possibility of using quantum entropy and other quantum information theoretical measures has been explored by the community of network scientists.
For classical complex networks, von Neumann's entropy has been applied over one decade ago [79].The combinatorial Laplacian matrix L, obtained from the adjacency matrix representing the network, is rescaled by the number of edges in the network.The normalization of the matrix L guarantees that the corresponding eigenvalues are non-negative and sum up to 1-in order to be interpreted as probabilities [80]-and some other properties which makes the resulting object similar to a quantum density matrix ρ.Network entropy is defined according to von Neumann quantum entropy as By exploiting the eigen-decomposition of the Laplacian matrix, it can be shown that this entropy corresponds to the Shannon entropy of the eigenvalue spectrum of ρ.This entropy has been generalized to the case of multilayer systems [81], composite networks where units exhibit different types of relationships that are generally modeled different layers (see [17,18,82] for a thorough review).It has been recently shown that the von Neumann entropy calculated from the rescaled Laplacian does not satisfy the sub-additivity property in some circumstances [13,40].This undesirable feature can be addressed by means of a more grounded definition [13], whose rationale is to measure the entropy of a network by exploiting how information diffuses through its topology (see Fig. 4A-B).Information diffusion in this context is governed by the equation ψi with ψ i (t) the amount of information in node i at time t.The solution of this diffusion equation is given, in vector notation, by ψ(t) = exp(−Lt)ψ(0), whose normalized propagator is used to define the density matrix as where time plays the role of a resolution parameter allowing one to probe entropy at different scales [13].A similar approach, involving a modified Laplacian matrix, has been recently used for revealing the mesoscale structure of complex directed networks [83,84].This quantum-inspired framework provides a powerful basis to develop an information theory of complex networks, with direct applications in classical network science, such as system comparison.

Comparing classical networks.
A known problem in network science is to compare two networks, without relying on a specific subset of indicators.Network information entropy allows one to introduce relative entropies such as the Kullback-Leibler divergence, to compare two networks with density matrices ρ and σ respectively: By exploiting the well-known classical result that the minimization of Kullback-Leibler divergence between a reference distribution and its parametric model corresponds to the maximization of the likelihood, it has been shown that in a network context this allows one to define the network log-likelihood by The introduction of network likelihood (see Fig. 4C) opens the door to a variety of applications in statistical inference and model selection, based on concepts such as the Fisher information matrix, Akaike and Bayesian information criteria, and minimum description length, to cite some of them [13].This new framework has been used to compare networks for several purposes.For instance, in the case of pairs of networks, the graphs are first merged by connecting each node from one network to any other node in the other network.Successively, continuous-time quantum walks are used to explore the composite system and the quantum Jensen-Shannon divergence between the evolution of two walks is calculated.This divergence, that is a measure of (dis)similarity, is shown to be maximum when the two original networks are isomorphic [57].
The square root of quantum Jensen-Shannon divergence has the nice property of defining a metric, allowing one to define a distance between networks.If ρ and σ are two density matrices corresponding to two networks with N nodes, their Jensen-Shannon divergence is defined by that is the difference between the entropy of the mixture µ = 1 2 (ρ + σ) and the semi-sum of the entropies of the original systems.In the context of multilayer systems, this measure has been used to quantify the distance between layers of a multiplex network, cluster and aggregate them appropriately in order to reduce its structural complexity [40].
These ideas have quickly found direct applications in biology.In genetic molecular systems, such as the ones described by gene-protein interactions, layers might encode different relationships among molecules-functional, e.g.additive, suppressive and other types of association, or physical, e.g.co-localization or direct interaction.The  [13].The propagator governing the dynamics of how information flows through a network, represented by a Laplacian matrix L, is used to define a "density matrix" ρ = e −Lτ /Tr(e −Lτ ) for the calculation of network entropy (A).It has been shown that this entropy has desirable properties, including sub-additivity (B) and can be used to devise an inference framework-based on the concept of network likelihood-to fit the model's parameters by Kullback-Leibler divergence minimization or likelihood maximization.An example for a Watts-Strogatz network is shown in (C), where the number of neighbors K and the rewiring probability prew generate likelihood regions (color coded) to fit a small-world network with K = 6 and prew = 0.2.The Jensen-Shannon distance, defined to compare networks, has been used to identify redundant information in multiplex networks [18,81] and reduce their structural complexity in biological, social and transportation systems [40].In (D) it is used to identify community-based associations in the human microbiome, where each site is represented by a layer in a multiplex network and revealed community of layers [13] are in perfect agreement with biological expectation (black squares) [85].
information-theoretic framework described here allowed to show that such systems exhibit a certain level of redundancy, larger than the one observed in man-made systems [40], suggesting the existence of biological mechanisms devoted to maximize diversity of interactions.
In computational neuroscience studies, the connectome of the nematode C. elegans-one of the most studied in the field because of its small size, with approximately 300 neuronal cells-has been mapped to a multiplex network where layers encode synaptic, gap junction, and neuromodulator interactions.Here, the analysis of reducibility revealed that the monoamine networks have a unique structure, with information complementary to that provided by neuropeptide networks [86].The same analysis, applied to a multilayer functional representation of the human brain revealed, quantitatively, the importance of not disregarding or aggregating connectivity information for clinical classification of healthy and schizophrenic subjects [87].
In another application, whose result is shown in Fig. 4D, the Jensen-Shannon distance between layers of a multiplex system has been used to identify community-based associations in the human microbiome [13], where the microbial network corresponding to each body site is represented by a layer in a multiplex network, in perfect agreement with biological expectation [85].
The (dis)similarity between networks has also been quantified by using a combination of the classical Jensen-Shannon distance and the concept of network node dispersion, measuring the heterogeneity of a graph in terms of connectivity distance among its nodes [41].
Degree distribution of quantum networks.Nodes in a complex network have different roles and their influence on system dynamics can vary widely depending on their topological characteristics.One of the simpler (and widely applied) characteristics is the degree centrality, defined as the number of edges incident on that node.Many real world networks have been found to follow a widely heterogeneous distribution of degree values [88].Several models, based on mechanisms like preferential attachment [15], fitness [89] or constrained random wiring [53], to mention some of them, were developed to reproduce degree distributions commonly observed in empirical systems.Despite the complexity of the linking pattern, the degree distribution of a network affects in a simple way the ongoing dynamics.In fact, it can be shown that the probability of finding a memoryless random walker-at a given node of a symmetric network at the stationary state-is just proportional to the degree of such node [65].In [7] the authors consider the relationship between the stochastic and the quantum version of such processes, with the ultimate goal of shedding light on the meaning of degree centrality in the case of quantum networks.They consider a stochastic evolution governed by the Laplacian matrix L S = LD −1 , the stochastic generator that characterizes classical random walk dynamics and leads to an occupation probability proportional to node degree.In the quantum version, an hermitian generator is required and the authors proposed the symmetric Laplacian matrix , generating a valid quantum walk that, however, does not lead to a stationary state, making difficult a direct comparison between classical and quantum versions of the dynamics.A common and useful workaround to this issue is to average the occupation probability over time (4).
The generators of the two dynamics are spectrally similar (see Fig. 3) and share the same eigenvalues, while the eigenvectors are related by the transformation φ As a consequence, if the system is in the ground state the average probability to find the walker on a node will be the same as in the classic case, which will depend solely on the degree of each node.For the cases in which the system is not in the ground state, it is possible to define a quantumness measure describing how far from the classical case the probability distribution of the quantum walker will be.In the case of uniformly distributed initial state ρ 0 , this provides a measure for the heterogeneity of the degree distribution of a quantum network.
Mesoscale organization of quantum systems.Community detection, and in general mesoscopic structure detection, has been widely studied in the literature of classical complex networks [90,91].While the definition of community "a subset of nodes tightly connected compared to what is expected" is in general ill defined and part of an ongoing debate, the number of proposed algorithms is incredibly high and still growing.
The cross-pollination of community detection with quantum mechanics is in two levels.On the one hand, chronologically, the first attempt was to borrow tools from quantum mechanics for applications to classical systems [83,[92][93][94][95]. On the other hand an algorithm to find communities in complex quantum systems was proposed in [76].
In [92,93] the authors propose a method for data clustering similar to kernel density estimators, in a quantum framework.The given data points are mapped to a Gaus-sian wave function and, supposing that the latter is an eigenstate for some time-independent Schrödingher equation: the minimization of the potential V (x) leads to the desired clustering.An extension to dynamical quantum systems has been introduced in [94].In this case the expectation values of the position operator evolves in time toward the closer minimum of the potential.This formulation can leverage the acceleration of graphics hardware.
A method based on continuous-time quantum walks was proposed in [95].Here a node affinity measure based on the response of node population density to link failure was given.If the population on two nodes changes in a similar manner after link removal, they are more likely to belong to the same community.
A magnetic Laplacian, where a magnetic field is expected to traverse all cycles in the network, was used in [83].With an approach similar to chiral walks, previously described, the symmetric Laplacian is amended with the original link directionality by a phase term e ±iθ , with θ being a parameter for the method, and used for community detection in directed networks, a longstanding problem in network science.
In the case of quantum systems, partitioning in modular units has been often carried out on the basis of ad hoc considerations.In an effort to extend community detection to the quantum mechanics realm, Faccin et al. [76] introduced several closeness matrices inspired by different quantum quantities.Given the Hamiltonian H = ij H ij j| of the quantum system of interest, the authors consider a continuous-time random walk on the system topology.The first quantity is energy transport, porting to the quantum realm the concept applied in several classical algorithms where communities are interpreted as traps for the dynamical process.In this framework, two nodes are considered to be close if, on average, their in-between transport is high.If this average is computed over a short time period (compared to evolution time scales), then the closeness values are proportional to the Hamiltonian terms |H ij |, providing a classical approach to community detection (see Fig. 5).A second quantity, also proposed as a closeness measure, is related to the average fidelity of the evolving process compared to the initial state.In this case the localization of eigenstates is the characteristic determining the closeness of two nodes.These methods augment current ad hoc approaches to partitioning nodes in quantum transport systems with enhanced methods based on community detection algorithms.

IV. OUTLOOK IN QUANTUM NETWORK SCIENCE
Generalization of complex network methods to the quantum setting represents a foundational advancement required to understand complexity in physical systems.These methods represent a change of paradigm which bring several road blocks that must be faced.Centrally, the application domain of complex network methods to quantum physics must be expanded, whereas studies in the other direction, i.e.where methods from the quantum domain have now been ported to network science.
From a foundational perspective, as networks necessarily represent physical systems, such systems are inherently governed by the laws of information physics.In fact, a research line is emerging now that seeks to quantify, in terms of implicit information processing capacity, networked systems, with several applications to social, technological and biological systems [13,40,41,86,87].Although this interesting direction seems promising, yet it is comparably in its infancy, whereas it is still not known how to generalize classical concepts of complexity science to the quantum domain.
Another relevant research direction, crucial for applications in classical network science, concerns the interplay between structure and dynamics, which is almost entirely unclear in the case of quantum networks.Although scalefree networks have been considered in the quantum setting [32], the result is an-albeit interesting-toy model with theoretical predictions to be verified experimentally.Therefore, further advancement along this track is of central interest, because it might play a fundamental role in quantum enhanced technology and could lead to experiments devoted to test cross-disciplinary ideas in quantum and complexity science [33].
The quest for a theoretical foundation for quantum complex networks might have a deep impact in information and communication technology.While information processing in classical systems is well controlled, it is also rather limited and quantum computing might overcome such limitations [96,97].However, given that such systems are more sensitive to interactions with the environment, they are also more exposed to errors than their classical counterparts.Quantum error-correcting codes allow us to store and manipulate quantum information in the presence of certain types of noise that, in this context, might perturb the quantum system causing effects similar to random failures in classical complex networks.The development of quantum error correction techniques that make quantum computing and quantum communication possible can not prescind from the study of 'system resilience', a topic that found uncountable applications in classical network science [26,28,31,98,99].Other types of perturbations that are natural for classical systems, such as targeted attacks of network hubs [26] or cascade-based attacks [100], still have no clear quantum counterpart and their study, from both theoretical and experimental perspectives, will play a key role in the development of a quantum Internet [5].In fact, it is tantalizing to think about how quantum hubs should be protected by the quantum counterpart of typical denial of service attacks.
Continued advances in the theory of complexity in net-worked quantum systems will help address the challenges faced as quantum technologies scale up to commercially feasible products.Work towards a quantum theory of complex networked systems is already opening up novel avenues when facing contemporary complexity challenges.

FIG. 4 .
FIG.4.Network information theory inspired by quantum thermodynamics and computing[13].The propagator governing the dynamics of how information flows through a network, represented by a Laplacian matrix L, is used to define a "density matrix" ρ = e −Lτ /Tr(e −Lτ ) for the calculation of network entropy (A).It has been shown that this entropy has desirable properties, including sub-additivity (B) and can be used to devise an inference framework-based on the concept of network likelihood-to fit the model's parameters by Kullback-Leibler divergence minimization or likelihood maximization.An example for a Watts-Strogatz network is shown in (C), where the number of neighbors K and the rewiring probability prew generate likelihood regions (color coded) to fit a small-world network with K = 6 and prew = 0.2.The Jensen-Shannon distance, defined to compare networks, has been used to identify redundant information in multiplex networks[18,81] and reduce their structural complexity in biological, social and transportation systems[40].In (D) it is used to identify community-based associations in the human microbiome, where each site is represented by a layer in a multiplex network and revealed community of layers[13] are in perfect agreement with biological expectation (black squares)[85].

FIG. 5 .
FIG.5.Community Detection in a light harvesting network (LHCII)[76].(top left) Monomeric subunit of the LHCII complex with pigments Chl-a (red) and Chl-b (green) packed in the protein matrix (gray).(top center) Schematic representation of Chl-a and Chl-b in the monomeric subunit, here the labeling follows the usual nomenclature (b601, a602. . .). (top right) Network representation of the pigments in circular layout, colors represent the typical partitioning of the pigments into communities.The widths of the links represent the strength of the couplings |Hij| between nodes.Here the labels maintain only the ordering (b601→1, a602→2,. . .). (a,b,c) Partition of LHCII applying the quantum community detection algorithm in[76].Transport (for long times) and Fidelity approaches give similar results while short time transport is closer to classical community detection.Link width denotes the pairwise closeness of the nodes.