Abstract
Topological characterization through graph-theoretical methods translates chemical and structural data into quantitative values that represent the molecular system. Our research explores the use of topological indices to study fractal structures. Molecular fractals are complex geometric configurations that exhibit self-similarity at different levels and systematically formed by repeating a fundamental unit. This study focuses on coronene-based molecular fractals, where coronene, a benzenoid molecule with a symmetrical graphite-like structure, finds applications in organic semiconductors, sensors, and molecular electronics, due to its unique electronic and optical properties. Additionally, information entropy is employed to evaluate and compare the structural complexities of coronene fractals. Spectra-based energetic properties such as total -electron energy, HOMO-LUMO energy gaps, spectral diameter, delocalization and resonance energies are calculated to assess their kinetic and thermodynamic stability. Furthermore, predictive models are provided for estimating spectral characteristics across higher-dimensional coronene fractal structures.
1 Introduction
Benzenoid hydrocarbons are a group of polycyclic compounds consisting of six-member linked rings, characterized by their aroma and unique physicochemical properties. These substances create powerful inter-molecular bonds by acting in single and double bonds alternatively (Hill et al., 2004). Higher-order structural co-ordination is indicated by larger -conjugated complexes. These characteristics primarily make them useful for applications in opto-electronic devices, nanomaterials, and natural semiconductors (Pisula et al., 2010; Pisula et al., 2011). Coronene, a planar molecule with seven peri-fused benzene rings, is well-known for having delocalized -electrons, extended conjugation, and extreme symmetry (Newman, 1940; Robertson and White, 1945; Popov and Boldyrev, 2012). It serves as a fundamental polycyclic aromatic hydrocarbon (PAH) model for studying larger PAHs, graphene quantum dots, and graphene nanoflakes. Coronene-based structures enable precise theoretical investigations and bridge PAHs with graphene materials (Santa Daría et al., 2024; Tachikawa and Lund, 2022). It has a well-defined structure, fluorescence, and electronic properties which makes it a benchmark in theoretical and experimental studies. Coronene fractals exhibit exceptional electronic, optical, and energy-related properties, with strong -electron delocalization enhancing charge transport and stable -conjugation improving the performance of capacitors and batteries (Demir and Üngördü, 2023; Sanyal et al., 2013; Dobrowolski et al., 2011). Molecular stacking and employers are further enhanced by its symmetrical and planar architecture (Fedotov et al., 2013).
Fractal geometry, which explores recurring patterns at different scales, has evolved from describing physical theories to serve various applications such as complexes in medical and molecular engineering, neural networks, and laptop graphics, etc (Kirkby, 1983). Extensive research has been carried out using fractal methods. These deterministic fractals arise by combining benzene with hierarchical structure sequences, making them a significant tool for advancing nanotechnology and biotechnology (Uahabi and Atounti, 2015). Clar aromatic sextet theory is a concept introduced by Erich Clar to describe the electronic structure of polycyclic aromatic hydrocarbons (PAHs). It is particularly useful for understanding resonance, stability, and reactivity in PAH systems (Hosoya, 2005). Fractal molecular architecture, often analyzed through Clar’s system and golden ratio measurements, exhibits scaling properties that demonstrate its adaptability and potential (Lee and Chang, 1996). Studies of coronene-based fractals have shown that they can serve as supports for advanced nanomaterials (Nisha and Senthil Kumar, 2020). Despite significant advances in theoretical research, the integration of these complex systems remains a challenge, requiring further research (Kumar et al., 2017). Recent work emphasis on the unique aromatic properties and scaling behavior of fractal benzenoids, emphasizes their importance in development and fabrication of high-performance nanomaterials for optical and electronic device applications (Duan et al., 2021).
In computational chemistry, topological indices considered as are important tools that provide information on the chemical and structural characteristics of molecules (Estrada and Uriarte, 2001; Kumar and Das, 2024). Among these, the degree-based Zagreb index and the distance-based Wiener index have been crucial in forecasting molecular characteristics, including stability and boiling points (Wiener, 1947; Gutman and Trinajstić, 1972). In this article, we utilize modified reverse degree-based indices that incorporate a variable parameter, “,” which potentially alters the graph’s degree sequence. Unlike traditional methods with fixed-degree sequences, this approach allows customization of the “,” value to better correlate with specific datasets and their properties. This method is not limited to specific indices and can be applied to all degree-based indices. Notably, as the “,” value increases, these modified indices exhibit a high correlation with the physicochemical characteristics of corona, blood cancer, and heart disease treatment drug molecules (Arockiaraj et al., 2023a; Arockiaraj et al., 2023b; Arockiaraj et al., 2024). In addition, they are used for stability analysis in advanced materials like carbon nanosheets, metal-organic frameworks, and pent-heptagonal nanostructures (Abul Kalaam and Berin Greeni, 2024). Further, employing hybrid models allows for more precise predictions of molecular activity (Arockiaraj et al., 2023c).
Entropy analysis is a fundamental method in the field of information theory, which offers special insights into the complexity and stability of molecules. Shannon’s entropy measures structural randomness (Dehmer, 2008; Shannon, 1948), while graph entropy is related to the vertices and edges of molecular graphs, which makes it easier to analyze a system using graph structures. Higher entropy of a structure constitutes more disorderness in the macrostructure, which reduces structural stability. However, high entropy materials, such as high-entropy alloys (HEAs), exhibit unique properties due to their high configurational entropy, which can result in the formation of stable disordered solid solutions. While high entropy promotes disorderness, it can also contribute to distinctive structural stabilities and desirable properties. For instance, HEAs are known for their high strength, ductility, and resistance to wear and corrosion.
Research articles focused on molecular fractals have explored various structural and topological aspects (Malik et al., 2023; Xu and Liu, 2025; Yogalakshmi and Easwaramoorthy, 2024). Recent studies on coronene fractals have examined degree and degree-sum properties, reverse degree-based indices, and coronene frameworks, as discussed in (Arockiaraj et al., 2022; Ullah et al., 2024; Khabyah et al., 2023). This study explores coronene fractal structures, analyzing their entropy levels and complexity through modified reverse degree-based indices. By delving into their structural and spectral features, it aims to deepen our understanding of their stability, complexity, and overall properties.
2 Methodology
In this study, we examine three configurations of coronene fractals modeled as two-dimensional molecular graph structure and is represented by , with and denote the number of vertices and edges, respectively. The degree of a vertex denoted as , indicates the number of vertices directly connected to . The maximum degree, , represents the highest connectivity among all vertices in the graph . Recent modification in reverse degree is done by introducing a parameter (with ) that enhance the graph degree sequence to closely predict properties (Arockiaraj et al., 2023a). The modified reverse degree, represented , is defined as follows:
Modified reverse degree-based topological indices, , are employed to characterize coronene fractals by evaluating atom connectivity and providing insights into their molecular structure. For a graph , the formulation of is given as:where represents the edge connecting vertices and . This formula provides an intensive evaluation by using the contributions of all edges inside the graph. The edge set is divided into equivalent subsets, such that . Each subset of , where and , groups edges based on vertex connectivity in . For any subset , the corresponding is calculated as:where represents the number of edges in subset , and evaluates the contribution of modified reverse degree for the connected vertices.
The total for graph is obtained by summing the contributions from all subsets :
The topological index functions based on the modified reverse degree are outlined below.
• Modified reverse first Zagreb index :
• Modified reverse second Zagreb index :
• Modified reverse forgotten index :
• Modified reverse Sombor index :
• Modified reverse geometric arithmetic index :
• Modified reverse hyper-Zagreb index :
• Modified reverse harmonic index :
• Modified reverse first redefined Zagreb index :
• Modified reverse second redefined Zagreb index :
• Modified reverse bi-Zagreb index :
• Modified reverse tri-Zagreb index :
• Modified reverse geometric bi-Zagreb index :
3 Evaluation of modified reverse degree indices
We explore three coronene fractal configurations: ZHCF, AHCF, and RCF, as illustrated in Figures 1, 2. The structural parameters for these configurations are given by: and ; for AHCF, these are and ; and for RCF, they are and . All configurations share a maximum vertex degree of 3. The modified reverse degree metrics for each vertex are as follows:
FIGURE 1

The degree based enumeration of coronene fractal structure.
FIGURE 2

Configurations of coronene fractals (a) ZHCF(3) (b) AHCF(2) (c) RCF(5,3).
To calculate the modified reverse topological indices, edge partitioning is used, as illustrated in Figure 1, for three configurations of coronene fractals based on their standard vertex degrees, as detailed in Table 1. Each index involves complex computations with varying parameters. For instance, the calculation of the first Zagreb-based index is demonstrated using ZHCF coronene structures for different values of . When the variable parameter is set to , the degree pairs (2,2), (2,3), and (3,3) are modified to (2,2), (2,1), and (1,1), respectively. Therefore,
TABLE 1
| Bond type | Coronene fractals | ||||
|---|---|---|---|---|---|
| C-C | 2 | 2 | |||
| C-C | 2 | 3 | |||
| C-C | 3 | 3 | |||
Degree based edge partition of three configurations of coronene fractals.
For the degree classes are modified into (3,3), (3,2), and (2,2). Therefore,
Similarly for the degree classes are modified into (1,1), (1,3), and (3,3). Therefore,
The modified reverse degree-based indices illustrated in Equations 1–12, combined with the edge partitioning present in Table 1, are employed to compute the for three configurations of coronene fractals. The results, corresponding to the variable parameters , and 3, are summarized in Tables 2–4.
TABLE 2
| Zigzag hexagonal coronene fractal structure | |||
|---|---|---|---|
Modified reverse degree indices of ZHCF structure for variable parameters = 1, 2, and 3.
TABLE 3
| Armchair hexagonal coronene fractal structures | |||
|---|---|---|---|
Modified reverse degree indices of AHCF structure for variable parameters = 1, 2, and 3.
TABLE 4
| Rectangular coronene fractal structures | |||
|---|---|---|---|
Modified reverse degree indices of RCF structure for variable parameters = 1, 2, and 3.
4 Evaluation of graph entropy
A mathematical foundation for assessing a system’s randomness and uncertainty is provided by Shannon’s concept of entropy, which quantifies the content of possibility distributions. For a discrete random variable with chances , thus Shannon’s entropy , is expressed as:where , represents the frequency of a specific outcome , and is the total number of outcomes (Shannon, 1948). The information obtained from measuring the system is captured by using the logarithm base-2 to validate the entropy values in bits. This equation has a significant analogy to thermodynamic entropy, which quantifies the randomness of states in a physical system (Mowshowitz and Dehmer, 2012; Sabirov and O-sawa, 2015). In Physics, thermodynamic entropy is used to assess microstates, while Shannon’s entropy is generally applicable to abstract systems, such as graphs, and allows for the analysis of their structural complexity using attributes like vertices and edges (Mowshowitz, 1968).
Based on this foundation, incorporating topological indices into the entropy framework appears as a robust approach to assess molecular complexity. This approach focuses on graph edges and uses topological indices , which are mathematical structural characterization of molecular graphs (Arockiaraj et al., 2023d). The probability given to each edge of a molecular graph , with edges is defined as , where is a modified reverse degree-based function and is the associated index. The graph entropy is expressed as:
Further the graph entropy equation simplifies as:
By employing specific topological indices, the simplified representation makes it easier to calculate graph entropy for a molecular graphs of coronene fractals. For example, the modified first Zagreb index applied to a ZHCF structure. The result of substituting into the entropy formula is:
Employing degree-based edge partitions presented in Table 1, the entropy of ZHCF when and is:
The entropy expressions for all configuration of coronene fractals are too extensive to display. Therefore, Tables 5–7 present the comparison of numerical values of modified reverse degree-based entropy levels for the fractal structures. For rectangular coronene fractals, we assume .
TABLE 5
| Zigzag hexagonal coronene fractal structure | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 9.3772 | 10.5436 | 11.3720 | 12.0148 | 9.4136 | 10.5797 | 11.4078 | 12.0505 | 9.3651 | 10.5332 | 11.3624 | 12.0057 | |
| 9.2082 | 10.3750 | 11.2037 | 11.8468 | 9.3605 | 10.5271 | 11.3556 | 11.9986 | 9.2105 | 10.3858 | 11.2186 | 11.8641 | |
| 9.1938 | 10.3591 | 11.1871 | 11.8297 | 9.3556 | 10.5221 | 11.3505 | 11.9934 | 9.2797 | 10.4518 | 11.2830 | 11.9276 | |
| 9.3734 | 10.5397 | 11.3680 | 12.0109 | 9.4128 | 10.5788 | 11.4070 | 12.0497 | 9.3724 | 10.5404 | 11.3695 | 12.0127 | |
| 9.4300 | 10.5958 | 11.4237 | 12.0663 | 9.4304 | 10.5961 | 11.4241 | 12.0667 | 9.4278 | 10.5937 | 11.4217 | 12.0643 | |
| 9.3988 | 10.5656 | 11.3941 | 12.0370 | 9.4181 | 10.5842 | 11.4124 | 12.0551 | 9.2860 | 10.4524 | 11.2808 | 11.9237 | |
| 9.3825 | 10.5491 | 11.3776 | 12.0205 | 9.4149 | 10.5810 | 11.4091 | 12.0518 | 9.3367 | 10.5058 | 11.3356 | 11.9792 | |
| 9.3911 | 10.5581 | 11.3867 | 12.0297 | 9.4165 | 10.5827 | 11.4108 | 12.0536 | 9.2991 | 10.4667 | 11.2957 | 11.9390 | |
| 9.2027 | 10.3687 | 11.1970 | 11.8399 | 9.3583 | 10.5248 | 11.3533 | 11.9962 | 9.2519 | 10.4253 | 11.2571 | 11.9020 | |
| 9.3240 | 10.4908 | 11.3194 | 11.9624 | 9.3896 | 10.5560 | 11.3843 | 12.0272 | 9.2934 | 10.4646 | 11.2953 | 11.9396 | |
| 9.2001 | 10.3658 | 11.1940 | 11.8368 | 9.3574 | 10.5240 | 11.3524 | 11.9953 | 9.2625 | 10.4353 | 11.2669 | 11.9117 | |
| 9.4230 | 10.5890 | 11.4171 | 12.0597 | 9.4262 | 10.5920 | 11.4201 | 12.0627 | 9.4056 | 10.5718 | 11.4001 | 12.0428 | |
Comparison of entropy levels for ZHCF at , , and .
TABLE 6
| Armchair hexagonal coronene fractal structures | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 10.1830 | 11.6202 | 12.5802 | 13.3006 | 10.2193 | 11.6560 | 12.6158 | 13.3362 | 10.1716 | 11.6105 | 12.5712 | 13.2921 | |
| 10.0142 | 11.4519 | 12.4122 | 13.1329 | 10.1664 | 11.6038 | 12.5640 | 13.2845 | 10.0201 | 11.4661 | 12.4301 | 13.1528 | |
| 9.9991 | 11.4353 | 12.3951 | 13.1154 | 10.1615 | 11.5987 | 12.5588 | 13.2793 | 10.0880 | 11.5308 | 12.4933 | 13.2152 | |
| 10.1792 | 11.6162 | 12.5762 | 13.2967 | 10.2185 | 11.6552 | 12.6150 | 13.3353 | 10.1789 | 11.6175 | 12.5782 | 13.2990 | |
| 10.2356 | 11.6720 | 12.6317 | 13.3519 | 10.2360 | 11.6724 | 12.6320 | 13.3523 | 10.2334 | 11.6699 | 12.6296 | 13.3499 | |
| 10.2048 | 11.6422 | 12.6024 | 13.3229 | 10.2238 | 11.6606 | 12.6204 | 13.3408 | 10.0918 | 11.5290 | 12.4891 | 13.2096 | |
| 10.1885 | 11.6258 | 12.5859 | 13.3063 | 10.2206 | 11.6573 | 12.6172 | 13.3375 | 10.1437 | 11.5835 | 12.5448 | 13.2659 | |
| 10.1972 | 11.6348 | 12.5951 | 13.3157 | 10.2222 | 11.6591 | 12.6189 | 13.3393 | 10.1054 | 11.5438 | 12.5045 | 13.2253 | |
| 10.0083 | 11.4452 | 12.4053 | 13.1258 | 10.1642 | 11.6015 | 12.5616 | 13.2821 | 10.0607 | 11.5047 | 12.4678 | 13.1900 | |
| 10.1300 | 11.5676 | 12.5279 | 13.2484 | 10.1955 | 11.6326 | 12.5925 | 13.3130 | 10.1012 | 11.5431 | 12.5053 | 13.2269 | |
| 10.0056 | 11.4422 | 12.4022 | 13.1226 | 10.1633 | 11.6006 | 12.5607 | 13.2812 | 10.0711 | 11.5146 | 12.4775 | 13.1995 | |
| 10.2287 | 11.6653 | 12.6251 | 13.3453 | 10.2318 | 11.6683 | 12.6280 | 13.3483 | 10.2114 | 11.6483 | 12.6082 | 13.3286 | |
Comparison of entropy levels for AHCF at , , and .
TABLE 7
| Rectangular coronene fractal structure | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 9.1164 | 10.1830 | 10.9586 | 11.5688 | 9.1529 | 10.2193 | 10.9947 | 11.6047 | 9.1031 | 10.1716 | 10.9482 | 11.5590 | |
| 8.9473 | 10.0142 | 10.7901 | 11.4005 | 9.0995 | 10.1664 | 10.9422 | 11.5524 | 8.9443 | 10.0201 | 10.8008 | 11.4143 | |
| 8.9337 | 9.9991 | 10.7741 | 11.3840 | 9.0947 | 10.1615 | 10.9371 | 11.5474 | 9.0154 | 10.0880 | 10.8669 | 11.4791 | |
| 9.1126 | 10.1792 | 10.9548 | 11.5649 | 9.1521 | 10.2185 | 10.9939 | 11.6039 | 9.1106 | 10.1789 | 10.9554 | 11.5661 | |
| 9.1695 | 10.2356 | 11.0108 | 11.6207 | 9.1699 | 10.2360 | 11.0112 | 11.6211 | 9.1673 | 10.2334 | 11.0087 | 11.6186 | |
| 9.1376 | 10.2048 | 10.9806 | 11.5909 | 9.1574 | 10.2238 | 10.9992 | 11.6093 | 9.0253 | 10.0918 | 10.8674 | 11.4776 | |
| 9.1216 | 10.1885 | 10.9642 | 11.5744 | 9.1542 | 10.2206 | 10.9960 | 11.6060 | 9.0742 | 10.1437 | 10.9209 | 11.5320 | |
| 9.1298 | 10.1972 | 10.9731 | 11.5835 | 9.1558 | 10.2222 | 10.9977 | 11.6078 | 9.0376 | 10.1054 | 10.8818 | 11.4924 | |
| 8.9422 | 10.0083 | 10.7838 | 11.3939 | 9.0973 | 10.1642 | 10.9399 | 11.5501 | 8.9869 | 10.0607 | 10.8403 | 11.4530 | |
| 9.0630 | 10.1300 | 10.9059 | 11.5162 | 9.1288 | 10.1955 | 10.9711 | 11.5812 | 9.0296 | 10.1012 | 10.8796 | 11.4915 | |
| 8.9397 | 10.0056 | 10.7809 | 11.3909 | 9.0965 | 10.1633 | 10.9390 | 11.5492 | 8.9977 | 10.0711 | 10.8504 | 11.4629 | |
| 9.1624 | 10.2287 | 11.0040 | 11.6140 | 9.1656 | 10.2318 | 11.0071 | 11.6170 | 9.1448 | 10.2114 | 10.9869 | 11.5970 | |
Comparison of entropy levels for RCF where at , , and .
The entropy stages provided in Tables 5–7 monitor dynamic variation throughout the three configurations of coronene fractals for , and 3. Notably, entropy values continually peaks at in comparison to and . The entropy measures differ slightly in their decimal values across all indices. Among the configurations, AHCF demonstrates slightly higher entropy values than the other coronene structures, while RCF exhibits lower entropy values, indicating greater structural stability. However, direct comparisons of complexity measures across these fractal structure are complicated by differences in the number of edges. We utilize relative measures, including structural information content (SIC) and bond information content (BIC), derived from the computed entropy values. These metrics provide a exact evaluation of the structural complexity and stability of the three coronene fractal configurations.
4.1 Relative complexity metrics
This subsection offers numerical and graphical estimation of complexity across the configurations of coronene fractals, emphasizing the importance of accounting for molecular size differences. Since graph entropy values are depending on the size of the molecular graph, the application of relative complexity measures has become essential for higher comparisons among molecular systems of varying dimensions (Dehmer, 2008). To address this, two normalized measures, namely structural information content (SIC) and bond information content (BIC), are introduced. Graph entropy alone may not adequately reflect structural complexity, especially for systems with differing dimensional sizes, highlighting the necessity of employing relative metrics (Bonchev and Trinajstić, 1982; Sabirov and Shepelevich, 2021). The maximum entropy concept is used to establish these metrics, where the limiting entropy value for is defined as (Junias et al., 2024). This leads to SIC, which quantifies molecular structure and the most useful information, as shown below:
Similarly, BIC includes a molecular graph where edges are counted to compute relative complexity. The formula for the BIC normalizes the entropy using the logarithmic scale of the total number of edges, as shown here:
From Equations 13, 14, we calculate the SIC and BIC measures for coronene fractals. The analysis focuses on the entropy values of the Zagreb index when , where providing insight into the relative complexity assessment between the fractals. For example the ZHCF(3) system with , the Zagreb index value and . The values calculated by equations for SIC and BIC are as follows.
The SIC and BIC measures for other coronene fractals across various vertex ranges, are presented in Table 8. These relative complexity measures offer a comparative analysis of complexity across different sizes, with values ranging from 0 to 1, where 1 indicates the highest complexity and 0 the lowest. The SIC and BIC measures, are shown in Table 8, with a graphical comparison in Figure 3.
TABLE 8
| Vertex ranges | Structures | SIC | BIC | |||
|---|---|---|---|---|---|---|
| 900–1152 | ZHCF(3) | 10.57965 | 12.75489 | 10.59619 | 0.829459 | 0.998439 |
| AHCF(2) | 10.2193 | 12.39874 | 10.23601 | 0.824221 | 0.998367 | |
| RCF(3,3) | 10.2193 | 12.39874 | 10.23601 | 0.824221 | 0.998367 | |
| 2040–2424 | ZHCF(4) | 11.40779 | 13.57932 | 12.06676 | 0.840086 | 0.94539 |
| AHCF(3) | 11.65602 | 13.82814 | 11.67243 | 0.842921 | 0.998595 | |
| RCF(9,3) | 11.65602 | 13.82814 | 11.67243 | 0.842921 | 0.998595 | |
| 4572–4872 | ZHCF(6) | 12.57579 | 14.74357 | 12.59199 | 0.852967 | 0.998713 |
| AHCF(4) | 12.61584 | 14.78463 | 12.63209 | 0.853307 | 0.998714 | |
| RCF(9,6) | 12.6542 | 14.82277 | 12.67043 | 0.853700 | 0.998718 |
Relative complexity measures of three classes of coronene structures.
FIGURE 3

Graphical comparison of complexity measures among zigzag, armchair and rectangular coronene fractals. (a) Comparison of SIC measures across the range of k values for coronene fractals. (b) Comparison of BIC measures across the range of k values for coronene fractals.
Table 8; Figure 3 show that RCF and AHCF exhibit similar complexity values at small scales. However, with increasing size, the rectangular fractals exhibits slightly higher complexity compared to armchair configuration, while the zigzag-based coronene exhibits lower complexity than all other configurations because BIC is evaluated based on number of bonds in molecular graph and SIC obtained from maximum entropy. These two analyzes facilitate better comparisons, and help to determine the most appropriate indicator of complexity measures for molecular system.
From Figure 3; Tables 5–8, greater entropy variations are observed among the three configurations for smaller structures, while for the largest structure, all configurations approach the 2D graphitic sheet, and their entropy values converge to a limit. However, two types of isentropic structures exist: AHCF(2) and RCF(3,3) have the same number of vertices (900) and edges (1206); similarly, AHCF(3) and RCF(9,3) share the same number of vertices (2424) and edges (3264). Thus, we use spectral properties for a more conclusive analysis of stability.
5 Analysis of spectral properties in coronene fractals
This section focuses on the spectral properties of coronene fractals, using metrics derived from their graph spectra. Since these structures are two-dimensional and satisfy the Coulson-Rushbrook theorem, this method is more effective for analysis. It is not practical to perform complete calculations for complex systems such as AHCF(3) and RCF(9,3) with 2424 vertices and spectral eigenvalues (Arockiaraj et al., 2022). Consequently, machine learning techniques are needed to efficiently estimate stability in large, fractal structures. Significant spectal and energy properties such as total -electron energy, spectral diameter, HOMO-LUMO energy gap, delocalization energy, and resonance energy, are determined by combinatorial analysis of graph spectra (Prabhu et al., 2024). These parameters provide valuable insights into the thermodynamic and kinetic stability of the coronene fractals under investigation.
The total -electron energy is a critical measure of electronic stability in conjugated systems (Gutman and Trinajstić, 1972; Gutman, 1978). For coronene fractals, including zigzag, armchair, and rectangular patterns, is calculated using the eigenvalues of the graph spectra of the molecular graph (Kalaam et al., 2024; Graovac et al., 1977). For a system with atoms:
The -electron distribution depends on whether is even or odd.
The HOMO-LUMO energy gaps, defined as the difference between the highest molecular orbital (HOMO) denoted and the lowest unoccupied molecular orbital (LUMO) denoted . It plays an important role in analyzing molecular reactivity and kinetic stability. These difference is calculated by subtracting the lowest positive eigenvalue from the highest negative eigenvalue from the graph spectrum, expressed as (Wu et al., 2018; Li et al., 2013). Larger HOMO-LUMO energy differences indicate increased kinetic stability and low chemical reactivity, as more energy is required to transfer an electron from HOMO to LUMO, thus decreasing the chemical reactivity however this difference does not directly reflect thermodynamic stability.
Thermodynamic stability is closely related to parameters such as delocalization and resonance energies, which generally increase with molecule size, increasing the stability. The delocalization energy , is calculated as . Kekul counts (KC), which reflect the number of Kekul resonance structures in coronene fractals, are used to compute resonance energies as coronene fractals are benzenoid systems and bipartite graphs. Thus (KC) is derived from the square root of the constant term of the characteristic polynomial (Balasubramanian, 2023). According to Herndon’s definition of resonance, (Herndon and Ellzey, 1974). The increase in size of the coronene fractals increases both delocalization and resonance energies. Because of the stabilization of the molecular orbitals (Mazouin et al., 2022), the HOMO-LUMO energy gap decreases with increasing molecular size. The spectral diameter is calculated as the difference between the maximum and minimum eigenvalues: . These graph spectra based energy properties were assessed using programs such as newGRAPH and MATLAB software (Stevanović et al., 2021; MATLAB, 2022). Table 9 displays the results, which are given in units.
TABLE 9
| Structure | gaps | |||||
|---|---|---|---|---|---|---|
| ZHCF(1) | 194.662 | 0.7638 | 1.475 | 0.474712 | 0.156614 | 5.5937 |
| ZHCF(2) | 765.926 | 0.65732 | 1.484 | 0.484353 | 0.159251 | 5.64926 |
| ZHCF(3) | 1713.792 | 0.62768 | 1.488 | 0.487666 | 0.160157 | 5.65958 |
| ZHCF(4) | 3038.259 | 0.61518 | 1.489 | 0.489343 | 0.160615 | 5.66336 |
| AHCF(1) | 194.662 | 0.7638 | 1.475 | 0.474712 | 0.156614 | 5.5937 |
| AHCF(2) | 1337.189 | 0.63774 | 1.486 | 0.485766 | 0.159636 | 5.65682 |
| AHCF(3) | 3609.522 | 0.61322 | 1.489 | 0.489077 | 0.160541 | 5.66408 |
| RCF(1) | 194.662 | 0.7638 | 1.475 | 0.474712 | 0.156614 | 5.5937 |
| RCF(2) | 640.3919 | 0.66974 | 1.482 | 0.482389 | 0.158714 | 5.64364 |
| RCF(3) | 1337.19 | 0.6379 | 1.486 | 0.485766 | 0.159637 | 5.65572 |
| RCF(4) | 2285.055 | 0.6229 | 1.488 | 0.487666 | 0.160156 | 5.66072 |
Energetic properties of three classes of polycyclic aromatic hydrocarbons.
The data present in Table 9; Figure 4 show that the HOMO-LUMO energy gaps decrease as the size of coronene structures increase. This suggests that larger structures have more electronic delocalization and resonance energy, which results in lower energy differences between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). Meanwhile, both delocalization and resonance energy show an increasing trend, reflecting enhanced stability and conjugation within these structures. Among the fractal configurations analyzed, the rectangular coronene fractals have the largest HOMO-LUMO energy gaps, suggesting high kinetic stability, lower reactivity, and the lowest delocalization and resonance energies. On the other hand, armchair coronene fractals display the smallest HOMO-LUMO energy gaps, indicating less kinetic stability, larger chemical reactivity, higher electron delocalization, and resonance energies, all of which lead to greater stability with efficient electron transfer. This study emphasizes the significance of structural configuration on stability and reactivity in coronene fractals.
FIGURE 4

Graphical representation of energetic properties across kekulene tessellations. (a) HOMO-LUMO energy gap (b) Energy per bond (c) Delocalization energy per bond (d) Resonance energy per bond.
6 Predictive models
The prediction of the spectral properties of chemical structure by graph-entropy measures utilizes structure-property models which play an important role in characterizing and prediction chemical properties using topological indices (Raza et al., 2024; Rauf et al., 2022) These models offer a cost-effective alternative to experimental studies, offering reliability, accuracy and robustness (Hayat et al., 2019). For coronene fractals, we examine the relationship between spectral features and entropy measurements obtained from the reverse degree-based indices. Our findings show that there is a better correlation between the entropy measures and spectral properties, except for the HOMO-LUMO energy gap, which exhibits negative correlation due to its decrease in energy gaps with increasing system size. As noted in the previous section, the first Zagreb index was employed to compare relative complexity measures among the structures. We found that entropy measures associated with demonstrate the strongest correlation with spectral characteristics. Linear regression analysis was used to develop predictive models for spectral properties. The linear regression equation given as , where is the spectral properties, is the regression coefficient, and is the regression constant. The statistical parameters such as , , -values, and are utilized to validate model’s performance.
The regression models optimized to predict spectral characteristics are given detailed in Table 10 and illustrated in Figure 5. The selection was based on their unique performance indicators, such as , adjusted , high -values, in addition to reduced error (SE) objectives. These metrics confirm the reliability and accuracy of the models. The developed models are particularly effective in estimating the energy value of high-aspect ratio coronene explosions. An efficient method based on linear regression was used to ensure accurate predictions while minimizing computational complexity.
TABLE 10
| P | Optimal regression equation | ||||
|---|---|---|---|---|---|
| 0.248–2.775 | 0.999 | 0.998 | 0.886 | 17582849.51 | |
| energy gaps | −0.038+1.036 | 0.951 | 0.946 | 0.015 | 175.129 |
| 0.004+1.448 | 0.973 | 0.970 | 0.001 | 329.292 | |
| 0.004+0.448 | 0.973 | 0.970 | 0.001 | 329.292 | |
| 0.001+0.149 | 0.973 | 0.970 | 0.006 | 327.889 | |
| 0.018+5.467 | 0.915 | 0.906 | 0.009 | 96.845 |
Statistically derived optimal regression models for predicting energetic properties.
FIGURE 5

Linear regression models for the energetic properties. (a) Total π electron energy (b) Homo-Lumo energy gap (c) Energy per bond (d) Delocalization energy (e) Resonance energy per bond (f) Spectral diameter.
7 Conclusion
In this paper, we develop topological expressions based on modified reverse degree-based indices for three configurations of two-dimensional coronene fractals. These indices capture structural complexities and are effective in predicting physico-chemical properties. The computed indices function as graph-based metrics for evaluating entropy levels and relative complexity. The resulting entropy values offer insights into the structural challenges of these fractal systems, providing a foundation for further investigation into their properties. When paired with graph spectra, these approaches form a comprehensive machine learning framework for efficiently and accurately computing the spectral and thermodynamic properties of fractals and other two-dimensional materials. By integrating graph-theoretic methods with advanced statistical techniques, this study contributes to the development of improved computational chemistry algorithms, particularly for QSAR and QSPR studies aimed at predicting the stability and characteristics of complex chemical systems.
Statements
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
AK: Conceptualization, Formal Analysis, Methodology, Writing – original draft. AG: Investigation, Methodology, Supervision, Validation, Writing – review and 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 no Generative AI was used in the creation of this manuscript.
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
Abul Kalaam A. R. Berin Greeni A. (2024). Comparative analysis of modified reverse degree topological indices for certain carbon nanosheets using entropy measures and multi criteria decision-making analysis, Int. J. Quantum Chem.124(1), pe.27326, 10.1002/qua.27326
2
Arockiaraj M. Greeni A. B. Kalaam A. A. (2023a). Linear versus cubic regression models for analyzing generalized reverse degree based topological indices of certain latest corona treatment drug molecules, Int. J. Quantum. Chem.123(10), pe.27136, 10.1002/qua.27136
3
Arockiaraj M. Greeni A. B. Kalaam A. A. (2023b). Comparative analysis of reverse degree and entropy topological indices for drug molecules in blood cancer treatment through QSPR regression models. Polycycl. Aromat. Compd., 1–18. 10.1080/10406638.2023.2271648
4
Arockiaraj M. Greeni A. B. Kalaam A. A. Aziz T. Alharbi M. (2024). Mathematical modeling for prediction of physicochemical characteristics of cardiovascular drugs via modified reverse degree topological indices. EPJE47 (8), 53. 10.1140/epje/s10189-024-00446-3
5
Arockiaraj M. Jency J. Abraham J. Ruth Julie Kavitha S. Balasubramanian K. (2022). Two-dimensional coronene fractal structures: topological entropy measures, energetics, NMR and ESR spectroscopic patterns and existence of isentropic structures. Mol. Phys.120 (11), e2079568. 10.1080/00268976.2022.2079568
6
Arockiaraj M. Jency J. Mushtaq S. Shalini A. J. Balasubramanian K. (2023d). Covalent organic frameworks: topological characterizations, spectral patterns and graph entropies. J. Math. Chem.61, 1633–1664. 10.1007/s10910-023-01477-5
7
Arockiaraj M. Paul D. Clement J. Tigga S. Jacob K. Balasubramanian K. (2023c). Novel molecular hybrid geometric-harmonic-Zagreb degree based descriptors and their efficacy in QSPR studies of polycyclic aromatic hydrocarbons. SAR QSAR Environ. Res.34, 569–589. 10.1080/1062936x.2023.2239149
8
Balasubramanian K. (2023). Topological indices, graph spectra, entropies, Laplacians, and matching polynomials of n-dimensional hypercubes. Symmetry15 (2), 557. 10.3390/sym15020557
9
Bonchev D. Trinajstić N. (1982). Chemical information theory: structural aspects. Int. J. Quantum Chem.22 (16), 463–480. 10.1002/qua.560220845
10
Dehmer M. (2008). Information processing in complex networks: graph entropy and information functionals. Appl Math Comput.201, 82–94. 10.1016/j.amc.2007.12.010
11
Demir C. Üngördü A. (2023). The design of push-pull substituted coronene molecules for optoelectronic applications. Mater. Chem. Phys.301, 127631. 10.1016/j.matchemphys.2023.127631
12
Dobrowolski M. A. Ciesielski A. Cyrański M. K. (2011). On the aromatic stabilization of corannulene and coronene. Phys. Chem. Chem. Phys.13 (46), 20557–20563. 10.1039/c1cp21994d
13
Duan Q. An J. Mao H. Liang D. Li H. Wang S. et al (2021). Review about the application of fractal theory in the research of packaging materials. Materials14 (4), 860. 10.3390/ma14040860
14
Estrada E. Uriarte E. (2001). Recent advances on the role of topological indices in drug discovery research. Curr. Med. Chem.8 (13), 1573–1588. 10.2174/0929867013371923
15
Fedotov P. V. Chernov A. I. Talyzin A. V. Anoshkin I. V. Nasibulin A. G. Kauppinen E. I. et al (2013). Optical study of nanotube and coronene composites. J. Nanoelectron. Optoelectron.8 (1), 16–22. 10.1166/jno.2013.1428
16
Graovac A. Gutman I. Trinajstić N. (1977). Topological approach to the chemistry of conjugated molecules. Berlin: Springer-Verlag.
17
Gutman I. (1978). The energy of a graph. Ber. Math. Stat. Sekt. Forschungszentrum Graz103, 1–22.
18
Gutman I. Trinajstić N. (1972). Graph theory and molecular orbitals. Total π-electron energy of alternant hydrocarbons. Chem. Phys. Lett.17 (4), 535–538. 10.1016/0009-2614(72)85099-1
19
Hayat S. Imran M. Liu J. B. (2019). Correlation between the Estrada index and π-electronic energies for benzenoid hydrocarbons with applications to boron nanotubes. Int. J. Quantum Chem.119 (23), e26016. 10.1002/qua.26016
20
Herndon W. C. Ellzey M. L. Jr. (1974). Resonance theory. V. Resonance energies of benzenoid and nonbenzenoid π systems. J. Am. Chem. Soc.96 (21), 6631–6642. 10.1021/ja00828a015
21
Hill J. P. Jin W. Kosaka A. Fukushima T. Ichihara H. Shimomura T. et al (2004). Self-assembled hexa-peri-hexabenzocoronene graphitic nanotube. Science304 (5676), 1481–1483. 10.1126/science.1097789
22
Hosoya H. (2005). “Clar’s aromatic sextet and sextet polynomial,” in Advances in the theory of benzenoid hydrocarbons (Springer Berlin Heidelberg), 255–272.
23
Junias J. S. Clement J. Rahul M. P. Arockiaraj M. (2024). Two-dimensional phthalocyanine frameworks: topological descriptors, predictive models for physical properties and comparative analysis of entropies with different computational methods. Comput. Mater. Sci.235, 112844. 10.1016/j.commatsci.2024.112844
24
Kalaam A. A. Greeni A. B. Arockiaraj M. (2024). Modified reverse degree descriptors for combined topological and entropy characterizations of 2D metal organic frameworks: applications in graph energy prediction. Front. Chem.12, 1470231. 10.3389/fchem.2024.1470231
25
Khabyah A. Ahmad A. Azeem M. Ahmad Y. Koam A. N. (2023). Reverse-degree-based topological indices of two-dimensional coronene fractal structures, Phys. Scr.99(1), pe.015216, 10.1088/1402-4896/ad10db
26
Kirkby M. J. (1983). The fractal geometry of nature. Benoit B. Mandelbrot. W. H. Freeman and co., San Francisco, 1982. No. of pages: 460. Price: £22.75 (hardback). Earth Surf. Process. Landforms8, 406. 10.1002/esp.3290080415
27
Kumar A. Duran M. Solá M. (2017). Is coronene better described by Clar’s aromatic π-sextet model or by the AdNDP representation?. J. Comput. Chem.38 (18), 1606–1611. 10.1002/jcc.24801
28
Kumar V. Das S. (2024). On structure sensitivity and chemical applicability of some novel degree-based topological indices. Comput. Chem.92 (1), 165–203. 10.46793/match.92-1.165k
29
Lee J. S. Chang K. S. (1996). Applications of chaos and fractals in process systems engineering. J. Process Control6 (2-3), 71–87. 10.1016/0959-1524(95)00051-8
30
Li X. Li Y. Shi Y. Gutman I. (2013). Note on the HOMO-LUMO index of graphs. MATCH Commun. Math. Comput. Chem.70 (1), 85–96.
31
Malik M. A. Imran M. Adeel M. (2023). On distance-based topological indices and co-indices of fractal-type molecular graphs and their respective graph entropies. PLoS One18 (11), e0290047. 10.1371/journal.pone.0290047
32
MATLAB (2022). MATLAB version: 9.13.0 (R2022b), Natick, Massachusetts: The MathWorks Inc.
33
Mazouin B. Schöpfer A. A. von Lilienfeld O. A. (2022). Selected machine learning of HOMO-LUMO gaps with improved data-efficiency. Mater. Adv.3 (22), 8306–8316. 10.1039/d2ma00742h
34
Mowshowitz A. (1968). Entropy and the complexity of graphs: I. An index of the relative complexity of a graph. Bull. Math. Biophys.30, 175–204. 10.1007/bf02476948
35
Mowshowitz A. Dehmer M. (2012). Entropy and the complexity of graphs revisited. Entropy14 (3), 559–570. 10.3390/e14030559
36
Newman M. S. (1940). A new synthesis of coronene. J. Am. Chem. Soc.62 (7), 1683–1687. 10.1021/ja01864a014
37
Nisha S. Senthil Kumar A. (2020). π-Self-assembly of a coronene on carbon nanomaterial-modified electrode and its symmetrical redox and H2O2 electrocatalytic reduction functionalities. ACS Omega5 (20), 11817–11828. 10.1021/acsomega.0c01258
38
Pisula W. Feng X. Müllen K. (2010). Tuning the columnar organization of discotic polycyclic aromatic hydrocarbons. Adv. Mater.22 (33), 3634–3649. 10.1002/adma.201000585
39
Pisula W. Feng X. Müllen K. (2011). Charge-carrier transporting graphene-type molecules. Chem. Mater.23 (3), 554–567. 10.1021/cm102252w
40
Popov I. A. Boldyrev A. I. (2012). Chemical bonding in coronene, isocoronene, and circumcoronene. Eur. J. Org. Chem.2012 (18), 3485–3491. 10.1002/ejoc.201200256
41
Prabhu S. Arulperumjothi M. Manimozhi V. Balasubramanian K. (2024). Topological characterizations on hexagonal and rectangular tessellations of antikekulenes and its computed spectral, nuclear magnetic resonance and electron spin resonance characterizations. Int. J. Quantum Chem.124 (7), 27365. 10.1002/qua.27365
42
Rauf R. A. Naeem M. Bukhari S. U. (2022). Quantitative structure–property relationship of Ev‐degree and Ve‐degree based topological indices: physico‐chemical properties of benzene derivatives. Int. J. Quantum Chem.122 (5), e26851. 10.1002/qua.26851
43
Raza Z. Arockiaraj M. Maaran A. Shalini A. J. (2024). A comparative study of topological entropy characterization and graph energy prediction for Marta variants of covalent organic frameworks. Front. Chem.12, 1511678. 10.3389/fchem.2024.1511678
44
Robertson J. M. White J. G. (1945). 164. The crystal structure of coronene: a quantitative X-ray investigation. J. Chem. Soc. (Resumed), 607–617. 10.1039/jr9450000607
45
Sabirov D. S. O-sawa E. (2015). Information entropy of fullerenes. J. Chem. Inf. Model.55 (8), 1576–1584. 10.1021/acs.jcim.5b00334
46
Sabirov D. S. Shepelevich I. S. (2021). Information entropy in chemistry: an overview. Entropy23 (10), 1240. 10.3390/e23101240
47
Santa Daría A. M. González-Sánchez L. Gómez S. (2024). Coronene: a model for ultrafast dynamics in graphene nanoflakes and PAHs. PCCP26 (1), 174–184. 10.1039/d3cp03656a
48
Sanyal S. Manna A. K. Pati S. K. (2013). Effect of imide functionalization on the electronic, optical, and charge transport properties of coronene: a theoretical study. J. Phys. Chem. C117 (2), 825–836. 10.1021/jp310362c
49
Shannon C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J.27 (3), 3–55. 10.1145/584091.584093
50
Stevanović L. Brankov V. Cvetković D. Simić S. (2021). newGRAPH: a fully integrated environment used for research process in graph theory. Available online at: http://www.mi.sanu.ac.rs/newgraph/index.html.
51
Tachikawa H. Lund A. (2022). Structures and electronic states of trimer radical cations of coronene: DFT-ESR simulation study. PCCP24 (17), 10318–10324. 10.1039/d1cp04638a
52
Uahabi K. L. Atounti M. (2015). Applications of fractals in medicine. Ann. Univ. Craiova-Math. Comput. Sci. Ser.42 (1), 167–174.
53
Ullah A. Nazir M. Zaman S. Hamed Y. S. Jabeen S. (2024). Fractal configurations of zigzag hexagonal type coronoid molecules: graph-theoretical modeling and its impact on physicochemical behavior, Phys. Scr.100(1), pe.015237, 10.1088/1402-4896/ad9a1f
54
Wiener H. (1947). Structural determination of paraffin boiling points. J. Am. Chem. Soc.69, 17–20. 10.1021/ja01193a005
55
Wu Y. Yang Y. Ye D. (2018). A note on median eigenvalues of bipartite graphs. MATCH Commun. Math. Comput. Chem.80 (3), 853–862.
56
Xu S. A. Liu J. B. (2025). Entropies and degree-based topological indices of coronene fractal structures. Fractal Fract.9 (3), 133. 10.3390/fractalfract9030133
57
Yogalakshmi K. Easwaramoorthy D. (2024). A new approach for generalized fractal dimensions based on topological indices for pyracyclene and pentahexoctite networks. IEEE Access12, 166176–166187. 10.1109/access.2024.3489218
Summary
Keywords
coronene fractals, modified reverse degree-based indices, information entropy, spectral properties, predictive models
Citation
Kalaam ARA and Greeni AB (2025) Two-dimensional coronene fractals: modified reverse degree indices, comparative analysis of information entropy and predictive modeling of spectral properties. Front. Chem. 13:1588942. doi: 10.3389/fchem.2025.1588942
Received
06 March 2025
Accepted
18 April 2025
Published
06 May 2025
Volume
13 - 2025
Edited by
Craig A. Bayse, Old Dominion University, United States
Reviewed by
Joseph Clement, VIT University, India
Muhammad Kashif Masood, Southeast University, China
Updates
Copyright
© 2025 Kalaam and Greeni.
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: A. Berin Greeni, beringreeni@gmail.com
Disclaimer
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.