Advances in Physical Chemistry

Volume 2015, Article ID 506894, 6 pages

http://dx.doi.org/10.1155/2015/506894

## Geometry, Energy, and Some Electronic Properties of Carbon Polyprismanes: *Ab Initio* and Tight-Binding Study

^{1}Department of Condensed Matter Physics, National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Kashirskoe Shosse 31, Moscow 115409, Russia^{2}Laboratory of Computational Design of Nanostructures, Nanodevices and Nanotechnologies, Research Institute for the Development of Scientific and Educational Potential of Youth, Aviatorov Street 14/55, Moscow 119620, Russia

Received 30 July 2015; Revised 7 October 2015; Accepted 11 October 2015

Academic Editor: Miquel Solà

Copyright © 2015 Konstantin P. Katin 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 report geometry, energy, and some electronic properties of [*n*,4]- and [*n*,5]prismanes (polyprismanes): a special type of carbon nanotubes constructed from dehydrogenated cycloalkane C_{4}- and C_{5}-rings, respectively. Binding energies, interatomic bonds, and the energy gaps between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) have been calculated using density functional approach and nonorthogonal tight-binding model for the systems up to thirty layers. It is found that polyprismanes become more thermodynamically stable as their effective length increases. Moreover, they may possess semiconducting properties in the bulk limit.

#### 1. Introduction

Carbon [*n*,*m*]prismanes (polyprismanes) can be regarded as stacked layers of dehydrogenated cycloalkane molecules, where is the number of vertices of the closed carbon ring and is the number of layers [1]. However, the polyprismanes’ ends are passivated by hydrogen atoms to avoid the dangling bonds. For large , polyprismanes (PP) represent single-walled carbon nanotube analogs with an extremely small cross-section as a regular polygon. Despite the fact that PP are of considerable fundamental and practical interest, their synthesis is very difficult task due to their high-strained carbon frame. Unlike carbon nanotubes, polyprismanes have rectangles on their surfaces instead of hexagons. So the angles between covalent C–C bonds are different from the value 109.5° usual for the sp^{3}-hybridized carbon orbitals. Only “simplest” PP representatives are obtained: [2,3]prismane C_{6}H_{6} (Ladenburg’s prismane) [2], [2,4]prismane C_{8}H_{8} (hydrocarbon cubane) [3], and [2,5]prismane C_{10}H_{10} [4]. Attempts to synthesize [2,*m*]prismanes starting from [2,6]prismane C_{12}H_{12} failed. Nevertheless, hypothetical ways of synthesis based on high-level density functional calculations are predicted for [3,4]prismanes and [4,4]prismanes [1]. Moreover, structural characteristics and electronic properties of [*n*,*m*]prismanes with , [5]; , [6]; , [7, 8] are actively studied in terms of the density functional theory. The evolution of aromaticity along the thermally forbidden and dimerizations of cyclobutadiene and benzene to give [2,4]prismane and [2,6]prismane, respectively, is also discussed in the literature [9]. The special place is occupied by investigations of mechanical properties of polyprismanes, including the analysis of the edge doping (by replacing hydrogen atoms with chemical functional groups or radicals) influence on the PP mechanical characteristics [10, 11]. It is established that polyprismanes are molecular auxetics, that is, compounds having negative Poisson’s ratio [12]. In addition, high kinetic stability of small prismanes (namely, [2,6]prismanes and [2,8]prismanes) is also confirmed via direct molecular dynamics simulations [13]. Long [*n*,*m*]prismanes with the large values of are of special interest. They are considered as molecular wires for nanoelectronic applications [14]. In addition, recent* ab initio* calculations predicted that “long” (up to ) and even infinite () [*n*,4]prismanes were stable and could exhibit the metallic behavior [15]. Moreover, the particular electronic properties of long PP make them potentially suitable in the biomedical sector as biological ion carriers [16–18]. So they can be used as drug delivery systems in the field of nanomedicine. But despite series of evident achievements in physical chemistry of long PP, their geometry, energy, and electronic properties are not well studied.

Nevertheless, PP are quite promising compounds. For example, polyprismanes can be functional nanowires with controlled electronic characteristics and can be used in computational logic elements. The incorporation of metastable carbon-nitrogen structures (nitrogen chains) inside the higher PP can stabilize the latter, which will make it possible to obtain an efficient high-energy material. The same principle of the formation of endohedral complexes (the encapsulation of therapeutic preparations inside polyprismanes) can be used for developing drug delivery systems. The use of PP as elements of measuring equipment, such as probes (tips) of the atomic-force microscopes, is also urgent. The absence of free covalent bonds makes polyprismanes insensitive to contaminants of the external medium (free radicals), which lowers the reaction ability of the tip making it possible to attain atomic resolutions in a microscope. All these possible applications pose quite concrete questions for the researchers: what are the limiting operating temperature modes of materials based on PP; how can we estimate lifetimes of the samples in extreme conditions (e.g., at cryogenic temperatures); how will electronic characteristics vary (e.g., conductivity) with an increase in the efficient sample length; how will edge doping affect their mechanical properties; and finally what compounds from the PP family are most thermodynamically and kinetically stable and the most probable candidates for the practical use?

In this paper, we present density functional theory and tight-binding calculations for the [*n*,4]- and [*n*,5]prismanes’ family with . The main purpose of this work is to obtain their binding energies, interatomic bonds, and HOMO-LUMO gaps and to analyze the behavior of these values in the bulk limit (). We hope that our study will stimulate the subsequent theoretical and experimental research of PP.

#### 2. Materials and Methods

In our study we use density functional theory as well as the nonorthogonal tight-binding model (NTBM) [20]. Density functional theory (DFT) with Becke’s three-parameter hybrid method and the Lee-Yang-Parr exchange-correlation energy functional (B3LYP) [21, 22] with the electron basis set of 6-311G(d) [23] were used to optimize the geometries and obtain the structural, energy, and electronic (namely, HOMO-LUMO gaps) characteristics of the polyprismanes. The geometries of all polyprismanes were optimized using the GAMESS program package [24]. In the tight-binding model, the total potential energy of the system is the sum of the quantum-mechanical electronic energy and the “classical” ionic repulsive energy. The former is defined as the sum of one-electron energies over the occupied states, the energy spectrum being determined from the stationary Schrödinger equation (2S, 2P_{x}, 2P_{y}, and 2P_{z} orbitals of carbon atoms and 1S orbitals of hydrogen atoms are considered). The electronic Hamiltonian is calculated from overlap matrix elements using extended Hückel approximation [25] with Anderson distant dependence for Wolfsberg-Helmholtz parameter [26]. The overlap integrals are calculated using standard Slater-Koster-Roothaan procedure [27, 28]. Therefore, the analytical expressions for the overlap integrals greatly simplify the calculation of Hamiltonian. The parameters’ fitting of tight-binding model is based on the criterion of the best correspondence between the computed and experimental (not derived from* ab initio* calculations) values of binding energies, bond lengths, and valence angles of several selected small hydrocarbon molecules (for parameterisation details, see [20, 29]). To validate the model developed, we computed the binding energies, bond lengths, and valence angles for various H–C–N–O compounds not used in the fitting procedure and compared them with the experimental data from NIST database [19]. The resulting tight-binding potential is well suited to modeling various H_{k}C_{l}N_{m}O_{n} molecules from small clusters to large systems as graphene, diamond, peptides, polyprismanes, and so forth [20]. For example, NTBM binding energies and structural properties for peptides such as bradykinin and colistin are in good agreement with the corresponding DFT calculations [20]. Unfortunately, we are not aware of experimental data derived for any polyprismane, except for cubane C_{8}H_{8}. For cubane, NTBM model is in excellent agreement with the existent experimental data. For example, experimental C–C and C–H bond lengths are 1.570 and 1.082 Å, respectively, and NTBM model gives 1.571 and 1.097 Å, respectively [30]. Binding energy obtained within the NTBM model (4.42 eV/atom) is also in good agreement with the experimental value (4.47 eV/atom) [30]. In addition, earlier we calculated the binding energies for hexa- and octaprismanes using NTBM and DFT/B3LYP/6-311G levels of theory. The binding energies of hexaprismane and octaprismane in terms of NTBM model were 4.56 and 4.47 eV/atom, respectively. These results agree well with our data found using the DFT/B3LYP/6-311G: 4.50 and 4.40 eV/atom for hexaprismane and octaprismane, respectively [13]. It should be noted that* ab initio* methods (say density functional theory) are very accurate but computer-resources demanding and hence applicable to small systems only. Calculation time within the DFT approach is estimated as , where is the number of atoms in molecular system. To check the convergence of the binding energies and the HOMO-LUMO gaps of “long” polyprismanes, we considered the molecular systems up to 500 carbon atoms. So the geometry and energy optimization of such systems in the frame of DFT is very computer resource-intensive or impossible. Therefore, the tight-binding approach is rather good compromise between the accuracy and computational costs and is commonly used for quantum chemistry modeling, including the calculation of binding energies, bond lengths, valence angles, and HOMO-LUMO gaps [31–33].

The samples studied are [*n*,*m*]prismanes with* m* = 4 and 5 (tetra- and pentaprismanes, resp.) and (see Figure 1). To obtain the equilibrium structures of polyprismanes, we used the method of structural relaxation so that the corresponding initial configuration relaxed to a state with the local or global energy minimum under the influence of intramolecular forces only. First of all, the forces acting on the all atoms were calculated using the Hellmann-Feynman theorem [34]. Next, atoms were shifted in the direction of the forces obtained proportional to the corresponding forces. Then, the relaxation step was repeated. The whole structure is relaxed up to residual atomic forces smaller than eV/Å.