The Annealed Entropy of Wiener Number on Random Double Hexagonal Chains

Abstract

We study a random planar honeycomb lattice model, namely the random double hexagonal chains. This is a lattice system with nonperiodic boundary condition. The Wiener number is an important molecular descriptor based on the distances, which was introduced by the chemist Harold Wiener in 1947. By applying probabilistic method and combinatorial techniques we obtained an explicit analytical expression for the expected value of Wiener number of a random double hexagonal chain, and the limiting behaviors on the annealed entropy of Wiener number when the random double hexagonal chain becomes infinite in length are analyzed.

Share and Cite:

Ren, H. and Su, X. (2017) The Annealed Entropy of Wiener Number on Random Double Hexagonal Chains. Applied Mathematics, 8, 1473-1480. doi: 10.4236/am.2017.810108.

1. Introduction

Topological indices (molecular structure descriptors) based on the distances between the vertices of a graph are widely used in theoretical chemistry to establish relations between the structure and the properties of molecules and provide correlations with physical, chemical, and thermodynamic parameters of chemical compounds [1] . Among the variety of these indices, the Wiener number, denoted by W ( G ) , is the best known one which was introduced by the chemist Harold Wiener in 1947 [2] as a simple parameter. Wiener number has been found to correlate with various physicochemical properties of a molecule (modeled by a graph): Boiling point, heat of vaporization, heat of isomerization, surface energy, specific dispersion, and sound velocity. In addition, the parameter also correlates with some π-electron characteristics of conjugated polymers; for example, the total π-electron energy and HOMO-LUMO (Highest Occupied Molecular Orbital-Lowest Unoccupied Molecular Orbital) separation, for details see [3] [4] and the references therein.

Let G be a connected graph with vertices labeled as 1,2, , n . The standard distance between two vertices i and j, denoted by d G ( i , j ) , is the length of a shortest path connecting them. In terms of graph theory, Wiener number [2] [3] is defined as the sum of distances between all pairs of vertices in G, i.e.,

W ( G ) = i < j d G ( i , j ) = 1 2 i V ( G ) j V ( G ) d G ( i , j ) = 1 2 i V ( G ) d ( i | G ) , (1)

where d ( i | G ) is the Wiener number of vertex i in G, defined by d ( i | G ) = j G i d ( i , j ) .

A hexagonal system is a 2-connected plane graph whose every interior face is bounded by a regular hexagon of unit length 1. Hexagonal systems are of great importance for theoretical chemistry because they are the natural graph representation of benzenoid hydrocarbons [5] [6] [7] [8] . A hexagonal system H is said to be catacondensed if all its vertices are on the outerface, otherwise H is said to be pericondensed. In [4] , Gutman et al. obtained an explicit analytical expression for the expected value of the Wiener number of a random benzenoid chain with n hexagons(a graph of unbranched catacondensed benzenoid-like structure). The random multiple chain was introduced in [9] , the generating procedure of which is inspired by the growth of single walled zigzag nanotubes [10] . Some results on double hexagonal chains (a special type of pericondensed hexagonal system with nonperiodic boundary condition, which is constructed by successively fusing a series of naphthalenes), can be found, for example, in [11] [12] [13] [14] [15] and the references therein. In statistical mechanics, entropy is related to the number of microscopic configurations that a thermodynamic system can have when in a state as specified by some macroscopic variables. In this paper, we study the annealed entropy of Wiener number on random double hexagonal chains.

The random double hexagonal chain R D 2 × n can be obtained from a naphthalene by stepwise triple-edge fusion of a new naphthalene. For convenience, we orient each naphthalene so that its interior edges are horizontal. There are two types of triple-edge fusion of two naphthalenes: α-type fusion and β-type fusion, as shown in Figure 1(a). At each step k ( k 3 ) , a random selection is taken from one of the two possible fusions: α-type fusion with probability p and β-type fusion with probability 1 p . In our model, we assume that the probability p is a constant, invariant to the step parameter k. That is, the process described is a zeroth-order Markov process. For n 2 , R D 2 × n is a pericondensed hexagonal system. Random (double) hexagonal chains offer a good model for a class of conjugated polymers [16] , many features of which have already been established [4] [11] [16] . By applying probabilistic method and combinatorial techniques an explicit analytical expression for the expected value W n ( = E ( W ( R D 2 × n ) ) ) of the Wiener number of a random double hexagonal chain with n naphthalenes is obtained. We note that the expression of W n is a

Figure 1. The triple-edge fusions. (a) α-type fusion and β-type fusion. (b) D 2 × n is obtained by two possible fusions.

polynomial in the variable n. Furthermore, the limiting behaviors on the annealed entropy of Wiener number when the random double hexagonal chain becomes infinite in length are analyzed.

2. Some Fundamental Recursion Relations

A double hexagonal chain D 2 × n with n naphthalenes can be constructed from D 2 × ( n 1 ) by attaching a new terminal naphthalene spanned by vertices { s n 1 , t n 1 , g n 1 , h n 1 , x 6 , x 5 , x 4 , x 3 , x 2 , x 1 } (or { r n 1 , s n 1 , t n 1 , g n 1 , x 1 , x 2 , x 3 , x 4 , x 5 , x 6 } ) (see Figure 1(b)).

Lemma 2.1. If D 2 × n is obtained from D 2 × ( n 1 ) by α-type fusion, then

d ( x 1 | D 2 × n ) = 6 n + 13 + d ( s n 1 | D 2 × ( n 1 ) ) , d ( x 2 | D 2 × n ) = 12 n + 1 + d ( s n 1 | D 2 × ( n 1 ) ) ,

d ( x 3 | D 2 × n ) = 6 n + 7 + d ( g n 1 | D 2 × ( n 1 ) ) , d ( x 4 | D 2 × n ) = 12 n + 5 + d ( g n 1 | D 2 × ( n 1 ) ) ,

d ( x 5 | D 2 × n ) = 12 n + 7 + d ( h n 1 | D 2 × ( n 1 ) ) , d ( x 6 | D 2 × n ) = 6 n + 13 + d ( h n 1 | D 2 × ( n 1 ) ) ,

and W ( D 2 × n ) = W ( D 2 × ( n 1 ) ) + 2 [ d ( s n 1 | D 2 × ( n 1 ) ) + d ( g n 1 | D 2 × ( n 1 ) ) + d ( h n 1 | D 2 × ( n 1 ) ) ] + 54 n + 11.

Similarly, if D 2 × n is obtained from D 2 × ( n 1 ) by β-type fusion, then

d ( x 1 | D 2 × n ) = 6 n + 13 + d ( g n 1 | D 2 × ( n 1 ) ) , d ( x 2 | D 2 × n ) = 12 n + 1 + d ( g n 1 | D 2 × ( n 1 ) ) ,

d ( x 3 | D 2 × n ) = 6 n + 7 + d ( s n 1 | D 2 × ( n 1 ) ) , d ( x 4 | D 2 × n ) = 12 n + 5 + d ( s n 1 | D 2 × ( n 1 ) ) ,

d ( x 5 | D 2 × n ) = 12 n + 7 + d ( r n 1 | D 2 × ( n 1 ) ) , d ( x 6 | D 2 × n ) = 6 n + 13 + d ( r n 1 | D 2 × ( n 1 ) ) ,

and W ( D 2 × n ) = W ( D 2 × ( n 1 ) ) + 2 [ d ( s n 1 | D 2 × ( n 1 ) ) + d ( g n 1 | D 2 × ( n 1 ) ) + d ( r n 1 | D 2 × ( n 1 ) ) ] + 54 n + 11 ,

where d ( s 1 | D 2 × 1 ) = d ( g 1 | D 2 × 1 ) = 21 , d ( h 1 | D 2 × 1 ) = d ( r 1 | D 2 × 1 ) = 25 .

Proof. By (1), we have d ( s 1 | D 2 × 1 ) = d ( g 1 | D 2 × 1 ) = 21 , d ( h 1 | D 2 × 1 ) = d ( r 1 | D 2 × 1 ) = 25 . Since a naphthalene is not a vertex rotation symmetry, we must distinguish between two different situations for the Wiener number of the vertex in the terminal naphthalene. In the following, we only consider the case of α-type fusion (the argument for the case of β-type fusion is analogous). Note that D 2 × ( n 1 ) has 6 n 2 vertices. If D 2 × n is obtained from D 2 × ( n 1 ) by α-type fusion then, for n 2 , we have the following relations:

d ( x 1 | D 2 × n ) = 6 n 2 + d ( s n 1 | D 2 × ( n 1 ) ) + 15 ,

d ( x 2 | D 2 × n ) = 2 ( 6 n 2 ) + Δ n 1 6 + 11 ,

d ( x 3 | D 2 × n ) = 6 n 2 + d ( g n 1 | D 2 × ( n 1 ) ) + 9 ,

d ( x 4 | D 2 × n ) = 2 ( 6 n 2 ) + d ( g n 1 | D 2 × ( n 1 ) ) + 9 ,

d ( x 5 | D 2 × n ) = 2 ( 6 n 2 ) + d ( h n 1 | D 2 × ( n 1 ) ) + 11 ,

d ( x 6 | D 2 × n ) = 6 n 2 + d ( h n 1 | D 2 × ( n 1 ) ) + 15.

and W ( D 2 × n ) = W ( D 2 × ( n 1 ) ) + d ( s n 1 | D 2 × ( n 1 ) ) + 2 [ d ( g n 1 | D 2 × ( n 1 ) ) + d ( h n 1 | D 2 × ( n 1 ) ) ] + Δ n 1 + 54 n + 11 ,

where Δ n 1 = min { d ( s n 1 | D 2 × ( n 1 ) ) , d ( g n 1 | D 2 × ( n 1 ) ) } . In this case, s n 1 and g n 1 coincide with the vertices x 2 and x 4 , respectively. Thus, d ( s n 1 | D 2 × ( n 1 ) ) d ( g n 1 | D 2 × ( n 1 ) ) = Δ n 2 4 d ( g n 2 | D 2 × ( n 1 ) ) < 0 . By induction on n, we know that Δ n 1 = d ( s n 1 | D 2 × ( n 1 ) ) . The proof is completed.

For a random double benzenoid chain R D 2 × n , there are two cases to be considered:

Case 1. R D 2 × n R D 2 × ( n + 1 ) by α-type fusion with probability p. In this case, r n , s n , t n , g n , h n (of R D 2 × n ) coincide with the vertices x 1 , x 2 , x 3 , x 4 , x 5 , respectively.

Case 2. R D 2 × n R D 2 × ( n + 1 ) by β-type fusion with probability 1 p . In this case, r n , s n , t n , g n , h n (of R D 2 × n ) coincide with the vertices x 5 , x 4 , x 3 , x 2 , x 1 , respectively.

The distance d ( r n | R D 2 × n ) , d ( s n | R D 2 × n ) , d ( g n | R D 2 × n ) , d ( h n | R D 2 × n ) and W ( R D 2 × ( n 1 ) ) are random variables. We denote their expected values by n = E ( d ( r n | R D 2 × n ) ) , S n = E ( d ( s n | R D 2 × n ) ) , G n = E ( d ( g n | R D 2 × n ) ) , n = E ( d ( h n | R D 2 × n ) ) and W n = E ( W ( R D 2 × n ) ) (or E ( W ( n ) ) in brief), respectively. Then, by Case 1, Case 2 and Lemma 2.1 we have

n = p ( 6 n + 13 + d ( s n 1 | R D 2 × ( n 1 ) ) ) + ( 1 p ) ( 12 n + 7 + d ( r n 1 | R D 2 × ( n 1 ) ) ) ;

S n = p ( 12 n + 1 + d ( s n 1 | R D 2 × ( n 1 ) ) ) + ( 1 p ) ( 12 n + 5 + d ( s n 1 | R D 2 × ( n 1 ) ) ) ;

G n = p ( 12 n + 5 + d ( g n 1 | R D 2 × ( n 1 ) ) ) + ( 1 p ) ( 12 n + 1 + d ( g n 1 | R D 2 × ( n 1 ) ) ) ;

n = p ( 12 n + 7 + d ( h n 1 | R D 2 × ( n 1 ) ) ) + ( 1 p ) ( 6 n + 13 + d ( g n 1 | R D 2 × ( n 1 ) ) )

and

W n = W ( R D 2 × ( n 1 ) ) + 2 ( d ( s n 1 | R D 2 × ( n 1 ) ) + d ( g n 1 | R D 2 × ( n 1 ) ) ) + 2 p d ( h n 1 | R D 2 × ( n 1 ) ) + 2 ( 1 p ) d ( r n 1 | R D 2 × ( n 1 ) ) + 54 n + 11.

Theorem 2.1. For n 2 , then

n = p ( 6 n + 13 + S n 1 ) + ( 1 p ) ( 12 n + 7 + n 1 ) ; S n = 12 n + 5 4 p + S n 1 ; G n = 12 n + 1 + 4 p + G n 1 ; n = p ( 12 n + 7 + n 1 ) + ( 1 p ) ( 6 n + 13 + G n 1 ) ; W n = W n 1 + 2 ( S n 1 + G n 1 ) + 2 p n 1 + 2 ( 1 p ) n 1 + 54 n + 11 (2)

with boundary conditions W 1 = 109 , S 1 = G 1 = 21 and 1 = 1 = 25 .

Proof. Noting that E ( n ) = n , E ( S n ) = S n , E ( G n ) = G n , E ( n ) = n and E ( W n ) = W n . Thus, Theorem 2.1 holds.

3. The Explicit Analytical Expression for W n

From (2), by successive subtraction method we have

S n = 6 n 2 + 11 n 4 p n + 4 p + 4 , G n = 6 n 2 + 7 n + 4 p n 4 p + 8. (3)

Therefore,

n = 12 n + 5 p 7 n p + 6 n 2 p + 8 p 2 4 n p 2 + 7 + ( 1 p ) n 1 , n = 20 + n + 6 n 2 21 p + 15 n p 6 n 2 p + 8 p 2 4 n p 2 + p n 1 , W n = W n 1 + 2 p n 1 + 2 ( 1 p ) n 1 + 23 + 42 n + 24 n 2 . (4)

Thus, we get the following result.

Theorem 3.1. For n 2 , if p = 0 then n = 6 n 2 + 13 n + 6 , S n = 6 n 2 + 11 n + 4 , G n = 6 n 2 + 7 n + 8 , n = 20 + n + 6 n 2 , W n = 12 n 3 + 40 n 2 + 49 n + 8 ; and if p = 1 then n = 6 n 2 + n + 20 , S n = 6 n 2 + 7 n + 8 , G n = 6 n 2 + 11 n + 4 , n = 6 n 2 + 13 n + 6 , W n = 12 n 3 + 40 n 2 + 49 n + 8 .

We now consider the case when 0 < p < 1 by using the method of generating functions. Let

( t ) = n 0 n t n , ( t ) = n 0 n t n , W ( t ) = n 0 W n t n .

Then, by (4) we get

( t ) 0 = ( 1 p ) t ( t ) + ( 19 + 4 p + 4 p 2 ) δ 3 + ( 12 + 5 p 4 p 2 ) δ 2 + 6 p δ 1 , ( t ) 0 = p t ( t ) + ( 27 12 p + 4 p 2 ) δ 3 + ( 13 + 3 p 4 p 2 ) δ 2 + ( 6 6 p ) δ 1 , W ( t ) W 0 = t W ( t ) + 2 p t ( t ) + 2 ( 1 p ) t ( t ) + 89 δ 3 + 90 δ 2 + 24 δ 1 , (5)

where δ 1 = n 0 n 2 t n + 1 , δ 2 = n 0 n t n + 1 and δ 3 = n 0 t n + 1 . Since ( 1 t ) j = n 0 C n j 1 t n j + 1 . Then ( 1 p t ) 1 = n 0 ( p t ) n and ( 1 t ) 1 = n 0 t n . Thus, by (5) we have

( t ) = n 0 ( ( 1 p ) t ) n [ 0 + ( 19 + 4 p + 4 p 2 ) m 1 t m + ( 12 + 5 p 4 p 2 ) m 1 ( m 1 ) t m + 6 p m 1 ( m 1 ) 2 t m ] = n 0 0 ( ( 1 p ) t ) n + k 0 n + m = k , m 1 ( 1 p ) n [ ( 19 + 4 p + 4 p 2 ) + ( 12 + 5 p 4 p 2 ) ( m 1 ) + 6 p ( m 1 ) 2 ] t k = n 0 0 ( ( 1 p ) t ) n + n 0 m = 1 n ( 1 p ) n m [ ( 19 + 4 p + 4 p 2 ) + ( 12 + 5 p 4 p 2 ) ( m 1 ) + 6 p ( m 1 ) 2 ] ] t n = n 0 n t n ,

where,

n = 0 ( 1 p ) n + m = 2 n ( 1 p ) n m [ 7 + 12 m + 5 p 7 m p + 6 m 2 p + 8 p 2 4 m p 2 ] ,

( t ) = n 0 ( p t ) n ( 0 + m 1 [ ( 27 12 p + 4 p 2 ) + ( 13 + 3 p 4 p 2 ) ( m 1 ) + ( 6 6 p ) ( m 1 ) 2 ] t m ) = n 0 0 ( p t ) n + k 0 n + m = k , m 1 p n ( 20 + m + 6 m 2 21 p + 15 m p 6 m 2 p + 8 p 2 4 m p 2 ) t k = n 0 0 ( p t ) n + n 0 m = 1 n p n m [ 20 + m + 6 m 2 21 p + 15 m p 6 m 2 p + 8 p 2 4 m p 2 ] t n = n 0 n t n ,

n = 0 p n + m = 1 n p n m [ 20 + m + 6 m 2 21 p + 15 m p 6 m 2 p + 8 p 2 4 m p 2 ] ,

W ( t ) = n 0 t n [ W 0 + 2 p m 1 m 1 t m + 2 ( 1 p ) m 1 m 1 t m + m 1 ( 23 + 42 m + 24 m 2 ) t m ] = n 0 [ W 0 + m = 1 n ( 2 p m 1 + 2 ( 1 p ) m 1 + 23 + 42 m + 24 m 2 ) t n ] = n 0 W n t n ,

W n = W 0 + m = 1 n ( 2 p m 1 + 2 ( 1 p ) m 1 + 23 + 42 m + 24 m 2 ) .

Thus, we reach the following result.

Theorem 3.2. If n 2 and 0 < p < 1 , then

n = 0 ( 1 p ) n + m = 2 n ( 1 p ) n m [ 7 + 12 m + 5 p 7 m p + 6 m 2 p + 8 p 2 4 m p 2 ] , S n = 6 n 2 + 11 n 4 p n + 4 p + 4 , G n = 6 n 2 + 7 n + 4 p n 4 p + 8 , n = 0 p n + m = 1 n p n m [ 20 + m + 6 m 2 21 p + 15 m p 6 m 2 p + 8 p 2 4 m p 2 ] , W n = W 0 + m = 1 n ( 2 p m 1 + 2 ( 1 p ) m 1 + 23 + 42 m + 24 m 2 ) . (6)

It is easy to verify that m = 0 n p n m = 1 p n + 1 1 p , m = 1 n m p n m = p n ( n + 1 ) p 1 + p n 1 ( p 1 ) 2 , m = 1 n m p m 1 = 1 ( n + 1 ) p n 1 p + p ( 1 p n ) ( 1 p ) 2 and

m = 1 n m 2 p n m = ( n 2 p + 2 n p + 2 n 2 p p 2 2 n p 2 n 2 p 2 + p 1 + n + p 2 + n ) / ( 1 + p ) 3 .

So by (4), we can calculate the values of 0 , 0 and W 0 , and by using Mathematica software(Mathematica 9.0) to (6) we obtain the solution of (6) as follows:

Theorem 3.3. For n 2 and 0 < p < 1 , then

n = 1 ( 1 p ) p ( 8 8 ( 1 p ) n + 5 n p + 6 n 2 p + 6 ( 1 p ) n p 4 p 2 9 n p 2 6 n 2 p 2 4 p 3 + 4 n p 3 ) ,

S n = 6 n 2 + 11 n 4 p n + 4 p + 4 ,

G n = 6 n 2 + 7 n + 4 p n 4 p + 8 ,

n = 1 ( 1 + p ) p ( 20 p n p 6 n 2 p + 16 p 2 3 n p 2 + 6 n 2 p 2 4 p 3 + 4 n p 3 + 2 p n + 6 p 1 + n ) ,

W n = 1 ( 1 + p ) 2 p 2 ( 16 + 16 ( 1 p ) n + 44 p + 16 n p 44 ( 1 p ) n p 32 p 2 + 13 n p 2 + 32 n 2 p 2 + 12 n 3 p 2 + 40 ( 1 p ) n p 2 40 p 3 34 n p 3 72 n 2 p 3 24 n 3 p 3 12 ( 1 p ) n p 3 + 60 p 4 43 n p 4 + 56 n 2 p 4 + 12 n 3 p 4 48 p 5 + 72 n p 5 24 n 2 p 5 + 16 p 6 24 n p 6 + 8 n 2 p 6 + 4 p 2 + n + 12 p 3 + n )

4. Conclusions and Suggestions

From Theorem 3.1 and Theorem 3.3, we know that W n ~ 12 n 3 as n tends to infinity. Note that W n is a polynomial in the variable n, which is different from the expected value of the Wiener number of the random benzenoid chain in the general case. And the limiting behaviors on the annealed entropy l o g W n when the random double hexagonal chain becomes infinite in length are

lim m 1 V log [ W n ] = 0 ,

where V is the number of the vertices in R D 2 × n , i.e., V = 6 n + 4 .

Open problem. For any p 0,1 and n 1 , is it true that the annealed entropy l o g W n of m-tuple random hexagonal chain is zero as m or m , n .

Acknowledgements

Sincere thanks to the members of JAMP for their professional performance, and special thanks to managing editor for a rare attitude of high quality. This research supported by NSFC(11551003) and the Natural Science Foundation of Qinghai(2015-ZJ-911).

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Gutman, I. and Furtula, B., Eds. (2012) Distance in Molecular Graphs-Theory (Distance in Molecular Graphs-Applications). University of Kragujevac, Kragujevac.
[2] Wiener, H. (1947) Structural Determination of Paraffin Boiling Points. Journal of the American Chemical Society, 69, 17-20. https://doi.org/10.1021/ja01193a005
[3] Dobrynin, A., Entringer, R. and Gutman, I. (2001) Wiener Index of Trees: Theory and Applications. Acta Applicandae Mathematicae, 66, 211-249.
https://doi.org/10.1023/A:1010767517079
[4] Gutman, I., Kennedy, J.W. and Quintas, L.V. (1990) Wiener Numbers of Random Benzenoid Chains. Chemical Physics Letters, 173, 403-408. https://doi.org/10.1016/0009-2614(90)85292-K
[5] Gutman, I. and Polansky, O.E. (1986) Mathematical Concepts in Organic Chemistry. Springer, Berlin. https://doi.org/10.1007/978-3-642-70982-1
[6] Gutman, I. and Cyvin, S.J., Eds. (1989) Introduction to the Theory of Benzenoid Hydrocarbons. Springer, Berlin. https://doi.org/10.1007/978-3-642-87143-6
[7] Gutman, I. and Cyvin, S.J., Eds. (1990) Advances in the Theory of Benzenoid Hydrocarbons. Topics in Current Chemistry, Vol. 153, Springer, Berlin. https://doi.org/10.1007/3-540-51505-4
[8] Gutman, I., Ed. (1992) Advances in the Theory of Benzenoid Hydrocarbons I. Topics in Current Chemistry, Vol. 162, Springer, Berlin.
[9] Ren, H.Z., Zhang, F.J. and Qian, J.G. (2012) Dimer Coverings on Random Multiple Chains of Planar Honeycomb Lattices. Journal of Statistical Mechanics: Theory and Experiment, 2012, P08002. https://doi.org/10.1088/1742-5468/2012/08/P08002
[10] Dresselhaus, M.S., Dresselhaus, G. and Avouris, P. (2001) Carbon Nanotubes: Synthesis, Structure, Properties and Applications. Springer, Berlin. https://doi.org/10.1007/3-540-39947-X
[11] Gutman I., Randic, M., Balaban, A., Furtula, B. and Vuckovic, V. (2005) π-Electron Contents of Rings in the Double-Hexagonal-Chain Homologous Series (Pyrene, Anthanthrene and Other Acenoacenes). Polycyclic Aromatic Compounds, 25, 215-226.
https://doi.org/10.1080/10406630591007080
[12] Ren, H.Z. and Zhang, F.J. (2007) Double Hexagonal Chains with Maximal Energy. International Journal of Quantum Chemistry, 107, 1437-1445. https://doi.org/10.1002/qua.21256
[13] Ren, H.Z. and Zhang, F.J. (2007) Double Hexagonal Chains with Minimal Total π-Electron Energy. Journal of Mathematical Chemistry, 42, 1041-1056. https://doi.org/10.1007/s10910-006-9159-9
[14] Ren, H.Z. and Zhang, F.J. (2007) Extremal Double Hexagonal Chains with Respectfully Yours, to k-Matchings and k-Independent Sets. Discrete Applied Mathematics, 155, 2269-2281. https://doi.org/10.1016/j.dam.2007.06.003
[15] Deng, H. (2008) The Anti-Forcing Number of Double Hexagonal Chains. MATCH Communications in Mathematical and in Computer Chemistry, 60, 183-192.
[16] Klein, D.J., ?ivkovi?, T.P. and Trinajsti?, N. (1987) Resonance in Random π-Network Polymers. Journal of Mathematical Chemistry, 1, 309-334. https://doi.org/10.1007/BF01179796

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.