# VDB Entropy Measures and Irregularity-Based Indices for the Rectangular Kekulene System

Authors: Weidong Zhao; K. Julietraja; P. Venugopal; Xiujun Zhang
Journal: Journal of Mathematics
 (2021)

---

## Abstract

Theoretical chemists are fascinated by polycyclic aromatic hydrocarbons (PAHs) because of their unique electromagnetic and other significant properties, such as superaromaticity. The study of PAHs has been steadily increasing because of their wide-ranging applications in several fields, like steel manufacturing, shale oil extraction, coal gasification, production of coke, tar distillation, and nanosciences. Topological indices (TIs) are numerical quantities that give a mathematical expression for the chemical structures. They are useful and cost-effective tools for predicting the properties of chemical compounds theoretically. Entropic network measures are a type of TIs with a broad array of applications, involving quantitative characterization of molecular structures and the investigation of some specific chemical properties of molecular graphs. Irregularity indices are numerical parameters that quantify the irregularity of a molecular graph and are used to predict some of the chemical properties, including boiling points, resistance, enthalpy of vaporization, entropy, melting points, and toxicity. This study aims to determine analytical expressions for the VDB entropy and irregularity-based indices in the rectangular Kekulene system.


---

## Body

## 1. Introduction

Chemical graph theory is one of the important branches of mathematics, which helps researchers to explore chemical compounds by reconstructing their underlying molecular structures as graphs. Amolecular graph is a simple graph that represents the skeleton of an organic chemical compound in which the vertices of the molecular graph represent carbon atoms, and its edges represent the links between those carbon atoms [1]. According to IUPAC, the molecular structural descriptor is described as “A molecular structural descriptor is a numerical quantity associated with the chemical constitution that can be used to correlate chemical structure with several physical characteristics, biological activity, or chemical reactivity” [2].Cycloarenes are conjugated macrocyclic chemical compounds composed of circularly fused benzene rings that enclose a cavity with inward-facing C-H bonds. Thus, cycloarenes are a subdivision of the Kekulene system. The magnetic susceptibility, aromaticity, and vibrational frequencies of cycloarenes have all been the subject of theoretical research [3, 4].The rectangular Kekulene systemRKm,n is a type of polycyclic aromatic compound that has a broad range of applications, from everyday materials to cutting-edge industries like nanotechnology. As a result, theoretical chemists have long been attracted to the study of their fundamental features, including carcinogenicity, toxicity, observed bioactivities, and other relevant features. The fascinating potential of increased stability with interesting magnetic and magnetocaloric effects is offered by the structure of the Kekulene molecule, which has 12 annulated benzene rings and a central cavity.Numerous theoretical and experimental studies on coronoids have been sparked by the synthesis of Kekulene and septulene [5–7]. Due to their intriguing property of superaromaticity, the study of Kekulene systems is also gaining momentum. Extra thermodynamic stability owing to macrocyclic conjugation in super-ring compounds like Kekulene is referred to as superaromaticity or macrocyclic aromaticity, and it accounts for a minor portion of global aromaticity. All of the coronoids that have been investigated so far are superaromatic in some way, with positive superaromatic stabilisation energies (SSEs). This gives rise to a more in-depth investigation of the features of Kekulene systems and their linkages to the underlying molecular structures. This research can be used in a variety of nanotechnology fields.The cavities, which are an intrinsic element of the Kekulene system, serve as prototypes for developing and synthesizing novel nanomaterials with applications in nanotechnology and biotechnology, as well as the burgeoning field of nanomedicine. They have also been used in the design and synthesis of porous and mesoporous materials based on calixarenes and mesoporous silica for the storage and complex formation of hazardous nuclear waste and other contaminants [8–10]. It has been demonstrated that by utilizing Kekulene system knowledge, one can accurately compute total pi-electron energy, resonance energy, and coronoid hydrocarbon enumeration [11].Graph-theoretical techniques are more efficient in acquiring the properties of the Kekulene system than computationally intensive quantum chemical calculations. Aihara [12] pointed out in a recent study that graph theory is beneficial not only for predicting topological resonance energies but also for exposing severe flaws in prior aromaticity theories. Recently, Julietraja and Venugopal computed the degree-based indices for coronoid structures [13]. Julietraja et al. analyzed the degree-based molecular descriptors using M-polynomial for certain classes of benzenoid systems [14, 15]. Julietraja et al. used entropy measures to investigate three prominent classes of PAHs [16].However, despite the several promising applications of the rectangular kekulene systems, there have been no attempts till now to study these systems from a structural perspective. Our research, detailed in this article, is an effort towards bridging this research gap. The objective of this article is to compute the VDB entropy measures and irregularity-based indices for the rectangular Kekulene system. Since we tackle two open problems in this article, the number of known TIs for these systems will be doubled. This increased number of known molecular descriptors will contribute to the current body of knowledge on these systems. The TIs will also form the basis on which new Kekulene structures can be designed and synthesised.The subsequent sections cover a literature survey on graph-theoretical concepts, degree-based entropy measures, and irregularity-based indices. There is a brief summary of the methods used for research, followed by the computation of the analytical expressions of both entropy measures and irregularity indices for rectangular Kekulene systems. The analytical expressions of the descriptors are used to compute the numerical values and observe their behaviour using graphical plots.


## 2. Graph-Theoretical Concepts

In this article,Γ1V,E denotes a connected graph, where V and E denote the vertex and edge sets, respectively. The degree of a vertex m1 in a graph Γ1 is the number of edges that are adjacent to that vertex m1 and is denoted by degΓ1m1 [17]. There is a large body of knowledge about degree-based topological descriptors, which have a strong co-relation with numerous physicochemical properties of PAHs such as Zagreb and its co-indices, ABC index, GA index, Randić index, sum-connectivity index, and SDD index. These indices are also used to generate degree-based entropy measures.


### 2.1. Edge Weight- and Degree-Based Entropy

Letn be the order of a graph of size m and φ be some meaningful information function. The Shannon’s entropy [18, 19] of a graph Γ1 is defined as(1)ENTφΓ1=−∑l=1nφdegΓ1vl∑m=1nφdegΓ1vmlogφdegΓ1vl∑m=1nφdegΓ1vm.Letvl∈VΓ1 and the degree of vl be represented by the information function φvl, that is, φvl=degΓ1vl. Then, equation (1) can be rewritten as(2)ENTφΓ1=−∑l=1ndegΓ1vl∑m=1ndegΓ1vmlogdegΓ1vl∑m=1ndegΓ1vm.According to the fundamental theorem of graph theory,∑m=1ndegΓ1vm=2E. Hence, equation (2) reduces to(3)ENTφΓ1=log2E−12Elog∑l=1ndegΓ1vldegΓvl.The entropy measure of an edge-weighted graph was first proposed by Chen et al. [20]. If Γ1=VΓ1;EΓ1;φl1m1 is an edge-weighted graph, where VΓ1,EΓ1, and φl1m1 represents vertex set, edge set, and edge weight of the edge l1m1 of Γ1, then we have(4)ENTφΓ1=−∑l1′m1′∈EΓφl1′m1′∑l1m1∈EΓ1φl1m1logφl1′m1′∑l1m1∈EΓ1φl1m1.By using the equations of (1), (3), and (4), we get the VDB entropy measures listed in Table 1.Table 1
Edge weight- and degree-based entropy.
S. no.Degree-based entropyMathematical expressions1First Zagreb entropyENTM1Γ1=logM1Γ1−1/M1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1degΓ1l1+degΓ1m12Second Zagreb entropyENTM2Γ1=logM2Γ1−1/M2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m13Second modified Zagreb entropyENTM2mΓ1=logM2mΓ1−1/M2mΓ1∏l1m1∈EΓ11/degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m14Reduced second Zagreb entropyENTRM2Γ1=logRM2Γ1−1/RM2Γ1log∏l1m1∈EΓ1degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−15Hyper Zagreb entropyENTHMΓ1=logHMΓ1−1/HMΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m12degΓ1l1+degΓ1m126Augmented Zagreb entropyENTAΓ1=logAΓ1−1/AΓ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−237Atom bond connectivity entropyENTABCΓ1=logABCΓ1−1/ABCΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m18Geometric arithmetic entropyENTGAΓ1=logGAΓ1−1/GAΓ1log∏l1m1∈EΓ12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m19Symmetric division deg (SDD) entropyENTSDDΓ1=logSDDΓ1−1/SDDΓlog∏l1m1∈EΓ1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l110Randić entropyENTRΓ1=logRΓ1−1/RΓ1log∏l1m1∈EΓ11/degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m111Forgotten entropyENTFΓ1=logFΓ1−1/FΓ1log∏l1m1∈EΓ1degΓ1l12+degΓ1m12degΓ1l12+degΓ1m1212Sum-connectivity entropyENTχΓ1=logχΓ1−1/χΓ1log∏l1m1∈EΓ11/degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m113First redefined Zagreb entropyENTReZG1Γ1=logReZG1Γ1−1/ReZG1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m114Second redefined Zagreb entropyENTReZG2Γ1=logReZG2Γ1−1/ReZG2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1


### 2.2. Irregularity-Based Indices for QSPR Analysis

Irregularity indices are numerical parameters that quantify the irregularity of a molecular graph. The study of irregular graphs was brought to the limelight by Paul Erdös in [21]. Erdös raised an open question on irregular graphs [22], at the Second Krakow Conference on Graph Theory.If a topological descriptor becomes zero for a normal graph and nonzero for a nonregular graph, it is known as an irregularity index. Bell proposed the first irregularity index in [23]. Irregularity indices have proved to be particularly useful in quantifying the topology of nonregular graphs. In QSPR and QSAR studies, irregularity of graphs is used to predict a range of physical and chemical properties, including melting and boiling points, enthalpy of vaporization, resistance, entropy, and toxicity. It has been proved that using regression analysis in [24], 4-octane isomer properties such as enthalpy of vaporization (HVAP), entropy, acentric factor (AcenFac), and standard enthalpy of vaporization (DHVAP) can be determined using irregularity indices with a correlation coefficient of magnitude higher than 0.9. Some commonly used irregularity indices and their analytical expressions are listed in Table 2.Table 2
Irregularity-based indices for QSPR analysis.
S. no.Irregularity-based indices1VARΓ1=∑l1∈VΓ1degΓ1l1−2m/n2=M1Γ/n−2m/n22IR1Γ1=∑l1∈VΓ1degΓ1l13−2m/n∑l1∈VΓ1degΓ1l12=FΓ−2m/nM1Γ3IR2Γ1=∑l1m1∈EΓ1degΓ1l1⋅degΓ1m1/m−2m/n=M2Γ/m−2m/n4IRDIFΓ1=∑l1m1∈EΓ1degΓ1l1/degΓ1m1−degΓ1m1/degΓ1m15ALΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m16IRLΓ1=∑l1m1∈EΓ1lndegΓ1l1−lndegΓ1m17IRLUΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m1/mindegΓ1l1⋅degΓ1m18IRLFΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m1/degΓ1l1⋅degΓ1m19IRFΓ1=∑l1m1∈EΓ1du−dv2=FΓ1−2M2Γ110IRLAΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m1/degΓ1l1+degΓ1m111IRD1Γ1=∑l1m1∈EΓ1ln1+degΓ1l1−degΓ1m112IRAΓ1=∑l1m1∈EΓ1degΓ1l1−1/2−degΓ1m1−1/2213IRBΓ1=∑l1m1∈EΓ1degΓ1l1−1/2−degΓ1m1−1/2214IRBΓ1=∑l1m1∈EΓ1degΓ1l1⋅degΓ1m1/m−2m/n=RRΓ1/m−2m/n15IRGAΓ1=∑l1m1∈EΓ1lndegΓ1l1+degΓ1m1/2degΓ1l1⋅degΓ1m116IRRtΓ1=1/2∑l1m1∈EΓ1degΓ1l1−degΓ1m1The recent work and development of irregularity indices can be seen more elaborately in [25–29].


## 2.1. Edge Weight- and Degree-Based Entropy

Letn be the order of a graph of size m and φ be some meaningful information function. The Shannon’s entropy [18, 19] of a graph Γ1 is defined as(1)ENTφΓ1=−∑l=1nφdegΓ1vl∑m=1nφdegΓ1vmlogφdegΓ1vl∑m=1nφdegΓ1vm.Letvl∈VΓ1 and the degree of vl be represented by the information function φvl, that is, φvl=degΓ1vl. Then, equation (1) can be rewritten as(2)ENTφΓ1=−∑l=1ndegΓ1vl∑m=1ndegΓ1vmlogdegΓ1vl∑m=1ndegΓ1vm.According to the fundamental theorem of graph theory,∑m=1ndegΓ1vm=2E. Hence, equation (2) reduces to(3)ENTφΓ1=log2E−12Elog∑l=1ndegΓ1vldegΓvl.The entropy measure of an edge-weighted graph was first proposed by Chen et al. [20]. If Γ1=VΓ1;EΓ1;φl1m1 is an edge-weighted graph, where VΓ1,EΓ1, and φl1m1 represents vertex set, edge set, and edge weight of the edge l1m1 of Γ1, then we have(4)ENTφΓ1=−∑l1′m1′∈EΓφl1′m1′∑l1m1∈EΓ1φl1m1logφl1′m1′∑l1m1∈EΓ1φl1m1.By using the equations of (1), (3), and (4), we get the VDB entropy measures listed in Table 1.Table 1
Edge weight- and degree-based entropy.
S. no.Degree-based entropyMathematical expressions1First Zagreb entropyENTM1Γ1=logM1Γ1−1/M1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1degΓ1l1+degΓ1m12Second Zagreb entropyENTM2Γ1=logM2Γ1−1/M2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m13Second modified Zagreb entropyENTM2mΓ1=logM2mΓ1−1/M2mΓ1∏l1m1∈EΓ11/degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m14Reduced second Zagreb entropyENTRM2Γ1=logRM2Γ1−1/RM2Γ1log∏l1m1∈EΓ1degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−15Hyper Zagreb entropyENTHMΓ1=logHMΓ1−1/HMΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m12degΓ1l1+degΓ1m126Augmented Zagreb entropyENTAΓ1=logAΓ1−1/AΓ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−237Atom bond connectivity entropyENTABCΓ1=logABCΓ1−1/ABCΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m18Geometric arithmetic entropyENTGAΓ1=logGAΓ1−1/GAΓ1log∏l1m1∈EΓ12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m19Symmetric division deg (SDD) entropyENTSDDΓ1=logSDDΓ1−1/SDDΓlog∏l1m1∈EΓ1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l110Randić entropyENTRΓ1=logRΓ1−1/RΓ1log∏l1m1∈EΓ11/degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m111Forgotten entropyENTFΓ1=logFΓ1−1/FΓ1log∏l1m1∈EΓ1degΓ1l12+degΓ1m12degΓ1l12+degΓ1m1212Sum-connectivity entropyENTχΓ1=logχΓ1−1/χΓ1log∏l1m1∈EΓ11/degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m113First redefined Zagreb entropyENTReZG1Γ1=logReZG1Γ1−1/ReZG1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m114Second redefined Zagreb entropyENTReZG2Γ1=logReZG2Γ1−1/ReZG2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1


## 2.2. Irregularity-Based Indices for QSPR Analysis

Irregularity indices are numerical parameters that quantify the irregularity of a molecular graph. The study of irregular graphs was brought to the limelight by Paul Erdös in [21]. Erdös raised an open question on irregular graphs [22], at the Second Krakow Conference on Graph Theory.If a topological descriptor becomes zero for a normal graph and nonzero for a nonregular graph, it is known as an irregularity index. Bell proposed the first irregularity index in [23]. Irregularity indices have proved to be particularly useful in quantifying the topology of nonregular graphs. In QSPR and QSAR studies, irregularity of graphs is used to predict a range of physical and chemical properties, including melting and boiling points, enthalpy of vaporization, resistance, entropy, and toxicity. It has been proved that using regression analysis in [24], 4-octane isomer properties such as enthalpy of vaporization (HVAP), entropy, acentric factor (AcenFac), and standard enthalpy of vaporization (DHVAP) can be determined using irregularity indices with a correlation coefficient of magnitude higher than 0.9. Some commonly used irregularity indices and their analytical expressions are listed in Table 2.Table 2
Irregularity-based indices for QSPR analysis.
S. no.Irregularity-based indices1VARΓ1=∑l1∈VΓ1degΓ1l1−2m/n2=M1Γ/n−2m/n22IR1Γ1=∑l1∈VΓ1degΓ1l13−2m/n∑l1∈VΓ1degΓ1l12=FΓ−2m/nM1Γ3IR2Γ1=∑l1m1∈EΓ1degΓ1l1⋅degΓ1m1/m−2m/n=M2Γ/m−2m/n4IRDIFΓ1=∑l1m1∈EΓ1degΓ1l1/degΓ1m1−degΓ1m1/degΓ1m15ALΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m16IRLΓ1=∑l1m1∈EΓ1lndegΓ1l1−lndegΓ1m17IRLUΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m1/mindegΓ1l1⋅degΓ1m18IRLFΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m1/degΓ1l1⋅degΓ1m19IRFΓ1=∑l1m1∈EΓ1du−dv2=FΓ1−2M2Γ110IRLAΓ1=∑l1m1∈EΓ1degΓ1l1−degΓ1m1/degΓ1l1+degΓ1m111IRD1Γ1=∑l1m1∈EΓ1ln1+degΓ1l1−degΓ1m112IRAΓ1=∑l1m1∈EΓ1degΓ1l1−1/2−degΓ1m1−1/2213IRBΓ1=∑l1m1∈EΓ1degΓ1l1−1/2−degΓ1m1−1/2214IRBΓ1=∑l1m1∈EΓ1degΓ1l1⋅degΓ1m1/m−2m/n=RRΓ1/m−2m/n15IRGAΓ1=∑l1m1∈EΓ1lndegΓ1l1+degΓ1m1/2degΓ1l1⋅degΓ1m116IRRtΓ1=1/2∑l1m1∈EΓ1degΓ1l1−degΓ1m1The recent work and development of irregularity indices can be seen more elaborately in [25–29].


## 3. Methods

In this paper, the computations are performed using graph-theoretical methods, the edge partition method, and analytical techniques. Maple is utilized in obtaining the analytical expressions for the degree-based entropy measures. Maple is used to calculate the numerical values of the VDB entropy measures and irregularity indices based on analytical expressions. The numerical results are represented visually using Origin, and Chem Draw Ultra is used for describing the molecular structures of the three PAHs.


### 3.1. Degree-Based Entropy for the Rectangular Kekulene System

LetΓ be a rectangular Kekulene system RKr,s. The number of vertices and edges of RKr,s is VΓ1=36rs+14s−2r and EΓ1=48rs+16s−4r. RKr,s consists of 12rs+10s+2r vertices of degree 2 and 24rs+4s−4r vertices of degree 3. It is picturised in Figure 1. The edge partition table of rectangular Kekulene system RKr,s is represented in Table 3.Figure 1
Rectangular Kekulene system.Table 3
Edge partition table ofRKr,s.
degΓ1l1,degΓ1m1Total number of edges2,24s+2r2,324rs+12s3,324rs−6r


### 3.2. First Zagreb Entropy

By using the edge partition in Table3,(5)M1Γ1=264s−28r+76s.Then, using the equation of first Zagreb entropy, the result is obtained as(6)ENTM1Γ1=logM1Γ1−1M1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1=logM1Γ−1M1Γlog∏l1m1∈E2,2degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1×∏l1m1∈E2,3degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1×∏l1m1∈E3,3degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1=log264s−28r+76s−1264s−28r+76slog4s+2r×44×24rs+12s×55×24rs−6r×66.


### 3.3. Second Zagreb Entropy

By using the edge partition in Table3,(7)M2Γ1=360s−46r+88s.Then, using the equation of second Zagreb entropy, the result is obtained as(8)ENTM2Γ1=logM2Γ1−1M2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1=logM2Γ1−1M2Γ1log∏l1m1∈E2,2degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1×∏l1m1∈E2,3degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1×∏l1m1∈E3,3degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1=log360s−46r+88s−1360s−46r+88slog4s+2r×44×24rs+12s×66×24rs−6r×99.


### 3.4. Second Modified Zagreb Entropy

By using the edge partition in Table3,(9)M2mΓ=1640s−1r+3s.Then, using the equation of the second modified Zagreb entropy, the result is obtained as(10)ENTM2mΓ1=logM2mΓ1−1M2mΓ1log∏l1m1∈EΓ11degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=logM2mΓ1−1M2mΓ1log∏l1m1∈E2,21degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=log1640s−1r+3s−11/640s−1r+3slog4s+2r×141/4×24rs+12s×161/6×24rs−6r×191/9.


### 3.5. Reduced Second Zagreb Entropy

By using the edge partition in Table3,(11)RM2Γ=144s−22r+28s.Then, using the equation of reduced second Zagreb entropy, the result is obtained as(12)ENTRM2Γ1=logRM2Γ1−1RM2Γ1log∏l1m1∈EΓ1degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1=logM2mΓ1−1M2mΓ1log∏l1m1∈E2,2degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1×∏l1m1∈E2,3degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1×∏l1m1∈E3,3degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1=log144s−22r+28−1144s−22r+28log4s+2r×11×24rs+12s×22×24rs−6r×44.


### 3.6. Hyper Zagreb Entropy

By using the edge partition in Table3,(13)HMΓ1=1464s−184r+364s.Then, using the equation of reduced second Zagreb entropy, the result is obtained as(14)ENTHMΓ1=logHMΓ1−1HMΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12=logHMΓ1−1HMΓ1log∏l1m1∈E2,2degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12×∏l1m1∈E2,3degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12×∏l1m1∈E3,3degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12=log1464s−184r+364s−11464s−184r+364slog4s+2r×1616×24rs+12s×2525×24rs−6r×3636.


### 3.7. Augmented Zagreb Entropy

By using the edge partition in Table3,(15)AΓ1=13214892s−1675r+128s.Then, using the equation of reduced second Zagreb entropy, the result is obtained as(16)ENTAΓ1=logAΓ1−1AΓ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23=logAΓ1−1AΓ1log∏l1m1∈E2,2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23×∏l1m1∈E2,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23×∏l1m1∈E3,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23=log13214892s−1675r+128s−11/3214892s−1675r+128slog4s+2r×88×24rs+12s×88×24rs−6r×72964729/64.


### 3.8. Atom Bond Connectivity Entropy

By using the edge partition in Table3,(17)ABCΓ1=12s+1r+8s2+16s−4r.Then, using the equation of atom bond connectivity entropy, the result is obtained as(18)ENTABCΓ1=logABCΓ1−1ABCΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1=logABCΓ1−1ABCΓ1log∏l1m1∈E2,2degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,3degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,3degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1=log12s+1r+8s2+16s−4r−112s+1r+8s2+16s−4rlog4s+2r×121/2×24rs+12s×121/2×24rs−6r×232/3.


### 3.9. Geometric Arithmetic Entropy

By using the edge partition in Table3,(19)GAΓ1=4s−4r+4856rs+2456s+24rs.Then, using the equation of geometric arithmetic entropy, the result is obtained as(20)ENTGAΓ1=logGAΓ1−1GAΓ1log∏l1m1∈EΓ12degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=logGAΓ1−1GAΓ1log∏l1m1∈E2,22degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E2,32degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E3,32degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=log4s−4r+4856rs+2456s+24rs−14s−4r+48/56rs+24/56s+24rslog4s+2r×11×24rs+12s×26526/5×24rs−6r×11.


### 3.10. Symmetric Division Deg (SDD) Entropy

By using the edge partition in Table3,(21)SDDΓ1=100rs−8r+34s.Then, using the equation of Symmetric division deg entropy, the result is obtained as(22)ENTSDDΓ1=logSDDΓ1−1SDDΓ1log∏l1m1∈EΓ1degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1=logSDDΓ1−1SDDΓ1log∏l1m1∈E2,2degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1×∏l1m1∈E2,3degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1×∏l1m1∈E3,3degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1=log100rs−8r+34s−1100rs−8r+34slog4s+2r×22×24rs+12s×13613/6×24rs−6r×22.


### 3.11. Randić Entropy

By using the edge partition in Table3,(23)RΓ1=46rs+26s+8rs−r+2s.Then, using the equation of Randić entropy, the result is obtained as(24)ENTRΓ1=logRΓ1−1RΓ1log∏l1m1∈EΓ11degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=logRΓ1−1RΓ1log∏l1m1∈E2,21degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=log46rs+26s+8rs−r+2s−146rs+26s+8rs−r+2slog4s+2r×121/2×24rs+12s×161/6×24rs−6r×131/3.


### 3.12. Forgotten Entropy

By using the edge partition in Table3,(25)FΓ1=744rs−92r+188s.Then, using the equation of forgotten entropy, the result is obtained as(26)ENTFΓ1=logFΓ1−1FΓ1log∏l1m1∈EΓ1degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12=logFΓ1−1FΓ1log∏l1m1∈E2,2degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12×∏l1m1∈E2,3degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12×∏l1m1∈E3,3degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12=log744rs−92r+188s−1744rs−92r+188slog4s+2r×88×24rs+12s×1313×24rs−6r×1818.


### 3.13. Sum-Connectivity Entropy

By using the edge partition in Table3,(27)χΓ1=2s+r+2455rs+1255s+46rs−6r.Then, using the equation of sum-connectivity entropy, the result is obtained as(28)ENTχΓ1=logχΓ1−1χΓ1log∏l1m1∈EΓ11degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1=logFΓ1−1FΓ1log∏l1m1∈E2,21degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1×∏l1m1∈E2,31degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1×∏l1m1∈E3,31degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1=log2s+r+2455rs+1255s+46rs−6r−12s+r+24/55rs+12/55s+46rs−6rlog4s+2r×121/2×24rs+12s×151/5×24rs−6r×161/6.


### 3.14. First Redefined Zagreb Entropy

By using the edge partition in Table3,(29)ReZG1Γ1=36rs−2r+14s.Then, using the equation of first redefined Zagreb entropy, the result is obtained as(30)ENTReZG1Γ1=logReZG1Γ1−1ReZG1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1=logReZG1Γ1−1ReZG1Γ1log∏l1m1∈E2,2degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,3degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,3degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1=log36rs−2r+14s−136rs−2r+14slog4s+2r×11×24rs+12s×565/6×24rs−6r×232/3.


### 3.15. Second Redefined Zagreb Entropy

By using the edge partition in Table3,(31)ReZG1Γ1=925s−7r+3245rs.Then, using the equation of the second redefined Zagreb entropy, the result is obtained as(32)ENTReZG2Γ1=logReZG2Γ1−1ReZG2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=logReZG2Γ1−1ReZG2Γ1log∏l1m1∈E2,2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E2,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E3,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=log925s−7r+3245rs−192/5s−7r+324/5rslog4s+2r×11×24rs+12s×656/5×24rs−6r×323/2.Theorem 1.
IfΓ1 is the RKr,s, then the irregularity-based descriptors are expressed as follows:(1)
VARΓ1=236r2s2+36rs2−r2−4rs+5s2/18rs−r+7s2(2)
IR1Γ1=4180r2s2+1434rs2−5r2−115rs+462s2/18rs−r+7s(3)
IR2Γ1=1/2182180rs−23r+44s/12rs−r+4srs−2180rs−23r+44s/12rs−r+4sr+72180rs−23r+44s/12rs−r+4ss−96rs+8r−32s/18rs−r+7s(4)
IRDIFΓ1=1/324rs+12s(5)
ALΓ1=24rs+12s(6)
IRLΓ1=0.4054524rs+12s(7)
IRLUΓ1=0.40545/224rs+12s(8)
IRLFΓ1=1/624rs+12s(9)
IRFΓ1=24rs+12s(10)
IRLAΓ1=2/524rs+12s(11)
IRD1Γ1=0.6931524rs+12s(12)
IRAΓ1=0.01683224rs+12s(13)
IRBΓ1=0.1010624rs+12s(14)
IRBΓ1=1/22166r2s2−126r2s+1926rs2−504r2s2−66rs+426s2+30r2s−444rs2−r2+11rs−100s2/12rs−r+4s18rs−r+7s(15)
IRGAΓ1=0.02039124rs+12s(16)
IRRtΓ1=1/224rs+12s


## 3.1. Degree-Based Entropy for the Rectangular Kekulene System

LetΓ be a rectangular Kekulene system RKr,s. The number of vertices and edges of RKr,s is VΓ1=36rs+14s−2r and EΓ1=48rs+16s−4r. RKr,s consists of 12rs+10s+2r vertices of degree 2 and 24rs+4s−4r vertices of degree 3. It is picturised in Figure 1. The edge partition table of rectangular Kekulene system RKr,s is represented in Table 3.Figure 1
Rectangular Kekulene system.Table 3
Edge partition table ofRKr,s.
degΓ1l1,degΓ1m1Total number of edges2,24s+2r2,324rs+12s3,324rs−6r


## 3.2. First Zagreb Entropy

By using the edge partition in Table3,(5)M1Γ1=264s−28r+76s.Then, using the equation of first Zagreb entropy, the result is obtained as(6)ENTM1Γ1=logM1Γ1−1M1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1=logM1Γ−1M1Γlog∏l1m1∈E2,2degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1×∏l1m1∈E2,3degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1×∏l1m1∈E3,3degΓ1l1+degΓ1m1degΓ1l1+degΓ1m1=log264s−28r+76s−1264s−28r+76slog4s+2r×44×24rs+12s×55×24rs−6r×66.


## 3.3. Second Zagreb Entropy

By using the edge partition in Table3,(7)M2Γ1=360s−46r+88s.Then, using the equation of second Zagreb entropy, the result is obtained as(8)ENTM2Γ1=logM2Γ1−1M2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1=logM2Γ1−1M2Γ1log∏l1m1∈E2,2degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1×∏l1m1∈E2,3degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1×∏l1m1∈E3,3degΓ1l1⋅degΓ1m1degΓ1l1⋅degΓ1m1=log360s−46r+88s−1360s−46r+88slog4s+2r×44×24rs+12s×66×24rs−6r×99.


## 3.4. Second Modified Zagreb Entropy

By using the edge partition in Table3,(9)M2mΓ=1640s−1r+3s.Then, using the equation of the second modified Zagreb entropy, the result is obtained as(10)ENTM2mΓ1=logM2mΓ1−1M2mΓ1log∏l1m1∈EΓ11degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=logM2mΓ1−1M2mΓ1log∏l1m1∈E2,21degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=log1640s−1r+3s−11/640s−1r+3slog4s+2r×141/4×24rs+12s×161/6×24rs−6r×191/9.


## 3.5. Reduced Second Zagreb Entropy

By using the edge partition in Table3,(11)RM2Γ=144s−22r+28s.Then, using the equation of reduced second Zagreb entropy, the result is obtained as(12)ENTRM2Γ1=logRM2Γ1−1RM2Γ1log∏l1m1∈EΓ1degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1=logM2mΓ1−1M2mΓ1log∏l1m1∈E2,2degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1×∏l1m1∈E2,3degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1×∏l1m1∈E3,3degΓ1l1−1⋅degΓ1m1−1degΓ1l1−1⋅degΓ1m1−1=log144s−22r+28−1144s−22r+28log4s+2r×11×24rs+12s×22×24rs−6r×44.


## 3.6. Hyper Zagreb Entropy

By using the edge partition in Table3,(13)HMΓ1=1464s−184r+364s.Then, using the equation of reduced second Zagreb entropy, the result is obtained as(14)ENTHMΓ1=logHMΓ1−1HMΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12=logHMΓ1−1HMΓ1log∏l1m1∈E2,2degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12×∏l1m1∈E2,3degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12×∏l1m1∈E3,3degΓ1l1+degΓ1m12degΓ1l1+degΓ1m12=log1464s−184r+364s−11464s−184r+364slog4s+2r×1616×24rs+12s×2525×24rs−6r×3636.


## 3.7. Augmented Zagreb Entropy

By using the edge partition in Table3,(15)AΓ1=13214892s−1675r+128s.Then, using the equation of reduced second Zagreb entropy, the result is obtained as(16)ENTAΓ1=logAΓ1−1AΓ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23=logAΓ1−1AΓ1log∏l1m1∈E2,2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23×∏l1m1∈E2,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23×∏l1m1∈E3,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−23degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1−23=log13214892s−1675r+128s−11/3214892s−1675r+128slog4s+2r×88×24rs+12s×88×24rs−6r×72964729/64.


## 3.8. Atom Bond Connectivity Entropy

By using the edge partition in Table3,(17)ABCΓ1=12s+1r+8s2+16s−4r.Then, using the equation of atom bond connectivity entropy, the result is obtained as(18)ENTABCΓ1=logABCΓ1−1ABCΓ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1=logABCΓ1−1ABCΓ1log∏l1m1∈E2,2degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,3degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,3degΓ1l1+degΓ1m1−2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1−2/degΓ1l1⋅degΓ1m1=log12s+1r+8s2+16s−4r−112s+1r+8s2+16s−4rlog4s+2r×121/2×24rs+12s×121/2×24rs−6r×232/3.


## 3.9. Geometric Arithmetic Entropy

By using the edge partition in Table3,(19)GAΓ1=4s−4r+4856rs+2456s+24rs.Then, using the equation of geometric arithmetic entropy, the result is obtained as(20)ENTGAΓ1=logGAΓ1−1GAΓ1log∏l1m1∈EΓ12degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=logGAΓ1−1GAΓ1log∏l1m1∈E2,22degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E2,32degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E3,32degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m12degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=log4s−4r+4856rs+2456s+24rs−14s−4r+48/56rs+24/56s+24rslog4s+2r×11×24rs+12s×26526/5×24rs−6r×11.


## 3.10. Symmetric Division Deg (SDD) Entropy

By using the edge partition in Table3,(21)SDDΓ1=100rs−8r+34s.Then, using the equation of Symmetric division deg entropy, the result is obtained as(22)ENTSDDΓ1=logSDDΓ1−1SDDΓ1log∏l1m1∈EΓ1degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1=logSDDΓ1−1SDDΓ1log∏l1m1∈E2,2degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1×∏l1m1∈E2,3degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1×∏l1m1∈E3,3degΓ1l1degΓ1m1+degΓ1m1degΓ1l1degΓ1l1/degΓ1m1+degΓ1m1/degΓ1l1=log100rs−8r+34s−1100rs−8r+34slog4s+2r×22×24rs+12s×13613/6×24rs−6r×22.


## 3.11. Randić Entropy

By using the edge partition in Table3,(23)RΓ1=46rs+26s+8rs−r+2s.Then, using the equation of Randić entropy, the result is obtained as(24)ENTRΓ1=logRΓ1−1RΓ1log∏l1m1∈EΓ11degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=logRΓ1−1RΓ1log∏l1m1∈E2,21degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,31degΓ1l1⋅degΓ1m11/degΓ1l1⋅degΓ1m1=log46rs+26s+8rs−r+2s−146rs+26s+8rs−r+2slog4s+2r×121/2×24rs+12s×161/6×24rs−6r×131/3.


## 3.12. Forgotten Entropy

By using the edge partition in Table3,(25)FΓ1=744rs−92r+188s.Then, using the equation of forgotten entropy, the result is obtained as(26)ENTFΓ1=logFΓ1−1FΓ1log∏l1m1∈EΓ1degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12=logFΓ1−1FΓ1log∏l1m1∈E2,2degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12×∏l1m1∈E2,3degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12×∏l1m1∈E3,3degΓ1l12+degΓ1m12degΓ1l12+degΓ1m12=log744rs−92r+188s−1744rs−92r+188slog4s+2r×88×24rs+12s×1313×24rs−6r×1818.


## 3.13. Sum-Connectivity Entropy

By using the edge partition in Table3,(27)χΓ1=2s+r+2455rs+1255s+46rs−6r.Then, using the equation of sum-connectivity entropy, the result is obtained as(28)ENTχΓ1=logχΓ1−1χΓ1log∏l1m1∈EΓ11degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1=logFΓ1−1FΓ1log∏l1m1∈E2,21degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1×∏l1m1∈E2,31degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1×∏l1m1∈E3,31degΓ1l1+degΓ1m11/degΓ1l1+degΓ1m1=log2s+r+2455rs+1255s+46rs−6r−12s+r+24/55rs+12/55s+46rs−6rlog4s+2r×121/2×24rs+12s×151/5×24rs−6r×161/6.


## 3.14. First Redefined Zagreb Entropy

By using the edge partition in Table3,(29)ReZG1Γ1=36rs−2r+14s.Then, using the equation of first redefined Zagreb entropy, the result is obtained as(30)ENTReZG1Γ1=logReZG1Γ1−1ReZG1Γ1log∏l1m1∈EΓ1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1=logReZG1Γ1−1ReZG1Γ1log∏l1m1∈E2,2degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1×∏l1m1∈E2,3degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1×∏l1m1∈E3,3degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1/degΓ1l1⋅degΓ1m1=log36rs−2r+14s−136rs−2r+14slog4s+2r×11×24rs+12s×565/6×24rs−6r×232/3.


## 3.15. Second Redefined Zagreb Entropy

By using the edge partition in Table3,(31)ReZG1Γ1=925s−7r+3245rs.Then, using the equation of the second redefined Zagreb entropy, the result is obtained as(32)ENTReZG2Γ1=logReZG2Γ1−1ReZG2Γ1log∏l1m1∈EΓ1degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=logReZG2Γ1−1ReZG2Γ1log∏l1m1∈E2,2degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E2,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1×∏l1m1∈E3,3degΓ1l1⋅degΓ1m1degΓ1l1+degΓ1m1degΓ1l1⋅degΓ1m1/degΓ1l1+degΓ1m1=log925s−7r+3245rs−192/5s−7r+324/5rslog4s+2r×11×24rs+12s×656/5×24rs−6r×323/2.Theorem 1.
IfΓ1 is the RKr,s, then the irregularity-based descriptors are expressed as follows:(1)
VARΓ1=236r2s2+36rs2−r2−4rs+5s2/18rs−r+7s2(2)
IR1Γ1=4180r2s2+1434rs2−5r2−115rs+462s2/18rs−r+7s(3)
IR2Γ1=1/2182180rs−23r+44s/12rs−r+4srs−2180rs−23r+44s/12rs−r+4sr+72180rs−23r+44s/12rs−r+4ss−96rs+8r−32s/18rs−r+7s(4)
IRDIFΓ1=1/324rs+12s(5)
ALΓ1=24rs+12s(6)
IRLΓ1=0.4054524rs+12s(7)
IRLUΓ1=0.40545/224rs+12s(8)
IRLFΓ1=1/624rs+12s(9)
IRFΓ1=24rs+12s(10)
IRLAΓ1=2/524rs+12s(11)
IRD1Γ1=0.6931524rs+12s(12)
IRAΓ1=0.01683224rs+12s(13)
IRBΓ1=0.1010624rs+12s(14)
IRBΓ1=1/22166r2s2−126r2s+1926rs2−504r2s2−66rs+426s2+30r2s−444rs2−r2+11rs−100s2/12rs−r+4s18rs−r+7s(15)
IRGAΓ1=0.02039124rs+12s(16)
IRRtΓ1=1/224rs+12s


## 4. Numerical and Graphical Analysis ofRKr,s

In this section, by changing the values of the variablesrands from 1 to 9, the numerical values of VDB entropy measures are obtained using the analytical expressions of the rectangular Kekulene system. Tables 4 and 5 show the results of these calculations. By looking at their numerical values, the comparison can be done easily to see whether there are any similarities or differences between the individual topological indices. Figures 2 and 3 show this behaviour as three-dimensional graphical representations. The differences between each topological index for a specific structure can be seen in these 3D plots.Table 4
Comparison table ofRKr,s.
RKr,sENTM1Γ1ENTM2Γ1ENTM2mΓ1ENTRM2Γ1ENTHMΓ1ENTAΓ1ENTABCΓ11,15.655.891.564.927.256.173.572,27.027.303.176.308.697.574.943,37.828.114.057.129.518.385.744,48.388.684.667.7110.088.956.315,58.839.125.128.1610.539.396.75Table 5
Comparison table ofRKr,s.
RKr,sENTGAΓ1ENTSDDΓ1ENTRΓ1ENTFΓ1ENTχΓ1ENTReZG1Γ1ENTReZG2Γ11,13.964.842.896.603.033.724.242,25.316.084.308.034.445.065.603,36.116.875.118.845.255.856.404,46.687.435.699.415.836.416.975,57.127.876.139.856.276.857.42Figure 2
Graphical visualisation for Table4.Figure 3
Graphical visualisation for Table5.As it is evident from the above calculations, it is quite convenient to compute topological indices. Since topological indices correlate well with physiochemical properties, these can be directly used to predict the behaviour of compounds. Such a theoretical analysis based on topology significantly reduces the time and efforts required for analysis of compounds, when compared to intensive quantum chemical calculations or laborious experiments. But, not all properties of a compound can be mathematically associated to topological indices. Some features would still require practical experiments to be studied. The prediction accuracy of QSAR/QSPR methods which use TIs can also be comparatively less than that of experimental observations. However, TIs have been widely accepted by the scientific community as predictive methods for a subset of physiochemical properties due to their ease of use and universal applicability.


## 5. Conclusion

In this paper, the mathematical expressions of VDB entropy measures and irregularity-based indices have been computed forRKr,s. Furthermore, the numerical values for these expressions of VDB entropy measures and irregularity-based indices are calculated using the obtained analytical expressions. The 3D graphical comparisons of these analytical expressions of VDB entropy measures have been represented visually to find the behavioural pattern of the molecular compound. The study of Wiener polarity-based TIs, eccentricity-based TIs, and the recent development of VDB indices for this structure can also be explored as they have never been studied earlier. It is still an area where further research could be done in the future.



---

*Source: 7404529-2021-12-14.xml*