Introduction
The understanding of chemical reactivity of buckminsterfullerene (C60),1 and other smaller and larger fullerenes would lead to the synthesis of a large number of fullerene derivatives whose properties and applications have been extensively investigated from many viewpoints.1-5One of the unique aspects of fullerenes in this field is the formation of interfullerene bondings, which makes a rich variety of nanoscale network structures such as dimers, oligomers, and one- and two-dimensional fullerene polymeric materials6-11 In general, it has been suggested that interaction between two fullerene cages can arise through one of the four possible approaches: C-C bond between C60 cages forming a [1+1] dimer (point-point mode); C=C bond forming a [2+2] dimer (side-side mode); forming a [5+5] dimer by face-face mode between two pentagons; and forming a [6+6] dimer by faceface mode between two hexagons.12 In recent years, more attention has been paid to fullerene dimers because their unique physical and chemical properties provide intriguing possibilities as model compounds for nano- and polymer science, and offer potential access to novel molecular electronic devices.13,14 In this respect, a series of fullerene dimers have been produced, such as C120,15 the carbon-bridged dimers: C121,16-19 C140, C131,20 and C122,15,19,21 and the heteroatom bridged dimer: C120O,22-24 which could be used as the basic units of fullerene chain structures. The simplest fullerene dimer C120, (C60)2, has been prepared by solid-state25 and by chemical methods. On the other hand, Osterodt and Vogtle,26 Fabre et al.,27 and Dragoe et al.28 isolated the C122 consisting of C60 fullerenes joined by a (C=C) bridge with sp2-hybridized C atoms which are added across C60 bonds shared by two hexagons (hexagon-hexagon bonds). Forman et al.29 reported the first experimental synthesis and characterization of five [2+2] structural isomers of fullerene dimers C140. The C131 is the first hybrid type of the dumbbell-like fullerene dimer that consists of two different sizes of cages, C60 and C70, with a central atom-bridge.20 Its formation, structure and properties may be more complicated and interesting than those of C121 or C141.16-20 Shvartsburg et al.30 used the chain of C2 units to design dimers of original fullerenes (C60 or C70). Finally, Manaa,31 and Anafcheh and Ghafouri32 proposed carbon and BN hexagons (benzene-like unit) as a building block for connecting between two C60 fullerene cages in order to yield unique electronic properties.
As mentioned above, it has been shown that fullerenes also react with themselves to generate fullerene dimers. Since carbon cages smaller than C60 violate the isolated pentagon rule (IPR), they have higher strain and reactivity due to the adjacent pentagons, thus they are good candidates to form dimers, polymers and solids.33 With this initial thought in mind, in this article we consider small fullerenes such as C24 and C36 as molecular bridges for joining higher fullerenes such as C60 and C80, for the first time, see Figure 1. Then we investigate their electronic and magnetic properties in comparison to those of their parents by calculating HOMO―LUMO gaps, binding energies, and NICS indices. In fact, the synthesis and characterization of such fullerene dimers are of main interest due to their unique structural, magnetic, superconducting and mechanical properties,7,9 which are considerably different from those of other carbon nanostructures such as carbon nanotubes and fullerene cages. Since the considered fullerenes in this study are borderless polycyclic conjugated systems with internal cavities, endohedral 3He NMR chemical shifts have proven to be a useful tool for characterizing them and their derivatives.34-36 Determining the chemical shift of encapsulated 3He nucleus into a fullerene cage and comparing with the 3He chemical shift outside gets a direct measure of the shielding of the magnetic field by the fullerene; such experiments are well known for fullerenes. Providing a good prediction for the endohedral 3He NMR chemical shift, nucleus independent chemical shift (NICS) was proposed by Schleyer et al. in 1996.37 Therefore, in order to probe the local effects of magnetic field inside each fullerene cage we employ grid distribution of NICS inside these molecular clusters and their parents. Magnetic field inside each cage is a consequence of diamagnetism and is related to the induced ring current in the fullerene molecular orbitals, which causes extra stabilization/destabilization in the case of aromatic/ antiaromatic compounds.38 Therefore, it can provide better insights of the electron delocalization, diamagnetic susceptibilities, molecular aromaticity and magnetic properties.
Figure 1.Schematic diagram of fullerene dimers together with the optimized geometries of C60-C36-C60 (C156) and C80-C24-C80 (C184).
Computational Methods. All density functional theory (DFT) quantum calculations are performed using Gaussian 98 program package.39 We consider small fullerenes C24 and C36 as molecular bridges for joining higher fullerenes C60 and C80, see Figure 1. Therefore, the structural geometries of six different configurations of the considered models are allowed to relax by all-atomic optimization. Because of the large sizes of the investigated systems optimization method is qualified step by step as follows: first C24, C36, C60 and C80 fullerene cages considered as the starting points for the design of these compounds are optimized at the B3LYP/6- 31G(d) level of theory.40 In the next step the geometries obtained in step 1 are used to create initial geometries of the fullerene dimers, 𝑖.𝑒., smaller fullerenes, C24 and C36, are located between two cages with approximate interlayer bond length of 1.6 Å (based on reported CC bond lengths for the fullerene dimers15,32); then optimization is first performed with 3-21G basis set for the resulted molecules and finally optimal geometries and normal mode frequencies for all the structures are obtained using standard 6-31G(d) basis set. The coordinates of all the optimized structures can be found in supplementary material. The standard 6-31G(d) basis set is employed due to being affordable and accurate enough for geometry optimization of even large molecules.40,41 Real frequencies obtained from frequency calculations confirm that all of them are minimum energy structures.
Figure 2.Ionization energy (eV) for the fullerene dimers.
As a stability criterion of different configurations, binding energies per atom have been calculated according to the following expression:
where ET is the total energy of the fullerene dimers. Systems with larger binding energies are more stable. To calculate the NICS, ghost atoms are placed along the principal axes of the considered fullerene dimers and their parents with a step size of 0.5Å. The zero point of the coordinate system is positioned at the bridge centers of the optimized structures of the considered fullerene dimers.
Results and Discussion
Geometrical and Stability Properties. We have chosen the C60 and C80 cages as parent molecules, which are fully optimized at the B3LYP/6-31G* level. The obtained structure of C60 is consistent with the literature; the prediction of bond lengths of the hexagon-hexagon (h-h) and hexagonpentagon (h-p) junctions (1.452 and 1.393 Å, respectively) is in excellent agreement with the experimental values (1.458 and 1.401 Å, respectively).42,43 The C36 and C24 fullerene cages are chosen as bridges between two fullerene cages. The structure of C24 can be regarded as a [12]trannulene44 capped with two hexagonal terminal caps. Its optimized structure is indeed affirmative of such a consideration with the uniform C-C bond lengths (1.420 Å) of the hexagonal terminal caps and the much longer C-C bond lengths of 1.530 Å between the hexagonal caps and the central [12]trannulenic ring. Moreover, the localized C=C bond lengths of 1.365 Å and C-C bond lengths of 1.462 Å in the central [12]trannulenic ring are in excellent agreement with the reported experimental values (1.365 and 1.463 Å, respectively). 44,45 The structure of C36 fullerene cage has been described in our previous work.46
In this work, we use the face-to-face linking model to create six types of fullerene dimers with C24 and C36 bridges, see Figure 1. When two neighboring carbon cages share their hexagonal caps with the face-to-face pattern, sp2 hybridization transfers to sp3 hybridization for the carbon atoms in the merged region. As the structures are being reported for the first time, their geometrical characteristics are discussed briefly with the aim of giving better interpretations of these clusters. Based on the optimized structures, the intercage bond lengths of six fullerene dimer are in the range 1.603- 1.622 Å, which can be compared with the intercage bond lengths of 1.60 Å reported by Fowler et al.47 for a linear chain of D6h-C36 cages, and the electron-diffraction pattern of C36-based solid which suggested an intercage distance shorter than 1.7 Å. Moreover, the calculated C-C bond lengths of 1.657-1.660 Å in the two hexagonal terminal caps of C24 bridge cages are slightly longer than the corresponding C-C bond lengths (1.522-1.604 Å) in the C36 bridge cages
Table 1.Total energy (ET), Binding energy (Ebin), and HOMO-LUMO energy gaps (Eg) in the fullerene dimers
To compare the obtained results with those available in the literature, binding energies per atom are calculated for the C36 and C60 fullerenes to be 8.880 and 7.820 eV/atom, which are in agreement with the previously reported values (8.55 and 7.72 eV/atom).2,45 The little difference observed can be due to the different computational methods used. We first note that binding energies for the fullerene dimers with C36 bridges are always lower than those of the fullerene dimers with C24 bridges. Secondly, the replacement of C60 cages with C80 ones (increasing the size of fullerene cage) leads to an increase in binding energy in these systems, see Table 1.
The energy difference between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), Eg, indicates that C36 and C60 are semiconductors with the Eg of 1.02 and 1.66 eV, respectively, which are in agreement with the reported values in the literature (0.8 and 1.8 eV).45,48 As can be seen in Table 1, in contrast to binding energies, the Eg values of the fullerene dimers with C36 bridges are larger than those of the fullerene dimers with C24 bridges. This occurrence may bring about a change in the related electrical conductivity since it is well known that the Eg (or band gap in the bulk materials) is a major factor that determines the electrical conductivity of the material. A classic relation between them is as follows49:
where σ is the electrical conductance and k is the Boltzmann’s constant. According to the equation, smaller Eg leads to higher conductance at a given temperature. Therefore, all of the considered models in this study are semiconductors with Eg values of 1.05-1.26 eV when C36 cage is located between two fullerene cages and 0.78-0.88 eV when C24 is sandwiched between two fullerene cages. Moreover, it is also found that band gaps increase with increasing the size, 𝑖.𝑒., replacing C60 with C80 leads to larger HOMO-LUMO gap for each group of fullerene dimers, see Table 1.
The first ionization potentials (IP) are calculated under the Koopmans’ theorem for closed-shell molecules, based on the frozen orbital approximations and the finite difference approach. In the other words, they are expressed in terms of the highest occupied molecular orbital (HOMO) energies, EHOMO: IP ≈ −EHOMO. At this point it is necessary to mention that the focus is not to find precise ionization potential values; instead, the primary purpose is to study the evolution and the trend of IP in the considered models, and this is just an approximate comparison. Figure 2 depicts the trend of ionization potential for the considered models. The IP plot indicates that the fullerene dimers with C36 bridges have higher ionization potentials compared to the fullerene dimers with C24 bridges, and are thus harder to oxidize (oxidation). Meanwhile, replacement of C60 cages with C80 ones (increasing the size of fullerene cage) leads to an increase in ionization potential in these systems.
NICS Characterization. Aromaticity describes molecules that benefit energetically from the presence of cyclic or spherical electron delocalization in closed circuits of mobile electrons. It is well known that the stability is not directly related to aromatic stabilization but with strain reduction. This is consistent with the high aromaticity of C36 compared to the antiaromaticity of C24, as it is measured by the magnetic aromaticity index of nucleus independent chemical shift (NICS) evaluated at the center of the cages.
Since the aromaticity is not an observable characteristic, there is no magnitude that defines it clearly, and so it is generally evaluated indirectly on the basis of energetic, geometric, or electronic criteria. Especially, it can be followed by obtaining information from the magnetic properties. In fact, the most important methods among several ones to evaluate aromaticity are based on NMR chemical shifts and diamagnetic susceptibilities. Compounds with considerably exalted diamagnetic susceptibility are considered as aromatic structures. The ring currents generated in such molecules by an external magnetic field result in special properties such as “exalted” magnetic susceptibilities and NMR chemical shifts displaced from their normal ranges.37,38 Such particular magnetic influences typically are especially large inside aromatic cyclic or cage-like molecules. To match the familiar NMR convention, NICS indices correspond to the negative of the magnetic shielding, a well-defined property of electronic systems, computed at chosen points designated using the Bq ghost atoms. Significantly negative NICS values in interior positions of cages (magnetically shielded) indicate the presence of induced diatropic ring currents or aromaticity.
On the other hand, antiaromatic cages are identified by their positive NICS values (magnetically deshielded), indicating paratropic ring currents.
Figure 3.Computed NICS values (ppm) along the principal axes of fullerene dimers connected through C24 and C36 bridge cages. The zero point of the coordinate system is positioned at the centers of bridge cages.
As mentioned above, C24 cage can be regarded as a [12]trannulene capped with two benzene rings at both sides. The local ring currents are diatropic within the six membered rings and, in sharp contrast, paratropic within the [12]trannulenic ring. Compensation of these two local effects results in NICS value of 37.89 ppm at the cage center,50 so C24 fullerene is antiaromatic at all. The C36 fullerene cage can be viewed as being made of a zigzag (6, 0) tubular belt, six-membered cyclic polyacene, joined to hexagonal terminal caps, which results in the high aromaticity of the C36 with the calculated NICS value of −26.52,50 in agreement with the experimental value of −28.8 ppm previously reported by Saunders et al..45
Figure 3 depicts variations of NICS values versus distances from the bridge center for the considered compounds and their parent cages (The calculated NICS values for the fullerene dimers have been shown in Table S1 of the Supplementary material). Dimerization of C60 and C80 fullerenes with C24 and C36 significantly change their aromatic characters, leading to a decrease in antiaromaticity of the C24 and C80 with paratropic characters, while an increase in the aromaticity of the aromatic C36 and C60 cages. In other words, NICS values reflect substantial differences in magnetic properties at the cage centers of fullerene dimers. For example, weakly diatrophic (aromatic) C60 has a moderate NICS value of −2.82 ppm,50 while it receives more aromaticity with NICS value of −7.35 ppm in the fullerene dimer C156 (C60-C36-C60). In this compound NICS in aromatic fullerene C36 changes from −26.52 ppm to the value of −29.76 ppm, 𝑖.𝑒., aromaticity increases. The inverse behavior is observed for C24 and C80 which are severely antiaromatic with high positive NICS values of 37.89 and 53.22 ppm, respectively. Compensation between diatropic and paratropic ring currents leads to a decrease in NICS values to 7.16 and 0.93 ppm at the cage centers of C24 and C80, respectively, in the fullerene dimer C184(C80-C24-C80). These trends reveal that fullerene dimerization causes major changes in the magnetic properties at the cage centers.
Now we are interested to make an attempt to correlate the stability of fullerene dimers with their aromaticity character. DFT results indicate that those fullerene dimers involving the initially harshly antiaromatic C24 or C80 cages are more energetically favorable configuration, with the binding energies of 8.14-9.03 eV/atom, than the fullerene dimers involving C36 and C60, with the binding energies of 3.06-3.96 eV/ atom. It is noted that the Ebin of the fullerene dimer C184 (C80- C24-C80) with three antiaromatic cage is larger than those of the other fullerene dimers. Hence, a change in aromaticity character, especially decrease of antiaromaticity, plays a major role in the stability of fullerene dimer. In fact, decrement of antiaromatic character (NICS) inside the joined antiaromatic fullerene cages leads to an increase in the ability of these compounds to sustain an induced ring current, which causes extra stabilization in the case of fullerene dimers of C80 fullerene with C24 bridge cages. Therefore, it seems that the results of this section and those of stability character in previous section mostly support each other.
Conclusion
We have performed a DFT theoretical description to evaluate the electronic and magnetic properties of fullerene dimers of C60 and C80 connected through C24 and C36 bridge cages with the face-to-face linking model. By comparing the results obtained in the present investigation, we emphasize the following points. First, binding energies for the fullerene dimers with C36 bridges are lower than those of the fullerene dimers with C24 bridges. Second, replacement of C60 cages with C80 ones always leads to the increase of binding energy in these systems. Third, HOMO-LUMO gaps, Eg, of the fullerene dimers with C36 bridges are larger than those of the fullerene dimers with C24 bridges and also replacement of C80 cage for C60 leads to larger Eg for the fullerene dimer. Fourth, variations of NICS values versus distances from the bridge center for the considered compounds and the parent cages indicate that dimerization of C60 and C80 fullerenes with C24 and C36 leads to a significant decrease in antiaromaticity of the antiaromatic cages C24 and C80, and an increase in the aromaticity of the aromatic cages C36 and C60. Finally, fullerene dimers involving the initially harshly antiaromatic C24 or C80 cages are more energetically favorable configurations than the fullerene dimers involving C36 and C60.
참고문헌
- Karaulova, E. N.; Bagrii, E. I. Russ. Chem. Rev. 1999, 68, 889. https://doi.org/10.1070/RC1999v068n11ABEH000499
- Dresselhaus, M. S.; Dresselhaus, G.; Eklunf, P. C. Science of Fullerenes and Carbon Nanotubes; Academic Press: New York, 1996.
- Hirsch, A.; Brettreich, M.; Wudl, F. Fullerenes: Chemistry and Reactions; Wiley-VCH: Weinhem, Germany, 2005.
- Chitta, R.; D'Souza, F. J. Mater. Chem. 2008, 18, 1440. https://doi.org/10.1039/b717502g
- Kharisov, B. I.; Kharissova, O. V.; Gomez, M. J.; Mendez, U. O. Ind. Eng. Chem. Res. 2009, 48, 545. https://doi.org/10.1021/ie800602j
- Stephens, P. W.; Bortel, G.; Faigel, G.; Tegze, M.; Janossy, A.; Pekker, S.; Oszalnyi, G.; Forro L. Nature 1994, 370, 636. https://doi.org/10.1038/370636a0
- Oszlanyi, G.; Bortel, G.; Faigel, G.; Granasy, L.; Bendele, G. M.; Stephens, P. W.; Forro, L. Phys. Rev. B 1996, 54, 11849. https://doi.org/10.1103/PhysRevB.54.11849
- Oszlanyi, G.; Baumgartner, G.; Faigel, G.; Forro, L. Phys. Rev. Lett. 1997, 78, 4438. https://doi.org/10.1103/PhysRevLett.78.4438
- Bendele, G. M.; Stephen, P. W.; Prassides, K.; Vavekis, K.; Kordatos, K.; Tanigaki, K. Phys. Rev. Lett. 1998, 80, 736. https://doi.org/10.1103/PhysRevLett.80.736
- Wang, G.-W.; Komatsu, K.; Murata, Y.; Shiro, M. Nature 1997, 387, 583. https://doi.org/10.1038/42439
- Chi, D. H.; Iwasa, Y.; Chen, X. H.; Takenobu, T.; Ito, T.; Mitani, T.; Nishibori, E.; Takata, M.; Sakata, M.; Kubozono, Y. Chem. Phys. Lett. 2002, 359, 177. https://doi.org/10.1016/S0009-2614(02)00297-X
- Ma, F.; Li, Z.-R.; Zhou, Z.-J.; Wu, D.; Li, Y.; Wang, Y.-F.; Li, Z.-S. J. Phys. Chem. C 2010, 114, 11242. https://doi.org/10.1021/jp9116479
- Segura, J. L.; Martin, N. Chem. Soc. Rev. 2000, 29, 13-25. https://doi.org/10.1039/a903716k
- Gao, X.; Zhao, Y.; Yuan, H.; Chen, Z.; Chai, Z. Chem. Phys. Lett. 2006, 418, 24. https://doi.org/10.1016/j.cplett.2005.10.092
- Komatsu, K.; Wang, G.-W.; Murata, Y.; Tanaka, T.; Fujiwara, K.; Yamamoto, K.; Saunders, M. J. Org. Chem. 1998, 63, 9358. https://doi.org/10.1021/jo981319t
- Miao, X.; Ren, T.; Sun, N.; Hu, J.; Zhu, Z.; Shao, Y.; Sun, B.; Zhao, Y.; Li, M. J. Electroanal. Chem. 2009, 629, 152. https://doi.org/10.1016/j.jelechem.2009.02.008
- Dragoe, N.; Shimotani, H.; Wang, J.; Iwaya, M.; De Bettencourt-Dias, A.; Balch, A. L.; Kitazawa, K. J. Am. Chem. Soc. 2001, 123, 1294. https://doi.org/10.1021/ja003350u
- Zhao, Y.; Chen, Z.; Yuan, H.; Gao, X.; Qu, L.; Chai, Z.; Xing, G.; Yoshimoto, S.; Tsutsumi, E.; Itaya, K. J. Am. Chem. Soc. 2004, 126, 11134. https://doi.org/10.1021/ja048232b
- Gao, X.; Yuan, H.; Chen, Z.; Zhao, Y. J. Comput. Chem. 2004, 25, 2023. https://doi.org/10.1002/jcc.20128
- Fowler, P. W.; Mitchell, D.; Taylor, R.; Seifert, G. J. Chem. Soc., Perkin Trans. 1997, 2, 1901.
- Dragoe, N.; Shimotani, H.; Hayashi, M.; Saigo, K.; De Bettencourt-Dias, A.; Balch, A. L.; Miyake, Y.; Achiba, Y.; Kitazawa, K. J. Org. Chem. 2000, 65, 3269. https://doi.org/10.1021/jo991676j
- Lebedkin, S.; Ballenweg, S.; Gross, J.; Taylor, R.; Kratschmer, W. Tetrahedron Lett. 1995, 36, 4971. https://doi.org/10.1016/0040-4039(95)00784-A
- Fujitsuka, M.; Takahashi, H.; Kudo, T.; Tohji, K.; Kasuya, A.; Ito, O. J. Phys. Chem. A 2001, 105, 675. https://doi.org/10.1021/jp002681q
- Balch, A. L.; Costa, D. A.; Fawcett, W. R.; Winkler, K. J. Phys. Chem. 1996, 100, 4823. https://doi.org/10.1021/jp953144m
- Lebedkin, S.; Gromov, A.; Giesa, S.; Gleiter, R.; Renker, B.; Rietschel, H.; Kratschmer, W. Chem. Phys. Lett. 1998, 285, 210. https://doi.org/10.1016/S0009-2614(98)00030-X
- Osterodt, J.; Vogtle, F. Chem. Commun. 1996, 547.
- Fabre, T. S.; Treleaven, W. D.; McCarley, T. D.; Newton, C. L.; Landry, R. M.; Saraiva, M. C.; Strongin, R. M. J. Org. Chem. 1998, 63, 3522. https://doi.org/10.1021/jo980219k
- Dragoe, N.; Tanibayashi, S.; Nakahara, K.; Nakao, S.; Shimotani, H.; Xiao, L.; Kitazawa, K.; Achiba, Y.; Kikuchi, K.; Nojima, K. Chem. Commun. 1999, 85.
- Forman, G. S.; Tagmatarchis, N.; Shinohara, H. J. Am. Chem. Soc. 2002, 124, 178. https://doi.org/10.1021/ja0168662
- Shvartsburg, A. A.; Hudgins, R. R.; Gutierrez, R.; Jungnickel, G.; Frauenheim, T.; Jackson, K. A.; Jarrold, M. F. J. Phys. Chem. A 1999, 103, 5275. https://doi.org/10.1021/jp9906379
- Manaa, M. R. J. Comput. Theor. Nanosci. 2009, 6, 397. https://doi.org/10.1166/jctn.2009.1049
- Anafcheh, M.; Ghafouri, R. Comput. Theor. Chem. 2012, 1000, 85. https://doi.org/10.1016/j.comptc.2012.09.026
- Fowler, P. W.; Mitchell, D.; Zerbetto, F. J. Am. Chem. Soc. 1999, 121, 3218. https://doi.org/10.1021/ja983853o
- Chen, Z.; Cioslowski, J.; Rao, N.; Moncrieff, D.; Buhl, M.; Hirsch, A.; Thiel, W. Theor. Chem. Acc. 2001, 106, 364. https://doi.org/10.1007/s002140100283
- Buehl, M.; Thiel, W.; Jiao, H.; Schleyer, P. V. R.; Saunders, M.; Anet, F. A. L. J. Am. Chem. Soc. 1994, 116, 6005. https://doi.org/10.1021/ja00092a076
- Chen, Z.; King, R. B. Chem. Rev. 2005, 105, 3613. https://doi.org/10.1021/cr0300892
- Schleyer, P. v. R.; Maerker, C.; Dransfeld, A.; Jiao, H.; van Eikema Hommes, N. J. R. J. Am. Chem. Soc. 1996, 118, 6317. https://doi.org/10.1021/ja960582d
- Buhl, M.; Hirsch, A. Chem. Rev. 2001, 101, 1153. https://doi.org/10.1021/cr990332q
- Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Zakrzewski Cheeseman, J. R. V. G.; Montgomery, J. A., Jr.; Stratmann, R. E.; Burant, J. C.; Dapprich, S.; Millam, J. M.; Daniels, A. D.; Kudin, K. N.; Strain, M. C.; Farkas, O.; Tomasi, J.; Barone, V.; Cossi, M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Clifford, S.; Ochterski, J.; Petersson, G. A.; Ayala, P. Y.; Cui, Q.; Morokuma, K.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Cioslowski, J.; Ortiz, J. V.; Baboul, A. G.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Gonzalez, C.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Andres, J. L.; Gonzalez, C.; Head-Gordon, M.; Replogle, E. S.; Pople, J. A., Gaussian, 98, Gaussian. Inc., Pittsburgh, PA, 1998.
- Hariharan, P. C.; Pople, J. A. Mol. Phys. 1974, 27, 209. https://doi.org/10.1080/00268977400100171
- Zhang, Y.; Wu, A.; Xu, X.; Yan, Y. J. Phys. Chem. A 2007, 111, 9431.
- Hale, P. D. J. Am. Chem. Soc. 1986, 108, 6087. https://doi.org/10.1021/ja00279a094
- Scuseria, G. E. Chem. Phys. Lett. 1991, 176, 423. https://doi.org/10.1016/0009-2614(91)90231-W
- Fokin, A. A.; Jiao, H.; Schleyer, P. v. R. J. Am. Chem. Soc. 1998, 120, 9364. https://doi.org/10.1021/ja9810437
- Lu, X.; Chen, Z. Chem. Rev. 2005, 105, 3643. https://doi.org/10.1021/cr030093d
- Anafcheh, M.; Ghafouri, R. Comput. Theor. Chem. 2013, 1017, 1. https://doi.org/10.1016/j.comptc.2013.04.018
- Fowler, P. W.; Heine, T.; Rogers, K. M.; Sandall, J. P. B.; Seifert, G.; Zerbetto, F. Chem. Phys. Lett. 1999, 300, 369. https://doi.org/10.1016/S0009-2614(98)01385-2
- Rivelino, R.; de Brito Mota, F. Nano Lett. 2007, 7, 1526. https://doi.org/10.1021/nl070308p
- Li, S. Semiconductor Physical Electronics, 2nd ed.; Springer, USA, 2006.
- Ghafouri, R.; Anafcheh, M. Physica E 2012, 44, 1386. https://doi.org/10.1016/j.physe.2012.02.023