Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Chemistry

formerly Central European Journal of Chemistry

IMPACT FACTOR 2018: 1.512
5-year IMPACT FACTOR: 1.599

CiteScore 2018: 1.58

SCImago Journal Rank (SJR) 2018: 0.345
Source Normalized Impact per Paper (SNIP) 2018: 0.684

ICV 2017: 165.27

Open Access
See all formats and pricing
More options …
Volume 16, Issue 1


Volume 13 (2015)

M-Polynomials And Topological Indices Of Zigzag And Rhombic Benzenoid Systems

Ashaq Ali / Waqas Nazeer / Mobeen Munir / Shin Min Kang
  • Corresponding author
  • Department of Mathematics and RINS, Gyeongsang National University, Jinju, 52828, Korea
  • Center for General Education, China Medical University, Taichung 40402, Taiwan
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2018-02-24 | DOI: https://doi.org/10.1515/chem-2018-0010


M-polynomial of different molecular structures helps to calculate many topological indices. This polynomial is a new idea and its beauty is the wealth of information it contains about the closed forms of degree-based topological indices of molecular graph G of the structure. It is a well-known fact that topological indices play significant role in determining properties of the chemical compound [1, 2, 3, 4]. In this article, we computed the closed form of M-polynomial of zigzag and rhombic benzenoid systemsbecause of their extensive usages in industry. Moreover we give graphs of M-polynomials and their relations with the parameters of structures.

Keywords: M-polynomial; Topological index; Benzenoid

1 Introduction

In mathematical chemistry, mathematical tools such as polynomials and numbers predict properties of compounds without using; quantum mechanics. These tools, in combination, capture information hidden in the symmetry of molecular graphs. Most commonly known invariants of such kinds are known as degree-based topological indices. These are the numerical values that correlate the structure with various physical properties, chemical reactivity and biological activities [5, 6, 7, 12]. It is an established fact that many properties such as heat of formation, boiling point, strain energy, rigidity and fracture toughness of a molecule are strongly connected to its graphical structure and this fact plays a synergic role in chemical graph theory. A graph G with vertex set V(G) and edge set E(G) is connected, if there exists a connection between any pair of vertices in G. The distance between two vertices u and v is denoted as d (u, v) and is the length of the shortest path between u and v in graph G. The number of vertices of G, adjacent to a given vertex v, is the “degree” of this vertex, and will be denoted by dv. For details on basics of graph theory, any standard text such as [13] can be of great help.

Several algebraic polynomials have useful applications in chemistry such as Hosoya Polynomial (also called Wiener polynomial) [8]. M-polynomial [14], introduced in 2015 helps in determining many degree-based topological indices. Benzenoid hydrocarbons play a vital role in our environment, and in the food and chemical industries. Benzenoid molecular graphs are systems with deleted hydrogens. It is a connected geometric figure obtained by arranging congruent regular hexagons in a plane, so that two hexagons are either disjoint or have a common edge. This figure divides the plane into one infinite (external) region and a number of finite (internal) regions. All internal regions must be regular hexagons. Benzenoid systems are of considerable importance in theoretical chemistry because they are the natural graph representation of benzenoid hydrocarbons. A vertex of a hexagonal system belongs to, at most, three hexagons. A vertex shared by three hexagons is called an internal vertex. Under this definition in [15] the figure under discussion is not a benzenoid system as one internal central region is a regular 14 sided polygon. These systems are planar consisting of regular hexagons as shown in Figure 1.

Example of benzenoid system left and a non – benzenoid system right.
Figure 1

Example of benzenoid system left and a non – benzenoid system right.

Definition 1

Let G be a simple connected graph. The M-polynomial of G is defined as: M(G,x,y)=δijmij(G)xiyj.

Where δ = Mind{dv | v ∈ V (G)}, Δ = Max{dv ∈ V (G)}, and mij(G) is the edge vuE(G) such that {dv, du} = {i, j}.

This polynomial has been one of the key areas of interest in computational aspects of materials. From this M-polynomial, we can calculate many topological indices. The topological index of a molecule structure can be considered as a non-empirical numerical quantity which quantifies the molecular structure and its branching pattern in many ways. M-polynomial of different molecular structures have been computed in [9, 10, 12, 16, 17,]. Essentially, the topological index can be regarded as a score function which maps each molecular structure to a real number and is used as a descriptor of the molecule under testing [1, 18, 19, 20]. Topologioal indices provide a good prediction of various physico-chemical properties of chemical compounds including boiling point, heat of evaporation, heat of formation, chromatographic retention times, surface tension, vapor pressure etc. Since the 1970s, two degree based graph invariants have been extensively studied. These are the first Zagreb index M1 and the second Zagreb index M2, introduced by Gutman and Trinajstic’ [2] and defined as: M1(G)=vV(G)(dv)2 and M2(G)=uvE(G)dudv.

Results obtained in the theory of Zagreb indices are summarized in the review [21].

Second modified Zagreb index is defined as: mM2(G)=uvE(G)1dudv.

In 1998, working independently, Bollobas and Erdos [22] and Amic et al. [3] proposed general Randic index. It has been extensively studied by both mathematicians and theoretical chemists

(See, for example, [23]). The Randic’ index is defined as: Rα(G)=uvE(G)(dudv)α,

where α is an arbitrary real number.

Where Dx=x(f(x,y)x,Dy=y(f(x,y)y,Sx=0xf(t,y)tdf,Sy==0yf(x,t)tdf,J(f(x,y))=f(x,x),Qα(f(x,y))=xαf(x,y).

Symmetric division index is defined as: SDD(G)=uvE(G)min(du,dv)max(du,dv)+max(du,dv)min(du,dv).

Another variant of Randic’ index is the harmonic index defined as: H(G)=vuE(G)2du+dv.

The Inverse sum index is defined as: I(G)=vuE(G)dudvdu+dv.

The augmented Zagreb index is defined as: A(G)=vuE(G)dudvdu+dv23,

and it is useful for computing heat of formation of alkanes [24, 25].

For detailed study about degree-based topological indices, we refer [26, 27, 28, 29, 30, 31, 32] and the references therein.

These topological indices can be recovered from M-polynomial [14], see following Table 1.

Table 1

Derivation of some degree-based topological indices from M-polynomial.

In this article, we compute the closed form of the M-polynomial for two famous benzenoid systems Zigzag benzenoid system and Rhombic benzenoid system. We also computed some degree-based topological indices.

2 Methodology

At first we obtain general pattern of vertex and edge partitions of 2D molecular graph connected to the Zigzag benzenoid system and Rhombic benzenoid system based on the degree of end vertices of edges. From this edge partition, using definition, we obtain M-polynomials of these systems. The 3D graph of M-polynomials are sketched by using maple 2015. Then using mathematical operators and table 1, we reach at the different degree-based indices.

Ethical approval: The conducted research is not related to either human or animals use.

3 Results and Discussions

In this part we give our main computational results in two sections.

3.1 Computational aspects of Zigzag benzenoid system

Let n be number of rows in graph of zigzag benzenoid system Zn with two hexagons in each row. Since first row contain two hexagons with twelve edges and one edge is common so we obtain total eleven edges in the first row and combining first and second row we obtain 24 total edges with three edges in common so we obtain 21 different edges. Continuing in the same way we obtain 10n + 1 edges and 8n + 2 vertices. We partition edges on the basis of degrees of endpoints of edges of the graph. All vertices are either of degree two or three. One can observe that at each row we have two edges of type {2,2}, one upside and one downside the chain except the endpoints where we have two more edges with end vertex having degree two.

Thus |E{2,2}|=|{e=uvE(Zn)|du=2,dv=2}=2n+4. Next we can see that four edges of type (2,3) are present in each row of the chain so |E{2,3}|=|{e=uvE(Zn)|du==2dv=3}=4n. Remaining are edges of type (3,3) given as . |E{3,3}|=|{e=uvE(Zn)|du=3,dv=3}|=(10n+1)(2n+4)4n=4n3.

Theorem 1

Consider the zigzag benzenoid system Zn, then its M-polynomial is M(Zn;x,y)=2(n+2)x2y2+4nx2y3+(4n3)x3y3.


Let Zn be the zigzag benzenoid system, then from the above decision |V (Zn)| = 8n + 2 and |E (Zn)| = 10n + 1. Also from the above decision, we can divide the edge set into the following three partitions: E1(Zn)={e=uvE(Zn):du=dv=2},E2(Zn)={e=uvE(Zn):du=2,dv=3},E3(Zn)={e=uvE(Zn):du=dv=3}.

In addition, |E1(Zn)|=2(n+2),|E2(Zn)|=4n,|E3(Zn)|=4n3.

Now by definition of M-polynomial, we have M(Zn;x,y)=ijmijxiyj=22m22x2y2+23m23x2y3+33m33x3y3=uvE1(Zn)m22x2y2+uvE2(Zn)m23x2y3+uvE3(Zn)m33x3y3=|E1(Zn)|x2y2+|E2(Zn)|x2y3+|E3(Zn)|x3y3=2(n+2)x2y2+4nx2y3+(4n3)x3y3.

Now we compute some degree-based topological indices of zigzag benzenoid from this M-polynomial.

Graph of zigzag benzenoid system Zn.
Figure 2

Graph of zigzag benzenoid system Zn.

Proposition 2

Consider the zigzag benzenoid system Zn, then

  1. mM2(Zn)=2918n+23.

  2. Rα(Zn)=4n9α+2n4α+4n6α − 3 × 9α+4 × 4α.

  3. Rα(Zn)=2n+44α+4n6α+4n39α.

  4. SSD(Zn)=623n+2.

  5. H(Zn)=313n+1.

  6. I(Zn)=645n12.

  7. A(Zn)=145716n45964.

The plot for the M-polynomial of Z1.
Figure 3

The plot for the M-polynomial of Z1.


Let M(Zn; x, y) = f (x, y) = 2(n + 2)x2 y2 + 4nx2 y3 + (4n - 3)x3 y3. Then Dxf(x,y)=4(n+2)x2y2+8nx2y3+3(4n3)x3y3,Dyf(x,y)=4(n+2)x2y2+12nx2y3+3(4n3)x3y3,DyDxf(x,y)=8(n+2)x2y2+24nx2y3+9(4n3)x3y3,Sy(f(x,y))=(n+2)x2y2+43nx2y3+13(4n3)x3y3,SxSy(f(x,y))=12(n+2)x2y2+23nx2y3+19(4n3)x3y3,Dyα(f(x,y))=2α+1(n+2)x2y2+3α4nx2y3+3α(4n3)x3y3,DxαDyα(f(x,y))=22α+1(n+2)x2y2+2α+23αnx2y3+32α(4n3)x3y3,Syα(f(x,y))=12α1(n+2)x2y2+43αnx2y3+13α(4n3)x3y3,SxαSyα(f(x,y))=122α1(n+2)x2y2+12α23αnx2y3+132α(4n3)x3y3,SyDx(f(x,y))=2(n+2)x2y2+8xnx2y3+(4n3)x3y3,SxDy(f(x,y))=2(n+2)x2y2+6nx2y3+(4n3)x3y3,.Jf(x,y)=2(n+2)x4+4nx5+(4n3)x6,SxJf(x,y)=12(n+2)x4+45nx5+16(4n3)x6,JDxDyf(x,y)=8(n+2)x4+24nx5+9(4n3)x6,SxJDxDyf(x,y)=2(n+2)x4+245nx5+32(4n3)x6,Dy3f(x(y)=16(n+2)x2y2+108nx2y3+27(4n3)x3y3,Dx3Dy3f(x,y)=108(n+2)x2y2+864nx2y3+729(4n3)x3y3,JDx3Dy3f(x,y)=108(n+2)x4+864nx5+729(4n3)x6,Q2JDx3Dy3f(x,y)=108(n+2)x2+864nx3+729(4n3)x4,Sx3Q2JDx3Dy3f(x,y)=272(n+2)x2+32nx3+72964(4n3)x4.

  1. mM2(Zn)=SxSy(f(x,y))|x=y=1=2918n+23.

  2. Rα(zn)=DxαDyα(f(x,y))|x=y=1=4n9α+2n4α+4n6α3×9α+4×4α.

  3. Rα(zn)=SxαSyα(f(x,y))|x=y=1=2n+44α+4n6α+6n39α.

  4. SSD(Zn)=(SyDx+SxDy)(f(x,y))|x=y=1=623n+2.

  5. H(Zn)=2SxJ(f(x,y))|x=1=313n+1.

  6. I(Zn)=SxJDxDy(f(x,y))x=1=645n12.

  7. A(Zn)=Sx3Q2JDx3Dy3(f(x,y))|x=1=145716n45964.

3.2 Computational aspects of the Rhombic benzenoid system

Take another benzenoid system in which hexagons are arranged to form a rhombic shape Rn, in which there are n rows of n hexagons as given in Figure 4. Then it has 2n(2+ 2) vertices and 3n2 + 4n − 1 edges. On the similar lines we partition edges on the basis of degrees of endpoints. Clearly, all vertices are either of 2nd or 3rd degree. appear only on the four corners of the rhomb, so |E{2,2}|=|{e=uvE(Rn)du=2,dv=2}|=6. Edges E{2, 3} appear in pairs along the boundary except the corners where they are single. So we obtain |E{2,3}|=|{e=uvE(Zn)|du=2,dv=3}|=8(n1). Remaining are edges of type (3,3) given as |E{3,3}|=|{e=uvE(Zn)|du=3,dv=3}|=3n2+4n1 − 6 − 8(n − 1) = )n2 − 4n + 1.

Graph of rhombic benzenoid system having n rows of n hexagons.
Figure 4

Graph of rhombic benzenoid system having n rows of n hexagons.

So we obtain the following result.

Theorem 3

Consider the rhombic benzenoid system Rn, then its M- polynomial is: M(Rn;x,y)=6x2y2+8(n1)x2y3+(n(3n4)+1)x3y3.


The proof is similar to theorem 1.

Now we compute some degree-based topological indices from this M-polynomial.

Proportion 4

Consider the rhombic benzenoid system Rn, then

  1. mM2(Rn)=518+43n+19n(3n4).

  2. Rα(Rn)=9αn(3n4)+8n6α+9α+6×4α8×6α.

  3. Rα(Rn)=64α+8n86α+n(3n4)+19α.

  4. SSD(Rn)=103+523n+2n(3n4).

  5. H(Rn)=215+165n+13n(3n4).

  6. I(Rn)=2110+485n+32n(3n4).

  7. A(Rn)=29564+64n+72964n(3n4).


The proof of this proposition is similar to proposition 2.

The plot for the M-polynomial of R1.
Figure 6

The plot for the M-polynomial of R1.

4 Conclusions and Discussion

In this article we computed M-polynomials and closed forms of degree-based topological indices of zigzag and rhombic benzenoid system. These indices play important role in determining properties of compound under investigations. We gave graphs of M-polynomials against the number of hexagons n in each structure. These graphs in fact determine the dependency of above discussed topological indices relating to n.


Authors are highly grateful to the referees for their careful reading and technical suggestions to improve the quality of the article. This research is supported by Gyeongsang National University, Jinju 52828, Korea


  • [1]

    Rücker G., Rücker C., On topological indices, boiling points, and cycloalkanes. Journal of chemical information and computer sciences. 1999, 39(5),788-802. CrossrefGoogle Scholar

  • [2]

    Gutman I,, Trinajstić N., Graph theory and molecular orbitals. Total φ-electron energy of alternant hydrocarbons. Chemical Physics Letters. 197, 17(4), 535-8. Google Scholar

  • [3]

    Amić D., Bešlo D., Lučlć B., Nikolić S., Trinajstić N., The vertex-connectivity index revisited. Journal of chemical information and computer sciences. 1998, 38(5), 819-822. CrossrefGoogle Scholar

  • [4]

    Kier L.B., Hall L.H., Molecular connectivity in structure-activity analysis. Research Studies, 1986. Google Scholar

  • [5]

    Gutman I., Molecular graphs with minimal and maximal Randić indices. Croatica chemica acta, 2002, 75(2), 357-369. Google Scholar

  • [6]

    Gutman I., Degree-based topological indices. Croatica Chemica Acta, 2013, 86(4), 351-361. CrossrefWeb of ScienceGoogle Scholar

  • [7]

    Vukičević, D., On the edge degrees of trees. Glas. Mat. Ser. III, 2009, 44(64), 259–266. CrossrefGoogle Scholar

  • [8]

    Dobrynin A.A., Entringer R., Gutman I., Wiener index of trees: theory and applications. Acta Applicandae Mathematicae, 2001, 66(3), 211-249. CrossrefGoogle Scholar

  • [9]

    Munir M., Nazeer W., Rafique S., Kang S.M., M-polynomial and degree-based topological indices of polyhex nanotubes. Symmetry, 2016, 8(12), 149. Web of ScienceCrossrefGoogle Scholar

  • [10]

    Munir M., Nazeer W., Rafique S., Kang S.M., M-polynomial and degree-based topological indices of polyhex nanotubes. Symmetry, 2016, 8(12),149. Web of ScienceCrossrefGoogle Scholar

  • [11]

    Ajmal M., Kang S.M., Nazeer W., Munir M., Jung C.Y., Some Topological Invariants of the Möbius Ladders. Global Journal of Pure and Applied Mathematics, 2016, 12(6), 5317-5327.Google Scholar

  • [12]

    Munir M., Nazeer W., Rafique S., Nizami A., Kang S. M., Some Computational Aspects of Triangular Boron Nanotubes, Symmetry, 2017, 9, 6, . CrossrefWeb of ScienceGoogle Scholar

  • [13]

    West D.B., Introduction to Graph Theory Prentice Hall Upper Saddle River. NJ Google Scholar, 1996. Google Scholar

  • [14]

    Deutsch E., Klavzar S., M-Polynomial and Degree-Based Topological Indices. Iranian Journal of Mathematical Chemistry, 2015, 6(2), 93-102. Google Scholar

  • [15]

    Wu R., Deng H., The general connectivity indices of benzenoid systems and phenylenes. MATCH Commun. Math. Comput. Chem., 2010, 64, 459-470. Google Scholar

  • [16]

    Javaid M., Jung C. Y., M-Polynomials and Topological Indices of Silicate and Oxide Networks, International Journal of Pure and Applied Mathematics, 2017, 115(1), 129-152,  CrossrefGoogle Scholar

  • [17]

    Kwun Y.C., Munir M., Nazeer W., Rafique S., Kang S.M., M-Polynomials and topological indices of V-Phenylenic Nanotubes and Nanotori. Scientific reports, 201, 7(1), 8756. Web of ScienceGoogle Scholar

  • [18]

    Klavžar S., Gutman I., A comparison of the Schultz molecular topological index with the Wiener index. Journal of chemical information and computer sciences. 1996 Sep 24;36(5):1001- 1003. CrossrefGoogle Scholar

  • [19]

    Deng H., Huang G., Jiang X., A unified linear-programming modeling of some topological indices. Journal of Combinatorial Optimization, 2015, 30(3), 826-837. CrossrefWeb of ScienceGoogle Scholar

  • [20]

    Deng H., Yang J., Xia F., A general modeling of some vertex- degree based topological indices in benzenoid systems and phenylenes. Computers & Mathematics with Applications, 2011, 61(10), 3017-23. Web of ScienceCrossrefGoogle Scholar

  • [21]

    Gutman I., Das K.C., The first Zagreb index 30 years after. MATCH Commun. Math. Comput. Chem., 2004, 50, 83-92. Google Scholar

  • [22]

    Bollobás B., Erdös P., Graphs of extremal weights. Ars Combinatoria, 1998, 50, 225-233. Google Scholar

  • [23]

    Kier L.B., Hall L.H., Molecular connectivity in structure-activity analysis. Research Studies, 1986. Google Scholar

  • [24]

    Huang Y., Liu B., Gan L., Augmented Zagreb index of connected graphs. Match-Communications in Mathematical and Computer Chemistry, 2012, 67(2), 483-494. Google Scholar

  • [25]

    Furtula B., Graovac A., Vukičević D., Augmented Zagreb index, J. Math. Chem., 2010, 48, 370–380. Web of ScienceCrossrefGoogle Scholar

  • [26]

    Farahani M.R., Rajesh Kanna M.R., Jamil M.K., Imran M., Computing the M-Polynomial of Benzenoid Molecular Graphs. Science International (Lahore), 2016, 28, 3251-3255. Google Scholar

  • [27]

    Sardar M. S., Zafar S., Farahani M.R., THE GENERALIZED ZAGREB INDEX OF CAPRA-DESIGNED PLANAR BENZENOID SERIES $Ca_k(C_6)$, Open J. Math. Sci., 2017, 1 44–51. CrossrefGoogle Scholar

  • [28]

    Mutee ur Rehman H., Sardar R., Raza A., Computing Topological Indices of Hex Board and its Line Graph, Open J. Math. Sci., 2017, 1, 62 - 71. CrossrefGoogle Scholar

  • [29]

    Sardar, M. S., Pan X.X., Gao W., Farahani M. R., Computing Sanskruti Index of Titania Nanotubes, Open J. Math. Sci., 2017, 126 - 131 Google Scholar

  • [30]

    Farahani M.R., Gao W., Baig A.Q., Khalid W., Molecular description of copper (II) oxide. Macedonian Journal of Chemistry and Chemical Engineering, 2017, 36(1), 93-99. Web of ScienceGoogle Scholar

  • [31]

    Gao W., Wang Y., Wang W., Shi L., The first multiplication atombond connectivity index of molecular structures in drugs. Saudi Pharmaceutical Journal, 2017, 25(4), 548-555. CrossrefGoogle Scholar

  • [32]

    Gao W., Farahani M.R., Wang S., Husin M.N., On the edge- version atom-bond connectivity and geometric arithmetic indices of certain graph operations. Applied Mathematics and Computation, 2017, 308, 11-17. Web of ScienceCrossrefGoogle Scholar

About the article

Received: 2017-08-22

Accepted: 2017-12-04

Published Online: 2018-02-24

Conflicts of Interest: The authors declare no conflict of interest.

Author Contributions: All authors contribute equally in writing of this paper.

Citation Information: Open Chemistry, Volume 16, Issue 1, Pages 73–78, ISSN (Online) 2391-5420, DOI: https://doi.org/10.1515/chem-2018-0010.

Export Citation

© 2018 Ashaq Ali et al., published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

Lian Chen, Abid Mehboob, Haseeb Ahmad, Waqas Nazeer, Muhammad Hussain, and M. Reza Farahani
Discrete Dynamics in Nature and Society, 2019, Volume 2019, Page 1
Zafar Hussain, Mobeen Munir, Muhammad Bilal, Alam Ameer, Shazia Rafique, and Shin Min Kang
Open Chemistry, 2019, Volume 17, Number 1, Page 177
Lina Zheng, Yiqiao Wang, and Wei Gao
Open Chemistry, 2019, Volume 17, Number 1, Page 81
Zeeshan Saleem Mufti, Muhammad Faisal Nadeem, Wei Gao, and Zaheer Ahmad
Open Chemistry, 2018, Volume 16, Number 1, Page 1200
Young Chel Kwun, Manzoor Ahmad Zahid, Waqas Nazeer, Ashaq Ali, Maqbool Ahmad, and Shin Min Kang
Open Physics, 2018, Volume 16, Number 1, Page 734
Zafar Hussain, Mobeen Munir, Shazia Rafique, and Shin Min Kang
Symmetry, 2018, Volume 10, Number 10, Page 478
Young Chel Kwun, Ashaq Ali, Waqas Nazeer, Maqbool Ahmad Chaudhary, and Shin Min Kang
Journal of Chemistry, 2018, Volume 2018, Page 1
Waqas Nazeer, Adeel Farooq, Muhammad Younas, Mobeen Munir, and Shin Kang
Biomolecules, 2018, Volume 8, Number 3, Page 92
Emeric Deutsch and Sandi Klavžar
Journal of Applied Mathematics and Computing, 2018

Comments (0)

Please log in or register to comment.
Log in