ORIGINAL RESEARCH article

Front. Chem., 06 May 2025

Sec. Theoretical and Computational Chemistry

Volume 13 - 2025 | https://doi.org/10.3389/fchem.2025.1588942

Two-dimensional coronene fractals: modified reverse degree indices, comparative analysis of information entropy and predictive modeling of spectral properties

  • School of Advanced Sciences, Vellore Institute of Technology, Chennai, India

Article metrics

View details

1

Citations

1,2k

Views

299

Downloads

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

FIGURE 1

The degree based enumeration of coronene fractal structure.

FIGURE 2

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 112, 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 24.

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 57 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 57 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

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 58, 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

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

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., 118. 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, 16331664. 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, 569589. 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), 463480. 10.1002/qua.560220845

  • 10

    Dehmer M. (2008). Information processing in complex networks: graph entropy and information functionals. Appl Math Comput.201, 8294. 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), 2055720563. 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), 15731588. 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), 1622. 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, 122.

  • 18

    Gutman I. Trinajstić N. (1972). Graph theory and molecular orbitals. Total π-electron energy of alternant hydrocarbons. Chem. Phys. Lett.17 (4), 535538. 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), 66316642. 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), 14811483. 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), 255272.

  • 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), 16061611. 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), 165203. 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), 7187. 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), 8596.

  • 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), 83068316. 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, 175204. 10.1007/bf02476948

  • 35

    Mowshowitz A. Dehmer M. (2012). Entropy and the complexity of graphs revisited. Entropy14 (3), 559570. 10.3390/e14030559

  • 36

    Newman M. S. (1940). A new synthesis of coronene. J. Am. Chem. Soc.62 (7), 16831687. 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), 1181711828. 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), 36343649. 10.1002/adma.201000585

  • 39

    Pisula W. Feng X. Müllen K. (2011). Charge-carrier transporting graphene-type molecules. Chem. Mater.23 (3), 554567. 10.1021/cm102252w

  • 40

    Popov I. A. Boldyrev A. I. (2012). Chemical bonding in coronene, isocoronene, and circumcoronene. Eur. J. Org. Chem.2012 (18), 34853491. 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), 607617. 10.1039/jr9450000607

  • 45

    Sabirov D. S. O-sawa E. (2015). Information entropy of fullerenes. J. Chem. Inf. Model.55 (8), 15761584. 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), 174184. 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), 825836. 10.1021/jp310362c

  • 49

    Shannon C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J.27 (3), 355. 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), 1031810324. 10.1039/d1cp04638a

  • 52

    Uahabi K. L. Atounti M. (2015). Applications of fractals in medicine. Ann. Univ. Craiova-Math. Comput. Sci. Ser.42 (1), 167174.

  • 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, 1720. 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), 853862.

  • 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, 166176166187. 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

*Correspondence: A. Berin Greeni,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics