Journal of Nanomaterials

Volume 2013, Article ID 451920, 8 pages

http://dx.doi.org/10.1155/2013/451920

## Thermodynamic and Theoretical Study of the Preparation of New Buckyballs from Corannulene, Coronene, and Circulene

^{1}Applied Chemical Division, Applied Science Department, University of Technology, Baghdad, Iraq^{2}Ministry of Sciences and Technology, Industrial Research & Development Directorate, Industrial Applications Center, Baghdad, Iraq^{3}Department of Chemical and Process Engineering, Faculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia

Received 29 July 2013; Revised 16 September 2013; Accepted 7 October 2013

Academic Editor: Mohammad Taghi Ahmadi

Copyright © 2013 Hasan R. Obayes et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We applied density functional theory (DFT) to study three polycyclic aromatic compounds (PAHs), corannulene, coronene, and circulene, for the preparation of twelve new buckyballs with molecular dimensions of less than a nanometer. The results showed that the corannulene molecule is bowl-shaped, the coronene molecule is planar, and the circulene molecule has a unique saddle-shaped structure. Cyclic polymerization of the three molecules can be used to prepare new buckyballs, and this process produces hydrogen molecules. The most symmetric buckyball is also the most stable based on the values of the HOMO energy levels and has the most efficient gap energy, making it potentially useful for solar cell applications.

#### 1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) are a class of unique compounds that consist of fused, conjugated aromatic rings that do not contain heteroatoms or carry substituents [1]. Alternate PAHs contain only six-membered rings, and certain alternant PAHs are called “benzenoid” PAHs. “Small” and “large” PAHs contain up to or more than six fused aromatic rings, respectively. Most research has focused on small PAHs due to their availability [2]. Corannulene has the shape of a bowl because it contains a five-membered ring, which inverts rapidly. In addition to its nonstandard geometry and dynamic behavior, this molecule has attracted considerable interest as an important building block for the synthesis of C_{53}. Additionally, the solid state packing behavior of corannulene is interesting [3]. As discussed in detail by Kawasa and Kurata, [4] not only bowl-shaped but also ball- and belt-shaped aromatic systems provide an exciting opportunity to study concave-convex interactions. Coronene (6-circulene) is an aromatic planar symmetric molecule, which has been studied, synthesized, and well characterized [5]. Both corannulene and coronene have interesting conductive properties due to their large electronic resonance. Circulene, the next member of this molecular family, consists of a central cyclooctatetraene fragment surrounded by phenyl rings. Despite the efforts of researchers [6], circulene has not yet been synthesized. Several theoretical studies [7], but no detailed studies, have been carried out on the structure of this molecule. Theoretical predictions of the existence of buckyball molecules appeared in the late 1960s and early 1970s, but these predictions were largely unnoticed. The discovery of buckyballs was unexpected because the scientists were producing carbon plasmas to replicate and characterize unidentified interstellar matter. Mass spectrometry analysis of their product indicated the formation of spheroidal carbon molecules [8].

#### 2. The Calculation Method

To calculate ground-state geometries, Gaussian 03, Revision C.01 [9] was optimized to a local minimum without symmetry restrictions using basis set 6-31 G [10, 11]. The combination of the Becke three-parameter hybrid (B3) [12, 13] exchange functional and the Lee-Yang-Parr (LYP) [14] correlation functional (B3LYP) [15, 16] was used for all geometry optimizations, thermodynamic functions (at 298.150 Kelvin and 1.0 atm), Highest Occupied Molecular Orbital Energies (E_{HOMO}), Lowest Unoccupied Molecular Orbital Energies (E_{LUMO}), and physical properties for the molecules in this study.

#### 3. Results and Discussion

Previous studies have shown that not all polycyclic aromatic hydrocarbons (PAHs) are flat molecules. We selected aromatic compounds known as circulenes for this study. The circulenes include 5-circulene (corannulene), 6-circulene (coronene), and 7-circulene (circulene). According to Density Function Theory (DFT) calculations, not all of these molecules are flat, as shown in Figure 1. Corannulene is bowl-shaped, coronene is planar, and circulene has a unique saddle-shaped structure, which are consistent with the literature.

##### 3.1. Density Function Theory (DFT)

A DFT calculation introduces an additional step to each major phase of a Hartree-Fock calculation. This additional step is a numerical integration of the functional (or various derivatives of the functional). Thus, in addition to the sources of numerical error in Hartree-Fock calculations (integral accuracy, SCF convergence, and CPHF convergence), the accuracy of DFT calculations also depends on the number of points used in the numerical integration. The “fine” integration grid is the default in Gaussian 03. This grid greatly enhances the calculation accuracy at minimal additional cost. We do not recommend using a smaller grid in production DFT calculations. It is important to use the same grid for all calculations when energies will be compared (e.g., computing energy differences, heats of formation, and so forth). Larger grids are available, for example, for tight optimization of certain systems. An alternate grid may be selected in the route section [17–20].

To prepare new buckyballs, a fundamental understanding of the polymerization of polycyclic aromatic hydrocarbons (PAHs) is important. The polymerization process for new buckyballs revealed the production of hydrogen molecules [21, 22], and the general reaction of formation of new buckyballs from three molecules of polycyclic aromatic hydrocarbons (PAHs) is as follows, where is number of carbon atoms, is , is 2 in corannulene and circulene, and is 0 in coronene.

###### 3.1.1. The Cyclic Polymerization of Three Corannulene Molecules

Scheme 1 shows all possible cyclic polymerization reactions of three molecules of corannulene. Reaction (1) produced a new buckyball C_{60}H_{2} (I) by forming five butagons, three pentagons, three hexagons, and two decagon cycles from three corannulene molecules. Reaction (2) produced a new buckyball C_{60}H_{2} (II) through the formation of three butagons, five pentagons, three hexagons, and two nonagon cycles. Reaction (3) produced a new buckyball C_{60}H_{2} (III) by forming three butagons, three pentagons, five hexagons, and two octagon cycles. Reaction (4) produced a new buckyball C_{60}H_{2} (IV) by forming eleven pentagons and two nonagon cycles. Thirteen different cycles were formed in these four reactions. All of these reactions are spontaneous and exothermic according to the values of the entropy change (), the Gibbs energy change (), and the enthalpy change (). The results for the E_{HOMO} (the Energy of High Occupied Molecular Orbital) and the total energy for the four reaction products in Table 1 reveal that the products are stable, that the new buckyball C_{60}H_{2} (IV) is the most stable among the four, and that the increase in the E_{HOMO} for C_{60}H_{2} (IV) is (−0.4530 eV), (−0.6177 eV), and (−0.6236 eV) relative to C_{60}H_{2} (I), C_{60}H_{2} (II), and C_{60}H_{2} (III) respectively. Additionally, the increase in the total energy is −0.1978 a.u or −124.121 KCalmol^{−1}, −0.1365 a.u or −85.655 KCalmol^{−1}, and −0.2260 a.u or −141.817 KCalmol^{−1} for C_{60}H_{2} (I), C_{60}H_{2} (II), and C_{60}H_{2} (III), respectively. The structures of the four new buckyballs are shown in Figure 2.

The molecular dimensions of the four new buckyballs are as follows: C_{60}H_{2} (I)—(0.85 -axis, 0.75 -axis, and 0.52 -axis) nm, C_{60}H_{2} (II)—(0.83 -axis, 0.73 -axis, and 0.49 -axis) nm, C_{60}H_{2} (III)—(0.84 -axis, 0.72 -axis, and 0.48 -axis) nm, and C_{60}H_{2} (IV)—(0.84 -axis, 0.73 -axis, and 0.48 -axis) nm.

###### 3.1.2. The Cyclic Polymerization of Three Molecules of Coronene

Scheme 2 shows all possible cyclic polymerization reactions of three molecules of coronene. Reaction (1) produced a new buckyball C_{72} (I) and forms nine butagons, six hexagons, and two nonagon cycles from the three coronene molecules. Reaction (2) produced a new buckyball C_{72} (II) by forming three butagons, ten pentagons, two hexagons, and two octagon cycles. Reaction (3) produced a new buckyball C_{72} (III) and formed two butagons, ten pentagons, three hexagons, and two heptagon cycles. Reaction (4) produced a new buckyball C_{72} (IV) and formed six butagons, and eleven hexagon cycles. Seventeen cycles were formed in each reaction. Every reaction is spontaneous and exothermic according to the values of , , and . The E_{HOMO} and the total energy for the four reaction products in Table 1 reveal that the products are stable, that the new buckyball C_{72} (IV) is the most stable among the four, and that the increase in E_{HOMO} for C_{72} (IV) is (−0.7102 eV), (−0.4599 eV), and (−0.6122 eV) relative to new buckyballs C_{72} (I), C_{72} (II) and C_{72} (III), respectively. Additionally, the increases in total energy are −0.4264 a.u or −267.570 KCalmol^{−1}, −0.0619 a.u or −38.843 KCalmol^{−1}, and −0.2572 a.u or −161.395 KCalmol^{−1} for new buckyballs C_{72} (I), C_{72} (II), and C_{72} (III), respectively. The structures of the four new buckyballs are shown in Figure 3.

The molecular dimensions of all four new C_{72} buckyballs are as follows: C_{72} (I)—(0.82 -axis, 0.80 -axis, and 0.52 -axis) nm, C_{72} (II)—(0.83 -axis, 0.78 -axis, and 0.55 -axis) nm, C_{72} (III)—(0.88 -axis, 0.77 -axis, and 0.55 -axis) nm, and C_{72} (IV)—(0.82 -axis, 0.79 -axis, and 0.56 -axis) nm.

###### 3.1.3. The Cyclic Polymerization for Three Molecules of Circulene

Scheme 3 shows all possible cyclic polymerization reactions of three molecules of circulene. Reaction (1) produced a new buckyball C_{84}H_{2} (I) by forming seven butagons, five pentagons, five hexagons, and two decagon cycles from three circulene molecules. Reaction (2) produced a new buckyball C_{84}H_{2} (II) by forming five butagons, seven pentagons, five hexagons, and two nonagon cycles. Reaction (3) produced a new buckyball C_{84}H_{2} (III) and formed five butagons, five pentagons, seven hexagons, and two octagon cycles. Reaction (4) produced a new buckyball C_{84}H_{2} (IV), seventeen pentagons, and two nonagon cycles. Nineteen different cycles were formed in each reaction. All of these reactions are spontaneous and exothermic according to the values of the change of entropy (), the change of Gibbs energy (), and the change of enthalpy (). However, the the Energy of the Highest Occupied Molecular Orbital (E_{HOMO}) and the total energy for the reaction products in Table 1 reveal that the product has more stability, in this state the product new buckyball C_{84}H_{2} (IV), most stable among the four new buckyballs, was C_{84}H_{2} (IV) with an increase in the E_{HOMO}, that is (−0.5197 eV), (−0.4968 eV), and (−0.6955 eV) relative to C_{84}H_{2} (I), C_{84}H_{2} (II), and C_{84}H_{2} (III), respectively. Additionally, the increases in total energy are −0.4272 a.u or −268.072 KCalmol^{−1}, −0.2805 a.u or −176.016 KCalmol^{−1}, and −0.2904 a.u or −182.228 KCalmol^{−1} relative to C_{84}H_{2} (I), C_{84}H_{2} (II) and C_{84}H_{2} (III), respectively. The structures of the four new C_{84}H_{2} buckyballs are shown in Figure 4. The molecular dimensions of all four new C_{84}H_{2} buckyballs are as follows: C_{84}H_{2} (I)—(0.86 -axis, 0.85 -axis, and 0.66 -axis) nm, C_{84}H_{2} (II)—(0.90 -axis, 0.88 -axis, and 0.61 -axis) nm, C_{84}H_{2} (III)—(0.89 -axis, 0.87 -axis, and 0.60 -axis) nm, and C_{84}H_{2} (IV)—(0.90 -axis, 0.87 -axis, and 0.60 -axis) nm.

##### 3.2. Energy Gap

The energy gap, which is also called the band gap, is an energy range in a solid where no electron states can exist. The gap energy generally refers to the energy difference (in electron volts) between the Lowest Unoccupied Molecular Orbital (LUMO) and the Highest Occupied Molecular Orbital (HOMO) in insulators and semiconductors. This gap energy is equivalent to the energy required to free an outer shell electron from its orbit about the nucleus to become a mobile charge carrier that moves freely within the solid material. The band gap is a major factor that determines the electrical conductivity of a solid. Substances with large gap energies are generally insulators, materials with smaller gap energies are semiconductors, and conducting materials have very small or no gaps energies. The Shockley-Queisser limit gives the maximum possible efficiency of single junction solar cells under unconcentrated sunlight as a function of the semiconductor band gap. If the band gap is too high, then the material cannot absorb most daylight photons; if the band gap is too low, then most photons have much more energy than is necessary to excite electrons across the band gap, and the rest is wasted. The semiconductors that are used commonly in commercial solar cells have band gaps near the peak of this curve shown in Figure 5.

In Table 1, the values of the energy gaps for all buckyballs in the (0.8335–2.3812) eV range, arranged by the increases in energy gap, are as follows:

#### 4. Conclusions

A quantum chemistry calculation is performed using the Density Function Theory (DFT) method to study the preparation of twelve new buckyballs from the cyclic polymerization of three polycyclic aromatic hydrocarbons (PAHs), corannulene, coronene, and circulene, and the production of hydrogen molecules. The results obtained for the new buckyballs show that the most symmetric buckyball is the most stable, depending on the values of E_{HOMO}. The molecular dimensions of all the new buckyballs are less than a nanometer, and the new buckyballs are characterized by the high efficiency of their energy gaps.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### References

- J. C. Fetzer,
*The Chemistry and Analysisof the Large Polycyclic Aromatic Hydrocarbons*, John Wiley, New York, NY, USA, 2000. - X. Feng, W. Pisula, and K. Müllen, “Large polycyclic aromatic hydrocarbons: synthesis and discotic organization,”
*Pure and Applied Chemistry*, vol. 81, no. 12, pp. 2203–2224, 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. T. Wu and J. S. Siegel, “Aromatic molecular-bowl hydrocarbons: synthetic derivatives, their structures, and physical properties,”
*Chemical Reviews*, vol. 106, pp. 4843–4867, 2006. View at Publisher · View at Google Scholar - T. Kawasa and H. Kurata, “Ball-, bowl-, and belt-shaped conjugated systems and their complexing abilities: exploration of the concave-Convex $\pi $-$\pi $ interaction,”
*Chemical Reviews*, vol. 106, pp. 5250–5273, 2006. View at Publisher · View at Google Scholar - J. C. Hanson and C. E. Nordman, “The crystal and molecular structure of corannulene, C
_{20}H_{10},”*Acta Crystallographica*, vol. B32, p. 1147, 1976. View at Google Scholar - K. Yano, M. Osatani, and K. Tani, “On atempted synthesis of 8-shaped circulene: synthesis and characterization of [2,2](4,${4}^{\prime}$)-biphenylo-(3,11)-dihenzo[c.l]chrysenophane and its cyclophadiene,”
*Bulletin of the Chemical Society of Japan*, vol. 73, pp. 185–189, 2000. View at Publisher · View at Google Scholar - T. Liljefors and O. Wennerstrom, “Molecular mechanics calculations on the structures and conformational properties of [8]circulene and some related phane compounds,”
*Tetrahedron*, vol. 33, no. 22, pp. 2999–3003, 1977. View at Google Scholar · View at Scopus - E. A. Katz, “Fullerene thin films as photovoltaic material,” in
*Nanostructured Materials for Solar Energy Conversion*, pp. 361–443, Elsevier, 2006. View at Google Scholar - J. A. Pople, M. J. Frisch, G. W. Trucks et al., “Gaussian Inc., 2004 Gaussian 03W (Revision C.01),” Gaussian Inc., Wallingford, Conn, USA, 2003.
- W. J. Pietro, M. M. Francl, W. J. Hehre, D. J. Defrees, J. A. Pople, and J. S. Binkley, “Self-consistent molecular orbital methods. 24. Supplemented small split-valence basis sets for second-row elements,”
*Journal of the American Chemical Society*, vol. 104, no. 19, pp. 5039–5048, 1982. View at Google Scholar · View at Scopus - K. D. Dobbs and W. J. Hehre, “Molecular orbital theory of the properties of inorganic and organometallic compounds 5. Extended basis sets for first-row transition metals,”
*Journal of Computational Chemistry*, pp. 861–879, 1987. View at Publisher · View at Google Scholar - A. D. Becke, “Density-functional exchange-energy approximation with correct asymptotic behavior,”
*Physical Review*, vol. 38, pp. 3098–3100, 1988. View at Publisher · View at Google Scholar - A. D. Becke, “Density-functional thermochemistry—III. The role of exact exchange,”
*The Journal of Chemical Physics*, vol. 98, no. 7, pp. 5648–5652, 1993. View at Google Scholar · View at Scopus - C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,”
*Physical Review*, vol. 37, pp. 785–789, 1988. View at Publisher · View at Google Scholar - A. A. Al-Amiery, H. D. Jaffar, H. R. Obayes et al., “Thermodynamic studies on 4-aminocoumarin tautomers,”
*International Journal of Electrochemical Science*, vol. 7, pp. 8468–8472, 2012. View at Google Scholar - J. H. Naama, G. H. Alwan, H. R. Obayes et al., “Curcuminoids as antioxidants and theoretical study of stability of curcumin isomers in gaseous state,”
*Research on Chemical Intermediates*, vol. 39, no. 9, pp. 4047–4059, 2013. View at Google Scholar - A. D. Becke, “Density-functional thermochemistry—V. Systematic optimization of exchange-correlation functionals,”
*Journal of Chemical Physics*, vol. 107, no. 20, pp. 8554–8560, 1997. View at Google Scholar · View at Scopus - S. W. Benson,
*Thermochemical Kinetics*, Wiley and Sons, New York, 1968. - R. Berger and M. Klessinger, “Algorithms for exact counting of energy levels of spectroscopic transitions at different temperatures,”
*Journal of Computational Chemistry*, vol. 18, no. 10, pp. 1312–1319, 1997. View at Google Scholar · View at Scopus - G. Zheng, H. A. Witek, P. Bobadova-Parvanova et al., “Parameter calibration of transition-metal elements for the spin-polartaed self-consistent-charge density-functional tight-binding (DFTB) Method: Sc, Ti, Fe, Co, and Ni,”
*Journal of Chemical Theory and Computation*, vol. 3, no. 4, pp. 1349–1367, 2007. View at Publisher · View at Google Scholar · View at Scopus - H. R. Obayes, “Theoretical study for the preparation of new buckyballs from the cyclic dimerization reaction of two molecules corannulene, coronene and circulene aromatic compounds,”
*Journal of Engineering & Technology*, vol. 31, no. 5, pp. 594–602, 2013. View at Google Scholar - H. R. Obayes, A. A. Al-Amiery, H. D. Jaffar et al., “Theoretical study for the preparation of sub-carbon nano tubes from the cyclic polymerization reaction of two molecules from corannulene, coronene and circulene aromatic compounds,”
*Journal of Computational and Theoretical Nanoscience*, vol. 10, pp. 2459–2463, 2013. View at Google Scholar