Spacing ratio statistics of multiplex directed networks
Tanu Raghav, Sarika Jalan

TL;DR
This paper investigates eigenvalue ratio statistics of multiplex directed networks using random matrix theory, revealing how multiplexing strength influences spectral properties and eigenvector delocalization.
Contribution
It introduces a numerical analysis of eigenvalue ratio statistics in multiplex directed networks represented by non-Hermitian and Hermitian matrices, highlighting the role of multiplexing strength.
Findings
Multiplexing strength governs eigenvalue ratio statistics.
Multiplexing strength influences eigenvector delocalization.
Results applicable to understanding dynamics in multilayer systems.
Abstract
Eigenvalues statistics of various many-body systems have been widely studied using the nearest neighbor spacing distribution under the random matrix theory framework. Here, we numerically analyze eigenvalue ratio statistics of multiplex networks consisting of directed Erdos-Renyi random networks layers represented as, first, weighted non-Hermitian random matrices and then weighted Hermitian random matrices. We report that the multiplexing strength rules the behavior of average spacing ratio statistics for multiplexing networks represented by the non-Hermitian and Hermitian matrices, respectively. Additionally, for both these representations of the directed multiplex networks, the multiplexing strength appears as a guiding parameter for the eigenvector delocalization of the entire system. These results could be important for driving dynamical processes in several real-world multilayer…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsSpectroscopy and Quantum Chemical Studies · Molecular spectroscopy and chirality · Complex Network Analysis Techniques
Spacing ratio statistics of multiplex directed networks
Tanu Raghav, and Sarika Jalan
Complex Systems Lab, Indian Institute of Technology Indore - Simrol, Indore - 453552, India
Abstract
Eigenvalues statistics of various many-body systems have been widely studied using the nearest neighbor spacing distribution under the random matrix theory framework. Here, we numerically analyze eigenvalue ratio statistics of multiplex networks consisting of directed Erdős-Rényi random networks layers represented as, first, weighted non-Hermitian random matrices and then weighted Hermitian random matrices. We report that the multiplexing strength rules the behavior of average spacing ratio statistics for multiplexing networks represented by the non-Hermitian and Hermitian matrices, respectively. Additionally, for both these representations of the directed multiplex networks, the multiplexing strength appears as a guiding parameter for the eigenvector delocalization of the entire system. These results could be important for driving dynamical processes in several real-world multilayer systems, particularly, understanding the significance of multiplexing in comprehending network properties.
Keywords: Eigenvalues, RMT, multiplex network
1 Introduction
Random matrix theory (RMT) has been successful in finding universality among spectral correlations of enormous class of many-body complex systems. Under this theory, instead of investigating a system’s Hamiltonian, an ensemble of random Hamiltonian matrices is investigated to outline the system specific properties [1]. RMT has been useful in many areas of physics and mathematics ranging from condensed matter to economical financial markets by successfully describing properties of spectral fluctuation of complex nuclei, atoms, and complex molecules [2]. One of the most commonly used measures in RMT is nearest neighbor spacing distribution (NNSD) [3]. NNSD corresponds to the probability density of the differences () between two consecutive levels, say with i=. For quantum systems with regular classical counterparts, NNSD coincides with the Poisson distribution , whereas for the Gaussian orthogonal ensemble (GOE) and Gaussian unitary ensemble (GUE), NNSD can be well approximated by the Wigner distribution [3] and , respectively.
To calculate NNSD, one has to eliminate the impact of level density by a process called unfolding so that the mean level spacing becomes one [4]. However, due to the unavailability of an analytical form of the mean level density, the spectrum is unfolded numerically. Unfolding the spectra can be ambiguous as it requires polynomial curve fitting by choosing the fitting parameter arbitrarily [5]. To avoid this perplexity, a measure has been introduced called the eigenvalue spacing ratio, defined as the ratio between the nearest and next nearest neighbor spacing of the eigenvalues [6]. Since average mean spacing becomes insignificant in calculating spacing ratios, no unfolding is required. Over the years, the spacing ratio has been widely used to study the spectral statistics of several many-body systems, including complex networks [7, 8, 9].
Furthermore, complex systems are widely studied using network theory by recognizing that several real-world systems around us comprise of interacting units or nodes [10]. In recent years, network theory has been successful in providing insights into structure and dynamics of physical, biological, social, and several other natural systems [11]. Recently, the realization that various types of interactions may exist between the same set of nodes has paved the way for multilayer networks which not only incorporate different types of connectivity in individual layers but also provide a representation of interconnected systems [12]. These interconnected systems conceptualized as multilayer networks give new insights into several structural and dynamical properties compared to their monolayer counterparts [13, 14]. An elemental feature of representing multilayer networks is assessing the interconnectivity between various types of connections. Investigating such interconnectivity or inter-layer connections is essential for studying the role of dynamical switching between layers in a multiplex system. These inter-layer connections weights, referred to as multiplexing strength, provide a quantitative measure of impacts of dynamical and structural properties of one layer on those of the another layer [15, 16]. Various studies on multilayer networks have established the emerging aspects of the interplay between the intra and inter-layer connections of the associated network [17, 18]. In recent times, a new emerging paradigm in network science encapsulating higher order interactions has gained an immense population. The studies associated with higher order interactions under the framework of hypergraphs and simplicial complexes have unveiled interesting phenomena about complex networks [19].
Next, on a fundamental level, interactions are asymmetrical in nature for a large class of networks corresponding to real-world complex systems. In a network, a connection between two nodes often has a particular direction mathematically represented by directed networks and corresponding asymmetric adjacency matrices. Well-known examples include the food web network [20], gene regulatory network [21], citation [22], and world wide web networks [23]. Also, the spectra of adjacency matrix of a network [24] relate to various structural properties [25, 26, 30], as well as dynamical processes on the underlying networks [28, 29]. Nevertheless, spectral properties of directed networks [32, 33, 34] have not been investigated in as much detail as for undirected networks [35, 36]. Regardless of the significance of directed complex networks, spectral properties of directed networks are yet to be explored more to understand the system’s dynamical processes and their relation with the structural properties. Because of the complex eigenvalues of the asymmetrical adjacency matrix, the analytical approach can be challenging to analyze spectra of directed networks. Investigations of complex eigenvalues corresponding to directed networks and several other many-body systems have gained momentum in recent years. For example, J. Baron analyzed the spectra of directed complex networks owning heterogeneity in the network, and deduced effect of heterogeneity on the stability of corresponding system [37]. Recently, Metz and Neri analytically obtained an exact expressions for the inverse participation ratio by relating it with the structural properties of the directed networks. They showed that as network connectivity increases, the largest eigenvalue and the eigenvalue at the boundary of the bulk undergo from a localized to a delocalized regime [38].
Further, T. Peron and colleagues analyzed the spectra of directed single-layer networks using the eigenvalue spacing ratio (). By doing so, the authors set out that the spectral properties can be relevant to dynamical processes on the directed networks, and deduced the network parameters that can control dynamical transitions in the corresponding system [9]. In this paper, we extend the spacing ratio analysis of directed single-layer networks to multiplex networks consisting of two layers. One of the questions which we focus on is that for multiplexing networks what could be governing parameters for any such transition which is observed in single-layer networks. Both the layers are constructed using directed Erdős-Rényi (ER) random networks, represented by (i) weighted non-Hermitian random adjacency matrix and (ii) Hermitian random adjacency matrix. Notably, we use the average spacing ratio as a measure of complexity in the system as a function of different structural parameters of the corresponding networks. Later, is used to identify the delocalization transition of directed multiplex networks. We show that multiplexing strength not only can change the minimum probability required to commence the transition, but also can suppresses the transition.
The paper is organized as follows. Section. 2 comprises the model definition and construction along with the definition of technique employed in this paper i.e., eigenvalue spacing ratio for real and complex eigenvalues. Section. 3.1 and section. 3.2 discuss the results on the interplay of the ratio statistics and various structural parameters for Non-Hermitian and Hermitian representation of multiplex network. Section. 4 associates the localization-delocalization transition using Shannon entropy with the average ratio transition from the Poisson ensemble to GUE/RGE where RGE stands for real Ginibre ensemble. Section. 5 concludes the paper with future applications and directions.
2 Models and quantities
A multiplex network can be defined as where , being the set of nodes in the layer , . If , the edge is an intra-layer connection; that is, it is on the nodes of the same layer. If , the edge is an inter-layer edge between the nodes of layer and layer . There exists a special type of multilayer network referred as multiplex network, wherein each layer there exists same set of nodes, that is, . In such networks, inter-layer connections connect each node of one layer to its mirror node in the other layers. Here, we consider bi-layer multiplex networks, with each layer having nodes. The adjacency matrix for such a -layer multiplex network can be written as
[TABLE]
where and are the adjacency matrices corresponding to the layer and layer , respectively, and is the multiplexing or inter-layer strength. Degree of the -th node in any layer () is given by = where is the -th entry of matrix and corresponds to the inter-layer connection. Average degree of a network in layer is then =. The eigenvalues of an adjacency matrix of a multiplex network are denoted by , , where is the size of the multiplex network and .
2.1 Multiplex network with random weighted non-Hermitian matrix
In a multiplex network, we consider ER random networks for both the layers to generate a directed random network [31]. The ER network of size is constructed by connecting each pair of the nodes by a directed edge with a probability . Next, the edges are weighted with statistically independent random variables , drawn from a normal distribution (mean and variance). If there exists a directed edge from to , we assign a weight and for , the weight is (Fig. 2). We denote the weighted adjacency matrix of both the layers by and define (for layers 1, 2) as
[TABLE]
Because of the directed nature of the network () the matrix can be non-Hermitian. Such matrices are identified and have been defined as a diluted version of the real Ginibre ensemble (dRGE) [9]. Real Ginibre ensemble (RGE) comprises of random matrices formed from standard Gaussian variables. Note that, for (complete network), is drawn from the RGE distribution. Also, for , i.e., a completely disconnected network, becomes a diagonal random matrix and follows Poisson ensemble (PE) statistics. Therefore, for the adjacency matrix of multiplex network represented by , the block matrices and are independently distributed as representing layer and layer , respectively, whereas off-diagonal block matrices () correspond to bi-directional inter-layer connections where is the identity matrix connecting mirror nodes of layers and .
2.2 Multiplex network with random weighted Hermitian matrix
To study the Hermitian representation of the adjacency matrix of directed random graphs, we again consider ER random model for both the layers. Next, in the extracted adjacency matrix, the edges are weighted with statistically independent random variable (Fig. 3) drawn from normal distribution (zero mean and variance one). We denote the weighted adjacency matrix as for the layers and , and define as
[TABLE]
Above defined method has been recently introduced in References [39, 40]. Lately, the adjacency operator defined in Refs. [39, 40] has been used in various studies such as semi-directed graphs [41], cluster hypergraphs of semi-directed graphs [42], quantum walk [43] etc. Note that is constructed such that = . Here, as is increased, a transition is observed from the PE (=[math]) to fully connected real symmetric matrices (=). At =, the ensemble is similar to the Gaussian orthogonal ensemble (GOE) of random matrix theory (RMT), with the difference that; in the GOE, the diagonal matrix elements have twice the variance as the off-diagonal ones. Such a matrix has been defined as a magnetic adjacency matrix () [9]. Here, in the adjacency matrix of a multiplex network represented by , and are represented by , constituting both layers of the multiplex network.
2.3 Measures used for spectral properties
To study the spectral properties of multiplex networks with both layers consisting of directed random networks, we consider eigenvalue ratio statistics. The eigenvalue ratio statistics has been extensively used to characterize the localization to de-localization transition and, thus, Poisson to GOE [47]. In the case of the non-Hermitian adjacency matrix, complex eigenvalues () are obtained, and thus, complex spacing ratio is computed in the following manner [46]. is defined as the ratio (-th) between the distance of the nearest neighbor eigenvalue to the distance of the next to nearest neighbor eigenvalue,
[TABLE]
where and are nearest and next nearest neighbor respectively, of . However, for real eigenvalues (), in case of the Hermitian adjacency matrix, the -th ratio
[TABLE]
Here, and . For , the probability distribution function of for the PE and GUE [44] is given by,
[TABLE]
[TABLE]
GUE consists of random Hermitian matrices whose complex entries (both real and imaginary) are Gaussian variables. For GOE [44],
[TABLE]
Eq. (8) to (10) provide good approximations for also. In the case of , is not known. To validate our results for , we refer to the work [9]. However, recently, Dusa and Wettig derived an approximate expression of the probability distribution for the GUE [45].
3 Results and Discussions
Here, we first present the spectral properties i.e, spacing ratio results of the bi-layer multiplex networks. We compute the average ratio and from the eigenvalues obtained from the ensemble of matrices which are numerically diagonalized. For numerical investigations, network size, connection probability, and the average degree are kept the same and are changed simultaneously in both the layers.
3.1 Diluted Real Ginibre Ensemble
Fig. 4 presents the average ratio as a function of connection probability for the multiplex networks () with both layers represented by and multiplexed by . For different network sizes, is computed by taking the average over an ensemble of , being the size of each multiplex network. Fig. 4(a)-(c) depict the dependence of on the network size . Fig. 4(a) presents the small-size effects for , which can be seen in Fig. 4(b) plotted for two extreme values of . For (Fig. 4(b)), as increases from (when all the nodes are disconnected) to (complete network), transits smoothly from the Poisson ensemble statistics to the RGE . Though versus curves follow the same trend for different , transition occurs at the lower values of the as increases. This result implies that the PE-RGE transition of in a multiplex network begins at lower values of as is increased in both layers. Next, for an in-depth understanding of system size dependence on spectral statistics of multiplex networks, the change in position of - curves with a change in is ascertained by defining a parameter . is the value of in both layers of a multiplex network at which reaches the midpoint between its smallest and largest values. Fig. 4(c) consists of the horizontal dashed line parallel to the axis marking half of the complete transition at . In the inset of Fig. 4(d), is plotted for different values on the log-log scale, which is shown to change linearly with , indicating,
[TABLE]
On fitting the data, we find and . Following, we divide by ,
[TABLE]
Fig. 4(d) shows of a multiplex network as a function of . The Poisson-Ginibre transition of by varying is observed to transpire on the same value for all the network sizes considered here. This indicates that for a fixed average degree, of a multiplex network is also fixed independent of the network size . Also, Fig. 4(d), marks the region of of PE statistics (), and for , (RGE). Whereas exhibit intermediate statistics as transits from PE to RGE statistics in this region. Thus, and characterize the boundaries between regime of different statistics regions. The result implies that density of connections are determining factor in deciding the behavior of .
To further substantiate the system size invariance, probability distribution function is presented in Fig. 5. For all the five values of (marked as vertical lines in Fig. 4(d)), the histogram () is plotted for various network size of multiplex network, which is observed to fall on a single curve (Fig. 5). This result verifies the outcome that is independent of system size. Moreover, the results observed here for multiplex networks are aligned with the ones observed in single-layer networks [9]. However, values of that mark the start and end of the PE-Ginibre transition slightly differ in the case of multiplex networks. This behavior is attributed to the additional degree by virtue of the connection between the layers (for instance in Fig. 4 and Fig. 5).
Next, we focus on inter-layer connections and the influence of on the spectral statistics of multiplex networks. We study the statistics of as a function of and then as is varied. Fig. 6 presents as a function of for different network sizes and the multiplexing strength. We observe that as multiplexing strength is decreased, the value required for transition to occur increases. Also, for lower values of , as tends to in both the layers, saturates at lower value and do not approach RGE, hence a partial transition to RGE statistics is observed. This behavior is observed to be the same for all the network sizes.
The spectral response might be different for different values of . To characterize this; we define three regimes of the multiplexing strength : weak multiplexing, , where the layers have a negligible interaction between them, : intermediate multiplexing strength, , where the layers are coupled, and inter-layer connections play prominent role and : strong multiplexing, , where layers are superposed with each other. Fig. 7 (left) plots (value of at which reaches the midpoint between its smallest and largest values) versus to measure the change in PE-Ginibre transition curves owing to the change in value. It can be seen that in the regimes and , the ratio statistics are not affected by the multiplexing strength, whereas, in regime , decreases as is increased, which reflect that with an increase in , the transition occurs with a lower connection probability (Fig. 7). Thus, in regime , the change in affects the PE to RGE transition such that even when both the layers are represented by a complete network at , by means of lowering only, the is suppressed to transit to the Ginibre ensemble statistics. Additionally, Fig. 7 (right) presents statistics as a function of . versus curves are observed to follow the same behavior irrespective of the network size as multiplexing strength is changed.
3.2 Hermitian representation of adjacency matrix
This section studies the spectral properties of multiplex networks consisting of the directed network represented by Hermitian matrix in both the layers. To analyze the ratio statistics of , we make use of (Eq. (7)).
Fig. 8 presents the statistics of multiplex networks () consisting of a magnetic adjacency matrix in both layers multiplexed by . Similar to Fig. 4 (in the case of ), here also, we observe small-size effects mainly for (Fig. 8(a)), which is further shown in Fig. 8(b) for limiting values of connection probability . For , as a function of is observed to respond in the same manner as observed for . So, we again choose the parameter and plot it with (inset of Fig. 8(d)). After fitting the data and finding the scaling parameter (Eq. (11)), is plotted as a function of (Fig. 8(d)). Again, versus curves for different network sizes fall on the same curve, demonstrating that for known in both the layers, of a whole multiplex network can be identified. Thus, if is fixed in both the layers, of a multiplex network is also fixed.
The vs. curves show a smooth transition on increasing the value from at to at . On further increasing the , starts decreasing, and at , for as observed in single layer networks also. Here, values observed at , and correspond to RMT ensembles (), GUE () and GOE (), respectively, outlined in references [9, 44]. Thus, of a multiplex network undergoes the triple transition as a function of from PE to GUE and then GUE to GOE statistics. As changes, the network in both layers and thus the corresponding adjacency matrix undergoes the transformation; at , the block matrices and have only real diagonal entries indicating PE statistics, . Further increase in originates random uni-directional links in the network, causing random off-diagonal imaginary entries in the matrix and . At this point, GUE like statistics is observed for , however, the matrix cannot be said to be alike GUE. As , all the nodes connect to all other nodes giving rise to bi-directional links and, thus, vanishing imaginary entries in the matrix. Subsequently, GOE like statistics is observed. Note that the results obtained for multiplex networks show similar behavior to single-layer networks [9].
However, the range of for which the transition occurs differs slightly from those obtained in single-layer networks, as observed in the case of also. The results imply that on multiplexing two directed random networks, spectral statistics of the whole multiplex network remain the same as that of single-layer networks.
Now, we again investigate the impact of multiplexing strength on the relationship between structural properties and spectral properties of the multiplex networks consisting of in both layers. In Fig. 10, we present the statistics as a function of for different values and network sizes. First, for lower values of multiplexing strength ( and ), the connection probability required to begin the transition from PE to GUE statistics and then GUE to GOE statistics increases. Second: for a lowly multiplexed network (), with an increase in , does not complete the transition to GUE/GOE statistics; rather ceases at which shows that the spectral transition can be subdued by decreasing the multiplexing strength (). The results, here, infer and propound the multiplexing strength as the regulating parameter for statistics as a function of . Also, when , as p tends to , there is an abrupt transition from GUE to GOE statistics. However, for , transits gradually from near GUE to near GOE statistics. On further decreasing to , increases slightly with an increase in ; however, it remains in near PE statistics regime. Moreover, the behavior of spectral response with the variation in multiplexing strength is observed to be the same for different network sizes considered. Additionally, for a particular , curves as a function of coincide for different values.
4 Delocalization Transition
Entropy has been widely studied as a measure of the degree of irregularity [48], disorder [49], complexity [50], and randomness [51] in various systems. Also, it has been used to study the localization-delocalization transition in several systems [52, 53]. In this section, we apply entropy analysis, first: to quantify localization-delocalization transition in multiplex networks consisting of directed networks in both layers. Secondly, to verify the PE-RGE transition of and PE-GUE transition of by a change in the corresponding eigenvectors from the localized regime (PE) to extended (RGE or GUE). In the case of multilayer networks, delocalization transition has been studied using different measures [54, 55, 56]. Here, we make use of Shannon entropy , which for eigenvector corresponding to eigenvalue () is defined as
[TABLE]
has been used as an exemplary index for eigenvector localization-delocalization transition [57] as it provides information about the number of principal components of the eigenvector . has also been used to calculate the localization length in several systems [58]. Here, to compute Shannon entropy of multiplex networks, we numerically diagonalize the corresponding matrices of a large ensemble. Then the average is taken over all eigenvectors such that for an ensemble size of , is computed over eigenvectors. In the case of single-layer networks, for both and [9], when p=[math], only diagonal entries are present in the matrix, due to which, in the corresponding eigenvector, there is only one component with a magnitude, resulting in . As is increased to , now, because of the presence of off-diagonal entries, the corresponding eigenvectors will have nearly equal components resulting in . However, in the case of multiplex networks, an additional entry is present in each row of the corresponding adjacency matrix (), owing to interlayer connection () resulting in even at .
In Fig. 11, versus is presented for different network sizes () and multiplexing strength () of multiplex networks. Fig. 11(a) and 11(b) show for multiplex networks represented by and , respectively. For a particular multiplexing strength, as increases from [math] to , is observed to manifest a smooth transition from =[math] (when only diagonal entries are present, and no nodes are connected) to its maximum value (when all the nodes are connected). Also, for different network sizes , shows the same behavior with respect to . Now, when =[math] in both the layers, eigenvectors of a multiplex network are localized as , and as tends to , reaches its maximum value, then the eigenvectors are delocalized. The results obtained are the same for both the multiplex networks consisting of networks corresponding to and in both layers. However, for and is different, so to get the normalized value of average entropy, is computed for all the values of . Also, subject to , is which is the Shannon entropy of eigenvectors of RGE matrices. Note that considered here is referred from the recent work [9] where it was computed numerically for different network sizes. In the case of a multiplex network , = and [57]. Now, as we decrease the multiplexing strength, the average degree required to begin the transition increases. Also, for lower , as is increased, is observed to saturate at lower values than . Moreover, for =, because of the additional inter-layer connection, even at =[math], whereas for = and , .
5 Conclusion
To conclude, we have numerically analyzed the spectral properties of multiplex networks consisting of directed random networks in both layers. In particular, by employing the spacing ratio measure, we investigated the average ratio as a function of various structural parameters of the network in layers. Without a requisition of the unfolding, the eigenvalue spacing ratio can determine the spectral statistics of various systems. Here, we study ER random networks by characterizing them in two different manners, () diluted real Ginibre ensemble (dRGE), which are sparse non-Hermitian random matrices and are analyzed using the complex spacing ratio () (Eq. (4)), () sparse Hermitian matrices for which real spacing ratio () is computed (Eq. (5)). Multiplexing two networks yields the spacing ratio statistics, as that of the corresponding single-layer networks [9]. The average ratio / is shown to transit smoothly from PE to RGE/GUE statistics as a function of . The value for which this transition occurs manifests a decrease with an increase in the network size. After scaling and varying as a function of , the transition curves for different network size are observed to coincide on a single universal curve. Further, in recent years, studies on multiplex networks have revealed additional information than those of the corresponding single-layer networks. To illustrate this, we analyzed the spacing ratio statistics for different multiplexing strengths. For a multiplex network represented by a diluted version of real Ginibre ensemble matrices in both layers, varying the multiplexing strength affect the ratio statistics in two ways. First, at =, shows a transition from PE to RGE statistics as tends to . However, for , as tends to , saturates at lower values, and a complete transition to RGE statistics is suppressed. Second, with a decrease in the multiplexing strength, more connection probability is required to commence the transition. Next, for multiplex networks consisting of sparse Hermitian matrices in both the layers, a triple transition PE-GUE-GOE is observed with an increase in the connection probability. We show that upon decreasing the multiplexing strength hinders the PE-GUE and GUE-GOE statistics transition. Also, for =, as tends to , a sudden transition occurs from GUE-GOE statistics which for lower values of is found to occur gradually.
We wish to note here that the results presented here are akin to the phase transitions concerning Raman scattering [59] and materials such as MoS2 nanoflakes [60]. The results achieved here using multilayer model can be helpful in better understanding of several underlying mechanisms and dynamics of such systems. Additionally, we analyze the localization-delocalization transition of the eigenvectors of multiplex networks using Shannon entropy. We show that by decreasing the multiplexing strength, even for = (when the network is complete), the eigenvectors of a multiplex network are not wholly delocalized, which again implies that can act as a suppressing parameter in studying several processes in multiplex networks. In recent years, eigenvalue spacing ratio statistics has gained enormous attention from researchers because of its wide application in various fields such as condensed matter, chaotic systems, etc. Further, spacing ratio technique has already been used to measure the localization-delocalization transition in different systems, including single-layer networks. We believe that spacing ratio technique can potentially be used to study various structural and dynamical properties of complex systems represented by multiplex networks having directed connections. For instance, eigenvalues of laplacian matrix can give insight into synchronization of corresponding complex network. This work has focused on adjacency matrix only, one can analyze the laplacian of multiplex directed networks.
SJ and Tanu acknowledge Govt. of India, BRNS Grant No. 37(3)/14/11/2018-BRNS/37131 for financial support and SRF fellowship, respectively. SJ thankfully acknowledges DST grant SPF/2021/000136. We thank Ranveer Singh and Stefano Boccaletti for useful suggestions and interesting discussions.
References
- [1] T. Guhr et al., Phys. Rep. 299, 189-425, (1998).
- [2] J.J.M. Verbaarschot, T. Wettig. Random Matrix Theory and Chiral Symmetry in QCD. Annual Review of Nuclear and Particle Science. 50(1), 343-410, (2000); B. I. Shklovskii, B. Shapiro, B. Sears, P. Lambrianides, and H. Shore, Phys. Rev. B 47, 11487, (1993); V. Plerou, P. Gopikrishnan, B. Rosenow, L. A. Nunes, Amaral, and H. E. Stanley, Phys. Rev. Lett. 83, 1471, (1999); C. E. Porter, Statistical theories of spectra: fluctuations, Tech. Rep. (1965).
- [3] M.L. Mehta, Random Matrices, Academic Press, (1991).
- [4] F. J. Dyson and M. L. Mehta, J. Math. Phys. 4, 701, (1963)
- [5] I. O. Morales, E. Landa, P. Stránsky, and A. Frank, Phys. Rev. E 84, 016203, (2011); G. Torres-Vargas, R. Fossion, and J. Mendez-Bermudez, Physica A 545, 123298, (2020)
- [6] V. Oganesyan and D. A. Huse, Phys. Rev. B 75, 155111, (2007)
- [7] T. Kanazawa and T. Wettig. New universality classes of the non-Hermitian Dirac operator in QCD-like theories. Phys. Rev. D. 104(1), 014509, (2021); H. H. Jen. Quantum correlations of localized atomic excitations in a disordered atomic chain. Phys. Rev. A. 105(2), 023717, (2022); O. Giraud, N. Macé, E. Vernier and F. Alet. Probing Symmetries of Quantum Many-Body Systems through Gap Ratio Statistics. Phys. Rev. X. 12(1), 011006, (2022); L. Sá, P. Ribeiro, and T. Prosen. Complex Spacing Ratios: A Signature of Dissipative Quantum Chaos. Phys. Rev. X. 10(2), 021019, (2020).
- [8] A. Mishra and T. Raghav and S. Jalan. Eigenvalue ratio statistics of complex networks: Disorder versus randomness. Phys. Rev. E. 105(6), 064307, (2022); C.T. Martínez-Martínez, J.A. Méndez-Bermúdez, T. Peron, Y. Moreno. Statistical properties of mutualistic-competitive random networks.Chaos, Solitons and Fractals. 153(2), 111504, (2021).
- [9] T. Peron, and B. M. F. de Resende and F. A. Rodrigues, and L. da F.Costa, and J. A. Méndez-Bermúdez. Spacing ratio characterization of the spectra of directed random networks. Phys. Rev. E. 102(6), 062305, (2020).
- [10] M.E.J. Newman. The Structure and Function of Complex Networks. SIAM Review. 45(2), (2003).
- [11] R. Albert, and A-L. Barabási. Statistical mechanics of complex networks. Rev. Mod. Phys. 74(1), 47-97, (2002).
- [12] M. Kivelä, A. Arenas, M. Barthelemy, J. P. Gleeson, Y. Moreno and M A Porter. Multilayer networks, Journal of Complex Networks, (3), 203–271, (2014).
- [13] M. De Domenico, C. Granell, M. A. Porter and A. Arenas. The physics of spreading processes in multilayer networks. Nature Phys 12, 901–906, (2016).
- [14] S. Boccaletti, G. Bianconi, R. Criado, C.I. del Genio, J. Gómez-Gardeñes, M. Romance, I. Sendiña-Nadal, Z. Wang and M. Zanin. The structure and dynamics of multilayer networks. Physics Reports. 544(1), 1-122, (2014).
- [15] T. Raghav and S. Jalan. Random matrix analysis of multiplex networks. Physica A: Statistical Mechanics and its Applications. 586, 126457, (2022).
- [16] S. Jalan and V. Rathore and AD Kachhvah and A. Yadav. Inhibition-induced explosive synchronization in multiplex networks. Phys. Rev. E. 99(6), 062305, (2019).
- [17] R.S. Escoboza1, I.S. Nadal, I. Leyva, R. Gutiérrez, J.M. Buldú, S. Boccaletti, Inter-layer synchronization in multiplex networks of identical layers, Chaos 26(6), 065304, (2016).
- [18] G. F.de Arruda et al J. Stat. Mech. 103405, (2020).
- [19] S. Majhi, M. Perc and D. Ghosh. Dynamics on higher-order networks: a review. J. R. Soc. Interface. 19 (2022).
- [20] J.A. Dunne, R.J. Williams, N.D. Martinez, Food-web structure and network theory: The role of connectance and size, Proc. Natl. Acad. Sci. USA 99(20), 12917–12922, (2002).
- [21] G. Weng, U.S. Bhalla, R. Iyengar, Complexity in biological signaling systems, Science 284(5411), 92–96, (1999).
- [22] S. Redner, How popular is your paper? An empirical study of the citation distribution, Eur. Phys. J. B 4, 131–134, (1998).
- [23] R. Albert, H. Jeong and AL. Barabási. Diameter of the World-Wide Web. Nature 401, 130–131, (1999).
- [24] D. M. Cvetković, M. Doob, and H. Sachs, Spectra of Graphs: Theory and Applications, 3rd revised ed. Academic Press. (1997).
- [25] I. J. Farkas, I. Derenyi, A. L. Barabasi, T. Vicsek. Spectra of real-world graphs: Beyond the semicircle law. Phys. Rev. E. 64(2), 026704, (20).
- [26] D. M. Cvetković, P. Rowlinson, and S. Simic. An Introduction to the Theory of Graph Spectra (London Mathematical Society Student Texts). Cambridge: Cambridge University Press. (2009).
- [27] J. G. Restrepo, E. Ott, and B. R. Hunt. Onset of synchronization in large networks of coupled oscillators. Phys. Rev. E. 71(3), 036151, (2005).
- [28] J. G. Restrepo, E. Ott, and B. R. Hunt. Characterizing the dynamical importance of network nodes and links. Phys. Rev. Lett. 97(9), 094102, (2006).
- [29] M. A. M. de Aguiar and Y. Bar-Yam. Spectral analysis and the dynamic response of complex networks. Phys. Rev. E. 71(1), 016106, (2005).
- [30] M. Zhou, J. Tan, H. Liaoa, Z. Wang, and R. Mao. Dismantling complex networks based on the principal eigenvalue of the adjacency matrix. Chaos. 30(8), 083118, (2020).
- [31] P. Erdös and A. Rényi. On Random Graphs I. Publicationes Mathematicae Debrecen. 6, 290-297, (1959). P. Erdös and A. Rényi. On the Evolution of Random Graphs. Publication of the Mathematical Institute of the Hungarian Academy of Sciences. 5, 17-61, (1960). P. Erdös and A. Rényi. On the evolution of random graphs. Bulletin of the Institute of International Statistics. 38(4), 343-347, (1961).
- [32] S. Jalan, G. Zhu and B. Li. Spectral properties of directed random networks with modular structure. Phys. Rev. E. 84(4), 046107, (2011).
- [33] B. Ye and L. Qiu and X. Wang and T. Guhr. Spectral statistics in directed complex networks and universality of the Ginibre ensemble. Communications in Nonlinear Science and Numerical Simulation. 20(3), 1026-1032, (2015); B. Ye, K. Zuo and J. Jia, Random matrix analysis of spectral properties in directed complex networks, The 26th Chinese Control and Decision Conference (2014 CCDC), 616-620, (2014).
- [34] C. Li and Z. Qu. Distributed estimation of algebraic connectivity of directed networks. Systems and Control Letters. 62(6), 517-524, (2013).
- [35] I. J. Farkas, I. Derényi, A. L. Barabási, and T. Vicsek. Spectra of real-world graphs: Beyond the semicircle law. Phys. Rev. E. 64(2), 026704, (2001).
- [36] C. Sarkar, and S. Jalan. Spectral properties of complex networks. Chaos 28, 102101, (2018).
- [37] J. W. Baron. Eigenvalue spectra and stability of directed complex networks. arXiv:2206.13436
- [38] F. L. Metz and I. Neri. Localization and Universality of Eigenvectors in Directed Random Graphs. Phys. Rev. Lett. 126(4), 040604, (2021).
- [39] K. Guo and B. Mohar, J. Graph Theory 85(217), (2017).
- [40] J. Liu and X. Li, Linear Algebra Appl. 466(182), (2015).
- [41] S. Samanta, M. Pal, R. Mahapatra, K. Das and R. S. Bhadoria. A Study on Semi-directed Graphs for Social Media Networks. International Journal of Computational Intelligence Systems. 14(1), (2021).
- [42] S. Samanta, G. Muhiuddin, A. M. Alanazi, and K. Das. A Mathematical Approach on Representation of Competitions: Competition Cluster Hypergraphs. Mathematical Problems in Engineering. (2020).
- [43] S. Kubota, E. Segawa and T. Taniguchi. Quantum walks defined by digraphs and generalized Hermitian adjacency matrices. Quantum Information Processing, 20(95), (2021).
- [44] Y. Y. Atas, E. Bogomolny, O. Giraud, and G. Roux. Distribution of the Ratio of Consecutive Level Spacings in Random Matrix Ensembles. Phys. Rev. Lett. 110(8), 084101, (2013).
- [45] I. G. Dusa and T. Wettig. Approximation formula for complex spacing ratios in the Ginibre ensemble. Phys. Rev. E. 105(4), 044144, (2022).
- [46] L. Sá, P. Ribeiro and T. Prosen. Complex Spacing Ratios: A Signature of Dissipative Quantum Chaos. Phys. Rev. X. 10(2), 021019, (2020).
- [47] C.T. Martínez-Martínez and J.A. Méndez-Bermúdez and Y. Moreno and J. J. Pineda-Pineda and J. M. Sigarreta. Spectral and localization properties of random bipartite graphs. Chaos, Solitons and Fractals: X. 3, 100021, (2019).
- [48] Y. Sun, H. Zhao, J. Liang and Ma X. Eigenvalue-based entropy in directed complex networks. PLoS One, 16(6), (2021).
- [49] C.G.S. Freitas, A.L.L. Aquino and H.S. Ramos et al. A detailed characterization of complex networks using Information Theory. Sci Rep 9, 16689, (2019).
- [50] K. Anand, and G. Bianconi. Entropy measures for networks: Toward an information theory of complex topologies. Phys. Rev. E. 80(4), 045102, (2009).
- [51] Y.M. Omar and P. A Plapper. Survey of Information Entropy Metrics for Complex Networks. Entropy, 22, 1417, (2020).
- [52] O. Farzadian and T. Oikonomou and M.R.R. Good and M.D. Niry. Entropic analysis of the localization–delocalization transition in a one-dimensional correlated lattice. Physica A: Statistical Mechanics and its Applications. 545, 123350, (2020).
- [53] C. T. Martínez-Martínez and J. A Méndez-Bermúdez. Information Entropy of Tight-Binding Random Networks with Losses and Gain: Scaling and Universality. Entropy. 21(1), (2019).
- [54] G. F. de Arruda, E. Cozzo, T. p. Peixoto, F. A. Rodrigues, and Y. Moreno. Disease Localization in Multilayer Networks. Phys. Rev. X. 7(1), 011014, (2017).
- [55] J. A. Méndez-Bermúdez, G. F. de Arruda and F. A. Rodrigues and Y. Moreno. Scaling properties of multilayer random networks. Phys. Rev. E. 96(1), 012307, (2017).
- [56] G. Bianconi. Statistical mechanics of multiplex networks: Entropy and overlap. Phys. Rev. E. 87(6), 062806, (2013).
- [57] B. Mirbach and H.J. Korsch. A Generalized Entropy Measuring Quantum Localization. Annals of Physics. 265(1), 80-97, (1998).
- [58] J. A. Méndez-Bermúdez, A. Alcazar-López and A. J. Martínez-Mendoza, F. A. Rodrigues and T. Peron. Universality in the spectral and eigenfunction properties of random networks. PHYSICAL REVIEW E 91(3), 032122, (2015) and the references within.
- [59] P. Yogi et. al. Quantifying the Short-Range Order in Amorphous Silicon by Raman Scattering. Anal. Chem. 90(13), 8123–8129, (2018).
- [60] M. Tanwar et. al. Fano-Type Wavelength-Dependent Asymmetric Raman Line Shapes from MoS2 Nanoflakes. ACS Phys. Chem Au. 2(5), 417–422, (2022).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] T. Guhr et al., Phys. Rep. 299 , 189-425, (1998).
- 2[2] J.J.M. Verbaarschot, T. Wettig. Random Matrix Theory and Chiral Symmetry in QCD. Annual Review of Nuclear and Particle Science. 50 (1), 343-410, (2000); B. I. Shklovskii, B. Shapiro, B. Sears, P. Lambrianides, and H. Shore, Phys. Rev. B 47 , 11487, (1993); V. Plerou, P. Gopikrishnan, B. Rosenow, L. A. Nunes, Amaral, and H. E. Stanley, Phys. Rev. Lett. 83 , 1471, (1999); C. E. Porter, Statistical theories of spectra: fluctuations, Tech. Rep. (1965).
- 3[3] M.L. Mehta, Random Matrices, Academic Press, (1991).
- 4[4] F. J. Dyson and M. L. Mehta, J. Math. Phys. 4 , 701, (1963)
- 5[5] I. O. Morales, E. Landa, P. Stránsky, and A. Frank, Phys. Rev. E 84 , 016203, (2011); G. Torres-Vargas, R. Fossion, and J. Mendez-Bermudez, Physica A 545 , 123298, (2020)
- 6[6] V. Oganesyan and D. A. Huse, Phys. Rev. B 75 , 155111, (2007)
- 7[7] T. Kanazawa and T. Wettig. New universality classes of the non-Hermitian Dirac operator in QCD-like theories. Phys. Rev. D. 104 (1), 014509, (2021); H. H. Jen. Quantum correlations of localized atomic excitations in a disordered atomic chain. Phys. Rev. A. 105 (2), 023717, (2022); O. Giraud, N. Macé, E. Vernier and F. Alet. Probing Symmetries of Quantum Many-Body Systems through Gap Ratio Statistics. Phys. Rev. X. 12 (1), 011006, (2022); L. Sá, P. Ribeiro, and T. Prosen. Complex Spacin
- 8[8] A. Mishra and T. Raghav and S. Jalan. Eigenvalue ratio statistics of complex networks: Disorder versus randomness. Phys. Rev. E. 105 (6), 064307, (2022); C.T. Martínez-Martínez, J.A. Méndez-Bermúdez, T. Peron, Y. Moreno. Statistical properties of mutualistic-competitive random networks.Chaos, Solitons and Fractals. 153 (2), 111504, (2021).
