Graphene clusters consisting of 24 to 150 carbon atoms and hydrogen termination at the zigzag boundary edges have been studied, as well as clusters disordered by vacancy(s).Density Function Theory and Gaussian03 software were used to calculate graphene relative stability, desorption energy, band gap, density of states, surface shape, dipole momentum and electrical polarization of all clusters by applying the hybrid exchange-correlation functional Beke-Lee-Yang-Parr. Furthermore, infrared frequencies were calculated for two of them.Different basis sets, 6-31g**, 6-31g* and 6-31g, depending on the sizes of clusters are considered to compromise the effect of this selection on the calculated results. We found that relative stability and the gap decreases according to the size increase of the graphene cluster.Mulliken charge variation increase with the size. For about 500 carbon atoms, a zero HOMO-LUMO gap amount is predicted. Vacancy generally reduces the stability and having vacancy affects the stability differently according to the location of vacancies. Surface geometry of each cluster depends on the number of vacancies and their locations. The energy gap changes as with the location of vacancies in each cluster. The dipole momentum is dependent on the location of vacancies with respect to one another. The carbon-carbon length changes according to each covalence band distance from the boundary and vacancies. Two basis sets, 2 6-31g* and 6-31g**, predict equal amount for energy, gap and surface structure, but charge distribution results are completely different.
IntroductionA Graphene sheet is a 2-dimensional honey-comb lattice with carbon atoms occupying the corners of hexagons [1], in which each atom is connected to three other carbon atoms. Due to the observed physical [2], electronic [3] and material science properties [4,5] of graphene, this material has become an extremely promising and an absolutely vital element for basic sciences and modern technology. The interest in graphene is wide, from quantum electrodynamics [4,6], condensed matter [7-9] to genomic applications [10,11]. Several research groups have studied various properties and structures of graphene by different methods, like tight binding [12], semi-empirical [13], molecular dynamics [14-16], Hubbard model [17] and Density Function Theories (DFT) [18, 19]. The use of DFT in the ab initio calculation of molecular properties has increased dramatically. This can be attributed to (i) the development of new and more accurate density functions, (ii) the increasing versatility, efficiency, and availability of DFT codes and most importantly, (iii) the superior ratio of accuracy to efforts exhibited by DFT computations relative to other ab initio methodologies [20]. As a result, DFT has been accepted as one of the most accurate methods of calculation for a wide range of studies about graphene [21-23].It is well understood that the edge effect of finite graphene has a crucial role in its electronic properties [24,25] and DFT method has been used to calculat...