COMPUTING HYPER WIENER AND SCHULTZ INDICES OF TUZC6[P,Q]
NANOTUBE BY GAP PROGRAM
ALI IRANMANESH*, YASER ALIZADEHa
Department of Mathematics, Tarbiat Modares University P. O. Box: 14115-137, Tehran, Iran
aDepartment of Mathematics, Tarbiat Modares University P. O. Box: 14115-137, Tehran, Iran
In this paper, we give an algorithm for computing the Wiener and Schultz indices of any graph and by this algorithm; we obtain the hyper Wiener and Schultz indices of TUZC6 [p,q] nanotube.
(Received August 30, 2009; accepted September 27, 2009)
Keywords: Hyper Wiener index, Schultz index, Nanotubes, Gap programming.
1. Introduction
Topological indices of nanotubes are numerical descriptors that are derived from graph of chemical compounds. Such indices based on the distances in 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 chemist Harold Wiener [1] introduced Wiener index to demonstrate correlations between physicochemical properties of organic compounds and the index of their molecular graphs. The hyper-Wiener index is one of the recently conceived distance-based graph invariants, used as a structure-descriptor for predicting physico–chemical properties of organic compounds. The hyper-Wiener index was invented by Randic´ (1993) and was eventually extensively studied [2, 3]. 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. The degree of a vertex is the number of vertices joining to i and denoted by . The ( , entry of the adjacency matrix of G is denoted by
( ) i ∈V G ( )
v i i j) ( , )
A i j . The hyper Wiener index of a graph G is denoted by WW (G) and defined as:
1 1 2
( ) ( , ) ( , )
2i j 2i j
WW G d i j d i j
< <
=
∑
+∑
(1) where d(i, j) stands for the distance between the vertices i and j in the graph G.Another topological index is Schultz index and denoted by MTI. This index was introduced by Schultz in 1989, as the molecular topological index [4], and it is defined by:
* Corresponding author: [email protected]
{ , } ( )
( )( ( , ) ( , ) ).
i j V G
MTI v i d i j A i j
⊆
=
∑
+ (2) In this paper, we give an algorithm for computing the hyper-Wiener and Schultz indices of any graph and by this algorithm we obtain the hyper-Wiener and Schultz indices of TUZC6 [p,q]nanotube.
2. An algorithm for computing of the hyper-Wiener and Schultz indices of a graph
In this section, we give an algorithm that enables us to compute the hyper-Wiener and Schultz indices of any graph. For this purpose, the following algorithm is presented:
We assign to any vertex one number. We determine all of adjacent vertices set of the vertex (G) and this set denoted by The set of vertices that their distance to vertex i is equal to is denoted by and consider
V i
i, ∈ N(i ).
(
t≥0t
)
Di, t Di,0 ={i}.Let e=ijbe an edge connecting the vertices i and j, then we have the following result:a) V =Ut≥0Di t, , i∈V (G .)
b) ,
( ) 1
( , ) i t ,
j V G t
d i j t D
∈ ≥
= ×
∑ ∑
2 2 ,( ) 1
( , ) i t , ( )
j V G t
d i j t D i V G
∈ ≥
= × ∀ ∈
∑ ∑
c) 2 ,
, 1
( ) 1 ( )
4i V t i t
W W G t t D
∈ ≥
=
∑
+ ×d)
{ , } ( )
( ) ( )( ( , ) ( , ) )
i j V G
MTI G v i d i j A i j
⊆
=
∑
+( ) ( )
( ) ( ( , ) ( , ))
i V G j V G
v i d i j A i j
∈ ∈
=
∑
×∑
+( ) ( ) ( )\ ( )
( ) 2 ( , )
i V G j N i j V G N i
v i d i j
∈ ∈ ∈
⎛ ⎞
= ×⎜ + ⎟
⎝ ⎠
∑ ∑ ∑
,
2
,
( ) 2
2 ( ) ( )
i t
i t
i V G j D t
v i v i t D
∈ ∈ ≥
⎛ ⎞
= ⎜⎜ + × × ⎟⎟
⎝ ⎠
∑ ∑
According to the above relations, by determining , we can obtain the Hyper- Wiener and Schultz indices of the graph G. In the continue we obtain for each vertex
.
1
, ,t ≥ Dit
, 1
, , t ≥
Dit i
(i) The distance between vertex i and its adjacent vertices is equal to 1, therefore ).
i
N( For each , , ≥1, the distance between each vertex of set ) and the vertex i is equal tot+1, thus we have:
Di,1 = j∈Dit t N(j)\(Di,tUDi,t−1
, 1 , ( ( ) \ ( , , 1) , 1.
i t j Di t i t i t
D + =U ∈ N j D UD − t ≥
According to the above equation we can obtain for each (G). The molecular topological index studied in many papers [5-8]. In a series of papers, the Wiener and Schultz indices of some nanotubes computed [9-13].
, 2
, ,t≥
Dit i∈V
3. Disscussion and conclusions
In this section, according to the above algorithm we write a program by GAP [14] to compute the hyper Wiener and Schultz indices of nanotube.
nanotube is a carbon nanotube formed from a graphite sheet that is rolled up so that it has a zigzag edge.
] ,
6[p q
TUZC TUZC6[p,q]
Fig. 1. NanotubeTUZC6[8,8].
We denote the number of hexagon in the first row by and the number of rows by . In each row, there are vertices and hence the number of vertices in this nanotube is equal to .
p q
p
2 2pq
The following program computes the hyper Wiener and Schultz indices of nanotube for arbitrary p and q.
] ,
6[p q TUZC p:=4; q:=5; #(for example)
n:=2*p*q; N:=[];
K1:=[1..2*p]; V1:=[2..2*p-1];
for i in V1 do
if i mod 2=0 then N[i]:=[i-1,i+1,i+2*p];
else N[i]:=[i-1,i+1]; fi; od;
N[1]:=[2,2*p]; N[2*p]:=[1,2*p-1,4*p];
K:=[2*p+1..n];
K2:=Filtered(K,i->i mod (4*p) in [1..2*p]);
for i in K2 do
if i mod 2 =0 then N[i]:=[i-1,i+1,i+2*p];
else N[i]:=[i-1,i+1,i-2*p]; fi;
if i mod (4*p)=1 then N[i]:=[i+1,i-2*p,i+2*p-1]; fi;
if i mod (4*p)=2*p then N[i]:=[i-1,i-2*p+1,i+2*p];
fi; od;
K3:=Filtered(K,i->i mod (4*p) in Union([2*p+1..4*p-1],[0]));
for i in K3 do
if i mod 2=0 then N[i]:=[i-1,i+1,i-2*p];
else N[i]:=[i-1,i+1,i+2*p]; fi;
if i mod (4*p)=2*p+1 then N[i]:=[i+1,i+2*p-1,i+2*p]; fi;
if i mod (4*p)=0 then N[i]:=[i-1,i-2*p,i-2*p+1]; fi; od;
for i in [n-2*p+1..n] do if q mod 2=1 then
if i mod 2 =0 then N[i]:=[i-1,i+1];
else N[i]:=[i-1,i+1,i-2*p]; fi;
if i mod (4*p)=1 then N[i]:=[i+1,i-2*p,i+2*p-1]; fi;
if i mod (4*p)=2*p then N[i]:=[i-1,i-2*p+1];fi;
else
if i mod 2=0 then N[i]:=[i-1,i+1,i-2*p];
else N[i]:=[i-1,i+1]; fi;
if i mod (4*p)=2*p+1 then N[i]:=[i+1,i+2*p-1]; fi;
if i mod (4*p)=0 then N[i]:=[i-1,i-2*p,i-2*p+1]; fi;
fi; od;
ww:=0; Sc:=0;
v:=[]; D:=[];
for i in [1..n] do D[i]:=[]; u:=[i];
D[i][1]:=N[i]; v[i]:=Size(N[i]);
u:=Union(u,D[i][1]);
ww:=ww+Size(D[i][1]);
Sc:=Sc+v[i]*2*Size(D[i][1]);
r:=1; t:=1;
while r<>0 do D[i][t+1]:=[];
for j in D[i][t] do
for m in Difference (N[j],u) do AddSet(D[i][t+1],m);
od; od;
u:=Union(u,D[i][t+1]);
ww:=ww+(((t+1)^2)+(t+1))*Size(D[i][t+1]);
Sc:=Sc+v[i]*(t+1)*Size(D[i][t+1]);
if D[i][t+1]=[] then r:=0;fi;
t:=t+1; od; od;
ww:=ww/4; #(The hyper Wiener index of the graph) Sc; #(The Schultz index of the graph)
Table 1. The hyper Wiener and Schultz indices of TUZC6 [p,q]
p q Hyper
Wiener
Schultz
3 2 273 822
3 3 900 2352
4 3 2080 4816
5 6 32475 51090
6 4 16296 28356
6 6 53046 79452
7 5 48930 73780
7 8 201488 246330
8 8 289920 340912
References
[1] H. Wiener. J.Am.Chem. Soc. 69, 17 (1947).
[2] M. Randic. Chem. Phys. Lett. 211, 478 (1993).
[3] M. Randi_, X. Gou, T. Oxley, H. Krishnapriyan, and L. Naylor. J. Chem. Inf. Comput. Sci.
34, 361 (1994).
[4] H.P. Schultz. J. Chem. Inf. Comput. Sci. 29, 227 (1989).
[5] H.P. Schultz, J. Chem. Inf. Comput. Sci. 40, 1158 (2000).
[6] Y. Alizadeh, A. Iranmanesh and S. Mirzaie, Digest Journal of Nanomaterials and Biostructures, 4 (2009)
[7] Y. Alizadeh, A. Iranmanesh , Digest Journal of Nanomaterials and Biostructures, 4 (2009) 67.
[8] A. Iranmanesh, O. Khormali, I. Najafi Khalilsaraee, B. Soleiman, Journal of Nanomaterials and Biostructures, 4 (2009) 167.
[9] M. Stefu and M.V Diudea, MATCH Commun. Math Comput. Chem. 50, 133 (2004).
[10] A. Iranmanesh, Y. Alizadeh. Am. J. Applied Sci. 5, 1754 (2008).
[11] M. Randic. Acta Chim. Slov. 49, 483 (2002).
[12] Ali Reza Ashrafi, Shahram Yousefi, Journal of Nanomaterials and
Biostructures, 4 (2009) 407.
[13] Ali Reza Ashrafi, Hosein Shabani, Journal of Nanomaterials and Biostructures, 4 (2009) 453..
[14] Schonert, M. et al. GAP, Groups, Algorithms and Programming, Lehrstuhl De fur Mathematik, RWTH, Achen, (1992).