- Department of Mathematics, University of the Punjab, Lahore, Pakistan
This paper presents the practical applications of Laplacian and signless Laplacian spectra across various fields including theoretical chemistry, computer science, electrical engineering, and complex network analysis. By focusing on the spectrum-based evaluation of generalized mesh network and ladder graphs, the research aims to uncover valuable relationships with the structural properties of real-world networks. The study not only explores the theoretical underpinnings but also applies these spectra to calculate essential network measures such as mean-first passage time, average path length, spanning trees, and spectral radius. These analyses offer a deeper understanding of how graph spectra influence network characteristics, enriching our ability to predict and analyze complex networks. This comprehensive approach enhances our knowledge across multiple scientific disciplines, facilitating more informed predictions about drugs infrastructure.
1 Introduction
Spectral graph theory is an important part of algebraic graph theory which mainly utilizes matrix theory, polynomial theory, and combinatorial methods to study the different spectra (or ranges of values) that come from graphs. Furthermore, spectral theory investigates how these spectra are connected to the structure and properties of graphs. It links the algebraic (math-based) aspects of graphs to their topological (shape-based) aspects. By studying how a graph's spectrum relates to its structure, we can not only better understand these graphs but also find useful applications in areas like improving networks, designing computer circuits, and solving operational problems. Some key areas of study in spectral theory include the adjacency spectrum, Laplacian spectrum, signless-Laplacian spectrum, and distance spectrum of graphs. Among these, the Laplacian spectrum is the most studied and produces the most results. Studying the Laplacian spectrum is not only valuable for theoretical knowledge but also has many uses in chemistry, physics, complex networks, and electronic engineering. Over the last few decades, significant attention has been given to examining the structure of graphs, along with their spectral, topological, and combinatorial characteristics. Laplacian eigenvalues have proven useful in identifying various graph invariants, including the Kirchhoff index, global mean-first passage time, and the count of spanning trees. Typically, the characteristic polynomial and spectrum of the graph matrix for certain graph operations, such as the complement, union, Cartesian product, direct product, and strong product, can be derived from the factor graphs. Since complex molecular graphs can be effectively described using graph operations, it is feasible to represent the properties of these complex molecular graphs through the invariants of their factor graphs.
From the perspective of spectral graph theory, numerous structural features and dynamic behaviors of graphs have been investigated. The literature on spectral graph theory covers diverse aspects of Laplacian matrices across different graph structures. Merris provides a comprehensive survey, discussing the Laplacian matrix's spectrum, algebraic connectivity, and applications in areas such as chemistry [1]. Hong and Zhang investigate Laplacian eigenvalues in simple and bipartite graphs, establishing key bounds and structural insights, particularly for trees and regular graphs [2]. Agaev and Chebotarev (2006) extended the study to weighted directed graphs, exploring the connections between Laplacian and stochastic matrices and their semiconvergent properties [3]. Rojo and Soto focused on unweighted rooted trees, analyzing the eigenvalues of adjacency and Laplacian matrices based on symmetric tridiagonal matrices [4]. Ding and Jiang investigate the spectral norms and eigenvalue distributions of random graph matrices, revealing key convergence behaviors aligned with Wigner's semi-circular law [5]. Wu extends the use of Laplacian matrices into quantum mechanics, exploring conditions for separability in weighted graphs [6]. Kaveh and Rahami (2006) focus on the eigenvalues and eigenvectors of graph products, presenting efficient methods for solving eigenproblems in structural mechanics, especially for Cartesian and lexicographic products [7]. Spielman emphasizes the significance of Laplacian matrices in algorithm design, highlighting their role in fast solutions for linear equations and their application to graph theory through innovations like graph sparsifiers and local clustering [8].
In 2012, Estrada introduced path Laplacian matrices as a new concept that generalized the combinatorial Laplacian and applied them to consensus analysis in networks, thus showing its potential for enhancing network synchronization and other applications [9]. In their 2013 publication, Krishnan et al. came up with a novel scheme of multi-level preconditioning for Laplacian matrices used in computer graphics which resulted in substantial performance benefits in applications such as image colorization and mesh processing [10]. Another similar work by Dong et al. aimed toward exploring laplacian matrix learning for smooth graph signal representation that contributed to advancement in graph signal processing [11]. Pirani and Sundaram analyzed the smallest eigenvalue properties of grounded Laplacian matrices giving out insights into spectral graph theory and its applications [12]. Efficient methods have been developed by Bergamaschi and Martínez to approximate the generalized inverse of Laplacian matrices which are important when solving large scale graph problems [13]. Recent research on Laplacian matrices has led to notable advancements. In 2016, Jog and Kotambari analyzed the spectra of coalesced complete graphs, studying the adjacency, Laplacian, and signless Laplacian energies to understand their spectral properties and applications [19]. Moving to 2018, Bandeira explored random Laplacian matrices, revealing that the largest eigenvalue often approximates the largest diagonal entry, with implications for convex relaxation techniques and Erdos-Rényi graph connectivity thresholds [14]. Li provided insights into the constrained Rayleigh quotient for eigen-balanced Laplacian matrices, which proved valuable for cooperative control problems and convergence rates in consensus protocols [16]. The work by Bergamaschi and Bozzo focused on comparing algorithms for computing the smallest eigenpairs of graph Laplacians, including the Implicitly Restarted Lanczos Method and Jacobi-Davidson method, particularly for large, sparse networks [18]. Zhou et al. introduced an optimal neighborhood multi-view spectral clustering algorithm, which enhances clustering performance by effectively combining first-order and high-order Laplacian matrices [15]. Hermann and Konigorski addressed the optimization of edge weights in directed graph Laplacians to achieve desired spectral properties [17].
In 2019, Alhevaz et al. explored the Brouwer-type conjecture related to the eigenvalues of the distance signless Laplacian matrix. Their findings provided bounds for the sums of the largest and smallest eigenvalues, applying these results to graphs with specific diameters and transmission properties [20]. Moving forward to 2022, Ganie and Shang investigated the spectral radius and energy of the signless Laplacian matrix of digraphs, proposing new lower bounds and characterizing extremal digraphs based on vertex degrees and walk lengths [21]. Also in 2022, Morbidi examined matrix functions of the Laplacian matrix and their applications to distributed formation control, showing how these functions can enhance performance and flexibility in consensus protocols [22]. Recent studies have significantly advanced our understanding of various Laplacian matrices and their applications. In 2021, Reinhart introduced the normalized distance Laplacian matrix, offering new insights into its spectral properties and connections with the normalized Laplacian matrix. The study showed that this matrix has fewer cospectral pairs compared to other matrices [23]. The same year, Chakrabarty et al. explored the spectral properties of adjacency and Laplacian matrices in inhomogeneous Erdõs-Rényi random graphs. Their work detailed the empirical spectral distributions and their convergence to deterministic limits [24].
The paper by Alazemi et al. [25] explores chain graphs, a specific class of bipartite graphs, with unique Laplacian eigenvalues. The authors provide structural insights, degree constraints, and analyze the eigenspaces of these graphs. Notably, they highlight conditions such as the absence of vertex triplets sharing identical neighborhoods and propose applications in Laplacian dynamics, including the controllability of multi-agent systems. Meanwhile, Anđelić et al. [26] introduce a family of tridiagonal matrices with eigenvalues as perfect squares, applying this result to analyze the Laplacian controllability of half graphs, a subclass of chain graphs, further advancing the understanding of spectral graph theory. The authors in [28] investigate the Laplacian controllability of graphs formed using standard graph products, including joins, Cartesian, tensor, and strong products. The study provides theoretical insights and introduces an iterative method to construct infinite families of controllable Laplacian pairs. Additionally, Anđelić et al. [29] focus on the Q-index, the largest eigenvalue of the signless Laplacian matrix, for connected graphs with fixed order and size. The authors derive spectral bounds for the Q-index of nested split graphs, offering both theoretical results and computational comparisons to improve understanding of spectral properties in graph theory.
More recently, in 2023, Bapat et al. extended the concept of bipartite matrices by examining the bipartite Laplacian matrix of nonsingular trees. They provided a combinatorial description of this matrix and established several key identities [27]. Additionally, Mallik expanded the Matrix Tree Theorem to signed graphs, introducing a new oriented incidence matrix and offering a combinatorial formula for the determinant of the signless net Laplacian matrix [30]. Raza et al. focused on generalized prism graphs and found that spectral analysis helps measure network features like passage time and path length [31]. Raza and Munir extended this by showing how Laplacian and signless Laplacian spectra can be used to understand network properties and predict behaviors in various fields [32]. In a later study, Raza et al. applied these methods to torus grid graphs, deriving key network measures and improving our knowledge of network structures [33].
We define a path of length αm ∈ ℕ as the graph Pαm that has vertex set V = {v ∈ ℕ : 0 ≤ v ≤ αm} and where two vertices determine an edge if and only if |vi−vj| = 1 for vi, vj ∈ V. Then, a mesh network graph , often known as two-dimensional lattice graph or grid graph, is defined as the Cartesian product , exhibits a total of 2mn−m−n = (m−1)n+(n−1)m edges, reflecting the combined count of horizontal and vertical edges. Simultaneously, boasts mn vertices, aligning with the Cartesian product of the vertex sets of Pm and Pn. The mentioned graph operation have gained significant attention in the field of graph theory and computer science. These graphs are widely used to model spatial relationships and connectivity in various applications, such as computer networks, image processing, and computational geometry. The study of grid graphs has evolved over the years, with researchers exploring their properties, algorithms, and applications. Harel and Sardashti [34] presented a comprehensive analysis of the structural characteristics of mesh network graph, highlighting their regularity and symmetry. Smith et al. [35] investigated efficient algorithms for computing shortest paths in mesh network graphs, providing valuable insights into optimizing navigation in grid-based environments. Additionally, Chen and Du [36] explored the application of in wireless sensor networks, showcasing their relevance in practical scenarios. Recent work by Hinz and Holz auf der Heide [35] delved into the dynamic aspects, addressing challenges related to real-time updates and adaptability. Furthermore, the survey by Kumar et al. [37] offers a holistic overview of mesh network graph applications and algorithmic advancements. Building on earlier research about mesh network graphs and their spectra, our study thoroughly examined the Adjacency et al. Laplacian spectra of . We didn't just calculate these spectra; we also applied them to real-world network analysis. Using our results, we computed important network measures such as graph energies, Kirchhoff index, mean-first passage time, path length, spanning trees, and spectral radius. This approach aimed to give a deeper insight into the properties of , such as its connectivity, resilience, and efficiency. In this section, we review key findings from earlier studies that relate to the solutions discussed in this paper. This study distinguishes itself by conducting a comprehensive spectral analysis of generalized mesh networks, focusing on the adjacency, Laplacian, and signless Laplacian spectra to derive explicit expressions for critical network parameters such as the Kirchhoff index, spectral radius, average path length, global mean first passage time, graph energies, and the number of spanning trees. While previous research has applied spectral methods to various network types, such as torus networks and categorical product networks, your work uniquely emphasizes generalized mesh networks, providing detailed spectral characterizations that enhance the understanding of their structural and dynamic properties. The complete structure of our article is presented hierarchically in Figure 1. By presenting results graphically, your study offers clear visualizations of how these parameters vary with network dimensions, facilitating deeper insights into their interplay and impact. This approach not only broadens the applicability of spectral methods but also offers a robust framework for exploring and optimizing complex real-world networks, thereby contributing valuable perspectives across multiple scientific disciplines.
2 Preliminaries
Before discussing graph-based matrices and the related lemmas associated with the Kronecker product, let's first revisit the notion of ψ-sum graphs. Let $ = (V($), E($)) be a simple undirected graph, where V($) represents its vertex set and E($) represents its edge set. The number of vertices in $ is denoted by |V($)|, and the number of edges is denoted by |E($)|. If an edge e connects two vertices u and v, the edge uv can also be referred to as e. For a given vertex v ∈ V($), its neighborhood in $, denoted by N$(v), is the set of vertices adjacent to v, specifically N$(v) = {u ∈ V($)∣uv ∈ E($)}. The degree of a vertex v, symbolized as d$(v), is the number of vertices in its neighborhood, i.e., d$(v) = |N$(v)|.
Definition 1. The diagonal matrix is defined as for i = j and the Laplacian matrix is defined by the subtraction of the adjacency matrix from the diagonal matrix of vertex degrees. Elaborating in matrix form, is defined as:
Lemma 1. Let E ∈ Mp, p(G), F ∈ Mp, q(G), G ∈ Mq, p(G), H ∈ Mq, q(G) with H being invertible, such that
Then,
Lemma 2. Let C = (cij) ∈ Mr, s(G), D ∈ Mt, u(G). Then the Kronecker product of C and D is defined as
Lemma 3. Let C and D be square matrices of order r and s, respectively, with eigenvalues λi (1 ≤ i ≤ r) and νj (1 ≤ j ≤ s). Then the eigenvalues of C ⊗ Is+Ir ⊗ D are λi+νj. Moreover, if Vi is an eigenvector of C corresponding to λi and Wj is an eigenvector of D corresponding to νj, then Vi ⊗ Wj is an eigenvector of C ⊗ Is+Ir ⊗ D corresponding to λi+νj.
Lemma 4. Let C ∈ Mr, s(G), D ∈ Mt, u(G), E ∈ Mr, t(G), F ∈ Ms, u(G), and β ∈ G. The following properties hold:
(a) (C ⊗ D)T = CT ⊗ DT.
(b) (C ⊗ D)(E ⊗ F) = (CE) ⊗ (DF).
(c) (C ⊗ D) ⊗ E = C ⊗ (D ⊗ E).
(d) β(C ⊗ D) = βC ⊗ D = C ⊗ βD.
(e) If C and D are invertible, then (C ⊗ D)−1 = C−1 ⊗ D−1.
Lemma 5. The eigenvalues of the Adjacency matrix, Laplacian matrix, and Signless Laplacian matrix for a path graph are expressed as , , and , respectively, where k = 0, 1, 2, …, n−1.
3 Methodologies and results
In this section, we have evaluated the exact values for the Adjacency, Laplacian and signless Laplacian spectrum of the generalized Mesh Network graphs utilizing the graph and algebra techniques. Theorem 1 provides expressions for the sum of reciprocals and the product of adjacency eigenvalues of the generalized mesh graph , which are fundamental in understanding network connectivity and robustness [38]. The sum of the reciprocals of the eigenvalues is often associated with resistance distance and other network invariants, while their product is related to graph determinant properties, which have applications in quantum networks and structural analysis [39]. Extending this analysis, Theorem 2 focuses on the Laplacian eigenvalues, which play a crucial role in describing network dynamics such as diffusion processes and synchronization [41]. The sum of the reciprocals of the Laplacian eigenvalues is connected to important network measures like Kirchhoff's index, influencing resistance-based properties, while their product is associated with the number of spanning trees, a key quantity in evaluating network reliability and resilience [42]. Furthermore, Theorem 3 explores the Signless Laplacian spectrum, which is particularly useful in applications involving directed flows and energy distribution in networks [40]. The sum of the reciprocals of these eigenvalues helps in analyzing clustering tendencies in complex networks, whereas their product provides a measure of structural stability and modular properties. These interpretations establish strong connections between spectral properties and real-world applications, enhancing the accessibility of the results for researchers in diverse fields such as physics, computer science, and engineering.
Theorem 1. Let the summation of the reciprocals and the product of the adjacency eigenvalues of the generalized mesh graph be denoted by and , respectively. Then:
Proof. The adjacency matrix of the mesh graph is:
By matrix addition, it can be expressed as:
Using Lemma 2, we have:
The matrix
is the adjacency matrix of Qm, a path graph with m vertices. Thus:
Now, assume two invertible matrices U and V related to the matrices Qn and Qm, such that:
and
Since, the eigenvalues of the Adjacency matrix for a path graph Qn are given by so the diagonal entries of the upper triangular matrices are:
Consequently:
and the diagonal entries of this upper triangular matrix are given by:
Thus, the adjacency eigenvalues for the generalized mesh graph are:
Using this result, we obtain:
and
Corollary 1. For a mesh graph with equal dimensions (n = m), the product and sum of the reciprocals of the adjacency eigenvalues are given by:
The proof follows directly from Theorem 1.
Theorem 2. Let and denote the sum of the reciprocals and the product of the Laplacian eigenvalues, respectively, for the generalized mesh network graph . Then, these quantities are given by:
and
Proof. The Laplacian matrix associated with the mesh network graph is expressed as:
By decomposing this matrix, it can be rewritten as:
Referring to Lemma 2, this is further simplified as:
The matrix
is the Laplacian matrix for the path graph Qn with n nodes. Therefore:
Suppose U and V are invertible matrices related to the matrices Qn and Qm, respectively. Then:
Since, The eigenvalues of the Laplacian matrix for a path graph Qn are given by so the diagonal entries of the upper triangular matrices are:
where i = 0, 1, 2, …, m−1 and j = 0, 1, 2, …, n−1.
Clearly:
with diagonal elements of the resulting matrix given by:
where i = 0, 1, 2, …, m−1 and j = 0, 1, 2, …, n−1.
Thus, the eigenvalues of the Laplacian matrix for the mesh network graph are:
where i = 0, 1, 2, …, m−1 and j = 0, 1, 2, …, n−1.
Finally, using the eigenvalues in Equation 1, we derive:
and
Corollary 2. For a mesh network graph of equal dimensions (n = m), the product and reciprocal of the sum of the Laplacian eigenvalues are given by:
The proof follows directly from Theorem 2.
Theorem 3. Let the sum of the reciprocals and the product of all Signless Laplacian eigenvalues of the generalized mesh graph be denoted by and , respectively. Then,
Proof. The Signless Laplacian matrix for the mesh graph is expressed as:
which can be broken down by matrix addition as follows:
According to Lemma 1.1, we have:
The matrix given by
is, in fact, the Signless Laplacian matrix of the path graph Pn with n vertices. Thus, we obtain:
Introducing two invertible matrices A and B that correspond to the matrices Pn and Pm, we have:
and
Since, the eigenvalues of the Signless Laplacian matrix for a path graph Qn are given by so the diagonal entries of the upper triangular matrices are:
where i = 0, 1, 2, …, m−1 and j = 0, 1, 2, …, n−1.
Clearly, the following holds:
where the diagonal elements of this matrix are given by:
where i = 0, 1, 2, …, m−1 and j = 0, 1, 2, …, n−1.
Thus, the adjacency eigenvalues for the mesh graph can be expressed as:
where i = 0, 1, 2, …, m−1 and j = 0, 1, 2, …, n−1.
From the results in Equation 3, we can derive:
and
Corollary 3. For regular dimension mesh graph (n = m), the products and reciprocals of the sums of Signless Laplacian eigenvalues are defined as
The proof is obvious by Theorem 3.
4 Laplacian spectra and implementations in networking
The framework developed in the previous section allows for the calculation of important network metrics, including graph energy, Kirchhoff index , spectral radius , average path length , global mean first passage time , and the number of spanning trees . To facilitate these computations, two key quantities, and , are introduced. The quantity is defined as the product of all non-zero eigenvalues, denoted as λi, of a given matrix, while is the sum of the reciprocals of these eigenvalues:
Here, λi represents the eigenvalues of the Laplacian matrix associated with the graph , where i ranges from 1 to n. These quantities serve as the basis for further analysis and provide a deeper understanding of various network properties.
4.1 Average path length
The networks with an extremely short mean path length, often referred to as "Small-world" networks, are common in real-world applications. This trait is frequently observed, and various parameters, such as the clustering coefficient, mean path length, and degree distribution, serve as strong indicators of the network's structure. Specifically, for a given mesh graph , the average path length, denoted by , is defined as the average number of steps along the shortest path dij. This metric is essential for measuring the efficiency of material transport or information exchange between all possible node pairs within the network. For the network , is given by:
In an electrical network modeled as a complete graph, there exists a notable connection between the shortest paths and the effective resistance , as detailed in reference [43]:
Here, n represents the order of the complete graph , which is the total number of vertices. By combining the equations above, a simplified expression is derived that reveals the relationship within the graph:
Corollary 4. For a ladder graph, denoted as , the average path length can be derived from the general formula for the mesh graph by setting n = 2. The expression for the average path length of the ladder graph is given by:
4.2 The number of spanning trees
The count of spanning trees () plays a crucial role in various complex network phenomena, including random walks, network reliability, resistor networks, transport systems, loop-erased random walks, and self-organized criticality, as explored in studies like [44–47]. Kirchhoff's Matrix-Tree Theorem, as detailed in [48, 49], reveals a fundamental link by showing that the product of all nonzero eigenvalues of a graph's Laplacian matrix equals the total number of spanning trees. This theorem is a powerful tool for accurately computing for a generalized mesh graph, denoted by . Essentially, this method provides an efficient way to decipher the intricate connections within the graph, greatly aiding in the precise determination of spanning trees across different network configurations:
Corollary 5. For a ladder graph, denoted as , the number of spanning trees can be obtained by setting n = 2 in the general formula for the number of spanning trees of the mesh graph. The expression for the ladder graph is:
4.3 Global mean-first passage time
In network analysis, the global mean-first passage time () is a key metric for assessing the speed of random walks in complex networks, offering insights into how rapidly information or entities travel through the network. It is calculated by averaging individual first passage times over all node pairs. The formula for is:
where Fi, j is the first passage time from node i to node j, and n is the total number of nodes. This average is normalized to include all unique pairs. The commuting time () between nodes i and j is given by:
where Ri, j is a graph-specific metric. For a generalized mesh graph , the global mean-first passage time is computed as:
where n = nm and . Thus, becomes:
Corollary 6. By setting n = 2 in the formula for the global mean-first passage time of a generalized mesh network, we obtain the corresponding result for the ladder graph:
4.4 Spectral radius
The spectral radius is a crucial metric in numerous fields, each leveraging it to gain insights into different systems. In vibration theory, it helps analyze the vibrational patterns of complex systems. Theoretical chemistry uses it to explore molecular structures and interactions, advancing chemical research. In combinatorial optimization, it supports improved decision-making and resource management. Communication networks rely on it to assess data transmission efficiency and reliability, while robustness analysis employs it to test system resilience. Electrical networks utilize the spectral radius to understand component stability and performance. Its versatility and broad application make it an invaluable tool across scientific and engineering disciplines [50, 51].
For adjacency matrices, the spectral radius, denoted as , represents the largest eigenvalue, reflecting the graph's connectivity and dynamics. This value is computed as:
In the context of a generalized mesh graph, can be determined for different types of matrices as follows:
Corollary 7. By setting n = 2 (as n represents the vertical dimension for the ladder graph) in the spectral radius formulas for the generalized mesh graph, we derive the corresponding results for the ladder graph:
4.5 Kirchoff network index
The concept of resistance distance, introduced by Randic and Klein, represents a significant innovation in network analysis. This approach models each edge as a unit resistor, effectively capturing the resistive properties of a network within a graph, denoted as H [52]. In electrical network theory, resistance distance, denoted by dij, measures the effective resistance between nodes i and j. This measurement is derived using Ohm's law. Another key metric is the Kirchhoff index, which sums the resistance distances for all pairs of vertices in the graph G. This index offers a comprehensive view of the network's overall resistance characteristics, providing insights into the electrical connectivity and flow patterns between nodes:
where n is the number of vertices in the graph. The Kirchhoff index can also be expressed in terms of the non-zero eigenvalues λi of the graph:
For a generalized mesh graph , the Kirchhoff index is calculated as:
Corollary 8. By setting n = 2 (as n represents the vertical dimension in the ladder graph) in the Kirchhoff index formula for the generalized mesh graph, we obtain the corresponding result for the ladder graph:
4.6 Graph Energies
Graph energies, such as Laplacian and Randić energy, are essential for understanding graph structures and dynamics. These energies have significant applications in various fields. In network science, they are used to predict robustness, as demonstrated by Li et al. [53, 54]. In molecular graph theory, Wang et al. linked graph energies to molecular stability and reactivity, providing insights into chemistry and drug discovery [55, 56]. In social networks, Chen and Zhang applied these energies to evaluate node importance and information flow [57, 58].
Consider the adjacency matrix of a graph G, denoted by , and let λi represent its eigenvalues derived from the characteristic polynomial. The Adjacency Energy (AE) is expressed as:
Similarly, the Laplacian Energy (LE) and Signless Laplacian Energy (QE) are defined as:
Using these definitions, the energies for a generalized mesh graph can be calculated as follows:
Corollary 9. By setting n = 2 (as the ladder graph has two vertical sides) in the energy formulas for the generalized mesh graph , we derive the corresponding graph energies for the ladder graph:
5 Results and discussions
In this section, we developed a MATLAB Algorithm 1 with a total run time of 0.182 seconds to produce Tables 1, 2. These tables provide exact values for several key metrics: Kirchhoff index , Spectral radius , Average path length , Global mean first passage time , Graph energies , and the number of spanning trees . The algorithm is designed for the generalized mesh network graph . In Table 1, k is set to 2, while p varies from 2 to 15. For Table 2, k is fixed at 3. Exact values for these metrics are calculated to provide a detailed understanding of the network's behavior across different dimensions. In addition to the tables, Figure 2 visually represents the relationships between network size and variations in , , , , , and , further improving the interpretation of the results.

Table 1. Assessment of network-related parameters for the generalized mesh network graph with m set to 2, where 2 ≤ n ≤ 15.

Table 2. Assessment of network-related parameters for the generalized mesh network graph with m set to 3, where 2 ≤ n ≤ 15.

Figure 2. Comparative representation of numeric values evaluated in Tables 1 and 3 for Generalized mesh network graph with m set to 2, where 2 ≤ n ≤ 15. (A) Average Path Length of . (B) Kirchoff Index of . (C) Global passage time of . (D) Graph energies of .
A key observation in the graphical representations is the clear trend indicating that as the network expands, several key metrics increase significantly. These visuals enhance the understanding of the network's behavior, complementing the numerical data and offering a more intuitive grasp of the dynamics within the generalized mesh network. The graphical depiction of the results provides a glimpse into the potential of our methodologies. Researchers are encouraged to utilize our carefully developed algorithm and analysis framework to explore the complexities of more sophisticated real-world networks. The flexibility of our approach offers a valuable toolset, enabling a deeper understanding of network behavior and performance in various contexts. This work paves the way for further studies, serving as a platform for future exploration of complex networks with improved precision and efficiency Algorithm 2).
Table 1 evaluates several network-related parameters for the Generalized Mesh Network Graph across values of n from 2 to 15. The parameters included are the Spectral Radius (), Average Edge Length (), Knot Number (), Average Path Length (), and Global Mean First Passage Time (). As n increases, a notable trend is observed across these parameters. The Spectral Radius () shows a slight increase from 4.013 to 5.963, reflecting a gradual growth in the network's connectivity as more nodes are added. The Average Edge Length () also increases consistently, indicating that as the network grows, the average distance between connected nodes becomes larger. The Knot Number (), which quantifies the number of key nodes in the network, increases significantly, suggesting that more nodes are becoming central as the network expands. The Average Path Length (APL) increases from 0.819 to 21.634, showing that the average distance between any two nodes grows with the size of the network. Finally, the Global Mean First Passage Time () increases from 0.8356 to 4.1385, reflecting that it takes more time on average for a random walker to reach a target node as the network becomes larger. Table 1 provides similar parameters for the Generalized Mesh Network Graph . The parameters assessed are the Spectral Radius (), Average Path Length (APL), Average Edge Length (), Global Mean First Passage Time (), and Knot Number (). Trends in these parameters show a clear pattern of growth and increase with respect to the network size.
The Spectral Radius () increases from 5.025 to 6.9582, indicating a growth in the connectivity strength as the network size increases. The Average Edge Length () increases with network size, which is consistent with the observed trend in Table 3, suggesting that longer edges become more prevalent in larger networks. The Average Path Length (APL) shows a notable increase from 1.45 to 24.458, similar to Table 3, indicating that as the network grows, the average distance between nodes becomes larger. The Global Mean First Passage Time () also increases from 1.1165 to 3.4895, indicating that it takes more time, on average, for a random walker to reach a target node in a larger network. The Knot Number () increases significantly from 14.5 to 2150.9, suggesting a rise in the centrality and importance of nodes within the network as its size expands. Both tables show consistent trends with increasing network size. For , the parameters reflect a gradual increase in spectral radius, average edge length, Kirchoff Index, and average path length, leading to a higher global mean first passage time. Similarly, for , there is a clear upward trend in the spectral radius, average edge length, average path length, and global mean first passage time, with a much more pronounced increase in the Kirchoff Index. These trends illustrate that as the network size increases, the network's complexity grows, resulting in longer paths and higher passage times. The increasing Kirchoff Index indicates more significant central nodes or hubs, which can be crucial for understanding the network's connectivity and efficiency.
6 Conclusion
In summary, this article presents a comprehensive investigation into the spectral properties of the generalized mesh network graph, focusing on adjacency, Laplacian, and signless Laplacian spectra. Through advanced algebraic techniques, we have effectively analyzed these spectral characteristics to derive critical network parameters, including the Kirchhoff index (), Spectral radius (), Average path length (APL), Global mean first passage time (), Graph energies (), and the number of spanning trees (). Our analysis highlights the utility of Laplacian spectra in calculating and understanding various aspects of network behavior. By presenting the results in graphical form, we have provided a clear visualization of how these parameters vary with network dimensions, enhancing our understanding of their interplay and impact. This work not only deepens our insight into the structural and dynamic properties of generalized mesh networks but also offers a robust framework applicable to more complex real-world networks. The methods demonstrated here can be adapted and extended to address specific research needs, facilitating the exploration and optimization of various network configurations. The versatility and precision of these techniques underscore their potential in advancing the study of network systems and their practical applications.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
AR: Conceptualization, Visualization, Writing – original draft. MM: Software, Supervision, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that Gen AI was used in the creation of this manuscript. During the preparation of this work the author(s) used AIL Models in order to improve the English writing as authors belong to non-English region. After using this tool/service, the author(s) reviewed and edited the content as needed and take(s) full responsibility for the content of the publication.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Merris R. Laplacian matrices of graphs: a survey. Linear Algebra Appl. (1994) 197–198:143–76. doi: 10.1016/0024-3795(94)90486-3
2. Hong Y, Zhang XD. Sharp upper and lower bounds for largest eigenvalue of the Laplacian matrices of trees. Discrete Mathemat. (2005) 296:187–97. doi: 10.1016/j.disc.2005.04.001
3. Agaev R, Chebotarev P. On the spectra of nonsymmetric Laplacian matrices. Linear Algebra Appl. (2005) 399:157–68. doi: 10.1016/j.laa.2004.09.003
4. Rojo O, Soto R. The spectra of the adjacency matrix and Laplacian matrix for some balanced trees. Linear Algebra Appl. (2005) 403:97–117. doi: 10.1016/j.laa.2005.01.011
5. Ding X, Jiang T. Spectral distributions of adjacency and Laplacian matrices of random graphs. Ann Appl Probab. (2010) 20:2086–117. doi: 10.1214/10-AAP677
6. Wu CW. Conditions for separability in generalized Laplacian matrices and diagonally dominant matrices as density matrices. Physics Letters A. (2006) 351:18–22. doi: 10.1016/j.physleta.2005.10.049
7. Kaveh A, Rahami H. Block diagonalization of adjacency and Laplacian matrices for graph product; applications in structural mechanics. Int J Numer Methods Eng. (2008) 68:33–63. doi: 10.1002/nme.1696
8. Spielman DA. Algorithms, graph theory, and linear equations in Laplacian matrices. In: Proceedings of the International Congress of Mathematicians (ICM 2010). New Delhi: Hindustan Book Agency (2010). p. 2698–2722.
9. Estrada E. Path Laplacian matrices: Introduction and application to the analysis of consensus in networks. Linear Algebra Appl. (2012) 436:3373–91. doi: 10.1016/j.laa.2011.11.032
10. Krishnan D, Fattal R, Szeliski R. Efficient preconditioning of Laplacian matrices for computer graphics. ACM Trans Graph. (2013) 32:1–15. doi: 10.1145/2461912.2461992
11. Dong X, Thanou D, Frossard P, Vandergheynst P. Laplacian matrix learning for smooth graph signal representation. In: 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP). Piscataway, NJ: IEEE (2015). p. 3736–3740.
12. Pirani M, Sundaram S. On the smallest eigenvalue of grounded Laplacian matrices. IEEE Trans Automat Contr. (2016) 61:509–14. doi: 10.1109/TAC.2015.2444191
13. Bergamaschi L., Martínez A. (2015). Efficiently preconditioned inexact Newton methods for large symmetric eigenvalue problems. Optimiz Methods Softw. 30:301–22. doi: 10.1080/10556788.2014.908878
14. Bandeira AS. Random Laplacian matrices and convex relaxations. Found Comput Mathem. (2018) 18:345–79. doi: 10.1007/s10208-016-9341-9
15. Zhou S, Liu X, Liu J, Guo X, Zhao Y, Zhu E, et al. Multi-view spectral clustering with optimal neighborhood laplacian matrix. In: Proceedings of the AAAI Conference on Artificial Intelligence. Palo Alto, CA: AAAI Press (2020). p. 6965–72.
16. Li W. Quotient with a general orthogonality constraint and an eigen-balanced Laplacian matrix: the greatest lower bound and applications in cooperative control problems. IEEE Trans Automat Contr. (2018) 63:4024–31. doi: 10.1109/TAC.2018.2815179
17. Hermann J, Konigorski U. Eigenvalue assignment for the laplacian matrix of directed graphs. In: Proceedings of the 2019 American Control Conference (ACC). Philadelphia, PA: ACC. (2019). p. 4036–4042.
18. Bergamaschi L, Bozzo E. Computing the smallest Eigenpairs of the graph Laplacian. SeMA J. (2018) 75:1–16. doi: 10.1007/s40324-017-0108-2
19. Jog SR, Kotambari R. On the adjacency, laplacian, and signless laplacian spectrum of coalescence of complete graphs. J Graph Theory. (2016) 82:5906801. doi: 10.1155/2016/5906801
20. Alhevaz A, Baghipur M, Ganie HA, Pirzada S. Brouwer type conjecture for the eigenvalues of distance signless Laplacian matrix of a graph. Linear Multilinea Algebra. (2019) 69:2423–40. doi: 10.1080/03081087.2019.1679074
21. Ganie HA, Shang Y. On the spectral radius and energy of signless Laplacian matrix of digraphs. Heliyon. (2022) 8:e09186. doi: 10.1016/j.heliyon.2022.e09186
22. Morbidi F. Functions of the Laplacian matrix with application to distributed formation control. IEEE Trans Cont Netw Syst. (2022) 9:1459–67. doi: 10.1109/TCNS.2021.3113263
23. Reinhart C. The normalized distance Laplacian. Special Matrices. (2021) 9:1–18. doi: 10.1515/spma-2020-0114
24. Chakrabarty A, Hazra RS, Hollander F, Sfragara M. Spectra of adjacency and Laplacian matrices of inhomogeneous Erdos-Rényi random graphs. Random Matric. (2021) 10:2150009. doi: 10.1142/S201032632150009X
25. Alazemi A, Anđelić M, Koledin T, Stanic Z. Chain graphs with simple Laplacian eigenvalues and their Laplacian dynamics. Computat Appl Mathem. (2022) 42:6. doi: 10.1007/s40314-022-02141-5
26. Anđelić M, da Fonseca CM, Kilic E, Stanic Z. A Sylvester-Kac matrix type and the Laplacian controllability of half graphs. Electron J Linear Algebra. (2022) 2022:38. doi: 10.13001/ela.2022.6947
27. Bapat RB, Jana R, Pati S. The bipartite Laplacian matrix of a nonsingular tree. Special Matrices. (2023) 11:20230102. doi: 10.1515/spma-2023-0102
28. Anđelić M, Brunetti M, Stanic Z. Laplacian controllability for graphs obtained by some standard products. Graphs Combinat. (2020) 36:1593–602. doi: 10.1007/s00373-020-02212-6
29. Anđelić M, da Fonseca CM, Simić SK, Toć DV. Connected graphs of fixed order and size with maximal Q-index: some spectral bounds. Discrete Appl Mathemat. (2012) 160:448–59. doi: 10.1016/j.dam.2011.11.001
30. Mallik S. Matrix tree theorem for the net Laplacian matrix of a signed graph. Linear Multilinear Algebra. (2024) 72:1138–52. doi: 10.1080/03081087.2023.2172544
31. Raza A, Munir M, Abbas T, Eldin SM, Khan I. Spectrum of prism graph and relation with network related quantities. AIMS Math. (2023) 8:2634–47. doi: 10.3934/math.2023137
32. Raza A, Munir MM. Insights into network properties: spectrum-based analysis with Laplacian and signless Laplacian spectra. Eur Physi J Plus. (2023) 138:802. doi: 10.1140/epjp/s13360-023-04441-z
33. Raza A, Munir M, Hussain M, Tolasa FT. A spectrum-based approach to network analysis utilizing Laplacian and signless laplacian spectra to torus networks. IEEE Access. (2024) 12:52016–52029. doi: 10.1109/ACCESS.2024.3384300
34. Harel D, Sardashti A. Structural characteristics of grid graphs. J Graph Theory. (1998) 27:133–48.
35. Hinz AM, Holz auf der Heide C. An efficient algorithm to determine all shortest paths in Sierpinski graphs. Discrete Appl Math. (2014) 177:111–20. doi: 10.1016/j.dam.2014.05.049
36. Yang Y, Liu H, Hou J, Zhang N, Wu T. Data collection of wireless sensor network with grid structure based on compressed sensing. In: 2021 IEEE 4th International Conference on Electronic Information and Communication Technology (ICEICT), Xi'an, China (2021) p. 363–8. doi: 10.1109/ICEICT53123.2021.9531288
37. Fang X, Misra S, Xue G, Yang D. A survey of smart grid architectures, applications, benefits and standardization. J Netw Comput Appl. (2012) 35:1950–62. doi: 10.1016/j.jnca.2012.07.002
38. Trinajstic N, Babic D, Nikolic S, Plavšic D, Amic D, Mihalic Z. The Laplacian matrix in chemistry. J Chem Inform Comput Sci. (1994) 34:368–76. doi: 10.1021/ci00018a023
39. Chu Z, Munir M, Yousaf A, Qureshi MI, Liu J. Laplacian and signless Laplacian spectra and energies of multi-step wheel networks. Mathem Biosci Eng. (2020) 17:3649–59. doi: 10.3934/mbe.2020206
40. Xu Y, Feng Z, Qi X. Signless-Laplacian eigenvector centrality: a novel vital nodes identification method for complex networks. Pattern Recognit. Letters. (2021) 148:7–14. doi: 10.1016/j.patrec.2021.04.018
41. Pecora LM, Sorrentino F, Hagerstrom AM, Murphy TE, Roy R. How do the eigenvalues of the Laplacian matrix affect route to synchronization patterns. Phys Lett A. (2014) 378:2590–5. doi: 10.1016/j.physleta.2014.07.019
42. Mohar B. Some applications of Laplace eigenvalues of graphs. In: Hahn G, Sabidussi G, , editors. Graph Symmetry. NATO ASI Series, Vol. 497. Dordrecht: Springer (1997). doi: 10.1007/978-94-015-8937-6_6
43. Lukovits I, Nikolić S, Trinajstić N. Resistance distance in regular graphs. Int J Quant Chem. (1999) 71:217–225.
44. Szabó GJ, Alava M, Kertész J. Geometry of minimum spanning trees on scale-free networks. Physica A. (2003) 330:31–6. doi: 10.1016/j.physa.2003.08.031
45. Wu ZH, Braunstein LA, Havlin S, Stanley HE. Transport in weighted networks: partition into superhighways and roads, Phys. Rev Lett. (2006) 96:148702. doi: 10.1103/PhysRevLett.96.148702
46. Dhar D. Theoretical studies of self-organized criticality. Physica A. (2006) 369:29–70. doi: 10.1016/j.physa.2006.04.004
47. Dhar D, Dhar A. Distribution of sizes of erased loops for loop-erased random walks, Phys. Rev E. (1997) 55:R2093. doi: 10.1103/PhysRevE.55.R2093
48. Zhang ZZ, Wu B, Comellas F. The number of spanning trees in Apollonian networks, Discrete Appl. Math. (2014) 169:206–13. doi: 10.1016/j.dam.2014.01.015
49. Godsil C, Royle G. Algebraic Graph Theory, Graduate Texts in Mathematics. New York: Springer. (2001). p. 207.
50. Liu P, Liu G, Lv H. Power function method for finding the spectral radius of weakly irreducible nonnegative tensors. Symmetry. (2022) 14:2157. doi: 10.3390/sym14102157
51. Jiang Z, Li L, Xu Y. Spectral Radius Optimization for Neural Networks. Neural Networks. (2020) 122:359–69.
53. Li S, Zhang Q, Wang Y, Li X. Graph energies and network robustness. J Complex Netw. (2022) 1–15.
54. Lee JR, Hussain A, Fahad A, Raza A, Qureshi MI, et al. On ev and ve-Degree Based Topological Indices of Silicon Carbides. CMES. (2022) 130:871–85. doi: 10.32604/cmes.2022.016836
55. Fowler PW. Energies of graphs and molecules. AIP Conf Proc. (2007) 963:517–20. doi: 10.1063/1.2836127
56. Zhang X, Raza A, Fahad A, Jamil MK. (2020). On face index of silicon carbides. Discrete Dyn Nat Soc. 2020:6048438. doi: 10.1155/2020/6048438
57. Truong QD, Truong QB, Dkaki T. Graph methods for social network analysis. In: Vinh P, Barolli L, , editors. Nature of Computation and Communication. ICTCC 2016. Lecture Notes of the Institute for Computer Sciences, Social Informatics and Telecommunications Engineering, Vol. 168. Cham: Springer (2016). doi: 10.1007/978-3-319-46909-6_25
58. Alghazzawi D, Raza A, Munir U, Ali S. Chemical applicability of newly introduced topological invariants and their relation with polycyclic compounds. J Mathemat. (2022) 2022:5867040. doi: 10.1155/2022/5867040
59. Ramana KVM, Shireesha R. Applications of spectral graph theory in machine learning and data science. Global J Eng Innov Interdiscipl Res. (2025) 5:1–6. doi: 10.33425/3066-1226.1070
60. Morzy M, Kajdanowicz T. Graph energies of egocentric networks and their correlation with vertex centrality measures. Entropy. (2018) 20:916. doi: 10.3390/e20120916
Keywords: Laplacian spectrum, spectral radius, Kirchhoff index, network stability, first passage time
Citation: Raza A and Mobeen Munir M (2025) Laplacian spectra and structural insights: applications in chemistry and network science. Front. Appl. Math. Stat. 11:1519577. doi: 10.3389/fams.2025.1519577
Received: 01 November 2024; Accepted: 05 May 2025;
Published: 13 June 2025.
Edited by:
Gang Ren, The Molecular Foundry, Berkeley Lab (DOE), United StatesReviewed by:
Ahmad Qazza, Zarqa University, JordanNadeem Ur Rehman, Aligarh Muslim University, India
Copyright © 2025 Raza and Mobeen Munir. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ali Raza, YWxsZWVyYXp6YTc4NkBnbWFpbC5jb20=