The bond incident degree (BID) indices can be written as a linear combination of the number of edges xi,j with end vertices of degree i and j. We introduce two transformations, namely, linearizing and unbranching, on catacondensed pentagonal systems and show that BID indices are monotone with respect to these transformations. We derive a general expression for calculating the BID indices of any catacondensed pentagonal system with a given number of pentagons, angular pentagons, and branched pentagons. Finally, we characterize the CPSs for which BID indices assume extremal values and compute their BID indices.
National Key R&D Program of China2019YFA0706402Natural Science Foundation of Guangdong Province2018A03031301151. Introduction
A pentagonal system is a connected geometric figure obtained by concatenating congruent regular pentagons side to side in a plane in such a way that the figure divides the plane into one infinite (external) region and a number of finite (internal) regions, and all internal regions must be congruent regular pentagons. In a pentagonal system, two pentagons are adjacent if they share an edge. In this work, we consider only catacondensed pentagonal systems (CPSs), i.e., pentagonal systems which have no internal vertices. This is equivalent to say that pentagonal system has only pentagons of type L1,L2,A2, and A3 (see Figure 1). We will denote by l1,l2,a2, and a3, respectively, the number of L1,L2,A2, and A3 pentagons the CPS has. We will denote the set of all CPSs with n pentagons by CHn. For more details on pentagonal systems and this theory, we refer the readers to [1]. We are interested in studying the behavior of general bond incident degree (BID) indices over pentagonal systems. For details about the pentagonal systems, see [2–7]. One of the much studied topological indices was put forward by the Randić index [8]. For a graph G, it is denoted by RG and is defined as(1)RG=∑uvdudv−1/2,where du denotes the degree of the vertex u and sum runs over all edges uv of G. There are many topological indices which appeared in the literature of mathematical chemistry (see for example [9–13]). Among the degree-based topological descriptors, the most studied are the first and second Zagreb indices [14–17], the sum-connectivity index [17–19], the atom-bond connectivity index [17, 20], the augmented Zagreb index [17, 21, 22], the geometric arithmetic index [23–25], and the harmonic index [19, 26, 27]. Each of these BID indices can be expressed [28] as(2)BIDG=∑1≤i≤j≤Δθi,jxi,j,Δ is the maximum degree of G,where xi,j is the number of edges in G with end vertices of degree i and j and θi,j is a finite sequence of nonnegative real numbers. For instance, if BID is the Randić index, then θi,j=1/ij. From now on, BID indices of any graph G induced by the sequence θi,j are defined by equation (2). In case of CPS, we have only vertices of degree 2, 3, and 4; therefore the general BID indices over CHn will be induced by a sequence θ22,θ23,θ24,θ33,θ34 of nonnegative real numbers:(3)BIDG=θ22x22+θ23x23+θ24x24+θ33x33+θ34x34,for every G∈CHn. Now, following [29], we define two transformations on a CPS, called linearizing and unbranching transformations. We will show that BID indices are monotone with respect to these transformations. This will give us a reduction procedure to compute the BID indices of any CPS. Finally, the results obtained will be applied to find the minimal and maximal values of BID indices over CHn.
Types of pentagons.
2. Linearizing and Unbranching a Catacondensed Pentagonal System
In this section, we define two transformations: the linearizing transformation T1 and the unbranching transformation T2, for the catacondensed pentagonal systems. In order to define these two transformations, we first recall the definition of coalescence of two graphs [30]. Suppose G1 and G2 are graphs with edges u1v1∈EG1 and u2v2∈EG2. The coalescence of G1 and G2 with respect to edges u1v1 and u2v2 denoted by G is formed by identifying the edges u1v1 and u2v2 in the following way:(4)VG=VG1−u1,v1∪VG2−u2,v2∪u,v,and the two edges in G are adjacent if
These are adjacent in G1 or G2, or
One is u and the other one is adjacent to u1 in G1 or to u2 in G2, or
One is v and the other one is adjacent to v1 in G1 or v2 in G2, or
First one is u and the second one is v
Now, we define the linearizing transformation on a catacondensed pentagonal system. Let P1 be catacondensed pentagonal system with angular pentagon P and sub-catacondensed pentagonal systems Q1 and Q2 as shown in Figure 2. Let P2 be the catacondensed pentagonal system obtained by applying coalescing operation on P1 that is moving the catacondensed pentagonal systems Q2 to the 2−2 edge of angular pentagon P. We call the pentagonal system P2 the linearization of pentagonal system P1 at the pentagon P.
The linearization transformation T1.
Let Θ1=−θ22+2θ23−2θ24+3θ33−2θ34; then, the difference of BID indices of the catacondensed pentagonal systems P1 and P2 is calculated in the next theorem.
Theorem 1.
Let P2 be a CPS obtained from P1 by linearizing at angular pentagon P. Then,(5)BIDP2−BIDP1=Θ1.
Proof.
In Figure 2, bold edges in P1 and P2 are the ones whose degree will change after the linearizing transformation. We collect all the bold edges in the set E from P1. Similarly, collect all the bold edges in the set F from P2. Note that the set of edges EP1∖E and EP2∖F is in one to one correspondence in such a way that for each edge in EP1∖E, the degree of end vertices is equal to those of corresponding edge in EP2∖F. Hence,(6)BIDP2−BIDP1=4θ23+3θ33−θ22+2θ23+2θ24+2θ34=−θ22+2θ23−2θ24+3θ33−2θ34=Θ1.
If Θ1<0, then from above theorem, we have BIDP2<BIDP1. Hence, by applying a linearizing transformation, we can construct a new catacondensed pentagonal system whose BID indices are greater than the BID indices of original catacondensed pentagonal system. This fact is reflected in the next example.
Example 1.
The sequence of CPS in Figure 3 satisfies(7)BIDP1>BIDP2>BIDP3>BIDP4>BIDP5>BIDP6,when Θ1<0.
BIDP1>BIDP2>BIDP3>BIDP4>BIDP5>BIDP6.
Next, we define the unbranching transformation on catacondensed pentagonal system. Let Z1 be CPS with branched pentagon R. Suppose that Q2 is a pentagonal chain in Z1. Note that such branched pentagon always exists in Z1. Let Z2 be a catacondensed pentagonal system obtained by applying coalescing operation on Z1 that is moving the sub-catacondensed pentagonal system Q3 to any of the 2-2 edge in the last pentagon of Q2. This process is shown in Figure 4. Let Θ2=θ22−2θ23+θ33; then, the difference of BID indices of the catacondensed pentagonal system Z1 and Z2 can be expressed as the difference of Θ1 and Θ2. This is proved in the next theorem.
The unbranching transformation T2.
Theorem 2.
Let Z2 be a catacondensed pentagonal system obtained from Z1 by unbranching at R; then,(8)BIDZ2−BIDZ1=Θ1−Θ2.
Proof.
In Figure 4, bold edges in Z1 and Z2 are the ones whose degree will change after the unbranching transformation. We collect all the bold edges in the set E from Z1. Similarly, collect all the bold edges in the set F from Z2. Note that the set of edges EZ1∖E and EZ2∖F is in one to one correspondence in such a way that for each edge in EZ1∖E, the degree of end vertices is equal to those of the corresponding edge in EZ2∖F. Hence,(9)BIDZ2−BIDZ1=5θ23+3θ33−2θ22+θ23+2θ34+2θ24+θ33=−2θ22+4θ23−2θ24+2θ33−2θ34=−θ22+2θ23−2θ24+3θ33−2θ34−θ22−2θ23+θ33=Θ1−Θ2.
Note that if Θ1−Θ2<0, then each of the BID indices defined on the catacondensed pentagonal system is monotone decreasing when we apply unbranching transformation. Next example depicts this fact.
Example 2.
The sequence of CPS in Figure 5 satisfies(10)BIDP1>BIDP2>BIDP3,when Θ1−Θ2<0.
BIDP1>BIDP2>BIDP3.
3. Method to Compute BID Indices of Catacondensed Pentagonal System
Examples 1 and 2 show that we can transform any catacondensed pentagonal system into linear pentagonal chain by successively applying linearizing and unbranching transformations. The number of steps depends on the number of angular pentagons a2 and number of branched pentagons a3. Theorems 1 and 2 show that we can find the exact value of topological indices after applying these transformations. Denote Ln by a linear pentagonal chain with n pentagons; then, the topological indices of Ln can be computed by the following formula:(11)BIDLn=4θ22+2nθ23+2n−3θ33.
Our next results show that once we know the number of angular pentagons a2, the number of branched pentagons a3, and the number of pentagons n in a catacondensed pentagonal system, we can compute its topological indices.
Lemma 1.
If K is a pentagonal chain with a2 angular pentagons, then(12)BIDK=BIDLn−a2Θ1.
Proof.
We will prove it by induction on a2. If a2=0, then K=Ln and the result follows. Suppose that the result holds valid for pentagonal chain with less than a2 pentagons. Let K be a pentagonal chain with a2>0 angular pentagons and let L be the pentagonal chain obtained by applying linearizing transformation on an angular pentagon K1 of K. Then, by Theorem 1,(13)BIDL−BIDK=Θ1.
Since the pentagonal chain L has a2−1 angular pentagons, then by induction hypothesis,(14)BIDL=BIDLn−a2−1Θ1.
Our next theorem generalizes Lemma 1 to any catacondensed pentagonal system.
Theorem 3.
Let M be a catacondensed pentagonal system with n pentagons, a2 angular pentagons, and a3 branched pentagons; then,(16)BIDM=BIDLn−a2Θ1−a3Θ1−Θ2.
Proof.
We will prove it by induction on a3. If a3=0, then by Lemma 1, the result holds. Suppose that the result is true for catacondensed pentagonal systems with less than a3>0 branched pentagons and let M be a catacondensed pentagonal system with a3 branched pentagons. By applying unbranching transformation at branch M1 of M, we obtain a pentagonal system N which has a3−1 branched pentagons and a2 angular pentagons. Hence, by Theorem 2,(17)BIDN−BIDM=Θ1−Θ2.
Next example details the computation of BID topological indices of any catacondensed pentagonal system using Theorem 3.
Example 3.
We compute the topological indices of catacondensed pentagonal system P1 shown in Figure 1. It contains 1 angular pentagon, 1 branched pentagon, and in total 10 pentagons. Hence, by Theorem 3, we have(20)BIDP1=BIDL10−Θ1−Θ1−Θ2=BIDL10−2Θ1+Θ2=4θ22+16θ23+13θ33−2−θ22+2θ23−2θ24+3θ33−2θ34+θ22−2θ23+θ33=7θ22+10θ23+4θ24+8θ33+4θ34.
For example, if the BID topological index is the Randić index, then θij=1/ij and so(21)RP1=712+1016+418+813+4112=14.7701.
If the BID index is the first Zagreb index M1P1, then θij=i+j, and(22)M1P1=7×4+10×5+4×8+8×6+4×7=186.
Corollary 1.
If Θ1<0 and Θ1−Θ2<0, then the minimal value of BID indices over CPn is attained in Ln.
If Θ1>0 and Θ1−Θ2>0, then the maximal value of BID indices over CPn is attained in Ln.
Proof.
Suppose that Θ1<0 and Θ1−Θ2<0; then, by Theorem 3, the minimum value of BID indices is attained when a2=0 and a3=0. Hence, our catacondensed pentagonal system is a linear pentagonal system. Similarly, when Θ1>0 and Θ1−Θ1>0, then the maximum value is attained when a2=0 and a3=0 which occurs in the linear pentagonal system.
From Theorem 3, it is clear that the maximum and the minimum values of BID indices over CPn depend on the values of a2, a3, Θ1, and Θ2. Also, for any catacondensed pentagonal system CPn with n≥2, we have l1≥2 and 0≤l2,a2,a3≤n−2. Let En and Fn denote the CPSs with n pentagons as shown in Figures 6 and 7.
Catacondensed pentagonal system En for n=7.
Catacondensed pentagonal system Fn for n = odd and n = even.
Corollary 2.
If Θ1>0 and Θ1−Θ2<0, then the minimal value of BID indices over CPn is attained in En.
If Θ1>0 and Θ1−Θ2<0, then the maximal value of BID indices over CPn is attained in Fn.
Proof.
(i) Suppose that Θ1>0 and Θ1−Θ2<0; then, by Theorem 3, the minimum value of BID indices is attained when a3=0. Also, the maximum value of a2 can be n−2 which occurs in En.
(ii) By Theorem 3, we get the maximum value of BID indices when a2=0.
In a catacondensed pentagonal system, we have(23)l1=a3+2,n=a2+l2+a3+l1.
We deduce(24)a3=12n−a2+l2+2=12n−l2+2.
If n is even (odd), then the maximum value of a3 will be obtained when l2=0 (l2=1) and this occurs in Fn.
Corollary 3.
If Θ1<0 and Θ1−Θ2>0, then the minimal value of BID indices over CPn is attained in Fn.
If Θ1<0 and Θ1−Θ2>0, then the maximal value of BID indices over CPn is attained in En.
Example 4.
For the second Zagreb index, we have θij = ij. Then, Θ1=−5<0 and Θ1−Θ2=−6<0. Hence, from Corollary 1, the minimum value of BID indices is attained in Ln.
If BID indices are among the Randić, geometric arithmetic, harmonic, and sum-connectivity indices, then Θ1>0 and Θ1−Θ2<0. Hence, from Corollary 2, the maximum value of BID indices is attained in Fn and the minimum value is attained in En.
For the atom-bond connectivity index, we have θij=i+j−2/i+j. In this case, Θ1>0 and Θ1−Θ2>0; hence, from Corollary 1, the maximum value of BID indices is attained in Ln.
In Table 1, we have computed the values of Θ1 and Θ2 for different BID indices. This helps us to find the CPS which has the maximum and minimum values of the following BID indices.
Values of Θ1 and Θ2 for different BID indices.
BID indices
i+j
ij
1/ij
2ij/i+j
2/i+j
1/i+j
ij3/i+j−23
i+j−2/i+j
Θ1
−2
−5
0.032
0.0944
0.0619
0.0467
−1.476
0.0019
Θ2
0
1
0.168
0.404
0.333
0.138
3.3906
−0.404
Maximal
En
En
En
En
Ln
Minimal
Ln
Ln
Fn
Fn
Fn
Fn
Ln
4. Concluding Remarks
Among the well-known topological indices, there are various bond incident degree (BID) indices. In this paper, we have studied these BID indices for the catacondensed pentagonal systems and derived a general expression for calculating the BID indices of any catacondensed pentagonal system with a given number of pentagons, angular pentagons, and branched pentagons. We have also characterized the systems having maximum/minimum well-known BID indices from the class of all catacondensed pentagonal systems with a fixed number of pentagons. The present study can be extended in several directions. One of such directions is to study the general pentagonal systems for the BID indices. Also, there are many BID indices (see, for example, Table 1) for which we have not been able to characterize the systems having maximum/minimum values from the class of all catacondensed pentagonal systems with a fixed number of pentagons.
Data Availability
No data were used to support the study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the National Key R&D Program of China (grant no. 2019YFA0706402) and the Natural Science Foundation of Guangdong Province (grant no. 2018A0303130115).
AliA.RazaZ.BhattiA. A.Extremal pentagonal chains with respect to bond incident degree indices2016941087087610.1139/cjc-2016-03082-s2.0-84991607074GutmanI.YanW.YehY.-N.YangB.-Y.Generalized Wiener indices of zigzagging pentachains200742210311710.1007/s10910-006-9078-92-s2.0-34547249219XiaoC.ChenH.LiuL.Perfect matchings in random pentagonal chains20175591878188610.1007/s10910-017-0767-32-s2.0-85021070632RouvrayD. H.KingR. B.RouvrayD. H.KingR. B.Harry in the limelight: the life and times of Harry Wiener2002111510.1533/9780857099617.1WangY.ZhangW. W.Kirchhoff index of linear pentagonal chains201811015941604WangH.-Y.QinJ.GutmanI.Wiener numbers of random pentagonal chains201345976XiaoC.ChenH.RaigorodskiiA. M.A connection between the Kekulé structures of pentagonal chains and the Hosoya index of caterpillar trees201723223023410.1016/j.dam.2017.07.0242-s2.0-85028833261RandicM.Characterization of molecular branching197597236609661510.1021/ja00856a0012-s2.0-8644280181HoroldagvaB.GutmanI.On some vertex degree based graph invariants201165723730HayatS.ImranM.LiuJ. B.Correlation between the Estrada index and p-electronic energies for benzenoid hydrocarbons with applications to boron nanotubes201912923ImranM.SiddiquiM. K.AhmadS.HanifM. F.MuhammadM. H.FarahaniM. R.Topological properties of benzenoid, phenylenes and nanostar dendrimers20192271229124810.1080/09720529.2019.1701267ZogićE.GlogicE.A note on the Laplacian resolvent energy, Kirchhoff index and their relations201923237ZhangL. Z.TianW. F.Extremal catacondensed benzenoids2003341111122GutmanI.TrinajstićN.Graph theory and molecular orbitals. Total φ-electron energy of alternant hydrocarbons197217453553810.1016/0009-2614(72)85099-12-s2.0-33845352210AliA.GutmanI.MilovanovićE.MilovanovićI.Sum of powers of the degrees of graphs: extremal results and bounds201880584BorovićaninB.DasK. C.FurtulaB.GutmanI.Bounds for Zagreb indices20177817100GutmanI.Degree-based topological indices201386435136110.5562/cca22942-s2.0-84897584934ZhouB.TrinajstićN.On a novel connectivity index20094641252127010.1007/s10910-008-9515-z2-s2.0-70350151674AliA.ZhongL.GutmanI.Harmonic index and its generalizations: extremal results and bounds201981249311EstradaE.TorresL.Rodrí GuezL.GutmanI.An atom bond connectivity index modeling the enthalpy of formation of alkanes199837-A849855FurtulaB.GraovacA.VukičevićD.Augmented Zagreb index201048237038010.1007/s10910-010-9677-32-s2.0-77954457136AliA.FurtulaB.GutmanI.VukičevićD.Augmented Zagreb index: extremal results and bounds2021852Vukic̃evićD.FurtulaB.Topological index based on the ratios of geometrical and arithmetical means of end vertex degrees of edges20094613691376DasK. C.GutmanI.FurtulaB.Survey on geometricarithmetic indices of graphs201165595644PortillaA.RodríguezJ. M.SigarretaJ. M.Recent lower bounds for geometric-arithmetic index201915982ZhongL.The harmonic index for graphs201225356156610.1016/j.aml.2011.09.0592-s2.0-80955166259FajtlowiczS.On conjectures of Graffiti II198760187197AliA.RazaZ.BhattiA. A.Bond incident degree (BID) indices of polyomino chains: a unified approach2016287-288283710.1016/j.amc.2016.04.0122-s2.0-84968756914RadaJ.CruzR.GutmanI.Vertex-degree-based topological indices of catacondensed hexagonal systems201357215415710.1016/j.cplett.2013.04.0322-s2.0-84878114855RadaJ.Energy ordering of catacondensed hexagonal systems2005145343744310.1016/j.dam.2004.03.0072-s2.0-9344260294