Computing the Szeged and PI Indices of VC5C7[p,q] and HC5C7[p,q] Nanotubes

In this paper we give a GAP program for computing the Szeged and the PI indices of any graph. Also we compute the Szeged and PI indices of VC5C7 [ p,q] and HC5C7 [ p,q] nanotubes by this program.


Introduction
One of the main distinctive characteristics of modern chemistry is the use of theoretical tools for the molecular modeling of physicochemical processes, chemical reaction, medicinal and toxicological events, etc., in which chemicals are involved. The success of the molecular modeling is judged by the insights that it offers on the nature of the processes studied, which permit better comprehension and a their rational modification. These properties, measured experimentally, are almost invariably expressed in quantitative terms, for instance boiling point, refraction index, transition state energy, percentage of inhibition of some enzymatic activity, lethal dose, and so forth. The paradigm for the modeling of such properties is the relationship that exists between them and the molecular structure of chemical. This fact highlights the first challenge for molecular modeling: the properties are expressed as number while the molecular structure is not. The way to solve this problem is by using molecular descriptors, which are numbers representing information about different molecular features, to describe quantitatively the properties under study. These models are known as quantitative structureproperty (QSPR) and quantitative structure-activity relationships (QSAR), depending on the physicochemical or biological nature of the properties studied, respectively.
Topological indices of nanotubes are numerical descriptors that are derived from graphs of chemical compounds. Such indices based on the distances in a graph are widely used for establishing relationships between the structure of nanotubes and their physicochemical properties. Usage of topological indices in biology and chemistry began in 1947, when the chemist Harold Wiener [1] introduced the so-called Wiener index to demonstrate correlations between physicochemical properties of organic compounds and the index of their molecular graphs. Wiener originally defined his index (W) on trees and studied its use for correlations of physicochemical properties of alkanes, alcohols, amines and their analogous compounds [2]. Starting from the middle of the 1970s, the Wiener index gained much popularity and, since then, new results related to it are constantly being reported. For a review, historical details and further bibliography on the chemical applications of the Wiener index see [3][4][5]. Another topological index was introduced by Gutman and called the Szeged index, abbreviated as Sz [2]. Let G be a connected graph. The vertex-set and edge-set of G denoted by V(G) and E(G) respectively. The distance between the vertices u and v, d(u,v), in a graph is the number of edges in a shortest path connecting them. Two graph vertices are adjacent if they are joined by a graph edge. Let e be an edge of a graph G connecting the vertices u and v. Two sets 1 ( N e G ) and 2 ( N e G) are defined as follows: The number of elements of 1 ( N e G ) and 2 ( N e G) are denoted by 1 ( n e G ) and 2 ( ) n e G respectively.
The Szeged index of the graph G is defined as For the reason of the coincidence of Wiener and Szeged indices in case of trees the authors in [6] and [7] introduced another Szeged/Wiener-like topological index and named it the Padmakar-Ivan index, abbreviated as PI. In fact the PI index of the graph G is denoted by PI(G) and defined as follows: Applications of the PI index to QSRP/QSAR were studied in [8]. The index was mostly compared with the Wiener and the Szeged indices. It turned out that the PI index has similar discriminating power as the other two indices and in many cases (for instance to model ¢max, the so called difference in doublet of deformation mode, of unbranched cycloalkanes) it gives better results. As we already mentioned, The Szeged index incorporates the distribution of vertices of a molecular graph, while the PI index does this job for the edges. Hence it seems that a combination of both could give good results in QSRP/QSAR studies. Indeed, the combination of the PI index and the Szeged index is the best for modeling polychlorinated biphenyls (PCBs) in environment among the three possible pairs of indices selected from the PI index, the Szeged index, and the Wiener index [8]. For the Wiener and the Szeged indices such studies were previously done in [9,10]. The Szeged and PI indices of some nanotubes were computed in [11][12][13][14].
The computation of Szeged and PI indices seems straightforward, but this is not entirely true. For computing the Szeged (or PI ) index of any graph, we must obtain 1 ( n e G ) and 2 ( n e G ) for any edge in the graph and this takes very long time. Up to now, many papers have been published in the international scientific literature which computed the Szeged index (or PI index) of some nanotubes by the above explanation [11][12][13][14][15][16][17][18][19][20][21][22]. In our paper, we obtain an algorithm which is faster than the direct implementation. Also, since up to now, not give any algorithm for computing these indexes, it is the best algorithm by GAP program and therefore, our paper is the first paper in this subject which give an algorithm for computation of Szeged and PI indices of any graph.

An algorithm for the computation of the Szeged and PI indices for an arbitrary graph.
In this section, we give an algorithm that enables us to compute the Szeged and PI indices of any graph. For this purpose, the following algorithm is presented: 1-We assign one number to any vertex.

2-We determine all of adjacent vertices set of the vertex
(G) and this set is denoted by The set of vertices that their distance to vertex i is equal to is denoted by and considering Let be an edge connecting the vertices and According to the above relations, by determining , we can obtain 1 , for each edge e and therefore the Szeged and PI indices of the graph G is computed. In the following section we obtain the for each vertex i . , 1 , , ≥ t D t i 3-The distance between vertex i and its adjacent vertices is equal to 1, therefore For each , the distance between each vertex of set and the vertex i is equal to , thus we have According to the above equation we can obtain for each 2  [23]. The name was chosen to reflect the aim of the system, which is group theoretical software for solving computational problems in group theory. The last years have seen a rapid spread of interest in the understanding, design and even implementation of group theoretical algorithms. GAP software was constructed by GAP's team in Aachen. We encourage the reader to consult Refs. [24] and [25] for background materials and computational techniques related to applications of GAP in solving some problems in chemistry and

Example
The Wiener index of tree dendrimers is computed in [26,27]. Since the Wiener and Szeged index coincide on trees [28,29], thus the Szeged index of T is equal to its Wiener index.
, , 1, The following results are obtained in [26,27]. For every the tree has order 3, and its Szeged index is equal to Wiener index , i.e. 2 3 For computation of the Szeged and PI indices of by the above program, at first we assign to any vertex one number (See figure 1); according to this numbering, the set of adjacent vertices to each vertex, 1≤i≤n, is obtained by the following program (part 1). In fact part 1 of the program is the presentation of the graph. We use part 2 for compute the Szeged and PI indices of the graph.  Add(T,i); od; sz;# (The value of sz is equal to Szeged index of the graph) pi; # (The value of pi is equal to PI index of the graph)

Computation of the Szeged and PI indices of VC 5 C 7 [p,q] nanotubes with the GAP program
A C C net is a trivalent decoration made by alternating C and C . It can cover either a cylinder or a torus. In this section we compute the Szeged and PI indices of VC C [p,q] nanotubes by GAP program  nanotube.
We denote the number of pentagons in the first row by p, in this nanotube the four first rows of vertices and edges are repeated alternatively, we denote the number of this repetition by . In each period there are vertices and vertices which are joined to the end of the graph and hence the number of vertices in this nanotube is equal to We partition the vertices of the graph to the following sets: 1 K : The vertices of the first row whose number is p 6 .
2 K : The vertices of the first row in each period except the first one whose number is ) 1 The vertices of the second rows in each period whose number is pq 2 .
4 K : The vertices of the third row in each period whose number is pq 6 .

5
K : The vertices of the fourth row in each period whose number is pq 2 .
6 K : The last vertices of the graph whose number is . 3 p Figure 3. Rows of m-th period.
We write a program to obtain the adjacent vertices set to each vertex in the sets i=1…6. We can obtain the adjacent vertices set to each vertex by the join of these programs. In this program, the value of x is the assign number of vertex in that row.  Add(T,i); od; sz; # (The value of sz is equal to Szeged index of the graph ) pi; # (The value of pi is equal to PI index of the graph ) 5 7 nanotubes with the GAP program.

Computation of the Szeged and PI indices of HC C [p,q]
In this section we compute the Szeged and PI indices of HC 5 C 7 [p,q] nanotubes similar to the previous section. HC 5 C 7 [p,q] nanotubes consists of heptagon and pentagon nets as seen below: We denote the number of heptagons in the first row by p. In this nanotube the four first rows of vertices and edges are repeated alternatively; we denote the number of this repetition by q. In each period there are 16p vertices and 2p vertices are joined to the end of the graph, and hence the number of vertices in this nanotube is equal to 16pq + 2p.
The following program is the same as the last program. In this program, value of x is the number of vertex in a row.  Add(T,i); od; sz;# (The value of sz is equal to Szeged index of the graph) pi;# (The value of pi is equal to PI index of the graph)

Conclusions
It takes a long time to compute the Szeged and PI indices of a graph theoretically, especially when the graph has many vertices. According to the algorithm presented in this paper and using the GAP program, we can write a program to compute these indices quickly. This method has been used here for the first time. We teste the algorithm to calculate the Szeged and PI indices of denderimers T k,d that were computed in [26,27]. In addition we compute these indices for and