Conformational Analysis and Electronic Structure of Chiral Carbon and Carbon Nitride Nanotubes

Geometry and electronic structure of chiral carbon and carbon nitride (CNx) nanotubes were investigated through quantum chemical methods. Finite nanotubes with diameters ranging from 5 to 10 Å and containing up to 500 atoms were considered. CNx structures were built through random substitution of carbon atoms by nitrogen. The molecules were fully optimized by semi-empirical quantum chemical method (PM3). Our results show that the energy associated with nitrogen incorporation depends strongly upon the tube helicity and diameter. The doping of nanotubes with nitrogen contributes to reduce the stress caused by the small diameter of the studied systems. Density of States (DOS) results for pure carbon and CNx nanostructures, obtained through DFT and Hartree-Fock calculations, were analyzed. The introduction of nitrogen in the tube produce states in the gap region which characterizes the metallic behavior, as expected for these systems after N-doping.


Introduction
Since their discovery by Sumio Iijima 1 in 1991, carbon nanotubes (CNT) have been intensively investigated, both theoretically and experimentally, due to their outstanding electronic and mechanical properties. These systems are usually treated as one-dimensional semiconductors or metals, depending on the geometry of the tubes 2 . The classification of single walled carbon nanotubes (SWNTs) is based on two chiral indices (n,m) which give the geometry of the basic graphene ribbon that is rolled to form a nanotube. According to the usual nomenclature, nanotubes are said to be achiral when one of the indices is zero (zig-zag) or when n = m (armchair) and all the others are chiral. It has been proposed that CNTs behave as 1-D conductors when the difference between the chiral indices is a multiple of 3: n -m = 3q, where q is an integer 2 .
Semiconducting nanotubes are of interest in the fabrication of electronic devices as they combine the outstanding mechanical properties of small band gap semiconductors, altogether in systems of nanoscopic dimensions. Metallic nanotubes are interesting as well since they are prototypes of mechanically robust molecular wires. The development of experimental techniques that precisely synthesize carbon nanotubes with uniform helicity and electronic properties is still a challenge. This fact might impose great limitations on the technological applications of these nanostructures. Theoretical 3 and experimental 4-6 studies have shown that is possible to modify the electronic properties of the nanotubes by replacing some of carbon atoms with heteroatoms 7 . Futhermore the incorporation of these heteroatoms also changes the nanotube estructure 8,9 , chemical reactivity 10 and mechanical properties 11 , presenting the possibility of controlling nanotube properties. In recent years, various synthesis methods to produce CNx nanotubes have been reported, including chemical vapor deposition (CVD) 12 In our previous work 15,28 , simulations of a random doping of CNTs showed that substitutional nitrogen on the hexagonal carbon network produces localized distortions on the tube walls. Analyzing nitrogen incorporation energy results for molecules of similar diameter, we conclude that carbon atoms are more easily substituted by nitrogen atoms in zigzag than in armchair nanotubes.
In the present work we report a quantum chemical study on chiral nanotubes with diameters varying from 5 to 10 Å. We analyze the role played by nitrogen doping in the stabilization of these molecular systems. In the case of chiral nanotubes our theoretical results showed that the energies of nitrogen incorporation are close to the calculated for armchair nanotubes regardless of its helicity. The modifications in the electronic structure due to random nitrogen substitution are also analyzed. A band associated with donor states emerges below the bottom of the conduction band after doping.

Computational Details
The geometry of tubular structures composed by carbon and nitrogen containing form 100 to 500 atoms were fully optimized through the semi-empirical quantum chemical method Parametric Method 3 (PM3) 29 . PM3 is a semi-empirical method derived from the Hartree-Fock theory. The advantages of semiempirical calculations are that they are much faster than ab initio calculations, and can be used for large organic molecules. The disadvantage of semiempirical calculations is that some properties cannot be predicted reliably. In the case of the properties analyzed in this study, PM3 semiempirical method is very reliable to predict molecular geometries and heats of formation of carbon materials. PM3 error in heats of formation is about 8.0 kcal.mol -1 [30] , with respect to the experimental values. Average error in bond length is 0.05 Å [30] .
Terminal bonds at tube ends were saturated with hydrogen atoms. Chiral nanotubes were then nitrogen-doped and the geometries were re-optimized. Nitrogen atoms were randomly placed substituting carbons at a given concentrations. For these substitutions, we adopted the following criteria: (i) adjacent atoms should not be substituted; (ii) the substitution of even number of atoms is preferable because a closed shell system is formed.
The energy associated to nitrogen incorporation was calculated as the difference in formation enthalpy of N-doped and pure carbon systems divided by the number of nitrogens. These calculations were performed within the quantum chemical package GAMESS 31 .
The electronic structure of the optimized tubular molecules was obtained through Hartree-Fock theory adopting CEP-4G basis set and the functional BLYP in the 6-31G basis set as implemented in GAUSSIAN03 package 32 . Becke's 33 exchange functional along with the correlation functional of Lee, Yang, and Parr 34 (BLYP) has been used successful in the electronic properties calculations of carbon nanotubes. However, this functional has been shown to be not allowed to large systems as some model molecules analyzed in this work. The Compact Effective Potencial (CEP) proposed by Stevens and co-workers in 1984 [35] is a pseudo-orbital basis set which consists of four Gaussian primitives using a common set of exponents for the s and p functions. Comparing with complete basis sets, these core potentials for the first-and second-row atoms 35 showed good results for electronic properties on large molecules. Density of States was convoluted using Gaussian line shapes with a half-width broadening of 0.2 eV. The DOS profiles were used for the interpretation of the structure evolution upon N doping of chiral carbon nanotubes.

Structural properties and enthalpy of formation
In Figure 1 the geometries of the chiral nanotubes considered in the present study are schematically shown. These corresponds to the optimized conformations of tubes (5.3), (6.3), (6.4), (7.4), (8.5), and (9.3), which have diameters from 5 to 9 Å.
The random substitution of carbon by nitrogen atoms produces geometry distortions around the nitrogen atoms and become more important as more nitrogen atoms are incorporated, transforming the straight tube walls into corrugated walls ( Figure 2). As can be observed in Figure 1, before doping the model molecules do not show irregularities along the tube walls. Similar results were already obtained in our early works 28,36 . In these works 28,36 we showed that N-doped zig-zag tubes are more stable (the formation enthalpy decreases) than armchair tubes of similar diameter. At the same time the stressed small diameter tubes are more easily doped by nitrogen than the large diameter tubes.
In Figure 3 we show the evolution of the enthalpy associated in nitrogen incorporation as a function of the [N]/[C] concentration for chiral and achiral carbon nitride nanotubes. The enthalpy of incorporation was calculated as the relative enthalpy per added nitrogen ([∆H f (C n-x N x ) -∆H f (C n )]/x), where ∆H f is the heat of formation obtained from PM3 calculations. Full square symbols in this picture indicate the results for armchair nanotubes, and full cycle symbols are for zig-zag nanotubes. Open symbols are for chiral nanotubes. Analyzing the energy of incorporation of nitrogen in the chiral nanotubes studied, we conclude that nanotubes with small diameters are more easily doped by nitrogen than the large ones.
Moreover, comparing zigzag, armchair, and chiral nanotubes our theoretical results showed that the enthalpy of formation for N-doped chiral nanotubes is similar to the armchair nanotubes, independently of their helicity.
Based on the results obtained for relative enthalpy of formation, we concluded that the energy associated to incorporation of nitrogen depends on the tubes diameter and helicity for both chiral and achiral carbon nitride nanotubes.   properties and chemical reactivity of carbon nanotubes 37,38 . Nitrogendoped carbon nanotubes render them n-type material regardless of tube chirality 39 , and nitrogen in doped carbon nanotubes acts as a donor. In particular, nitrogen-doped carbon nanotubes are effective field emitters [40][41][42][43] . CNx nanotubes have been demonstrated to show lower turn-on and higher current density than undoped carbon nanotubes 42 . Incorporation of nitrogen atoms on carbon nanotubes changes their tubular structures. Simulations of a random doping of carbon nanotubes showed that substitutional nitrogen on the hexagonal carbon network produces localized distortions on the tube walls. However, the corrugation does not produce the bends needed to interpret the bamboo-like structure that are experimentally observed. Y-junction can be created introducing defects in the hexagonal network of bonds 35,36 . Theoretical studies have proposed different types of defects in CNx nanostructures as, for example: (i) odd-membered rings [44][45][46] , and (ii) vacancies, as pyridine-like 47,48 , and porphyin-like 49 structures. Like unbranched nanotubes, the Y-junction 50 CNx nanotube is also considered to be one of the very promising materials in the upcoming field of nanoelectronics. Although a number of theoretical and experimental studies have been carried out on the electronic and structural properties of carbon nitride nanotubes, some issues are not clear, such as the electronic properties dependence as a function of the concentration and distribution of nitrogen in N-doped nanostructures.
In this study, PM3 optimized geometries of carbon and CNx nanotubes were used as input for electronic properties calculations at the HF/CEP-4G and BLYP/6-31G level of theory. Results of density of states (DOS) obtained by BLYP functional are displayed in Figure 4a for chiral nanotube (6.3). Calculations of electronic properties on DFT level of theory showed to be not allowed to large molecules analyzed, as (5.3), (6.4), (8.5), (9.3) and (9.6) nanotubes. In this sense, it was recalculated the electronic properties of small diameter nanotubes using HF/CEP-4G level of theory and compared with previous results in DFT level theory. CEP-4G basis gave

Electronic Properties
The doping of carbon nanotubes with nitrogen atoms is a practical and effective way to tailor their electronic properties. Incorporation of nitrogen atoms on carbon nanotubes changes the electrical electronic properties in qualitative good agreement with the results of more accurate method for small diameter nanotubes. In this study, HF/CEP-4G basis was used to determinate the electronic properties of large model molecules.
Results of density of states calculations are shown in Figure 4a for carbon nanotubes and carbon nitride nanotubes. In Figure 4a, the curve below represents the DOS of chiral carbon nanotube. The other curves are related to carbon nitride nanotubes with some selected nitrogen (N)/(C) concentration. The introduction of nitrogen in chiral tubes results: (i) a new peak in the top of the valence band. This structure is associated to the nitrogen incorporation which contributes with a pair of electrons; (ii) the main peak is continuously shifted toward lower energies as more nitrogen content is increased. In Figure 4b, we show an inset of the density of states in the gap region for (6.3) chiral nanotube. This picture shows the last highest occupied molecular orbitals (HOMO) and the first lowest unoccupied molecular orbitals (LUMO) for (6.3) carbon and carbon nitride nanotube. Results for frontier orbitals energy using HF or DFT methods, to the all chiral and achiral nanotubes studied, resulted in the decreasing of HOMO-LUMO gap after nitrogen substitution. This indicates that the introduction of nitrogen in the tube produces states in the gap region which characterizes the metallic behavior, as expected for these nanostructures after N-doping.

Conclusions
In summary, we have presented structural and electronic properties calculations of nanotubes of differente chirallity. Results of enthalpy of formation for chiral nanotubes showed that small diameter structures are more easily doped by nitrogen than large diameter tubes. Structural stress caused due to effects of curvature is reduced by the nitrogen incorporation, which contributes to the energy stabilization of small diameter tubes. Comparing the results of heat of formation for zigzag, armchair and chiral carbon nitride nanotubes we concluded that the energy associated to incorporation of nitrogen depends of the tubes diameter and chirality.
Calculations based on HF and DFT methods demonstrate that the introduction of nitrogen generates new states in the gap region characterizing the metallic behavior, as expected for CNx nanotubes.