Potential energy of complex networks: a quantum mechanical perspective

We propose a characterization of complex networks, based on the potential of an associated Schrödinger equation. The potential is designed so that the energy spectrum of the Schrödinger equation coincides with the graph spectrum of the normalized Laplacian. Crucial information is retained in the reconstructed potential, which provides a compact representation of the properties of the network structure. The median potential over several random network realizations, which we call ensemble potential, is fitted via a Landau-like function, and its length scale is found to diverge as the critical connection probability is approached from above. The ruggedness of the ensemble potential profile is quantified by using the Higuchi fractal dimension, which displays a maximum at the critical connection probability. This demonstrates that this technique can be successfully employed in the study of random networks, as an alternative indicator of the percolation phase transition. We apply the proposed approach to the investigation of real-world networks describing infrastructures (US power grid). Curiously, although no notion of phase transition can be given for such networks, the fractality of the ensemble potential displays signatures of criticality. We also show that standard techniques (such as the scaling features of the largest connected component) do not detect any signature or remnant of criticality.

related to the emergence of a giant component, at a critical value of the connection probability between pairs of nodes 29 . We will show that our analysis tools, based on the reconstructed potentials, are able to capture the singular behavior of the network close to the transition. Specifically, three indicators of such criticality will be identified: the depth, the width and the Higuchi Fractal dimension (HFD) 30 of the pointwise median potential, computed on several realizations of the ER network with the same size and connection probability. According to the mathematical tool one adopts, many strategies are possible in order to characterize a potential that represents the whole ensemble (over many realizations of the network). We shall opt for the pointwise median of the potentials over different network realizations, since it is best suited to unearth the sought criticality features. As we shall discuss in the following, other mathematical averaging procedures can be employed, but turn out to be less effective. We shall henceforth refer to the pointwise median potential over many network realizations at fixed connection probability as the ensemble potential. Finally, to check the validity of the proposed approach in real-world systems, we shall reconstruct potentials from the graph spectra of publicly available complex networks describing infrastructures (US power grid) 31 . These are real networks, for which no notion of phase transition can be defined. However, interestingly, we shall find that the fractality of the ensemble potential displays signatures of criticality.
The content of this Article is organized as follows: in "Setting up the problem: from graph spectra to reconstructed potentials" section we present the properties of the Laplacian spectrum and discuss the method for reconstructing the associated potentials via dressing transformations; in "Erdös-Rényi phase transition in the reconstructed potential framework" section we examine the description of the ER network provided by the reconstructed potential framework, focusing on the critical behavior of its length scale, depth and HFD at the phase transition; in "Analysis of a real-world network" section we use this approach to investigate a real-world network, the US power grid. The details of the dressing transformation for the potential reconstruction, some additional properties of reconstructed potential, as well as a comparison among different averaging procedures are presented in the Supplementary Information.

Setting up the problem: from graph spectra to reconstructed potentials
Algebraic graph theory is a branch of graph theory devoted to the investigation of graph connectivity properties using results and methods from algebra. In this respect, one of the most interesting results concerns the possibility of infering a connectivity measure by looking at particular eigenvalues of the graph Laplacian 32 . However, more general properties of a graph could be unearthed by examining all the eigenvalues and the graph energy, defined as their sum 33,34 . This approach could result in an alternative perspective on the problem of measuring the robustness of real-world complex networks. In the attempt to give a compact representation of the information contained in the whole graph spectrum, we apply a methodology, based on non-linear equations, to retrieve a one-dimensional potential given a set of energy levels (see Fig. 1 for a schematic overview).
Accordingly, starting from a generic (unweighted and undirected) network, we first compute its Laplacian and eigenvalues, and then reconstruct the unique Schrödinger potential associated with the spectrum. The potential profile will yield a snapshot of the network configuration, providing a quantum mechanical perspective on the study of its connectivity.
Laplacian spectrum. We introduce here the formalism that is necessary to analyze complex networks in the framework of graph theory 32 and recall a few basic notions that will be useful in our analysis. A graph G = (N, E ) is defined through a set N of N nodes and a set E of edges connecting them. The adjacency matrix A of G is a matrix whose elements a ij are nonvanishing only if a connection between node i and node j exists. In general, graphs can be built by assigning a weight and an orientation to each edge. In the present work, we shall focus on undirected and unweighted graphs, whose adjacency matrices are symmetric and binary (i.e., consisting only of 0 and 1 elements). Moreover, for the sake of simplicity, the networks considered in this work will not include loops, namely links connecting a node to itself. For each graph G , given its adjacency matrix A = {a ij } , the number of connections of each node, namely the node degree d i , is simply calculated by summing column-or row-wise the adjacency matrix. Accordingly, one defines the degree matrix D as the diagonal matrix with D ii = d i . The Laplacian L = {L ij } of G is defined as the difference D − A , so that:

Figure 1.
Conceptual workflow of the proposed methodology: starting from a given network, one first computes the Laplacian and its spectrum; the potential of a 1D Schrödinger equation, yielding the given set of eigenvalues, is then obtained and used to characterize the whole network and investigate its connectivity.
Scientific Reports | (2020) 10:18387 | https://doi.org/10.1038/s41598-020-75147-w www.nature.com/scientificreports/ The normalized Laplacian L = {L ij } of G is defined as the matrix with elements which can be expressed as L = D −1/2 LD 1/2 , with the convention D −1/2 i,i = 0 for d i = 0 , namely if the node associated with index i is isolated. If the complex network features components that are disconnected from each other, both the L and L matrices can be recast into block-diagonal forms, with each block corresponding to a specific component. The spectrum of the normalized Laplacian, which is also called the spectrum of the graph, provides comprehensive information on the structure of the network, with regard to the number of its connected components and their dimensions. The eigenvalues 1 ≤ 2 ≤ · · · ≤ N of L satisfy 0 ≤ i ≤ 2 . In particular, 1 = 0 is always an eigenvalue, whose multiplicity coincides with the number of connected components in the network. The eigenvalue 2 , called algebraic connectivity or Fiedler eigenvalue, is therefore nonvanishing if and only if the network is connected 10 . To give an idea of the spectra of L corresponding to peculiar structures, it is worth reviewing a few notable cases which will be relevant for our analysis 35 : If the network includes components K n , P n , C n with n < N , the eigenvalues of the n × n Laplacians associated with such components will contribute to the spectrum of the whole network Laplacian (2).

Reconstructing potentials through dressing transformations.
We now associate to a given network a potential and a 1D Schrödinger equation, whose energy levels coincide with the spectrum of the normalized Laplacian of the network. For this purpose, we shall apply a method based on the dressing transformation, proposed in Ref. 20 and employed e.g. in Refs. 19,23 . For a given network, we compute the eigenvalues 1 ≤ 2 ≤ · · · ≤ N of the normalized Laplacian L , and consider the shifted spectrum whose values are in [−2, 0] . Hence, we derive the potential V(x) such that all the N < N nonvanishing shifted eigenvalues are energy levels of the Schrödinger equation ( ℏ = 1, m = 1/2) with ψ(x) a normalizable wavefunction. We obtain the potential V by applying a method, whose details are outlined in the Supplementary Information, that is based on the iteration of a two-step procedure: we first solve the Riccati equation and then update the potential according to The iteration starts at n = N , with V N+1 = 0 and proceeds until n = 1 , yielding the final potential V 1 that can be identified with V in Eq. (4). Notice that the potentials obtained by the above procedure are even by construction.
The sequence of equations (5) can be integrated with any standard numerical method. Following the approach proposed in Ref. 23 , we adopted a fourth-order Runge-Kutta (RK) method with step size h = 10 −5 to guarantee a high reconstruction accuracy. It is worth remarking that the dressing transformation is not the only technique to solve the inverse problem of finding the potential corresponding to a given set of levels. Another possible choice is the variational method described in Ref. 22 , which however requires a search in a very large parameter space. Additional details are given in the Supplementary Information.

Erdös-Rényi phase transition in the reconstructed potential framework
A random network is a set of N nodes, with edges randomly placed among them. The original model of a random network was introduced by Erdös and Rényi (ER) in their seminal work 24 , where they defined an ensemble G(N, L) of graphs having a fixed number L of randomly placed edges. Here we will focus on a different formulation, proposed by Gilbert 36 and characterized by the fact that the number of edges is not fixed. This model is based on an ensemble of networks G(N, p), in which each node pair is connected with probability p. A parameter to quantify the connectivity for networks in this ensemble is the average degree Erdös and Rényi predicted the existence of a phase transition (in the thermodynamic limit) at �k� = 1 , corresponding to the critical probability p c = 1/(N − 1) . This transition consists in the appearance of percolation and Reconstructed potentials and connection probabilities. We consider different realizations of the ER complex network, with N = 500 nodes and different values of the connection probability p. To get a qualitative understanding of the relation between the network structure, the graph spectrum and the reconstructed potential, let us first look at three particular cases: p = 10 −4 , p = 0.9 and the critical probability The statistical ensemble for p = 10 −4 is made of almost disconnected networks. In each realization, a few pairs of connected nodes are present, and connected components with more than two nodes are extremely rare. Thus, the typical spectrum n ( n = 1, . . . , N ) of L consists of the highly degenerate eigenvalue 0, whose multiplicity coincides with the number of connected components, and the eigenvalue 2, related to the presence of paths P 2 of two nodes (see the examples in "Laplacian spectrum" section), which are on average pN(N − 1)/2 ≃ 12.5 . These eigenvalues correspond to the values −2 and 0 of the shifted spectrum (3), respectively. This structure of the spectrum tends to appear in the vast majority of realizations, one of which is displayed in Fig. 2. We observe that the number of wells in the reconstructed potential tends to coincide with the ground state degeneracy, which, in this case, is equal to 488, while the excited state has degeneracy 12. The similarity of the reconstructed potentials is reflected in the regularity of the ensemble potential V m (x) (computed as the pointwise median over M = 100 network realizations), shown in the right panel of Fig. 2.
For p = 0.9 , in basically all network realizations the graph consists of a single component, with a high density of links. In such configuration, the spectrum n of L consists of a nondegenerate eigenvalue 0, separated by a gap from the other eigenvalues, that concentrate around 1. Indeed, when the connection probability is close to p = 1 , the spectrum approximates that of the complete graph K N . Also in this case the structure of the spectrum is weakly dependent on the specific network realization, a feature which is again reflected in the similarity of the reconstructed potentials. The typical potential profile is characterized by a single minimum around x = 0 , and rapidly increases approaching an almost constant value. The statistical variability of the reconstructed potentials is extremely low, leading to a very smooth median V m (x) (computed over M = 100 network realizations), manifest in Fig. 3.
At the critical probability p = p c = 1/(N − 1) , (almost) each realization of the network is characterized by the emergence of a giant component, involving a number of nodes of order N 2/324, 29,36 . The remaining nodes are either connected into smaller size components or isolated. The eigenvalues n of the graph spectrum are distributed in the full interval [0, 2], with a larger concentration around the endpoints. This behavior is mainly due to the presence of the giant component, which can be approximately described as a path, that contains also small ramifications and cycles. For this reason the contribution of the giant component to the graph spectrum is qualitatively similar to the spectrum of a path P n with n ∼ N 2/3 . Although the giant component displays general (and common) features in all the network realizations at critical p, the details of its nontrivial structure can hardly be reproduced. Therefore, since giant components corresponding to different realizations are generally characterized by very different patterns and micro-structures, and even different dimensions, the shifted eigenvalues E n (positions of the associated spectral lines in the central panel of Fig. 4) will be wildly fluctuating. Smaller components and isolated nodes in the network will contribute to the spectrum of L with sparse and degenerate eigenvalues, as in the case of low p. In Fig. 4 we show the reconstructed potential for a single realization, whose eigenvalues E n correspond to the shifted spectrum (3) lying between −2 and 0. The shape of this potential, characterized by irregular oscillations around a constant value, with no appreciable increase in the considered x range, is actually www.nature.com/scientificreports/ rather similar for all realizations at critical p. However, the features of these oscillations wildly differ for different realizations, since they are subject to the same variability that characterizes the spectrum associated with the giant component. As a result, the profile of the median V m (x) is very irregular. This was observed to be true also for a surprisingly small number of realizations. This observation is central: we shall argue that the irregularity of the average potential at criticality is but a manifestation of the emergence of fractality.
Length scale of the reconstructed potential. The ensemble potentials V m (x) obtained from M = 100 realizations of the ER network with N = 500 nodes, for different values of p, are shown in Fig. 5. One observes a very different behavior, depending on whether p is smaller or larger than p c . For p < p c , V m features irregular oscillations superposed to an increasing trend, that becomes less and less steep as criticality is approached. For p > p c , the curves become smoother (differentiable), and a minimum in the origin gradually appears, followed by a rise to an almost constant plateau. As we shall see, this behavior is well approximated by a Landau potential V (x) = a sech 2 (x/b) , with −a and b denoting the depth and width of the trapping potential, respectively 40 . The reason for the dependence of V m (x) on p(> p c ) becomes evident when one compares the shifted Laplacian spectra in the central panels of Figs. 3 and 4: by approaching the percolation transition from above, the reconstructed potentials tend to become steeper, in order to accommodate a larger number of nondegenerate eigenvalues in the upper part of the spectrum. This observation suggests a fit of the ensemble potential V m (x) with an expression that interpolates between the large-p and critical-p regimes: we consider an ansatz of the following type in the attempt to fit the slower increase of V m (x) towards its asymptotic (large-x) value in the vicinity of the phase transition. Figure 6 displays the comparison between V m (x) and V fit (x, p) at different values of p > p c = 2 × 10 −3 ,  www.nature.com/scientificreports/ and the dependence of the fit parameters a and b on the connection probability p (at N = 500 ). The relative errors on the fit parameters a and b are always smaller than 0.02% for p ≤ 5 × 10 −2 , and increase as the critical point is approached, reaching 1.16% for a and 1.65% for b at p = 3 × 10 −3 , where the oscillations of the ensemble potential profile reduce the fit accuracy. The ensemble potential captures the critical behavior of the system, as it emerges from the divergence of the fit parameters as p approaches p c from above. In particular, the scales a and b, which set respectively the depth and width of the trapping potential, diverge like at the percolation transition, with critical exponents α ≃ 1.60 ± 0.22 and β ≃ 1.35 ± 0.13 . This shows that the very structure of the potential of the 1D Schrödinger equation, as well as the parameters that characterize it, detect the percolation phase transition of the associated network. A few additional properties of the potential V fit are discussed in the Supplementary Information. In the next subsection, we shall explore in detail the association between the shape of the ensemble potential and the onset of the percolation phase transition in the collection of ER networks from which the pointwise median is computed. Actually, the roughness of the V m profile represents a signature of fractality, a property that emerges around p c as a result of the median operation, but is not detected in single-realization reconstructed potentials, characterized instead by a smooth profile at all connection probabilities.
Before moving on to the quantitative description of the fractality of the median potential V m , it is worth observing that a possible alternative approach to the study of the ER network ensemble would be based on analysing the pointwise mean potential V µ . However, as we detail in the Supplementary Information, although V µ follows qualitatively similar trends to those of V m , depicted in Fig. 5, its profile remains smooth, even close to the percolation phase transition, and shows no sign of fractality.

Transition characterization through Higuchi Fractal Dimension.
We observed (Fig. 5) that the ensemble potential profile is less smooth at the critical probability p c than in the low-and high-p regimes, due to the spectral fluctuations of the associated networks. Figure 7 shows, for different values of N, the ensemble potentials over M = 100 realizations of the network, always at the critical connection probabilities p c = 1/(N − 1) ; it can be noticed that the increasing trend in x of the ensemble potential is more manifest for networks of smaller size N. At N = 500 , one observes oscillations but no overall increasing trend with x. We now introduce a measure to quantify the lack of smoothness of V m . In Ref. 19 , the ruggedness of the potentials reconstructed from given sets of energy levels was characterized in terms of a qualitative measure of fractal dimension. Here, we shall use the Higuchi Fractal Dimension (HFD) 30 , that, having been originally introduced to characterize time series, is especially suited to describe the profile of functions of one variable. From a given sequence {F 1 , F 2 , . . . , F ν } one extracts the subsets of values corresponding to indices separated by a positive integer k < ν and starting from i. The quantities represent properly normalized measures of the mean distance between neighboring values in (10). The terms L i (k) corresponding to the same spacing k can then be averaged over all possible initial points to obtain If �L(k)� ∼ k −D , the exponent D is called the HFD of the sequence {F 1 , F 2 , . . . , F ν } . In practice, the above dependence holds only in a certain range of k; in our analysis, we have considered 2 ≤ k ≤ 800. Figure 8 shows the HFD values of the ensemble potentials computed on M realizations of ER networks with N = 500 nodes at fixed p, with M = 100 (see Fig. 5) and M = 1000 . From both plots, we observe that the fractal dimension peaks around the critical probability, while it drops to a value close to one in the low-and high-p regimes. The inset in Fig. 8 displays a (close) zoom of the peaks around p = p c = 1/ (N − 1) , showing that the maximum value is attained at the critical probability, for the explored values of p and N.
The maximum of the HFD can be therefore be considered as an indicator of criticality. We observe that the plotted values at each connection probability p depend on the specific ensemble of M randomly sampled realizations of the network. We actually checked that plots from independent sets, each containing M = 100 realizations, generally fluctuate around a mean curve. However, despite this variability, a general trend emerges, characterized by the presence of a peak around the critical probability, in all the considered sets of M = 100 realizations.
Finally, a comparison between the two plots displayed in Fig. 8 reveals that the values of HFD obtained for different numbers M of realizations are sensibly different when the connection probability approaches the critical value. This result motivates us to systematically investigate the behavior of the fractal dimensions at p = p c = 1/(N − 1) for networks with varying numbers of nodes N and realizations M concurring to the (10) . In all the plots the HFD displays an initial increase, due to the fact that fractality emerges when the ensemble potential is computed over a relevant number of realizations, while the single potentials are not fractal (HFD≃ 1 ). The position of the maximum is largely independent of M, and the HFD decreases for large M. We observe that such a decrease becomes less significant for larger values of N. It would be tempting to assume that, as N increases further, the HFD reaches a plateau after the maximum, as a function of M. However, our data do not enable us to safely draw this conclusion (even by a tentative analysis of finite-size scaling), nor to determine the precise value of the maximum for N, M → ∞ . Additional details on the numerical procedure are outlined in the Supplementary Information.

Analysis of a real-world network
The formalism based on the reconstructed potentials and the HFD, developed in the previous sections and used to characterize ensembles of ideal ER networks, will now be applied to the description of a real-world situation.
It is important to stress that a real-world network is not a random graph in the sense of ER, and does not enjoy the same idealized features. It is not obvious, for example, that a real-world network percolates, certainly not in the same way as a random graph. For instance, the very presence of a unique giant component and its scaling features cannot be taken for granted. We shall now look at a specific example and consider a number of typical quantities that characterize it. The data analyzed in this section are taken from a public database 31 and represent the power grid of the Western States of the USA 41 ). This system can be modeled as an undirected unweighted complex network with N = 4941 nodes and L = 6594 edges; a node indicates a generator or a transformator or a substation, and each edge represents a power supply line. We shall find it convenient to work here with the average degree k defined in Eq. (7) (and we recall that the percolation phase transition of an ER random network would take place at �k� = 1 ). The network, displayed in the top left panel of Fig. 10, is made up of a single connected component and is characterized by an average degree �k� = 2.67 ; its fill, defined as the ratio between the number of edges L and the maximum number of edges N(N − 1)/2 in an undirected network without loops, takes the value 5.40 · 10 −4 .
The US power grid network is deterministic, with the edges corresponding to an organized structure, constrained by infrastructural requirements. In order to enable a comparison with a random graph, we shall introduce randomness in the system, by sampling out a subset of N ′ < N nodes to form an subgraph, in which the www.nature.com/scientificreports/ edges connecting these nodes are inherited from the original network. The sub-networks corresponding to different numbers N ′ of sub-sampled nodes are reported in Fig. 10. The properties of the sampled subnetworks will fluctuate according to the specific subset of N ′ nodes. We performed a statistical analysis in the intermediate range 370 ≤ N ′ ≤ 4113 , in which the variability of the sampled subnetwork is expected to be larger, as 0.20 ≤ �k� ≤ 2.22 would be close to criticality ( �k� = 1 ) if the network were random. At fixed N ′ , we generated M = 100 subgraph realizations and computed, for each of them, the average degree k , the size N ′ LCC of the largest connected component (LCC) and the reconstructed potential. The left panel of Fig. 11 displays the dependence of the LCC relative size N ′ LCC /N ′ vs the average degree k of such distributions, with the error bars representing their standard deviations over M = 100 subgraph realizations (at fixed N ′ ) (the inset shows for comparison the same quantity in the case of the ER network ensemble 1 ). The dependence of the size of the LCC on k appears to be smooth and no signature of a phase transition is observed: a fit of the LCC fraction for k > 1 , displayed in the right panel of Fig. 11, yields N ′ LCC /N ′ − (N ′ LCC /N ′ )| �k�=1 ∼ (�k� − 1) 7.0±0.1 (for comparison, for an ER network at criticality, one would get the exponent 1 1 ). We stress again that the graph is not random and a phase transition is not properly defined in this case. Although this is not central to our analysis, we also observe that the number of subnetworks scales like N N ′ , so that the 100 realizations are practically independent and very different from each other.
We now test the performance of the reconstructed potential formalism as a tool of network analysis. In Fig. 12 (left) we show the profile of the ensemble potential V m for some values of N ′ ; the corresponding mean  www.nature.com/scientificreports/ values of k are also reported to facilitate comparison with the analogous plots for the ER networks. Finally, we computed the HFD of the ensemble potentials, represented in Fig. 12 (right) as a function of the average degree in the sampled subnetworks. Notice that it is impossible to investigate the behavior of the HFD for larger values of k , as the US power grid network, displayed in the top left panel of Fig. 10, has an average degree �k� ≃ 2.67 , that represents an upper bound for our analysis. The range of k values here analyzed roughly corresponds to probabilities 4 × 10 −4 ≤ p ≤ 4.4 × 10 −3 in Fig. 8  (random ER networks). The HFD displays a clear increasing trend up to �k� ≃ 0.75 , and then reaches a plateau for larger values of k , characterized by irregular (a few percent) fluctuations. Notice that, as emphasized before, the real US power grid network is characterized by an average degree �k� ≃ 2.67 , a value that cannot be reached in our analysis. Comparison with a bona fide random ER network would lead us to expect that the plateau ends at k 3.5 ( p 7 × 10 −3 ). On a wider scale the plateau would then appear as a (sharp) maximum, see Fig. 8. The behavior of the HFD detects therefore the occurrence of a significant amount of fractality in the ensemble of sampled subnetworks for the corresponding values of N ′ and k . We can interpret this behavior as a remnant of the criticality of the corresponding ER graph. Observe that such conclusion does not emerge from the quantity displayed in Fig. 11. The fractality of the potential, detected by the HFD, appears therefore to be a better indicator of a structure that is reminiscent of the ER phase transition. It is remarkable that we are able to reach this conclusion although the sampled subnetworks originate from a given, real-world deterministic graph and cannot, strictly speaking, be considered "random", as those of an ER graph.
A few additional comments are in order. We evaluated the HFD for 2 ≤ k ≤ 800 , k being the scale index introduced in Eq. (12). This choice was motivated by consistency with the analysis of "Transition characterization through Higuchi Fractal Dimension" section, pertaining to an ER random network. An investigation with 2 ≤ k ≤ 35 would yield the same qualitative features, with a lower (about 30%) fractality index, but a very similar plateau, starting again at �k� ≃ 0.75.

Conclusions and outlook
We have proposed a quantum-inspired approach to investigate complex networks: by using the mathematical framework provided by dressing transformations, we have developed a technique to uniquely associate a Schrödinger-like potential to the graph spectrum of a given network. We have verified that such potential probes and detects nontrivial phenomena in complex network dynamics, such as the phase transition in the Erdös-Rényi model, related to the emergence of a giant component. In particular, such transition can be investigated by monitoring the length scale and the Higuchi Fractal Dimension of the ensemble potential computed from several network realizations having the same size and connection probability.
We have applied this technique to the study of a real-world network, showing that the fractality of the ensemble potential displays a behavior that is reminescent of criticality, although no bona fide phase transition can be properly defined in such a case. We also observed that standard techniques (such as the scaling feaures of the LCC) are unable to detect any signature or remnant of criticality. Future work will extend this analysis to more complex models of artificial networks, which include other formation mechanisms, such as link rewiring and preferential attachment.
The reconstructed potential provides a snapshot representation of the structure of a network, yielding information on its connectivity and on the number of disconnected components. We observe that fractality has been related to network topology and self-similarity [42][43][44] ; to establish a relation between such fractality, arising from connectivity, and the one outlined in our work could foster interesting developments. For example, this aspect could be considered when examining real-world networks, since the reconstructed potential and its fractality could be used to test their robustness, shed light about hidden geometries and even diagnose possible weaknesses.
For this reason, we plan to explore further the characterization of real-world networks through reconstructed potentials, trying in particular to understand whether the typical patterns found in this analysis are an intrinsic feature of the specified domains. Moreover, we will investigate the possibility to improve our framework by combining the reconstructed potentials with other approaches to complex networks, based on entropy 11,45,46 and machine learning 47 .

Data availability
The datasets generated and analysed during the current study are available from the corresponding author on reasonable request.

Scientific Reports
| (2020) 10:18387 | https://doi.org/10.1038/s41598-020-75147-w www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.