Next Article in Journal
Stress–Strength Reliability Analysis for Different Distributions Using Progressive Type-II Censoring with Binomial Removal
Next Article in Special Issue
Eigenvalue −1 of the Vertex Quadrangulation of a 4-Regular Graph
Previous Article in Journal
Exotic Particle Dynamics Using Novel Hermitian Spin Matrices
Previous Article in Special Issue
Conversion of Unweighted Graphs to Weighted Graphs Satisfying Properties R and SR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Multivariable Zhang–Zhang Polynomial of Phenylenes

1
Faculty of Natural Sciences and Mathematics, University of Maribor, Koroška Cesta 160, 2000 Maribor, Slovenia
2
Institute of Mathematics, Physics and Mechanics, Jadranska 19, 1000 Ljubljana, Slovenia
Axioms 2023, 12(11), 1053; https://doi.org/10.3390/axioms12111053
Submission received: 6 October 2023 / Revised: 10 November 2023 / Accepted: 14 November 2023 / Published: 15 November 2023
(This article belongs to the Special Issue Spectral Graph Theory, Molecular Graph Theory and Their Applications)

Abstract

:
The Zhang–Zhang polynomial of a benzenoid system is a well-known counting polynomial that was introduced in 1996. It was designed to enumerate Clar covers, which are spanning subgraphs with only hexagons and edges as connected components. In 2018, the generalized Zhang–Zhang polynomial of two variables was defined such that it also takes into account 10-cycles of a benzenoid system. The aim of this paper is to introduce and study a new variation of the Zhang–Zhang polynomial for phenylenes, which are important molecular graphs composed of 6-membered and 4-membered rings. In our case, Clar covers can contain 4-cycles, 6-cycles, 8-cycles, and edges. Since this new polynomial has three variables, we call it the multivariable Zhang–Zhang (MZZ) polynomial. In the main part of the paper, some recursive formulas for calculating the MZZ polynomial from subgraphs of a given phenylene are developed and an algorithm for phenylene chains is deduced. Interestingly, computing the MZZ polynomial of a phenylene chain requires some techniques that are different to those used to calculate the (generalized) Zhang–Zhang polynomial of benzenoid chains. Finally, we prove a result that enables us to find the MZZ polynomial of a phenylene with branched hexagons.
MSC:
05C92; 05C31; 05C70; 05C30; 05C10

1. Introduction

Kekulé structures are often employed to provide an insight into the  π -electron structure of polycyclic conjugated molecules (note that a Kekulé structure of a polycyclic conjugated hydrocarbon is represented by a structural formula with double bonds between certain pairs of carbon atoms, such that each carbon atom is adjacent to exactly one double bond). When a molecule is modelled by the molecular graph, Kekulé structures are in one-to-one correspondence with the perfect matchings of the obtained graph [1,2]. Since Kekulé structures can be used to predict various properties and chemical behaviour of molecules [2,3], the investigation and enumeration of these structures is a classical topic in theoretical and mathematical chemistry.
A theory closely related to the concept of Kekulé structures is the Clar’s aromatic sextet theory [4]. In accordance with this theory, H. Zhang and F. Zhang introduced the concept of Clar covers [5]. Note that a Clar cover of a benzenoid system (i.e., a benzenoid graph that can be embedded into the regular hexagonal lattice) consists of some hexagons (which represent the aromatic sextets) and of double bonds such that all vertices of the graph are covered. The maximum possible number of hexagons among all Clar covers in a given benzenoid graph is known as the Clar number. By Clar’s theory, the most important Kekulé structures are those in which the number of aromatic sextets equals the Clar number. Some well-investigated concepts that are closely related to Clar covers are resonant sets, alternating sets, and the Fries number (for some research on these topics see [6,7,8,9]).
In [5], the authors also introduced the so called Clar covering polynomial that was later named the Zhang–Zhang polynomial. For a benzenoid system G, this polynomial is defined as
Z Z ( G ; x ) = p 0 z ( G ; p ) x p ,
where  z ( G ; p )  is the number of Clar covers of G with exactly p hexagons. It is interesting that the Clar number, the number of Kekulé structures, and the first Herndon number can be easily calculated by using this polynomial. For some recent papers on the Zhang–Zhang polynomial, see [10,11,12].
Another topic closely related to the above-mentioned concepts is the study of resonance graphs. These graphs are used to model interactions among different Kekulé structures; see [13] for a survey on this topic. It is worth mentioning that the relation between the Zhang–Zhang polynomial of a molecular graph and the cube polynomial of its resonance graph was considered for different families of graphs [12,14]. Some recent research on resonance graphs can be found, for example, in [15].
In [16], the generalized Zhang–Zhang polynomial (GZZ polynomial) was introduced in order to increase the sensibility of the standard Zhang–Zhang polynomial. This polynomial contains two variables and counts so-called generalized Clar covers, which can contain edges, hexagons, and also 10-cycles (representing two adjacent hexagons). Therefore, the GZZ polynomial gives explicit information on  π -electron cyclic conjugation within 10-membered rings. It was shown in the same paper that for some families of molecular graphs, the GZZ polynomial of a given graph is equal to the generalized cube polynomial of the corresponding resonance graph. Later, several recursive formulas for calculating the GZZ polynomial of benzenoid systems were derived and numerical results were obtained in order to test the relation between the GZZ polynomial and three energy-based quantities [17]. Furthermore, it was shown that the molecular vibrational energies can be related to Clar-structure-based parameters by using the GZZ polynomial [18].
The aim of this paper is to introduce another variant of the Zhang–Zhang polynomial for phenylenes. Note that phenylenes are conjugated systems composed of 6-membered rings and 4-membered rings and, therefore, the  π -electron properties of these molecules are of great interest for theoretical chemistry. For some results related to the Kekulé structures of phenylenes, see [19,20,21], and more information on the chemistry of phenylenes can be found in the references stated in [20]. In order to take into account 4-membered rings, 6-membered rings, and also combinations of these rings, the corresponding Clar covers should not contain only 6-cycles (hexagons) and edges, but also 4-cycles (quadrilaterals) and 8-cycles. Hence, our new polynomial has three variables; therefore, we call it the multivariable Zhang–Zhang polynomial (in short, the MZZ polynomial). It is clear that this polynomial provides some information on the distribution of  π -electrons in phenylenes. Moreover, as we already mentioned, such approach could have interesting chemical applications [17,18].
In the next section, we state some basic definitions and formally introduce the MZZ polynomial. Next, in Section 3, we prove several recursive formulas that enable use to calculate the MZZ polynomial of a phenylene by using the MZZ polynomials of its subgraphs. An efficient algorithm for computing the MZZ polynomial of any phenylene chain is then deduced in Section 4. Finally, in Section 5 we discuss how to calculate the MZZ polynomial of a phenylene with branched hexagons.

2. Preliminaries

In this paper, we consider only finite and simple graphs. For a graph G, we denote by  V ( G )  the set of vertices and by  E ( G )  the set of edges of G. The degree of a vertex  v V ( G )  is the number of vertices adjacent to v. If  k 3 , then a cycle of length k in G, denoted as  C k , will be called a k-cycle of G.
Let G be a plane graph. Two distinct faces of G are adjacent if they share a common edge. We denote the set of edges lying on some face f of G by  E ( f ) . Also, the subgraph induced by the edges in  E ( f )  is the boundary of f. An inner face of G whose boundary is a 6-cycle  C 6  is called a hexagon of G. Similarly, an inner face of G whose boundary is a 4-cycle  C 4  is called a quadrilateral of G. Furthermore, the vertices of G that belong to the outer face are known as boundary vertices and all the other vertices of G are interior vertices. Similarly, the edges lying on the outer face will be called boundary edges. An outerplane graph is a plane graph in which all vertices are boundary vertices.
A benzenoid graph is a bipartite 2-connected plane graph in which all interior vertices have degree 3, all boundary vertices have degree 2 or 3, and all inner faces are hexagons. It is worth mentioning that these graphs are sometimes referred to as fusenes [22]. Note also that by our definition, a benzenoid graph is not always a subgraph of the regular hexagonal lattice. Moreover, a benzenoid graph which does not contain any interior vertices (i.e., an outerplane benzenoid graph) is called a catacondensed benzenoid graph.
Let B be a catacondensed benzenoid graph. If we add a quadrilateral between any two adjacent hexagons of B, then the obtained graph is called a phenylene. In Figure 1, we can see a catacondensed benzenoid graph B and the corresponding phenylene G.
To deduce our main results of the paper, we need to consider a wider family of graphs and not just phenylenes. Therefore, throughout the paper a generalized phenylene will be an outerplane bipartite graph in which no two distinct quadrilaterals are adjacent. It is easy to check that in such a graph every 4-cycle is the boundary of some quadrilateral, every 6-cycle is the boundary of some hexagon, and in the interior of every 8-cycle there is either one inner face or one quadrilateral and one hexagon (with the common edge).
Let h be a hexagon of a generalized phenylene G such that h contains exactly two vertices of degree 2, which will be denoted as x and y. We say that h is angular if  x y  is an edge of G, and on the other hand, h is linear if x and y are not adjacent in G. Moreover, if G is a phenylene, then a hexagon of G is called branched if it is adjacent to exactly three quadrilaterals of G. It is easy to see that phenylene G from Figure 1 contains two angular hexagons, one linear hexagon, and no branched hexagons. On the other hand, a phenylene from Figure 2 contains one angular hexagon, two linear hexagons, and one branched hexagon.
Let G be a phenylene. If every hexagon of G is adjacent to at most two quadrilaterals, we say that G is a phenylene chain. In addition, a linear phenylene chain is a phenylene chain that does not contain any angular hexagons. Obviously, phenylene G from Figure 1 is a phenylene chain, but G is not a linear phenylene chain.
A subset of edges M of a graph G is called a perfect matching of G if every vertex of G is an end vertex of exactly one edge from M. It is well known that in chemistry, perfect matchings are usually referred to as Kekulé structures.
A spanning subgraph C of a generalized phenylene G is a (4,6,8)-Clar cover of G if every connected component of C is a 4-cycle  C 4  (a quadrilateral), a 6-cycle  C 6  (a hexagon), an 8-cycle  C 8 , or an edge  K 2 . Figure 2 shows a phenylene with a (4,6,8)-Clar cover composed of one 8-cycle, two hexagons, two quadrilaterals, and seven edges (the bold edges in the figure represent connected components of this (4,6,8)-Clar cover).
Finally, we introduce the multivariable Zhang–Zhang polynomial (or shortly MZZ polynomial M Z Z ( G ; x , y , z )  of a generalized phenylene G by setting
M Z Z ( G ; x , y , z ) = p , q , r 0 m z ( G ; p , q , r ) x p y q z r ,
where  m z ( G ; p , q , r )  represents the number of (4,6,8)-Clar covers of G that contain exactly p quadrilaterals, q hexagons, and r 8-cycles. Moreover, if G has no vertices, we consider the set ∅ as the unique (4,6,8)-Clar cover of G, so in this case we define  M Z Z ( G ; x , y , z ) = 1 . Furthermore, the MZZ polynomial of G will be often denoted simply as  M Z Z ( G ) , i.e.,  M Z Z ( G ) : = M Z Z ( G ; x , y , z ) .
Some simple observations about the  M Z Z  polynomial of a generalized phenylene G can now be stated:
  ( i )
The number of perfect matchings of G is exactly  m z ( G ; 0 , 0 , 0 ) , which is equal to  M Z Z ( G ; 0 , 0 , 0 ) . Therefore, the number of perfect matchings can be obtained by taking  x = y = z = 0  into the MZZ polynomial.
( i i )
If G does not have any perfect matching, then the set of all (4,6,8)-Clar covers of G is empty, so  M Z Z ( G ) = 0 .
( i i i )
If G has only one perfect matching, then G has exactly one (4,6,8)-Clar cover which contains only edges from the unique perfect matching, so  M Z Z ( G ) = 1 .
( i v )
It can be observed that if G is a phenylene, then the set of all hexagons of G is always a (4,6,8)-Clar cover of G (which contains only hexagons as connected components). Therefore, the MZZ polynomial of a phenylene with  h ( G )  hexagons always contains the term  y h ( G ) .
Suppose that G is a generalized phenylene and let  CC ( G )  be the set of all (4,6,8)-Clar covers of G. In addition, for any  C CC ( G )  let  q ( C )  be the number of quadrilaterals of C h ( C )  the number of hexagons of C, and  o ( C )  as the number of 8-cycles of C. It is not difficult to see that we can express the  M Z Z  polynomial in the following way.
Proposition 1.
If G is a generalized phenylene, then
M Z Z ( G ) = C CC ( G ) x q ( C ) y h ( C ) z o ( C ) .
Note that if  CC ( G ) =  in the above proposition, then the index set of the corresponding sum is empty, so  M Z Z ( G ) = 0 .
We finish this section with one additional notation. If G is a generalized phenylene, C a (4,6,8)-Clar cover of G, and H a quadrilateral, a hexagon, an 8-cycle, or an edge of G, then by writing  H C  we mean that H is a connected component of C.

3. Computing the MZZ Polynomial from Smaller Graphs

In this section, we provide several results that enable us to calculate the MZZ polynomial of a generalized phenylene by using the MZZ polynomials of subgraphs of the original graph. We should mention that the results stated in this section are analogous to those from [5,17], but several additional insights are needed in our case.
The following notation will be used for any generalized phenylene G:
  • If  e = u v  is an edge of G, then  G e = G u v  denotes the graph obtain from G by deleting edge e;
  • If  e = u v  is an edge of G, then  G u v  denotes the graph obtain from G by deleting vertices u and v, together with all the edges incident to u or v;
  • If f is an inner face of G, then  G f  denotes the graph obtain from G by deleting all the vertices of f, together with all the edges incident to the vertices of f;
  • If f and  f  are two adjacent inner faces of G, then  G f f  denotes the graph obtain from G by deleting all the vertices of f and all the vertices of  f , together with all the edges incident to these vertices;
  • If f is a hexagon of G and  f  a quadrilateral of G such that f and  f  are adjacent, then  f f  denotes the unique 8-cycle of G induced by f and  f .
In the first theorem we show that the MZZ polynomial of a generalized phenylene G is equal to the product of MZZ polynomials of connected components of G.
Theorem 1.
Let G be a generalized phenylene with connected components  G 1 , G 2 , , G p . Then it holds
M Z Z ( G ) = j = 1 p M Z Z ( G j ) .
Proof. 
If  CC ( G ) = , the result is obvious. Otherwise, we can write the set  CC ( G j )  as  CC ( G j ) = { C | G j : C CC ( G ) } , where  C | G j  denotes the restriction of C to  G j  for any  j { 1 , , p } . Therefore, one can obtain
M Z Z ( G ) = C CC ( G ) x q ( C ) y h ( C ) z o ( C ) = C CC ( G ) j = 1 p x q ( C | G j ) y h ( C | G j ) z o ( C | G j ) = C 1 CC ( G 1 ) C 2 CC ( G 2 ) C p CC ( G p ) j = 1 p x q ( C j ) y h ( C j ) z o ( C j ) = j = 1 p C j CC ( G j ) x q ( C j ) y h ( C j ) z o ( C j ) = j = 1 p M Z Z ( G j ) ,
which is the desired result.    □
The next theorem can be applied when one has a generalized phenylene composed of several generalized phenylenes joined by edges.
Theorem 2.
If G is a generalized phenylene and  e = u v  is an edge of G not belonging to any quadrilateral, hexagon or 8-cycle of G, then it holds
M Z Z ( G ) = M Z Z ( G u v ) + M Z Z ( G e ) .
Proof. 
If  CC ( G ) = , the result obviously follows. Therefore, we can suppose  CC ( G ) . Let  CC 1 = { C CC ( G ) : e C }  and  CC 2 = { C CC ( G ) : e C } . Therefore, we obtain
M Z Z ( G ) = C CC ( G ) x q ( C ) y h ( C ) z o ( C ) = C CC 1 x q ( C ) y h ( C ) z o ( C ) + C CC 2 x q ( C ) y h ( C ) z o ( C ) = C CC 1 x q ( C e ) y h ( C e ) z o ( C e ) + C CC 2 x q ( C ) y h ( C ) z o ( C ) .
We can easily see that
C CC 1 x q ( C e ) y h ( C e ) z o ( C e ) = C CC ( G u v ) x q ( C ) y h ( C ) z o ( C ) .
On the other hand, we know that e does not belong to any quadrilateral, hexagon or 8-cycle of G, which implies
C CC 2 x q ( C ) y h ( C ) z o ( C ) = C CC ( G e ) x q ( C ) y h ( C ) z o ( C ) .
Therefore, we obtain
M Z Z ( G ) = C CC ( G u v ) x q ( C ) y h ( C ) z o ( C ) + C CC ( G e ) x q ( C ) y h ( C ) z o ( C ) = M Z Z ( G u v ) + M Z Z ( G e ) ,
which completes the proof.    □
The following theorem will be essentially used in the proofs of some other results.
Theorem 3.
Suppose that G is a generalized phenylene. Let  e = u v  be a boundary edge of G such that e also lies on some hexagon f of G (see Figure 3). Also, let  p 5  denotes the number of quadrilaterals adjacent to f. If  p = 0 , then
M Z Z ( G ) = y M Z Z ( G f ) + M Z Z ( G u v ) + M Z Z ( G e ) .
If  p 0 , let  f 1 , , f p  be the quadrilaterals adjacent to f. Then
M Z Z ( G ) = y M Z Z ( G f ) + M Z Z ( G u v ) + z j = 1 p M Z Z ( G f f j ) + M Z Z ( G e ) .
Proof. 
If  CC ( G ) = , the result obviously follows. Therefore, assume  CC ( G ) . Note that the details of the proof will be stated only for  p 0 , since the case  p = 0  is easier. Firstly, we introduce the following sets:
CC 1 = { C CC ( G ) : f C } , CC 2 = { C CC ( G ) : e C } , CC 3 j = { C CC ( G ) : f f j C } for j { 1 , , p } , CC 4 = CC ( G ) CC 1 CC 2 j = 1 p CC 3 j .
Since the sets stated above are pairwise disjoint, we obtain
M Z Z ( G ) = C CC ( G ) x q ( C ) y h ( C ) z o ( C ) = j = 1 2 C CC j x q ( C ) y h ( C ) z o ( C ) + j = 1 p C CC 3 j x q ( C ) y h ( C ) z o ( C ) + C CC 4 x q ( C ) y h ( C ) z o ( C ) = C CC 1 x q ( C f ) y h ( C f ) + 1 z o ( C f ) + C CC 2 x q ( C e ) y h ( C e ) z o ( C e ) + j = 1 p C CC 3 j x q ( C f f j ) y h ( C f f j ) z o ( C f f j ) + 1 + C CC 4 x q ( C ) y h ( C ) z o ( C ) = y C CC ( G f ) x q ( C ) y h ( C ) z o ( C ) + C CC ( G u v ) x q ( C ) y h ( C ) z o ( C ) + z j = 1 p C CC ( G f f j ) x q ( C ) y h ( C ) z o ( C ) + C CC ( G e ) x q ( C ) y h ( C ) z o ( C ) = y M Z Z ( G f ) + M Z Z ( G u v ) + z j = 1 p M Z Z ( G f f j ) + M Z Z ( G e ) ,
which means we are done.    □
A statement similar to that in Theorem 3 can be obtained if f is a quadrilateral.
Theorem 4.
Suppose that G is a generalized phenylene. Let  e = u v  be a boundary edge of G such that e also lies on some quadrilateral f of G. Also, let  p 3  denotes the number of hexagons adjacent to f. If  p = 0 , then
M Z Z ( G ) = x M Z Z ( G f ) + M Z Z ( G u v ) + M Z Z ( G e ) .
If  p 0 , let  f 1 , , f p  be the hexagons adjacent to f. Then
M Z Z ( G ) = x M Z Z ( G f ) + M Z Z ( G u v ) + z j = 1 p M Z Z ( G f f j ) + M Z Z ( G e ) .
Proof. 
The proof is very similar to the proof of Theorem 3. Therefore, we skip the details.    □
The following proposition will be needed in several parts of the paper.
Proposition 2.
Suppose G is a generalized phenylene and  e = u v  is an edge of G. If e does not belong to any perfect matching of G, then
M Z Z ( G ) = M Z Z ( G e ) .
On the other hand, if e belongs to all perfect matchings of G (or if G does not have a perfect matching), then
M Z Z ( G ) = M Z Z ( G u v ) .
Proof. 
It is easy to see that both formulas hold if G does not have a perfect matching (in such a case  M Z Z ( G ) = M Z Z ( G e ) = M Z Z ( G u v ) = 0 ). Therefore, suppose that G has a perfect matching.
If e does not belong to any perfect matching or if it belongs to all perfect matchings, then obviously e can not be an edge of a 4-cycle, a 6-cycle, or an 8-cycle that is contained in some (4,6,8)-Clar cover of G. If e does not belong to any perfect matching, then the set  CC ( G )  is the same as the set  CC ( G e ) . On the other hand, if e belongs to all perfect matchings of G, we can find a bijection between the set  CC ( G )  and the set  CC ( G u v ) , which preserves the number of 4-cycles, 6-cycles, and 8-cycles in corresponding (4,6,8)-Clar covers. The result now follows directly from the definition of the MZZ polynomial.    □
To state the next theorem, we need the following assumption and some additional notation.
Assumption 1.
Let  G 1  be a generalized phenylene with a perfect matching and let  e 1 = u 1 v 1  be a boundary edge of  G 1  with both end vertices of degree 2 such that  e 1  also lies on some hexagon  f 1  of  G 1 . Similarly, let  G 2  be a generalized phenylene with a perfect matching and let  e 2 = u 2 v 2  be a boundary edge of  G 2  such that  e 2  also lies on some quadrilateral  f 2  of  G 2 .
By identifying edges  e 1  and  e 2  of graphs  G 1  and  G 2  satisfying Assumption 1, we obtain another generalized phenylene which will be denoted as  G 1 · G 2 , see Figure 4. Moreover, by  e = u v  we denote the new edge obtained by the identification of edges  e 1  and  e 2  (sometimes we also write e instead of  e 1  or  e 2 ). In addition, for  i { 1 , 2 } , denote by  G i  the graph  G i u i v i  and by  G i  the graph  G i f i .
In the following theorem, we prove a formula for computing the MZZ polynomial of  G 1 · G 2 .
Theorem 5.
If  G 1  and  G 2  are two generalized phenylenes satisfying Assumption 1, then the MZZ polynomial of  G 1 · G 2  can be computed as
M Z Z ( G 1 · G 2 ) = M Z Z ( G 1 ) M Z Z ( G 2 ) + M Z Z ( G 1 ) M Z Z ( G 2 ) + z M Z Z ( G 1 ) M Z Z ( G 2 ) M Z Z ( G 1 ) M Z Z ( G 2 ) .
Proof. 
Let  s 1 , , s p p 3 , be the quadrilaterals of  G 1  adjacent to  f 1 . We can assume that  p 1  (otherwise the proof is immediate). Let all the notation be as in Figure 4. Firstly, we apply Theorem 3 to edge  u u  of  G 1 · G 2 . Therefore,
M Z Z ( G 1 · G 2 ) = y M Z Z ( G 1 · G 2 f 1 ) + M Z Z ( G 1 · G 2 u u ) + z j = 1 p M Z Z ( G 1 · G 2 f 1 s j ) + z M Z Z ( G 1 · G 2 f 1 f 2 ) + M Z Z ( G 1 · G 2 u u ) .
By Theorem 1, we have
M Z Z ( G 1 · G 2 f 1 ) = M Z Z ( G 1 f 1 ) M Z Z ( G 2 ) ,
M Z Z ( G 1 · G 2 f 1 f 2 ) = M Z Z ( G 1 ) M Z Z ( G 2 ) ,
and
M Z Z ( G 1 · G 2 f 1 s j ) = M Z Z ( G 1 f 1 s j ) M Z Z ( G 2 )
for  j { 1 , , p } .
It is easy to see that the edge  v v  does not belong to any perfect matching of  G 1 · G 2 u u . This is obvious if  G 1 · G 2 u u  does not have a perfect matching. Otherwise, suppose that  v v  belongs to some perfect matching of  G 1 · G 2 u u . Then it follows that the graph  G 2 v  has a perfect matching, which is a contradiction because  G 2  has a perfect matching. Consequently, by Proposition 2 and Theorem 1 we obtain
M Z Z ( G 1 · G 2 u u ) = M Z Z ( G 1 ) M Z Z ( G 2 ) .
In a similar way we can see that if  G 1 · G 2 u u  has a perfect matching, then the edge  v v  belongs to any perfect matching of  G 1 · G 2 u u  (otherwise  G 2 u  has a perfect matching) and therefore, by Proposition 2 and Theorem 1 we obtain
M Z Z ( G 1 · G 2 u u ) = M Z Z ( G 1 u v ) M Z Z ( G 2 ) .
Hence, it follows
M Z Z ( G 1 · G 2 ) = M Z Z ( G 1 ) M Z Z ( G 2 ) + z M Z Z ( G 1 ) M Z Z ( G 2 ) + M Z Z ( G 2 ) ( y M Z Z ( G 1 f 1 ) + M Z Z ( G 1 u v ) + z j = 1 p M Z Z ( G 1 f 1 s j ) ) .
Next, we apply Theorem 3 to edge  e = u v  (or  e 1 = u 1 v 1 ) of  G 1 . Consequently,
M Z Z ( G 1 ) = y M Z Z ( G 1 f 1 ) + M Z Z ( G 1 u v ) + z j = 1 p M Z Z ( G 1 f 1 s j ) + M Z Z ( G 1 e ) .
Obviously, the edges  u u  and  v v  belong to all perfect matchings of  G 1 e  (if this graph has a perfect matching), so Proposition 2 implies
M Z Z ( G 1 e ) = M Z Z ( G 1 u v u v ) = M Z Z ( G 1 u v ) .
Therefore,
M Z Z ( G 1 ) M Z Z ( G 1 ) = y M Z Z ( G 1 f 1 ) + M Z Z ( G 1 u v ) + z j = 1 p M Z Z ( G 1 f 1 s j )
and we finally obtain
M Z Z ( G 1 · G 2 ) = M Z Z ( G 1 ) M Z Z ( G 2 ) + z M Z Z ( G 1 ) M Z Z ( G 2 ) + M Z Z ( G 1 ) M Z Z ( G 2 ) M Z Z ( G 1 ) M Z Z ( G 2 ) ,
which completes the proof.    □
Remark 1.
It can happen that  G 1 · G 2  does not have a perfect matching although each of  G 1  and  G 2  has a perfect matching. In this case, edge  e 1  should not belong to any perfect matching of  G 1  and  e 2  should not belong to any perfect matching of  G 2 . Then we have  M Z Z ( G 1 ) = M Z Z ( G 2 ) = M Z Z ( G 1 ) = M Z Z ( G 2 ) = 0 , so both sides of the formula stated in Theorem 5 are equal to 0.
The basic compound of a phenylene is the graph obtained from a quadrilateral and a hexagon by identifying two edges, see Figure 5a. We now apply the previous result in the case where  G 2  is the basic compound of a phenylene. Therefore, we obtain the MZZ polynomial of  G 1 · P , see Figure 5b.
Corollary 1.
Let P be the basic compound of a phenylene and let e be an edge of P shown in Figure 5a. If  G 1  is a generalized phenylene satisfying Assumption 1, then
M Z Z ( G 1 · P ) = ( y + 2 ) M Z Z ( G 1 ) + z M Z Z ( G 1 ) + ( x + z + 1 ) M Z Z ( G 1 ) .
Proof. 
Obviously,  M Z Z ( P ) = x + y + z + 3 M Z Z ( P ) = y + 2 , and  M Z Z ( P ) = 1 . Hence, by Theorem 5 we obtain
M Z Z ( G 1 · P ) = M Z Z ( G 1 ) M Z Z ( P ) + M Z Z ( G 1 ) M Z Z ( P ) + z M Z Z ( G 1 ) M Z Z ( P ) M Z Z ( G 1 ) M Z Z ( P ) = M Z Z ( G 1 ) ( y + 2 ) + M Z Z ( G 1 ) ( x + y + z + 3 ) + z M Z Z ( G 1 ) M Z Z ( G 1 ) ( y + 2 ) = ( y + 2 ) M Z Z ( G 1 ) + z M Z Z ( G 1 ) + ( x + z + 1 ) M Z Z ( G 1 ) ,
which finishes the proof.    □
In a similar way, we can also prove the next corollary.
Corollary 2.
Let Q be the quadrilateral and let e be an edge of Q. If  G 1  is a generalized phenylene satisfying Assumption 1, then
M Z Z ( G 1 · Q ) = M Z Z ( G 1 ) + z M Z Z ( G 1 ) + ( x + 1 ) M Z Z ( G 1 ) .

4. Phenylene Chains

In this section, we provide some techniques for computing the MZZ polynomial of phenylene chains. Interestingly, it turns our that this task requires a different approach than the one used to calculate the (generalized) Zhang–Zhang polynomial of benzenoid chains [5,17]. In fact, even for linear phenylene chains, the computation is not trivial.
Let  P h 0  be the graph with no vertices and let  P h 1  be the hexagon  h 1 . Moreover, for  n > 1  we denote by  P h n  a phenylene chain with exactly n hexagons obtained by adding the basic compound of a phenylene to the hexagon  h n 1  of  P h n 1  (to the edge with both end vertices of degree 2 in  P h n 1 ). The new hexagon is denoted by  h n . As a special case, denote by  P L n  the linear phenylene chain with n hexagons, where  n 0 . The linear phenylene chain  P L 6  is depicted in Figure 6.
Computing the (generalized) Zhang–Zhang polynomial of a linear benzenoid chain is straightforward; see [5,17]. On the other hand, obtaining the MZZ polynomial of a linear phenylene chain is more complicated. However, this polynomial can be computed by the recurrence relation described in the following proposition.
Proposition 3.
If  P L n  is the linear phenylene chain with n hexagons, then  M Z Z ( P L 0 ) = 1 M Z Z ( P L 1 ) = y + 2 , and for any  n 1  it holds
M Z Z ( P L n + 1 ) = ( y + 2 ) M Z Z ( P L n ) + ( x + 2 z + 1 ) M Z Z ( P L n 1 ) .
Proof. 
We obtain by definition that  M Z Z ( P L 0 ) = 1  and  M Z Z ( P L 1 ) = y + 2 . If  n 1 , then the graph  P L n + 1  can be obtained as  P L n · P , where P is the basic compound of a phenylene, see Figure 5. Obviously,  P L n = P L n 1  and therefore,  M Z Z ( P L n ) = M Z Z ( P L n 1 ) . Moreover, by Proposition 2 we have  M Z Z ( P L n ) = M Z Z ( P L n 1 ) .  Consequently, by Corollary 1 we obtain
M Z Z ( P L n + 1 ) = ( y + 2 ) M Z Z ( P L n ) + z M Z Z ( P L n ) + ( x + z + 1 ) M Z Z ( P L n ) = ( y + 2 ) M Z Z ( P L n ) + z M Z Z ( P L n 1 ) + ( x + z + 1 ) M Z Z ( P L n 1 ) = ( y + 2 ) M Z Z ( P L n ) + ( x + 2 z + 1 ) M Z Z ( P L n 1 ) ,
which is the desired result.    □
By using the above result, we immediately get  M Z Z ( P L 2 ) = y 2 + x + 4 y + 2 z + 5  and  M Z Z ( P L 3 ) = y 3 + 2 x y + 6 y 2 + 4 y z + 4 x + 14 y + 8 z + 12 .
As a simple consequence of the above proposition, we can calculate the number of perfect matchings (Kekulé structures) of the linear phenylene chain  P L n . This number will be denoted as  K n . Since  K n = M Z Z ( P L n ; 0 , 0 , 0 ) , it follows by Proposition 3 that  K 0 = 1 K 1 = 2 , and for any  n 1  it holds
K n + 1 = 2 K n + K n 1 .
By solving this recurrence relation, we obtain that the number of perfect matchings of  P L n  is equal to
K n = 2 + 2 4 1 + 2 n + 2 2 4 1 2 n .
However, we should point out that the number  K n  was already calculated in [19] and that this sequence is closely related to the sequence A000129 (the Pell numbers) in the On-line encyclopedia of integer sequences. Moreover, different methods for computing the number of perfect matching of various molecular graphs are presented in  [1,2].
Next, we describe an algorithm that calculates the MZZ polynomial of any phenylene chain. Some auxiliary definitions and results are needed for this purpose.
Let  P h n  be a phenylene chain with n hexagons, where  n 2 . Obviously, hexagon  h n  contains exactly four vertices of degree 2. By deleting these four vertices from  P h n , we obtain the corresponding reduced phenylene chain  P h n * . Additionally, let  P h 1 *  be an edge  K 2 . In Figure 7, we can see a phenylene chain  P h 4  and the corresponding reduced phenylene chain  P h 4 * .
First, we state the following proposition, which generalizes Proposition 3.
Proposition 4.
Let  P h n + 1  be a phenylene chain with  n + 1  hexagons, where  n 2 . If the hexagon  h n  is linear, then it holds
M Z Z ( P h n + 1 ) = ( y + 2 ) M Z Z ( P h n ) + ( x + 2 z + 1 ) M Z Z ( P h n 1 ) .
On the other hand, if the hexagon  h n  is angular, then it holds
M Z Z ( P h n + 1 ) = ( y + 2 ) M Z Z ( P h n ) + z M Z Z ( P h n 1 ) + ( x + z + 1 ) M Z Z ( P h n * ) .
Proof. 
As in the proof of Proposition 3,  P h n + 1  can be obtained as  P h n · P , where P is the basic compound of a phenylene, see Figure 5. Again, it holds  P h n = P h n 1 . In the first case, by Proposition 2 we have  M Z Z ( P h n ) = M Z Z ( P h n 1 ) . On the other hand, if  h n  is angular, we obtain  M Z Z ( P h n ) = M Z Z ( P h n * ) . The result now follows by Corollary 1.    □
A similar proposition can be stated also for reduced phenylene chains.
Proposition 5.
Let  P h n + 1  be a phenylene chain with  n + 1  hexagons, where  n 2 , and let  P h n + 1 *  be the corresponding reduced phenylene chain. If the hexagon  h n  is linear, then it holds
M Z Z ( P h n + 1 * ) = M Z Z ( P h n ) + ( x + z + 1 ) M Z Z ( P h n 1 ) .
On the other hand, if the hexagon  h n  is angular, then it holds
M Z Z ( P h n + 1 * ) = M Z Z ( P h n ) + z M Z Z ( P h n 1 ) + ( x + 1 ) M Z Z ( P h n * ) .
Proof. 
Obviously,  P h n + 1 *  can be obtained as  P h n · Q , where Q is the quadrilateral. The proof is now similar to the proof of Proposition 4, but we use Corollary 2.    □
To a phenylene chain  P h n  with n hexagons we assign the vector of length n, denoted as  ( s 1 , s 2 , , s n ) , such that  s 1 = s n = 0  and if  n > 2 , then for every  i { 2 , , n 1 }  we define  s i { 1 , 2 }  in the following way:  s i = 1  if the hexagon  h i  is linear; and  s i = 2  if the hexagon  h i  is angular. For example, the phenylene chain from Figure 1 has the corresponding vector  ( 0 , 1 , 2 , 2 , 0 ) .
By using this notation, we can now present Algorithm 1. In the algorithm, we use the above propositions and the following initial values:  M Z Z ( P h 1 ) = y + 2 M Z Z ( P h 1 * ) = 1 M Z Z ( P h 2 ) = y 2 + x + 4 y + 2 z + 5 , and  M Z Z ( P h 2 * ) = x + y + z + 3 .
By using our implementation of Algorithm 1 in SageMath, for phenylene G from Figure 1 we immediately obtain
M Z Z ( G ) = y 5 + x 3 y + 2 x 2 y 2 + 4 x y 3 + 10 y 4 + 2 x 2 y z + 4 x y 2 z + 8 y 3 z + x y z 2 + 2 y 2 z 2 + 3 x 3 + 14 x 2 y + 28 x y 2 + 44 y 3 + 8 x 2 z + 32 x y z + 52 y 2 z + 7 x z 2 + 21 y z 2 + 2 z 3 + 23 x 2 + 73 x y + 106 y 2 + 56 x z + 126 y z + 39 z 2 + 69 x + 140 y + 112 z + 81 .
We can also state the following theorem.
Theorem 6.
Algorithm 1 correctly computes the MZZ polynomial of a phenylene chain  P h n  with n hexagons and the corresponding reduced phenylene chain  P h n * . Moreover, it can be implemented in  O ( n )  time, in the model where addition and multiplication of polynomials of three variables can be performed in constant time.
Proof. 
By Propositions 4 and 5, the algorithm correctly computes the MZZ polynomial of  P h n  and  P h n * . Moreover, the algorithm contains one “for” loop with  n 2  steps and all the remaining commands can be calculated in constant time. Hence, it follows that Algorithm 1 can be implemented in linear time with respect to the number of hexagons.    □
Algorithm 1: The MZZ polynomial of a (reduced) phenylene chain
Axioms 12 01053 i001
As we will see in the next section, it turns out useful to consider also chains that start with a quadrilateral. In this paper, such chains will be called modified phenylene chains. Let  M P h 0  be an edge  K 2  and let  M P h 1  be the basic compound of a phenylene with hexagon  h 1 . Moreover, for  n > 1 , we denote by  M P h n  a modified phenylene chain with exactly n hexagons obtained by adding the basic compound of a phenylene to hexagon  h n 1  of  M P h n 1  (to the edge with both end vertices of degree 2 in  M P h n 1 ). The new hexagon is denoted by  h n .
Similar to before, the reduced modified phenylene chain can be also defined (by deleting the four vertices of degree 2 in hexagon  h n ). If it is obtained from  M P h n  with  n 1 , we will denote it as  M P h n * . Figure 8 shows a modified phenylene chain  M P h 4  and the corresponding reduced modified phenylene chain  M P h 4 * . Note that any modified phenylene chain is also a reduced phenylene chain, but on the other hand, reduced modified phenylene chains form another class of graphs.
The next proposition can be proved in the same way as Propositions 4 and 5. We notice that the same recurrence relations apply for modified phenylene chains (with different initial conditions).
Proposition 6.
Let  M P h n + 1  be a modified phenylene chain with  n + 1  hexagons, where  n 1 , and let  M P h n + 1 *  be the reduced modified phenylene chain. If the hexagon  h n  is linear, then it holds
M Z Z ( M P h n + 1 ) = ( y + 2 ) M Z Z ( M P h n ) + ( x + 2 z + 1 ) M Z Z ( M P h n 1 ) , M Z Z ( M P h n + 1 * ) = M Z Z ( M P h n ) + ( x + z + 1 ) M Z Z ( M P h n 1 ) .
On the other hand, if the hexagon  h n  is angular, then it holds
M Z Z ( M P h n + 1 ) = ( y + 2 ) M Z Z ( M P h n ) + z M Z Z ( M P h n 1 ) + ( x + z + 1 ) M Z Z ( M P h n * ) , M Z Z ( M P h n + 1 * ) = M Z Z ( M P h n ) + z M Z Z ( M P h n 1 ) + ( x + 1 ) M Z Z ( M P h n * ) .
Similarly as before, to a modified phenylene chain  M P h n  with n hexagons we assign the vector of length n ( s 1 , s 2 , , s n ) , such that  s n = 0  and if  n > 1 , then for every  i { 1 , , n 1 }  we define  s i { 1 , 2 }  in the following way:  s i = 1  if the hexagon  h i  is linear; and  s i = 2  if the hexagon  h i  is angular. By using the stated proposition and initial values  M Z Z ( M P h 0 ) = 1 M Z Z ( M P h 1 ) = x + y + z + 3 M Z Z ( M P h 1 * ) = x + 2 , we obtain Algorithm 2 which computes the  M Z Z  polynomial for any (reduced) modified phenylene chain. It can be noticed that this algorithm is a slight modification of Algorithm 1.
Algorithm 2: The MZZ polynomial of a (reduced) modified phenylene chain
Axioms 12 01053 i002
We can also check that Algorithm 2 correctly calculates the MZZ polynomial of a (reduced) modified phenylene chain in linear time with respect to the number of hexagons (in the model where addition and multiplication of polynomials of three variables can be performed in constant time).

5. Phenylenes with Branched Hexagons

In this section we show how the MZZ polynomial can be calculated for phenylenes that contain some branched hexagons. Note that the main theorem of the section is similar to analogous results from [5,17], but some technical differences are needed. In the theorem, the following assumption will be required.
Assumption 2.
Let  G 1 G 2 , and  G 3  be generalized phenylenes such that each of them has a perfect matching. Moreover, for any  i { 1 , 2 , 3 }  let  e i = u i v i  be a boundary edge of  G i  such that  e i  also lies on some quadrilateral  f i  of  G i .
By adding edges  u 1 v 3 v 1 u 2 , and  v 2 u 3  to the disjoint union of graphs  G 1 G 2 , and  G 3  satisfying Assumption 2, we obtain another generalized phenylene and denote it as  G 1 G 2 G 3 , see Figure 9 for an example. Moreover, let f be the hexagon of  G 1 G 2 G 3  that contains edges  u 1 v 3 v 1 u 2 , and  v 2 u 3 . It is easy to see that the graph  G 1 G 2 G 3  also has a perfect matching. Similarly as in Section 3, for any  i { 1 , 2 , 3 }  denote by  G i  the graph  G i u i v i  and by  G i  the graph  G i f i .
In the next theorem, we calculate the MZZ polynomial of  G 1 G 2 G 3 .
Theorem 7.
If  G 1 , G 2 , and  G 3  are three generalized phenylenes satisfying Assumption 2, then the MZZ polynomial of  G = G 1 G 2 G 3  can be computed as
M Z Z ( G ) = j = 1 3 M Z Z ( G j ) + ( y + 1 ) j = 1 3 M Z Z ( G j ) + z M Z Z ( G 1 ) M Z Z ( G 2 ) M Z Z ( G 3 ) + z M Z Z ( G 1 ) M Z Z ( G 2 ) M Z Z ( G 3 ) + z M Z Z ( G 1 ) M Z Z ( G 2 ) M Z Z ( G 3 ) .
Proof. 
Suppose that we have the notation from Assumption 2 (see also Figure 9). By using Theorem 3 on the edge  u 1 v 3  of G one can obtain
M Z Z ( G ) = M Z Z ( G u 1 v 3 ) + M Z Z ( G u 1 v 3 ) + y M Z Z ( G f ) + z M Z Z ( G f f 1 ) + z M Z Z ( G f f 2 ) + z M Z Z ( G f f 3 ) .
It is easy to see that edges  v 1 u 2  and  v 2 u 3  do not belong to any perfect matching of  G u 1 v 3  and, therefore, by Proposition 2 and Theorem 1 it follows
M Z Z ( G u 1 v 3 ) = M Z Z ( G u 1 v 3 v 1 u 2 v 2 u 3 ) = j = 1 3 M Z Z ( G j ) .
In a similar way we can conclude that edges  v 1 u 2  and  v 2 u 3  belong to every perfect matching of  G u 1 v 3  (if this graph has a perfect matching) and, consequently, by Proposition 2 and Theorem 1 we obtain
M Z Z ( G u 1 v 3 ) = M Z Z ( G f ) = j = 1 3 M Z Z ( G j ) .
Furthermore, Theorem 1 implies
M Z Z ( G f f 1 ) = M Z Z ( G 1 ) M Z Z ( G 2 ) M Z Z ( G 3 )
and analogous formulas can be obtained for  M Z Z ( G f f 2 )  and  M Z Z ( G f f 3 ) . By taking the obtained equalities into the first formula of this proof, we finally obtain the desired result. □
By using Theorem 7, we can reduce the problem of calculating the MZZ polynomial of any phenylene to the problem of calculating the MZZ polynomials of (reduced) phenylene chains and (reduced) modified phenylene chains, which can be easily solved by using Algorithms 1 and 2.
To show an example, let  G = G 1 G 2 G 3  be a phenylene from Figure 10.
We can easily calculate that
M Z Z ( G 1 ) = x 2 + x y + y 2 + x z + y z + 5 x + 5 y + 5 z + 8 , M Z Z ( G 1 ) = y 2 + x + 4 y + 2 z + 5 , M Z Z ( G 1 ) = x + y + z + 3
and
M Z Z ( G 2 ) = M Z Z ( G 3 ) = x + y + z + 3 , M Z Z ( G 2 ) = M Z Z ( G 3 ) = y + 2 , M Z Z ( G 2 ) = M Z Z ( G 3 ) = 1 .
Therefore, Theorem 7 finally gives
M Z Z ( G ) = y 5 + x 4 + 3 x 3 y + 4 x 2 y 2 + 4 x y 3 + 10 y 4 + 3 x 3 z + 7 x 2 y z + 8 x y 2 z + 8 y 3 z + 3 x 2 z 2 + 5 x y z 2 + 4 y 2 z 2 + x z 3 + y z 3 + 11 x 3 + 27 x 2 y + 32 x y 2 + 44 y 3 + 27 x 2 z + 54 x y z + 56 y 2 z + 21 x z 2 + 29 y z 2 + 5 z 3 + 47 x 2 + 93 x y + 108 y 2 + 93 x z + 143 y z + 50 z 2 + 97 x + 149 y + 133 z + 92 .
A more complicated case is obtained if we have a situation from Figure 11, where G is a phenylene containing two branched hexagons with one quadrilateral between them. To calculate  M Z Z ( G 1 ) , we can firstly apply Corollary 2 for the face  f 1 . Note that in this case Theorem 7 and Corollary 2 both require to calculate  M Z Z ( G 1 ) = M Z Z ( G 1 f 1 ) . However, we notice that by Theorem 1 it holds  M Z Z ( G 1 ) = M Z Z ( G 1 1 ) M Z Z ( G 1 2 )  since the edge  e  can not belong to any perfect matching of  G 1 = G 1 f 1 . Here  G 1 1  and  G 1 2  are the two connected components of the graph  G 1 e  (see Figure 11).

6. Conclusions

In the paper, we have introduced the MZZ polynomial for any generalized phenylene. By definition, this polynomial counts the number of so called (4,6,8)-Clar covers with specific numbers of 4-cycles, 6-cycles, and 8-cycles. We provided several results that enable us to calculate the MZZ polynomial of a generalized phenylene by using the MZZ polynomials of subgraphs of the original graph. Then, we focused on phenylene chains and developed an efficient algorithm that calculates the MZZ polynomial for any (reduced) phenylene chain in linear time with respect to the number of hexagons (assuming that addition and multiplication of polynomials of three variables can be calculated in constant time). Furthermore, the main result of Section 5 can be used to calculate the MZZ polynomial of a phenylene with branched hexagons.
Regarding the future work, it would be interesting to investigate the chemical applicability of the MZZ polynomial. In particular, relations between the MZZ polynomial and some energy-based quantities of phenylenes can be tested. Moreover, it would be nice to obtain new mathematical properties of the MZZ polynomial and investigate this polynomial for some special classes of phenylenes.

Funding

The author acknowledges the financial support from the Slovenian Research and Innovation Agency (research programme No. P1-0297 and project No. N1-0285).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Cyvin, S.J.; Gutman, I. Kekulé Structures in Benzenoid Hydrocarbons; Springer: Berlin, Germany, 1988. [Google Scholar]
  2. Gutman, I.; Cyvin, S.J. Introduction to the Theory of Benzenoid Hydrocarbons; Springer: Berlin, Germany, 1989. [Google Scholar]
  3. Gojak, S.; Gutman, I.; Radenković, S.; Vodopivec, A. Relating resonance energy with Zhang-Zhang polynomial. J. Serb. Chem. Soc. 2007, 72, 673–679. [Google Scholar] [CrossRef]
  4. Clar, E. The Aromatic Sextet; John Wiley & Sons: London, UK, 1972. [Google Scholar]
  5. Zhang, H.; Zhang, F. The Clar covering polynomial of hexagonal systems I. Discrete Appl. Math. 1996, 69, 147–167. [Google Scholar] [CrossRef]
  6. Fowler, P.W.; Myrvold, W.; Vandenberg, R.L.; Hartung, E.J.; Graver, J.E. Clar and Fries structures for fullerenes. Art Discret. Appl. Math. 2022. [Google Scholar] [CrossRef]
  7. Graver, J.E.; Hartung, E.J.; Souid, A.Y. Clar and Fries numbers for benzenoids. J. Math. Chem. 2013, 51, 1981–1989. [Google Scholar] [CrossRef]
  8. Klavžar, S.; Salem, K.; Taranenko, A. Maximum cardinality resonant sets and maximal alternating sets of hexagonal systems. Comput. Math. Appl. 2010, 59, 506–513. [Google Scholar] [CrossRef]
  9. Zhai, S.; Alrowaili, D.; Ye, D. Clar structures vs Fries structures in hexagonal systems. Appl. Math. Comput. 2018, 329, 384–394. [Google Scholar] [CrossRef]
  10. Langner, J.; Witek, H.A. Interface theory of benzenoids. MATCH Commun. Math. Comput. Chem. 2020, 84, 143–176. [Google Scholar]
  11. Li, G.; Pei, Y.; Wang, Y. Clar Covering polynomials with only real zeros. MATCH Commun. Math. Comput. Chem. 2020, 84, 217–228. [Google Scholar]
  12. Tratnik, N.; Ye, D. Resonance graphs on perfect matchings of graphs on surfaces. Graphs Comb. 2023, 39, 68. [Google Scholar] [CrossRef]
  13. Zhang, H. Z-transformation graphs of perfect matchings of plane bipartite graphs: A survey. MATCH Commun. Math. Comput. Chem. 2006, 56, 457–476. [Google Scholar] [CrossRef]
  14. Zhang, H.; Shiu, W.C.; Sun, P.K. A relation between Clar covering polynomial and cube polynomial. MATCH Commun. Math. Comput. Chem. 2013, 70, 477–492. [Google Scholar]
  15. Che, Z. Peripheral convex expansions of resonance graphs. Order 2021, 38, 365–376. [Google Scholar] [CrossRef]
  16. Žigert Pleteršek, P. Equivalence of the generalized Zhang-Zhang polynomial and the generalized cube polynomial. MATCH Commun. Math. Comput. Chem. 2018, 80, 215–226. [Google Scholar]
  17. Furtula, B.; Radenković, S.; Redžepović, I.; Tratnik, N.; Žigert Pleteršek, P. The generalized Zhang-Zhang polynomial of benzenoid systems—Theory and applications. Appl. Math. Comput. 2022, 418, 126822. [Google Scholar] [CrossRef]
  18. Radenković, S.; Redžepović, I.; Furtula, B.; Đorđević, S.; Tratnik, N.; Žigert Pleteršek, P. Relating vibrational energy with Kekulé- and Clar-structure-based parameters. Int. J. Quantum Chem. 2022, 122, e26867. [Google Scholar] [CrossRef]
  19. Gutman, I. Algebraic structure count of linear phenylenes. Indian J. Chem. 1993, 32A, 281–284. [Google Scholar]
  20. Gutman, I. Algebraic structure count of linear phenylenes and their congeners. J. Serb. Chem. Soc. 2003, 68, 391–399. [Google Scholar] [CrossRef]
  21. Gutman, I.; Furtula, B. A Kekulé structure basis for phenylenes. J. Mol. Struct. Theochem. 2006, 770, 67–71. [Google Scholar] [CrossRef]
  22. Brinkmann, G.; Caporossi, G.; Hansen, P. A constructive enumeration of fusenes and benzenoids. J. Algorithms 2002, 45, 155–166. [Google Scholar] [CrossRef]
Figure 1. A benzenoid graph B and the corresponding phenylene G.
Figure 1. A benzenoid graph B and the corresponding phenylene G.
Axioms 12 01053 g001
Figure 2. A (4,6,8)-Clar cover of a phenylene.
Figure 2. A (4,6,8)-Clar cover of a phenylene.
Axioms 12 01053 g002
Figure 3. Generalized phenylene G in Theorem 3.
Figure 3. Generalized phenylene G in Theorem 3.
Axioms 12 01053 g003
Figure 4. Graph  G 1 · G 2  in Theorem 5.
Figure 4. Graph  G 1 · G 2  in Theorem 5.
Axioms 12 01053 g004
Figure 5. (a) The basic compound of a phenylene P and (b) the graph  G 1 · P .
Figure 5. (a) The basic compound of a phenylene P and (b) the graph  G 1 · P .
Axioms 12 01053 g005
Figure 6. Linear phenylene chain  P L 6 .
Figure 6. Linear phenylene chain  P L 6 .
Axioms 12 01053 g006
Figure 7. A phenylene chain  P h 4  and the corresponding reduced phenylene chain  P h 4 * .
Figure 7. A phenylene chain  P h 4  and the corresponding reduced phenylene chain  P h 4 * .
Axioms 12 01053 g007
Figure 8. A modified phenylene chain  M P h 4  and the corresponding reduced modified phenylene chain  M P h 4 * .
Figure 8. A modified phenylene chain  M P h 4  and the corresponding reduced modified phenylene chain  M P h 4 * .
Axioms 12 01053 g008
Figure 9. A generalized phenylene  G = G 1 G 2 G 3 .
Figure 9. A generalized phenylene  G = G 1 G 2 G 3 .
Axioms 12 01053 g009
Figure 10. A phenylene  G = G 1 G 2 G 3 .
Figure 10. A phenylene  G = G 1 G 2 G 3 .
Axioms 12 01053 g010
Figure 11. A phenylene  G = G 1 G 2 G 3 , where  G 1 1  and  G 1 2  are connected components of  G 1 e .
Figure 11. A phenylene  G = G 1 G 2 G 3 , where  G 1 1  and  G 1 2  are connected components of  G 1 e .
Axioms 12 01053 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tratnik, N. The Multivariable Zhang–Zhang Polynomial of Phenylenes. Axioms 2023, 12, 1053. https://doi.org/10.3390/axioms12111053

AMA Style

Tratnik N. The Multivariable Zhang–Zhang Polynomial of Phenylenes. Axioms. 2023; 12(11):1053. https://doi.org/10.3390/axioms12111053

Chicago/Turabian Style

Tratnik, Niko. 2023. "The Multivariable Zhang–Zhang Polynomial of Phenylenes" Axioms 12, no. 11: 1053. https://doi.org/10.3390/axioms12111053

APA Style

Tratnik, N. (2023). The Multivariable Zhang–Zhang Polynomial of Phenylenes. Axioms, 12(11), 1053. https://doi.org/10.3390/axioms12111053

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop